Solving groundwater depletion in India while achieving food security.

Significant groundwater depletion in regions where grains are procured for public distribution is a primary sustainability challenge in India. We identify specific changes in the Indian Government's Procurement & Distribution System as a primary solution lever. Irrigation, using groundwater, facilitated by subsidized electricity, is seen as vital for meeting India's food security goals. Using over a century of daily climate data and recent spatially detailed economic, crop yield, and related parameters, we use an optimization model to show that by shifting the geographies where crops are procured from and grown, the government's procurement targets could be met on average even without irrigation, while increasing net farm income and arresting groundwater depletion. Allowing irrigation increases the average net farm income by 30%. The associated reduction in electricity subsidies in areas with significant groundwater depletion can help offset the needed spatial re-distribution of farm income, a key political obstacle to changes in the procurement system.

Explaining the trends and variability in the United States tornado records using climate teleconnections and shifts in observational practices.

The annual frequency of tornadoes during 1950-2018 across the major tornado-impacted states were examined and modeled using anthropogenic and large-scale climate covariates in a hierarchical Bayesian inference framework. Anthropogenic factors include increases in population density and better detection systems since the mid-1990s. Large-scale climate variables include El Niño Southern Oscillation (ENSO), Southern Oscillation Index (SOI), North Atlantic Oscillation (NAO), Pacific Decadal Oscillation (PDO), Arctic Oscillation (AO), and Atlantic Multi-decadal Oscillation (AMO). The model provides a robust way of estimating the response coefficients by considering pooling of information across groups of states that belong to Tornado Alley, Dixie Alley, and Other States, thereby reducing their uncertainty. The influence of the anthropogenic factors and the large-scale climate variables are modeled in a nested framework to unravel secular trend from cyclical variability. Population density explains the long-term trend in Dixie Alley. The step-increase induced due to the installation of the Doppler Radar systems explains the long-term trend in Tornado Alley. NAO and the interplay between NAO and ENSO explained the interannual to multi-decadal variability in Tornado Alley. PDO and AMO are also contributing to this multi-time scale variability. SOI and AO explain the cyclical variability in Dixie Alley. This improved understanding of the variability and trends in tornadoes should be of immense value to public planners, businesses, and insurance-based risk management agencies.

Crop switching reduces agricultural losses from climate change in the United States by half under RCP 8.5.

A key strategy for agriculture to adapt to climate change is by switching crops and relocating crop production. We develop an approach to estimate the economic potential of crop reallocation using a Bayesian hierarchical model of yields. We apply the model to six crops in the United States, and show that it outperforms traditional empirical models under cross-validation. The fitted model parameters provide evidence of considerable existing climate adaptation across counties. If crop locations are held constant in the future, total agriculture profits for the six crops will drop by 31% for the temperature patterns of 2070 under RCP 8.5. When crop lands are reallocated to avoid yield decreases and take advantage of yield increases, half of these losses are avoided (16% loss), but 57% of counties are allocated crops different from those currently planted. Our results provide a framework for identifying crop adaptation opportunities, but suggest limits to their potential.

Integrating the social, hydrological and ecological dimensions of freshwater health: The Freshwater Health Index.

Degradation of freshwater ecosystems and the services they provide is a primary cause of increasing water insecurity, raising the need for integrated solutions to freshwater management. While methods for characterizing the multi-faceted challenges of managing freshwater ecosystems abound, they tend to emphasize either social or ecological dimensions and fall short of being truly integrative. This paper suggests that management for sustainability of freshwater systems needs to consider the linkages between human water uses, freshwater ecosystems and governance. We present a conceptualization of freshwater resources as part of an integrated social-ecological system and propose a set of corresponding indicators to monitor freshwater ecosystem health and to highlight priorities for management. We demonstrate an application of this new framework -the Freshwater Health Index (FHI) - in the Dongjiang River Basin in southern China, where stakeholders are addressing multiple and conflicting freshwater demands. By combining empirical and modeled datasets with surveys to gauge stakeholders' preferences and elicit expert information about governance mechanisms, the FHI helps stakeholders understand the status of freshwater ecosystems in their basin, how ecosystems are being manipulated to enhance or decrease water-related services, and how well the existing water resource management regime is equipped to govern these dynamics over time. This framework helps to operationalize a truly integrated approach to water resource management by recognizing the interplay between governance, stakeholders, freshwater ecosystems and the services they provide.

Six centuries of Upper Indus Basin streamflow variability and its climatic drivers.

Our understanding of the full range of natural variability in streamflow, including how modern flow compares to the past, is poorly understood for the Upper Indus Basin (UIB) because of short instrumental gauge records. To help address this challenge, we use Hierarchical Bayesian Regression (HBR) with partial pooling to develop six centuries long (1394-2008 C.E.) streamflow reconstructions at three UIB gauges (Doyian, Gilgit, and Kachora), concurrently demonstrating that HBR can be used to reconstruct short records with interspersed missing data. At one gauge (Partab Bridge), with a longer instrumental record (47 years), we develop reconstructions using both Bayesian Regression (BR) and the more conventionally used Principal Components Regression (PCR). The reconstructions produced by PCR and BR at Partab Bridge are nearly identical and yield comparable reconstruction skill statistics, highlighting that the resulting tree-ring reconstruction of streamflow is not dependent on the choice of statistical method. Reconstructions at all four reconstructions indicate flow levels in the 1990s were higher than mean flow for the past six centuries. While streamflow appears most sensitive to accumulated winter (January-March) precipitation and summer (MJJAS) temperature, with warm summers contributing to high flow through increased melt of snow and glaciers, shifts in winter precipitation and summer temperatures cannot explain the anomalously high flow during the 1990s. Regardless, the sensitivity of streamflow to summer temperatures suggests that projected warming may increase streamflow in coming decades, though long-term water risk will additionally depend on changes in snowfall and glacial mass balance.

An Empirical, Nonparametric Simulator for Multivariate Random Variables with Differing Marginal Densities and Nonlinear Dependence with Hydroclimatic Applications.

Multivariate simulations of a set of random variables are often needed for risk analysis. Given a historical data set, the goal is to develop simulations that reproduce the dependence structure in that data set so that the risk of potentially correlated factors can be evaluated. A nonparametric, copula-based simulation approach is developed and exemplified. It can be applied to multiple variables or to spatial fields with arbitrary dependence structures and marginal densities. The nonparametric simulator uses logspline density estimation in the univariate setting, together with a sampling strategy to reproduce dependence across variables or spatial instances, through a nonparametric numerical approximation of the underlying copula function. The multivariate data vectors are assumed to be independent and identically distributed. A synthetic example is provided to illustrate the method, followed by an application to the risk of livestock losses in Mongolia.

Scaling of extreme rainfall areas at a planetary scale.

Event magnitude and area scaling relationships for rainfall over different regions of the world have been presented in the literature for relatively short durations and over relatively small areas. In this paper, we present the first ever results on a global analysis of the scaling characteristics of extreme rainfall areas for durations ranging from 1 to 30 days. Broken power law models are fit in each case. The past work has been focused largely on the time and space scales associated with local and regional convection. The work presented here suggests that power law scaling may also apply to planetary scale phenomenon, such as frontal and monsoonal systems, and their interaction with local moisture recycling. Such features may have persistence over large areas corresponding to extreme rain and regional flood events. As a result, they lead to considerable hazard exposure. A caveat is that methods used for empirical power law identification have difficulties with edge effects due to finite domains. This leads to problems with robust model identification and interpretability of the underlying relationships. We use recent algorithms that aim to address some of these issues in a principled way. Theoretical research that could explain why such results may emerge across the world, as analyzed for the first time in this paper, is needed.

Multimodel ensembles of streamflow forecasts: Role of predictor state in developing optimal combinations.

A new approach for developing multimodel streamflow forecasts is presented. The methodology combines streamflow forecasts from individual models by evaluating their skill, represented by rank probability score (RPS), contingent on the predictor state. Using average RPS estimated over the chosen neighbors in the predictor state space, the methodology assigns higher weights for a model that has better predictability under similar predictor conditions. We assess the performance of the proposed algorithm by developing multimodel streamflow forecasts for Falls Lake Reservoir in the Neuse River Basin, North Carolina (NC), by combining streamflow forecasts developed from two low-dimensional statistical models that use sea-surface temperature conditions as underlying predictors. To evaluate the proposed scheme thoroughly, we consider a total of seven multimodels that include existing multimodel combination techniques such as combining based on long-term predictability of individual models and by simple pooling of ensembles. Detailed nonparametric hypothesis tests comparing the performance of seven multimodels with two individual models show that the reduced RPS from multimodel forecasts developed using the proposed algorithm is statistically significant from the RPSs of individual models and from the RPSs of existing multimodel techniques. The study also shows that adding climatological ensembles improves the multimodel performance resulting in reduced average RPS. Contingency analyses on categorical (tercile) forecasts show that the proposed multimodel combination technique reduces average Brier score and total number of false alarms, resulting in improved reliability of forecasts. However, adding multiple models with climatology also increases the number of missed targets (in comparison to individual models' forecasts) which primarily results from the reduction of increased resolution that is exhibited in the individual models' forecasts under various forecast probabilities.