Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radars”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

Comparing Satellite- and Ground-Based Observations of Cloud Occurrence Over High Southern Latitudes

The 2B-CLDCLASS-LIDAR R05 (2BCL5) and the raDAR/liDAR (DARDAR) satellite retrievals of cloud occurrence are compared as a function of altitude and latitude. The largest disparities are observed at low altitudes over high southern latitudes. These data sets are cross referenced to ground-based measurements from the Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE) campaign at McMurdo Station, Antarctica. Compared with AWARE observations, both 2BCL5 and DARDAR underestimate cloud occurrence below 1.5 km, with 2BCL5 and DARDAR distinguishing roughly one third of cloud occurrences observed by AWARE at 0.5 km. While DARDAR identifies greater cloud occurrences than 2BCL5 below 1.5 km, cloud occurrence values for the two data sets have similar differences relative to ground-based measurements. Therefore, the DARDAR retrievals of greater cloud occurrence at low altitudes are likely due to a larger quantity of false positives associated with radar ground clutter or attenuated lidar retrievals. DARDAR cloud occurrences match better with AWARE than 2BCL5 above 5 km. However, the likely underestimation of ground-based measurements at higher altitudes suggests DARDAR may underestimate high level cloud occurrence. Finally, both data sets indicate the presence of liquid containing clouds at temperatures within the homogeneous freezing regime, despite the fact that the ECMWF-AUX data set implemented in their processing clearly indicates temperatures below -38°C. Finally, using AWARE radiosonde (ECMWF-AUX) temperature data, we find that 2BCL5 detects 13.3% (13.8%) of mixed phase clouds below -38°C, while DARDAR detects 5.7% (6.6%) of mixed phase and 1.1% (1.3%) of liquid phase clouds below -38°C.

54 ENVIRONMENTAL SCIENCES↗

A Global High‐Resolution Mesoscale Convective System Database Using Satellite‐Derived Cloud Tops, Surface Precipitation, and Tracking

Abstract A new methodology is developed to construct a global (60°S–60°N) long‐term (2000–2019) high‐resolution (∼10‐km h) mesoscale convective system (MCS) database by tracking MCS jointly using geostationary satellite infrared brightness temperature ( T b ) and precipitation feature (PF) characteristics from the Integrated Multi‐satellitE Retrievals for GPM precipitation data sets. Independent validation shows that the satellite‐based MCS data set is able to reproduce important MCS statistics derived from ground‐based radar network observations in the United States and China. We show that by carefully considering key PF characteristics in addition to T b signatures, the new method significantly improves upon previous T b ‐only methods in detecting MCSs in the midlatitudes for all seasons. Results show that MCSs account for over 50% of annual total rainfall across most of the tropical belt and in selected regions of the midlatitudes, with a strong seasonality over many regions of the globe. The tracking database allows Lagrangian aspects such as MCS lifetime and translational speed and direction to be analyzed. The longest‐lived MCSs preferentially occur over the subtropical oceans. The land MCSs have higher cloud‐tops associated with more intense convection, and oceanic MCSs have much higher rainfall production. While MCSs are observed in many regions of the globe, there are fundamental differences in their dynamic and thermodynamic structures that warrant a better understanding of processes that control their evolution. This global database provides significant opportunities for observational and modeling studies of MCSs, their characteristics, and roles in regional and global water and energy cycles, as well as their hydrologic and other impacts.

54 ENVIRONMENTAL SCIENCES↗

Simulation of Continental Shallow Cumulus Populations Using an Observation-Constrained Cloud-System Resolving Model

Continental shallow cumulus (ShCu) clouds observed on 30 August 2016 during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign are simulated by using an observation-constrained cloud-system resolving model. On this day, ShCu forms over Oklahoma and southern Kansas and some of these clouds transition to deeper, precipitating convection during the afternoon. We apply a four-dimensional ensemble-variational (4DEnVar) hybrid technique in the Community Gridpoint Statistical Interpolation (GSI) system to assimilate operational data sets and unique boundary layer measurements including a Raman lidar, radar wind profilers, radiosondes, and surface stations collected by the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) atmospheric observatory into the Weather Research and Forecasting (WRF) model to ascertain how improved environmental conditions can influence forecasts of ShCu populations and the transition to deeper convection. Independent observations from aircraft, satellite, as well as ARM's remote sensors are used to evaluate model performance in different aspects. Several model experiments are conducted to identify the impact of data assimilation (DA) on the prediction of clouds evolution. The analyses indicate that ShCu populations are more accurately reproduced after DA in terms of cloud initiation time and cloud base height, which can be attributed to an improved representation of the ambient meteorological conditions and the convective boundary layer. Extending the assimilation to 18 UTC (local noon) also improved the simulation of shallow-to-deep transitions of convective clouds.

54 ENVIRONMENTAL SCIENCES↗

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate

Characterizing the properties of groundwater aquifers is essential for predicting aquifer response and managing groundwater resources. In this work, we develop a high-dimensional scalable Bayesian inversion framework governed by a three-dimensional quasi-static linear poroelastic model to characterize lateral permeability variations in groundwater aquifers. In this work, we determine the maximum a posteriori (MAP) point of the posterior permeability distribution from centimeter-level surface deformation measurements obtained from Interferometric Synthetic Aperture Radar (InSAR). The scalability of our method to high parameter dimension is achieved through the use of adjoint-based derivatives, inexact Newton methods to determine the MAP point, and a Mat´ern class sparse prior precision operator. Together, these guarantee that the MAP point is found at a cost, measured in number of forward/adjoint poroelasticity solves, that is independent of the parameter dimension. We apply our methodology to a test case for a municipal well in Mesquite, Nevada, in which InSAR and GPS surface deformation data are available. We solve problems with up to 320,824 state variable degrees of freedom (DOFs) and 16,896 parameter DOFs. A consistent treatment of noise level is employed so that the aquifer characterization result does not depend on the pixel spacing of surface deformation data. Our results show that the use of InSAR data significantly improves characterization of lateral aquifer heterogeneity, and the InSAR-based aquifer characterization recovers complex lateral displacement trends observed by independent daily GPS measurements.

54 ENVIRONMENTAL SCIENCES↗

What Determines the Number and the Timing of Pulses in Afternoon Precipitation in the Green Ocean Amazon (GoAmazon) Observations?

Using GoAmazon observations over the central Amazon, we investigate double-pulse events of afternoon precipitation for the first time. Relative humidity is the dominant factor in determining the onset timing and the number of pulses of precipitation. A moister free troposphere on single-pulse days features heavy precipitation in the early afternoon, which stabilizes the atmosphere and prohibits further convection. On early–onset double-pulse days, weaker precipitation from shallower convection in the first pulse leads to less consumption of atmospheric instability, making the environment favorable for a second pulse of stronger precipitation. A cluster tracking algorithm using the scanning radar reflectivity is developed to examine the relationship between precipitation pulses. Through explicitly tracking convective clusters that exist from the first to the second pulse, convective aggregation plays a more dominant role to the survival of clusters with larger final size, and on average around 20% of cluster survivals are due to natural growth.

54 ENVIRONMENTAL SCIENCES↗

Ensemble‐Based Assimilation of Satellite All‐Sky Microwave Radiances Improves Intensity and Rainfall Predictions for Hurricane Harvey (2017)

Abstract Ensemble‐based data assimilation of radar observations across inner‐core regions of tropical cyclones (TCs) in tandem with satellite all‐sky infrared (IR) radiances across the TC domain improves TC track and intensity forecasts. This study further investigates potential enhancements in TC track, intensity, and rainfall forecasts via assimilation of all‐sky microwave (MW) radiances using Hurricane Harvey (2017) as an example. Assimilating Global Precipitation Measurement constellation all‐sky MW radiances in addition to GOES‐16 all‐sky IR radiances reduces the forecast errors in the TC track, rapid intensification (RI), and peak intensity compared to assimilating all‐sky IR radiances alone, including a 24‐hr increase in forecast lead‐time for RI. Assimilating all‐sky MW radiances also improves Harvey's hydrometeor fields, which leads to improved forecasts of rainfall after Harvey's landfall. This study indicates that avenues exist for producing more accurate forecasts for TCs using available yet underutilized data, leading to better warnings of and preparedness for TC‐associated hazards in the future.

54 ENVIRONMENTAL SCIENCES↗

On the Importance of Using Event-Specific Wave Diffusion Rates in Modeling Diffuse Electron Precipitation

A few to tens of keV electron precipitation that carries substantial energy source down to the upper atmosphere to create aurora is manifested as an important magnetosphere-ionosphere coupling process. The precipitation is usually caused by scattering processes associated with plasma waves in the magnetosphere. The scattering process is often quantified by wave diffusion rates that indicate how fast an electron is scattered. Global models commonly use diffusion coefficients that are derived from statistical wave models. However, due to the statistical nature, many localized, transient features could be smeared out. In this study, we investigate electron precipitation using event-specific diffusion coefficients that are obtained based on simultaneous in-situ measured/inferred, rather than statistical, chorus wave dynamics. We find that the application of the event-specific diffusion coefficients associated with a more dynamic and intense chorus wave model leads more electrons, particularly at several to tens of keV in the dawn-to-noon sector at L > 3, to precipitate than using statistical coefficients. Here, the new simulation roughly captures both the intensity and variability of the precipitating flux as detected by the NOAA/POES satellites. Ionospheric electron density in the lower E region (100–120 km) observed by the mid-latitude Millstone Hill radar is also much better reproduced, while the case using statistical diffusion coefficients underestimates the ionization rate. This study implies the importance of using event-specific diffusion rates in simulating the diffuse electron precipitation and understanding the magnetosphere-ionosphere coupling.

79 ASTRONOMY AND ASTROPHYSICS↗

Adjoint-Based Inversion of Geodetic Data for Sources of Deformation and Strain

An adjoint-based formulation leads to a particularly efficient approach for inverting geodetic measurements for the source of the deformation. Specifically, the quantities necessary to iteratively improve the fit to the observations can be computed with just three forward calculations, one to obtain the current residuals, another to solve the adjoint problem, and a third to compute the step length. An inversion algorithm utilizing the adjoint-based gradient is applied to a set of Interferometric Synthetic Aperture Radar (InSAR) data gathered between 2016 and 2018 over the Tulare Basin in California's Central Valley. Because the measured deformation is due to groundwater withdrawal, a penalty function is included in the inversion to avoid placing aquifer volume change in locations that are far from any documented wells. The solution of the inverse problem provides estimates of aquifer compaction that provide a match to the observed range changes while honoring the well data. The solution indicates an average aquifer volume loss of 2.17 km 3 /year over the two year period from January 2016 to January 2018, encompassing one drought year (2016) and one wet year (2017). Finally, this magnitude of lost volume is compatible with the 3.1 km 3 /year decrease in water volume for the entire Central Valley, estimated from GRACE satellite gravity data.

58 GEOSCIENCES↗

Shallow Aseismic Slip in the Delaware Basin Determined by Sentinel-1 InSAR

The Delaware Basin, Texas is currently a hot-spot of induced seismicity and ground deformation due to fluid extraction and injection associated with horizontal drilling techniques; however, the driving mechanism behind the seismicity and deformation remains under debate. Here, using vertical and east-west horizontal surface deformation measurements derived from Sentinel-1 interferometric synthetic aperture radar (InSAR), we show that the subsurface responds differently to oil and gas activity in the northern and southeastern portions of the basin. In the north, where there is little seismicity, deformation patterns display long-wavelengths and equidimensional patterns. In contrast, the southeast region hosts most of the seismicity and displays spatial deformation patterns with narrow linear features that strike parallel to the maximum principal horizontal stress and to trends in seismicity, suggesting movement along normal faults. We model a linear deformation feature using edge dislocations and show that the InSAR observations can be reproduced by slip on normal faults contained within the Delaware Mountain Group (DMG), the formation that hosts local wastewater injection and the majority of earthquakes. Our model consists of three parallel, high-angle normal faults, with two dipping toward one another in a graben structure. Slip magnitudes reach up to 25 cm and are spatially correlated with injection wells. Measured seismicity can only explain ~2% of the fault motion predicted by our fault model, suggesting that slip leading to the deformation is predominantly aseismic. We conclude that seismic and aseismic fault motion in the southeastern Delaware Basin is likely driven by wastewater injection near critically-stressed normal faults within the DMG.

58 GEOSCIENCES↗

Evaluation of the CAM6 Climate Model Using Cloud Observations at McMurdo Station, Antarctica

A comparative analysis between observational data from McMurdo Station, Antarctica and the Community Atmosphere Model version 6 (CAM6) simulation is performed focusing on cloud characteristics and their thermodynamic conditions. Ka-band Zenith Radar (KAZR) and High Spectral Resolution Lidar (HSRL) retrievals are used as the basis of cloud fraction and cloud phase identifications. Radiosondes released at 12-hour increments provide atmospheric profiles for evaluating the simulated thermodynamic conditions. Our findings show that the CAM6 simulation consistently overestimates (underestimates) cloud fraction above (below) 3 km in four seasons of a year. Normalized by total in-cloud samples, ice and mixed phase occurrence frequencies are underestimated and liquid phase frequency is overestimated by the model at cloud fractions above 0.6, while at cloud fractions below 0.6 ice phase frequency is overestimated and liquid-containing phase frequency is underestimated by the model. The cloud fraction biases are closely associated with concurrent biases in relative humidity (RH), i.e., high (low) RH biases above (below) 2 km. Frequencies of correctly simulating ice and liquid-containing phase increase when the absolute biases of RH decrease. Cloud fraction biases also show a positive correlation with RH biases. Water vapor mixing ratio biases are the primary contributor to RH biases, and hence, likely a key factor controlling the cloud biases. Finally, this diagnosis of the evident shortfalls of representations of cloud characteristics in CAM6 simulation at McMurdo Station brings new insight in improving the governing model physics therein.

54 ENVIRONMENTAL SCIENCES↗

All-sky aerosol direct radiative effects at the ARM SGP site

All-sky aerosol direct radiative effect (DRE) was estimated for the first time at the Atmospheric Radiation Measurement Southern Great Plains site using multi-year ground-based observations. The NASA Langley Fu-Liou radiation model was employed. Observed inputs for the radiation model include aerosol and cloud vertical extinction profile from Raman lidar; spectral aerosol optical depth, single-scattering albedo, and asymmetry factor from Aerosol Robotic Network; cloud water content profiles from radars; temperature and water vapor profiles from radiosondes; and surface shortwave spectral albedo from radiometers. A cloudy-sky radiative closure experiment was performed. The relative mean differences between modeled and observed surface downwelling shortwave total fluxes were 6% (7%) for transparent (opaque) cloudy-skies. The estimated annual mean all-sky aerosol DRE is -2.13±0.54W m -2 at the top of atmosphere (TOA) and -5.95±0.87 W m -2 at the surface, compared to -3.00±0.58 W m -2 AND -6.85±1.00 W m -2 , respectively, under clear-sky conditions. The seasonal cycle of all-sky aerosol DRE is similar to that of the clear-sky, except with secondary influences of the clouds: The cloud radiative effect is strongest (most negative) in the spring, which reduces the all-sky aerosol DRE. The relative uncertainties in all-sky aerosol DRE due to measurement errors are generally comparable to those in clear-sky conditions except for the aerosol single-scattering albedo. Finally, the TOA all-sky aerosol DRE relative uncertainty due to aerosol single-scattering albedo uncertainty is larger than that in clear-sky, leading to a larger total relative uncertainty. The measurement errors in cloud properties have small effects on the all-sky aerosol DRE.

54 ENVIRONMENTAL SCIENCES↗

The Behaviors of Intraseasonal Cloud Organization During DYNAMO/AMIE

This study investigates the organization of tropical convection associated with the Madden-Julian Oscillation (MJO) during the Dynamics of the MJO/Atmospheric Radiation Measurement MJO Investigation Experiment field campaign. While it is known that tropical clouds can organize and impact the large-scale environment, how this occurs, and its underlying mechanism are not fully understood. Application of several existing cloud organization indices showed inconsistent evolutions in the measured degree of organization with the MJO. The inconsistency arises from the varying definitions and assumptions of cloud organization behaviors that are applied to each index. While these indices often combine different properties of clouds, such as their number density, size, and distance between them, the analysis of these properties separately provided further understanding of how clouds organize with the MJO. Using the rainfall clusters identified from the S-Polka radar, we find that deep convective rainfall clusters begin to increase their number density before the arrival of MJO enhanced convective center, which is accompanied by increased proximity (shorter distance to each other) and followed by growth in their size. However, the nonrandomness in the spatial distribution of rainfall clusters maximizes as MJO convection decays. Deep convective clusters become the least randomly distributed as the clusters decay because of the suppression and decay of isolated deep convective cells, while clustered deep convective cells exist longer. This evolution of cloud organization is analogous to mesoscale convective systems, indicating that the duration and frequency of their organization stages are altered by the large-scale environmental perturbations associated with the MJO.

54 ENVIRONMENTAL SCIENCES↗

Using TRMM Latent Heat as a Source to Estimate Convection Induced Gravity Wave Momentum Flux in the Lower Stratosphere

Due to a lack of observations, it is a big challenge to quantify the momentum flux transported from the troposphere into the lower stratosphere by gravity waves generated by convection. This limits our understanding of the dynamics of general circulation in the stratosphere and the realism of simulations of the Quasi Biennial Oscillation. In the past decade, some general circulation models have linked the momentum flux to the latent heating from subgrid scale convective precipitation with some success. However, there is still a large uncertainty in the sources of subgrid scale convection in these models. This study applies the links between the momentum flux and convective latent heating directly to the convective precipitation derived from 16 years of TRMM precipitation radar observations in the tropics and subtropics. The total and directional momentum flux at 100 hPa are derived at individual convective pixels by using the TRMM Spectrum Latent Heating product and large-scale wind profiles from the ERA-Interim dataset. For the first time, we are able to estimate the geographical distribution of momentum flux at 100 hPa in tropics and subtropics from the observed convective sources. The diurnal, seasonal, and interannual variations of the derived momentum flux are presented and discussed. These results could provide a reference to validate the dynamic coupling between the tropical troposphere and the stratosphere in the general circulation models.

54 ENVIRONMENTAL SCIENCES↗

Achievements and Prospects of Global Broadband Seismographic Networks After 30 Years of Continuous Geophysical Observations

Global seismographic networks (GSNs) emerged during the late nineteenth and early twentieth centuries, facilitated by seminal international developments in theory, technology, instrumentation, and data exchange. The mid- to late-twentieth century saw the creation of the World-Wide Standardized Seismographic Network (1961) and International Deployment of Accelerometers (1976), which advanced global geographic coverage as seismometer bandwidth increased greatly allowing for the recording of the Earth's principal seismic spectrum. The modern era of global observations and rapid data access began during the 1980s, and notably included the inception of the GEOSCOPE initiative (1982) and GSN (1988). Through continual improvements, GEOSCOPE and the GSN have realized near-real time recording of ground motion with state-of-art data quality, dynamic range, and timing precision to encompass 180 seismic stations, many in very remote locations. Data from GSNs are increasingly integrated with other geophysical data (e.g., space geodesy, infrasound and Interferometric Synthetic Aperture Radar). Globally distributed seismic data are critical to resolving crust, mantle, and core structure; illuminating features of the plate tectonic and mantle convection system; rapid characterization of earthquakes; identification of potential tsunamis; global nuclear test verification; and provide sensitive proxies for environmental changes. As the global geosciences community continues to advance our understanding of Earth structure and processes controlling elastic wave propagation, GSN infrastructure offers a springboard to realize increasingly multi-instrument geophysical observatories. Here, we review the historical, scientific, and monitoring heritage of GSNs, summarize key discoveries, and discuss future associated opportunities for Earth Science.

58 GEOSCIENCES↗

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. 2. Quantifying the Uncertainty

Uncertainty quantification of groundwater (GW) aquifer parameters is critical for efficient management and sustainable extraction of GW resources. These uncertainties are introduced by the data, model, and prior information on the parameters. We develop a Bayesian inversion framework that uses Interferometric Synthetic Aperture Radar (InSAR) surface deformation data to infer the laterally heterogeneous permeability of a transient linear poroelastic model of a confined GW aquifer. The Bayesian solution of this inverse problem takes the form of a posterior probability density of the permeability. Exploring this posterior using classical Markov chain Monte Carlo (MCMC) methods is computationally prohibitive due to the large dimension of the discretized permeability field and the expense of solving the poroelastic forward problem. However, in many partial differential equation (PDE)-based Bayesian inversion problems, the data are only informative in a few directions in parameter space. For the poroelasticity problem, we prove this property theoretically for a one-dimensional problem and demonstrate it numerically for a three-dimensional aquifer model. Here we design a generalized preconditioned Crank-Nicolson (gpCN) MCMC method that exploits this intrinsic low dimensionality by using a low-rank-based Laplace approximation of the posterior as a proposal, which we build scalably. The feasibility of our approach is demonstrated through a real GW aquifer test in Nevada. The inherently two-dimensional nature of InSAR surface deformation data informs a sufficient number of modes of the permeability field to allow detection of major structures within the aquifer, significantly reducing the uncertainty in the pressure and the displacement quantities of interest.

54 ENVIRONMENTAL SCIENCES↗

Thunderstorm Activity Under Intermediate and Extreme Climate Change Scenarios

Assessing how increasing greenhouse gas concentrations may modify regional climates is an ongoing challenge. Relatively little work has examined how climate change may influence processes related to regional thunderstorm activity. This is important to consider in areas where thunderstorms are important for maintaining regional hydroclimates. Here, this study focuses on three convection-allowing climate simulations—namely, a retrospective simulation (1990–2005) and two possible climate change scenarios (2085–2100)—with domains encompassing the conterminous United States. Regional and seasonal variability is noted in the response of thunderstorm activity as measured by thresholds in simulated radar reflectivity factor in the two climate change scenarios. A decrease in thunderstorm activity is projected for the Southern Plains, whereas the Southeast and Midwest experience an increase in thunderstorm activity. An examination of environmental parameters related to thunderstorm activity reveals an overall increase in convective instability but spatially varying changes in convective inhibition.

54 ENVIRONMENTAL SCIENCES↗

Understanding the Microphysical Control and Spatial‐Temporal Variability of Warm Rain Probability Using CloudSat and MODIS Observations

Abstract By combining measurements from MODIS and the CloudSat radar, we develop a parameterization scheme to quantify the combined microphysical controls by liquid water path (LWP) and cloud droplet number concentration (CDNC) of the probability of precipitation (PoP) in marine low cloud over tropical oceans. We demonstrate that the spatial‐temporal variation of grid‐mean in‐cloud can be largely explained by the variation of the joint probability density function of LWP and CDNC in the phase space specified by the bivariate PoP (LWP and CDNC) function. Through a series of sensitivity tests guided by this understanding, we find that in the Southeastern Pacific and Atlantic the stratocumulus to cumulus transition of the is mainly due to the variation of CDNC while the annual cycle is mainly due to the variation of LWP. The results of this study provide a viable way to diagnose the root cause of warm rain problems in global climate models.

54 ENVIRONMENTAL SCIENCES↗

Groundwater Volume Loss in Mexico City Constrained by InSAR and GRACE Observations and Mechanical Models

Abstract Groundwater withdrawal can cause localized and rapid poroelastic subsidence, spatially broad elastic uplift of low amplitude, and changes in the gravity field. Constraining groundwater loss in Mexico City, we analyze data from the Gravity Recovery and Climate Experiment and its follow‐on mission (GRACE/FO) and Synthetic Aperture Radar (SAR) Sentinel‐1A/B images between 2014 and 2021. GRACE/FO observations yield a groundwater loss of 0.85–3.87 km 3 /yr for a region of ∼300 × 600 km surrounding Mexico City. Using the high‐resolution interferometric SAR data set, we measure >35 cm/yr subsidence within the city and up to 2 cm/yr of uplift in nearby areas. Attributing the long‐term subsidence to poroelastic aquifer compaction and the long‐term uplift to elastic unloading, we apply respective models informed by local geology, yielding groundwater loss of 0.86–12.57 km 3 /yr. Our results suggest Mexico City aquifers have been depleting at faster rates since 2015, exacerbating the socioeconomic and health impacts of long‐term groundwater overdrafts.

58 GEOSCIENCES↗