Advanced Quantitative Precipitation Information: Improving Monitoring and Forecasts of Precipitation, Streamflow, and Coastal Flooding in the San Francisco Bay Area

American Meteorological Society (AMS) | October 26th, 2022

Advanced Quantitative Precipitation Information (AQPI) is a synergistic project that combines observations and models to improve monitoring and forecasts of precipitation

Advanced Quantitative Precipitation Information (AQPI) is a synergistic project that combines observations and models to improve monitoring and forecasts of precipitation, streamflow, and coastal flooding in the San Francisco Bay area. As an experimental system, AQPI leverages more than a decade of research, innovation, and implementation of a statewide, state-of-the-art network of observations, and development of the next generation of weather and coastal forecast models. AQPI was developed as a prototype in response to requests from the water management community for improved information on precipitation, riverine, and coastal conditions to inform their decision making processes. Observation of precipitation in the complex Bay Area landscape of California’s coastal mountain ranges is known to be a challenging problem. But, with new advanced radar network techniques, AQPI is helping fill an important observational gap for this highly populated and vulnerable metropolitan area. The prototype AQPI system consists of improved weather radar data for precipitation estimation; additional surface measurements of precipitation, streamflow and soil moisture; and a suite of integrated forecast modeling systems to improve situational awareness about current and future water conditions from sky to sea. Together these tools will help improve emergency preparedness and public response to prevent loss of life and destruction of property during extreme storms accompanied by heavy precipitation and high coastal water levels - especially high-moisture laden atmospheric rivers. The Bay Area AQPI system could potentially be replicated in other urban regions in California, the United States, and world-wide.

Advancing Hydro-Economic Optimization to Identify Vulnerabilities and Adaptation Opportunities in California’s Water System

California Natural Resources Agency (CNRA) | August 27th, 2018

Long-term shifts in the timing and magnitude of reservoir inflows will affect water supply reliability in California. Hydro-economic models can help explore climate chang

Long-term shifts in the timing and magnitude of reservoir inflows will affect water supply reliability in California. Hydro-economic models can help explore climate change concerns by identifying system vulnerabilities and adaptation strategies for statewide water operations. This work contributes a new open-source limited foresight implementation of the CALVIN model, a hydro-economic model including roughly 90% of California’s urban and agricultural water demands. The model includes the ability to determine water allocations on an annual basis without knowledge of future availability, and to efficiently evaluate ensembles of streamflow projections representing a range of possible future climates. We then assess the vulnerability of the statewide system to changes in total annual runoff and the fraction of runoff occurring during winter, which primarily depends on temperature. Results are analyzed with a focus on adaptation strategies, aided by the economic representation of water demand in the model. These strategies include changes to reservoir operating policies, and conveyance and storage expansion. As water availability decreases, model results show quadratic increases in shortage cost, and corresponding increases in the marginal costs of adaptation strategies and environmental flow constraints. Reservoirs adapt to warmer climates by increasing average storage levels in winter and routing excess runoff to downstream reservoirs with available capacity. Both small and large changes to reservoir operations were observed compared to historical hydrology, showing that no single operating strategy achieves optimality for all reservoirs. Increasing the fraction of winter flow causes small increases in total shortage cost, indicating the ability to manage a changing hydrologic regime with adaptive reservoir operations. The results of this project improve estimates of the cost of climate change to California’s water system under a range of future conditions and highlight adaptation strategies that minimize costs of increased hydrologic variability. 

Agricultural Damages and Losses from ARkStorm Scenario Flooding in California

American Society of Civil Engineers (ASCE) | November 17th, 2016

Scientists designed the ARkStorm scenario to challenge the preparedness of California communities for widespread flooding with a historical precedence and increased l

Scientists designed the ARkStorm scenario to challenge the preparedness of California communities for widespread flooding with a historical precedence and increased likelihood under climate change. California is an important provider of vegetables, fruits, nuts, and other agricultural products to the nation. This study analyzes the agricultural damages and losses pertaining to annual crops, perennial crops, and livestock in California exposed to ARkStorm flooding. Statewide, flood damage is incurred on approximately 23% of annual crop acreage, 5% of perennial crop acreage, and 5% of livestock, e.g., dairy, feedlot, and poultry, acreage. The sum of field repair costs, forgone income, and product replacement costs span $3.7 and $7.1 billion (2009) for a range of inundation durations. Perennial crop loss estimates dominate, and the vulnerability of orchards and vineyards has likely increased with recent expansion. Crop reestablishment delays from levee repair and dewatering more than double annual crop losses in the delta islands, assuming the fragile system does not remain permanently flooded. The exposure of almost 200,000 dairy cows to ARkStorm flooding poses livestock evacuation challenges.

American and Sacramento Rivers, California, Erodibility Measurements and Model

U.S. Geological Survey (USGS) | August 11th, 2021

A previous report by the authors described sediment sampling and drilling by the U.S. Geological Survey (USGS) beside the American and Sacramento Rivers near Sacramento,

A previous report by the authors described sediment sampling and drilling by the U.S. Geological Survey (USGS) beside the American and Sacramento Rivers near Sacramento, California, in support of a U.S. Army Corps of Engineers project focused on regional flood control. The drilling was performed to define lithology, extract samples for laboratory testing, and perform borehole erosion tests (BETs). The U.S. Department of Agriculture (USDA) performed jet erodibility tests (JETs) near each drilling site, and a team from Texas A&M University performed laboratory tests with an erosion function apparatus (EFA). Collectively, the effort was intended to reveal spatial variations in sediment erodibility and provide data for use in a model to simulate morphological response to a major flood. The data collected by the USGS are available in a public data release. This report, developed in cooperation with the U.S. Army Corps of Engineers, provides comparisons of the three types of measurements of the erodibility of riverbed sediments. The BET is performed in the field and reveals erodibility of sediments below the bed surface. The JET is likewise performed in the field but reveals only erodibility of exposed sediments. The EFA test is done in the laboratory and was performed on soils extracted from different depths beneath the bed surface, in many cases reconstituted for laboratory testing. Tests were performed at nominally similar locations but differed by meters to tens of meters in horizontal locations. The comparison was undertaken to investigate differences among results obtained by the individual measurement approaches and to elucidate pros and cons of each method. The critical shear stress to initiate erosion and the rate of change of erosion rate per unit increase of excess shear stress, sometimes referred to as the erosion coefficient, served as the primary basis for comparison. The three test methods in some cases resulted in order of magnitude differences in estimates of these parameters. Some differences could be attributed to variances in site location or result from testing surface sediment versus a deeper layer, but systematic differences are also evident in the results. The tests performed in the laboratory using the EFA resulted in much lower values of critical shear stress and much higher values of the erosion coefficient compared to the JET tests performed by the USDA team on surface sediments. Critical shear stress was poorly resolved in the BET results because of the limited number of results per site, but the erosion coefficients derived from BET results were systematically lower than those obtained using the EFA. A new, simplified approach is also proposed to estimate the increase in channel cross-sectional area during a large flood, given data describing the initial river cross section, riverbed erodibility parameters, and peak flood discharge and duration. The model runs until the cross section erodes to an equilibrium condition or the flood ends. Output describes the area of the cross section at the end of the simulation and the time required to reach equilibrium if it was reached within the simulated period. The model assumes unique, constant values for both the critical shear stress and the erosion coefficient and represents the fluid mechanics in a simplified way, making it of limited value for quantitative predictions. It does, however, provide an indication of which cross sections are most likely to undergo the greatest change in the design event and can be used to investigate sensitivity of erosion predictions to variability in sediment erodibility measurements.

An assessment of future tidal marsh resilience in the San Francisco Estuary through modeling and quantifiable metrics of sustainability

Frontiers in Environmental Science (Frontiers) | November 25th, 2022

Quantitative, broadly applicable metrics of resilience are needed to effectively manage tidal marshes into the future. Here we quantified three metrics of temporal marsh

Quantitative, broadly applicable metrics of resilience are needed to effectively manage tidal marshes into the future. Here we quantified three metrics of temporal marsh resilience: time to marsh drowning, time to marsh tipping point, and the probability of a regime shift, defined as the conditional probability of a transition to an alternative super-optimal, suboptimal, or drowned state. We used organic matter content (loss on ignition, LOI) and peat age combined with the Coastal Wetland Equilibrium Model (CWEM) to track wetland development and resilience under different sea-level rise scenarios in the Sacramento-San Joaquin Delta (Delta) of California. A 100-year hindcast of the model showed excellent agreement (R2 = 0.96) between observed (2.86 mm/year) and predicted vertical accretion rates (2.98 mm/year) and correctly predicted a recovery in LOI (R2 = 0.76) after the California Gold Rush. Vertical accretion in the tidal freshwater marshes of the Delta is dominated by organic production. The large elevation range of the vegetation combined with high relative marsh elevation provides Delta marshes with resilience and elevation capital sufficiently great to tolerate centenary sea-level rise (CLSR) as high as 200 cm. The initial relative elevation of a marsh was a strong determinant of marsh survival time and tipping point. For a Delta marsh of average elevation, the tipping point at which vertical accretion no longer keeps up with the rate of sea-level rise is 50 years or more. Simulated, triennial additions of 6 mm of sediment via episodic atmospheric rivers increased the proportion of marshes surviving from 51% to 72% and decreased the proportion drowning from 49% to 28%. Our temporal metrics provide critical time frames for adaptively managing marshes, restoring marshes with the best chance of survival, and seizing opportunities for establishing migration corridors, which are all essential for safeguarding future habitats for sensitive species.

An Assessment of Potential Severe Droughts in the Colorado River Basin

American Water Resources Association (AWRA) | September 15th, 2022

An evaluation of atmospheric rivers over the North Pacific in CMIP5 and their response to warming under RCP 8.5

American Geophysical Union (AGU) | November 13th, 2015

Landfalling atmospheric rivers over the North Pacific are evaluated in historical (1980–2005) simulations from 28 models participating in the fifth phase of the Couple

Landfalling atmospheric rivers over the North Pacific are evaluated in historical (1980–2005) simulations from 28 models participating in the fifth phase of the Coupled Model Intercomparison Project (CMIP5) and compared to the Modern-Era Retrospective Analysis for Research and Applications and ERA-Interim reanalyses. Landfalling dates are identified as having spatially elongated high wind and moisture features in physical proximity to the coastline. Model performance relative to reanalysis is found to be quite variable. The majority can resolve the spatial structure of landfalling events, but few correctly resolve the frequency distribution, interannual variability in number and amplitude of moisture flux, and median landfalling latitude. The response of a subset of high-performing models to projected warming at the end of the 21st century is investigated using Representative Concentration Pathway (RCP) 8.5 (2070–2100) projections. Selected models show a broadening of the frequency distribution, with the largest increase in frequency occurring equatorward of peak historical frequency. While there is a robust increase in AR-related moisture transport compared to the historical period, this increase is not linked to more extreme AR events within the context of the projected future climate. The occurrence of years with many AR dates is greater at the end of the 21st century than during the historical period. The equatorward increase in peak historical frequency is colocated with increases in the 850 and 250 hPa zonal winds. The response in moisture flux to warming is mostly thermodynamic, but equatorward of its peak distribution, it is dominated by a dynamic response.

An Integrated Surface Water-Groundwater Interaction Model for the Carmel River

California State University, Monterey Bay | April 30th, 2010

The Carmel River supplies fresh water to the residents of the Monterey Peninsula within the water district served by the California American Water company (Cal-Am). The S

The Carmel River supplies fresh water to the residents of the Monterey Peninsula within the water district served by the California American Water company (Cal-Am). The State Water Resources Control Board, in 1995, ordered that Cal-Am reduce their annual diversions from the Carmel River Watershed to within legal right (3,376 acre-feet per year). An alternative water source to the Carmel River has yet to be determined, and as a result, Cal-Am must continue to divert more water than their legal entitlement. The current magnitude of diversions has a negative impact on the spawning and migrating habitat of Steelhead (Oncorhynchus mykiss). Existing software in the Tarsier Environmental Modeling Framework was used to model the spatial distribution of surface water along the Carmel River. The model simulated the flow of water downstream from catchment area to the Pacific Ocean. The river channel was represented by a network data set comprised of links, representing individual reaches of the river, connected by nodes. Prior to my work, systematic error in the model was thought to be partially the result of the model lacking a simulation of the interactions between the surface water and the underlying aquifer. A groundwater sub-model was developed to correct for the systematic error. The groundwater sub-model simulated the movement of water between the river channel and the aquifer. Stock variables representing a shallow and deep aquifer were added to each link of the network data set. Simulated water in the surface water stock of each link percolates to these aquifer stocks until the groundwater reaches aquifer capacity, allowing surface water to continue flowing downstream. The model also allows the lateral flow of groundwater according to Darcy’s Law. Quantitative and qualitative analysis of model output compared to observed data showed an increase in model accuracy. Quantitatively, a Nash- Sutcliff Coefficient was improved from 0.88 to 0.97 with the addition of the groundwater model. Qualitatively, a visualization of the longitudinal profile of the river system showed the simulated aquifer controlling the surface flow. A hypothetical application of this model is presented where reducing the pumping rate from the aquifer allowed the wetted river channel to increase by 2.5 km. Future work on the model should include a reservoir sub-model and accounting for spatial variability in the precipitation throughout the catchment area. With these improvements it will be possible to improve predictions of the spatial distribution of surface water along the Carmel River given hypothetical scenarios such as the rate and spatial distribution of pumping from the Carmel Aquifer. Using these simulations to inform decisions of river management could benefit all stakeholders of the river.

An Introduction to EcoAtlas: Applied Aquatic Science

San Francisco Estuary Institute (SFEI) | September 14th, 2016

The term “EcoAtlas” has been in use for over 20 years in association with the study, assessment, and reporting of aquatic resources in California. What began as “W

The term “EcoAtlas” has been in use for over 20 years in association with the study, assessment, and reporting of aquatic resources in California. What began as “Wetland Tracker” in 2000, with its goal of recording essential information about wetland restoration projects in the greater Bay Area, later blossomed into a broad technological ecosystem of interrelated tools, each focused on delivering specialized, program-focused features, as defined by key stakeholder groups. Having evolved over time, EcoAtlas comprises a diverse toolset, but the collection shares common libraries, common approaches, and common development methodologies. Though each emerged from different funding opportunities and requirement drivers, the toolset collectively addresses a synthesized “whole watershed approach.” Whether estimating the ideal riparian buffer width for a given stream or assessing the health of a wetland at the edge of the estuary, the EcoAtlas tools allow practitioners to deploy the right tool for the job across the entire watershed, thereby producing a complete picture through composite outputs.

An overview of multi-dimensional models of the Sacramento-San Joaquin Delta

University of California, Davis (UC Davis) | December 1st, 2016

Over the past fifteen years, the development and application of models in San Francisco Bay and the Sacramento-San Joaquin Delta has transformed our ability to analyze an

Over the past fifteen years, the development and application of models in San Francisco Bay and the Sacramento-San Joaquin Delta has transformed our ability to analyze and understand the underlying physics of the system. Initial applications of three-dimensional models focused primarily on salt intrusion, and provided a valuable resource for investigating how sea level rise and levee failures in the Delta could influence water quality in the Delta in the future. However, multi-dimensional models have also provided significant insights into some of the fundamental biological relationships that have shaped our thinking about the system by exploring the relationship between X2, flow, fish abundance, and the low salinity zone. Through modeling, it has been possible to move beyond salinity to understand how large-scale changes to the system are likely to affect sediment dynamics, and to assess the potential effects on species that rely on turbidity for habitat. Lastly, the coupling of multi-dimensional hydrodynamic models with particle tracking models has led to advances in our thinking about residence time, the retention of food organisms, the effect of south Delta exports on larval entrainment, and the pathways and behaviors of salmonids that travel through the Delta. This chapter provides an overview of these recent advances and how they have increased our understanding of the distribution and movement of fish and food organisms. The applications presented serve as a guide to the current state of the science of Delta modeling and provide examples of how we can use multidimensional models to predict how future Delta conditions will affect both fish and water supply.