This paper investigates the applicability of model emulation to speed up simulation time of CPU intensive environmental models. Polynomial chaos expansion (PCE) emulators are constructed for three ...case studies of increasing complexity. The level of emulator training and the order of polynomial necessary to sufficiently build accurate emulators for each model are investigated. Although the PCE emulators shown here do not approximate well the outputs of parameter rich models (80 + parameters), results demonstrate that the emulators mimic closely outputs of relatively simple, low dimensional, simulation models (15 parameters or less). Furthermore, the PCE emulators are tested with applications such as Global Sensitivity Analysis (GSA). Results illustrate the advantages and drawbacks of using classical PCE emulators for treating computational limitation of complex environmental models.
•PCE emulators are constructed with sparse grid sampling for 3 environmental models.•Accuracy, efficiency, and fidelity of emulators were tested.•For models that have low to medium complexity, the emulators worked well.•For a hyper-complex model with 87 parameters, the emulators were relatively weak.•Dimensionality reduction of the hyper-complex model allowed for better emulation.
Accurate and detailed knowledge of California's groundwater is of paramount importance for statewide water resources planning and management, and to sustain a multi-billion-dollar agriculture ...industry during prolonged droughts. In this study, we use water supply and demand information from California's Department of Water Resources to develop an aggregate groundwater storage model for California's Central Valley. The model is evaluated against 34 years of historic estimates of changes in groundwater storage derived from the United States Geological Survey's Central Valley Hydrologic Model (USGS CVHM) and NASA's Gravity Recovery and Climate Experiment (NASA GRACE) satellites. The calibrated model is then applied to predict future changes in groundwater storage for the years 2015-2050 under various precipitation scenarios from downscaled climate projections. We also discuss and project potential management strategies across different annual supply and demand variables and how they affect changes in groundwater storage. All simulations support the need for collective statewide management intervention to prevent continued depletion of groundwater availability.
The daunting complexity of ecosystems has led ecologists to use mathematical modelling to gain understanding of ecological relationships, processes and dynamics. In pursuit of mathematical ...tractability, these models use simplified descriptions of key patterns, processes and relationships observed in nature. In contrast, ecological data are often complex, scale‐dependent, space‐time correlated, and governed by nonlinear relations between organisms and their environment. This disparity in complexity between ecosystem models and data has created a large gap in ecology between model and data‐driven approaches. Here, we explore data assimilation (DA) with the Ensemble Kalman filter to fuse a two‐predator‐two‐prey model with abundance data from a 2600+ day experiment of a plankton community. We analyse how frequently we must assimilate measured abundances to predict accurately population dynamics, and benchmark our population model's forecast horizon against a simple null model. Results demonstrate that DA enhances the predictability and forecast horizon of complex community dynamics.
Oscillations in global modes of variability (MoV) form global teleconnections that affect regional climate variability and modify the potential for severe and damaging weather conditions. ...Un-derstanding the link between certain MoVs and regional climate can improve the ability to more accurately predict environmental conditions that impact human life and health. In this study, we explore the connection between different MoVs, including the Arctic Oscillation (AO), Eurasian teleconnection, Indian Ocean Dipole (IOD), North Atlantic Oscillation (NAO), and El Niño Southern Oscillation (Nino34), with winter and summer climate in the High Mountain Asia (HMA) region, including geopotential height at 250 hPa (z250), 2-m air temperature (T2M), total precipitation (PRECTOT), and fractional snow cover area (fSCA). Relationships are explored for the same monthly period between the MoVs and the climate variables, and also using a lagged correlation analysis to investigate whether any relationship exists at different time lags. We find that T2M has a negative correlation with the Eurasian teleconnection in the Inner Tibetan Plateau and Central China in both winter and summer and a positive correlation in Western China in summer. PRECTOT has a positive correlation with all MoV in most regions in winter, especially with the IOD, and a negative correlation in summer, especially with the Eurasian teleconnection. Snow cover in winter is positively correlated with most indices throughout many regions in HMA, likely due to wintertime precipitation also being positively correlated with most indices. Generally, the AO and NAO show similar correlation patterns with all climate variables, espe-cially in the winter, possibly due to their oscillations being so similar. Furthermore, the AO and NAO are shown to be less significant in explaining the variation in HMA climate compared to other MoVs such as the Eurasian teleconnection. Overall, our results identify different time-windows and specific regions within HMA that exhibit high correlation between climate and MoVs, which might offer additional predictability of the MoVs as well as of climate and weather patterns in HMA and throughout the globe.
Using climate model ensembles containing members that exhibit very high climate sensitivities to increasing CO
concentrations can result in biased projections. Various methods have been proposed to ...ameliorate this 'hot model' problem, such as model emulators or model culling. Here, we utilize Bayesian Model Averaging as a framework to address this problem without resorting to outright rejection of models from the ensemble. Taking advantage of multiple lines of evidence used to construct the best estimate of the earth's climate sensitivity, the Bayesian Model Averaging framework produces an unbiased posterior probability distribution of model weights. The updated multi-model ensemble projects end-of-century global mean surface temperature increases of 2
C for a low emissions scenario (SSP1-2.6) and 5
C for a high emissions scenario (SSP5-8.5). These estimates are lower than those produced using a simple multi-model mean for the CMIP6 ensemble. The results are also similar to results from a model culling approach, but retain some weight on low-probability models, allowing for consideration of the possibility that the true value could lie at the extremes of the assessed distribution. Our results showcase Bayesian Model Averaging as a path forward to project future climate change that is commensurate with the available scientific evidence.
Regions with high productivity of agriculture, such as the Beqaa Plain, Lebanon, often rely on groundwater supplies for irrigation demand. Recent reports have indicated that groundwater consumption ...in this region has been unsustainable, and quantifying rates of groundwater depletion has remained a challenge. Here, we utilize 15 years of data (June 2002–April 2017) from the Gravity Recovery and Climate Experiment (GRACE) satellite mission to show Total Water Storage (TWS) changes in Lebanon’s Beqaa Plain. We then obtain complimentary information on various hydrologic cycle variables, such as soil moisture storage, snow water equivalent, and canopy water storage from the Global Land Data Assimilation System (GLDAS) model, and surface water data from the largest body of water in this region, the Qaraaoun Reservoir, to disentangle the TWS signal and calculate groundwater storage changes. After combining the information from the remaining hydrologic cycle variables, we determine that the majority of the losses in TWS are due to groundwater depletion in the Beqaa Plain. Results show that the rate of groundwater storage change in the West Beqaa is nearly +0.08 cm/year, in the Rashaya District is −0.01 cm/year, and in the Zahle District the level of depletion is roughly −1.10 cm/year. Results are confirmed using Sentinel-1 interferometric synthetic aperture radar (InSAR) data, which provide high-precision measurements of land subsidence changes caused by intense groundwater usage. Furthermore, data from local monitoring wells are utilized to further showcase the significant drop in groundwater level that is occurring through much of the region. For monitoring groundwater storage changes, our recommendation is to combine various data sources, and in areas where groundwater measurements are lacking, we especially recommend the use of data from remote sensing.
The seasonal-to-decadal terrestrial water balance on river basin scales depends on several well-characterized but uncertain soil physical processes, including soil moisture, plant available water, ...rooting depth, and recharge to lower soil layers. Reducing uncertainties in these quantities using observations is a key step toward improving the data
fidelity and skill of land surface models. In this study, we quantitatively
characterize the capability of Gravity Recovery and Climate Experiment (NASA-GRACE) measurements – a key constraint on total water storage (TWS) – to inform and constrain these processes. We use a reduced-complexity physically based model capable of simulating the hydrologic cycle, and we apply Bayesian inference on the model parameters using a Markov chain Monte Carlo algorithm, to minimize mismatches between model-simulated and GRACE-observed TWS anomalies. Based on the prior and posterior model parameter distributions, we further quantify information gain with regard to terrestrial water states, associated fluxes, and time-invariant process parameters. We show that the data-constrained terrestrial water storage model can capture basic physics of the hydrologic cycle for a watershed in the western Amazon during the period January 2003 through December 2012, with an r2 of 0.98 and root mean square error of 30.99 mm between observed and simulated TWS. Furthermore, we show a reduction of uncertainty in many of the parameters and state variables, ranging from a 2 % reduction in uncertainty for the porosity parameter to an 85 % reduction for the rooting depth parameter. The annual and interannual variability of the system are also simulated accurately, with the model simulations capturing the impacts of the 2005–2006 and 2010–2011 South American droughts. The results shown here suggest the potential of using gravimetric observations of TWS to identify and constrain key parameters in soil hydrologic models.
Given the increasing use of climate projections and multi-model ensemble weighting for a diverse array of applications, this project assesses the sensitivities of climate model weighting strategies ...and their resulting ensemble means to multiple components, such as the weighting schemes, climate variables, or spatial domains of interest. The purpose of this study is to assess the sensitivities associated with multi-model weighting strategies. The analysis makes use of global climate models from the Coupled Model Intercomparison Project Phase 5 (CMIP5) and their statistically downscaled counterparts created with the localized constructed analogs (LOCA) method. This work focuses on historical and projected future mean precipitation and daily high temperatures of the south-central United States. Results suggest that the model weights and the corresponding weighted model means can be sensitive to the weighting strategy that is applied. For instance, when estimating model weights based on Louisiana precipitation, the weighted projections show a wetter and cooler south-central domain in the future compared to other weighting strategies. Alternatively, for example, when estimating model weights based on New Mexico temperature, the weighted projections show a drier and warmer south-central domain in the future. However, when considering the entire south-central domain in estimating the model weights, the weighted future projections show a compromise in the precipitation and temperature estimates. As for uncertainty, our matrix of results provided a more certain picture of future climate compared to the spread in the original model ensemble. If future impact assessments utilize weighting strategies, then our findings suggest that how the specific weighting strategy is used with climate projections may depend on the needs of an impact assessment or adaptation plan.
Recently, assessments of global climate model (GCM) ensembles have transitioned from using unweighted means to weighted means designed to account for skill and interdependence among models. Although ...ensemble-weighting schemes are typically derived using a GCM ensemble, statistically downscaled projections are used in climate change assessments. This study applies four ensemble-weighting schemes for model averaging to precipitation projections in the south-central United States. The weighting schemes are applied to (1) a 26-member GCM ensemble and (2) those 26 members downscaled using Localized Canonical Analogs (LOCA). This study is distinct from prior research because it compares the interactions of ensemble-weighting schemes with GCMs and statistical downscaling to produce summarized climate projection products. The analysis indicates that statistical downscaling improves the ensemble accuracy (LOCA average root mean square error is 100 mm less than the CMIP5 average root mean square error) and reduces the uncertainty of the projected ensemble-mean change. Furthermore, averaging the LOCA ensemble using Bayesian Model Averaging reduces the uncertainty beyond any other combination of weighting schemes and ensemble (standard deviation of the mean projected change in the domain is reduced by 40–50 mm). The results also indicate that it is inappropriate to assume that a weighting scheme derived from a GCM ensemble matches the same weights derived using a downscaled ensemble.
Seasonal variability of the global hydrologic cycle directly impacts human activities, including hazard assessment and mitigation, agricultural decisions, and water resources management. This is ...particularly true across the High Mountain Asia (HMA) region, where availability of water resources can change depending on local seasonality of the hydrologic cycle. Forecasting the atmospheric states and surface conditions, including hydrometeorological relevant variables, at subseasonal-to-seasonal (S2S) lead times of weeks-to-months is an area of active research and development. NASA’s 15 Goddard Earth Observing System (GEOS) S2S prediction system has been developed with this research goal in mind. Here, we benchmark the forecast skill of GEOS-S2S (version 2) hydrometeorological forecasts at 1-3 month lead times in the HMA region, including a portion of the Indian Subcontinent, during the retrospective forecast period, 1981-2016. To assess forecast skill, we evaluate 2-m air temperature, total precipitation, fractional snow cover, snow water equivalent, surface soil moisture, and terrestrial water storage forecasts against the Modern-Era Retrospective analysis for Research and 20 Applications, Version 2 (MERRA-2) and independent reanalysis data, satellite observations, and data fusion products. Anomaly correlation is highest when the forecasts are evaluated against MERRA-2 and particularly in variables with long memory in the climate system, likely due to similar initial conditions and model architecture used in GEOS-S2S and MERRA-2. When compared to MERRA-2, results for the 1-month forecast skill range from anomaly correlation of Ranom=0.18 for precipitation to Ranom=0.62 for soil moisture. Anomaly correlations are consistently lower when forecasts are 25 evaluated against independent observations; results for the 1-month forecast skill range from Ranom=0.13 for snow water equivalent to Ranom=0.24 for fractional snow cover. We find that, generally, hydrometeorological forecast skill is dependent on the forecast lead time, the memory of the variable within the physical system, and the validation dataset used. Overall, these results benchmark the GEOS-S2S system’s ability to forecast HMA hydrometeorology.