Terrestrial groundwater travels through subterranean estuaries before reaching the sea. Groundwater‐derived nutrients drive coastal water quality, primary production, and eutrophication. We determined how dissolved inorganic nitrogen (DIN), dissolved inorganic phosphorus (DIP), and dissolved organic nitrogen (DON) are transformed within subterranean estuaries and estimated submarine groundwater discharge (SGD) nutrient loads compiling > 10,000 groundwater samples from 216 sites worldwide. Nutrients exhibited complex, nonconservative behavior in subterranean estuaries. Fresh groundwater DIN and DIP are usually produced, and DON is consumed during transport. Median total SGD (saline and fresh) fluxes globally were 5.4, 2.6, and 0.18 Tmol yr⁻¹ for DIN, DON, and DIP, respectively. Despite large natural variability, total SGD fluxes likely exceed global riverine nutrient export. Fresh SGD is a small source of new nutrients, but saline SGD is an important source of mostly recycled nutrients. Nutrients exported via SGD via subterranean estuaries are critical to coastal biogeochemistry and a significant nutrient source to the oceans.

Temporal gaps within the Gravity Recovery and Climate Experiment (GRACE) (gap: 20 months), between GRACE and GRACE Follow-On (GRACE-FO) missions (gap: 11 months), and within GRACE-FO record (gap: 2 months) make it difficult to analyze and interpret spatiotemporal variability in GRACE- and GRACE-FO-derived terrestrial water storage (TWS_GRACE) time series. In this study, an overview of data and approaches used to fill these gaps and reconstruct the TWS_GRACE record at the global scale is provided. In addition, the study provides an innovative approach that integrates three machine learning techniques (deep-learning neural networks [DNN], generalized linear model [GLM], and gradient boosting machine [GBM]) and eight climatic and hydrological input variables to fill these gaps and reconstruct the TWS_GRACE data record at both global grid and basin scales. For each basin and grid cell, the model performance was assessed using Nash–Sutcliffe efficiency coefficient (NSE), correlation coefficient (CC), and normalized root-mean-square error (NRMSE), a leader model was selected based on the model performance, and variables that significantly control leader model outputs were defined. Results indicate that (1) the leader model reconstructed the TWS_GRACE with high accuracy over both grid and local scales, particularly in wet and low anthropogenically active regions (grid scale: NSE = 0.65 ± 0.20, CC = 0.81 ± 0.13, and NSE = 0.56 ± 0.16; basin scale: NSE = 0.78 ± 0.14, CC = 0.89 ± 0.07, and NRMSE = 0.43 ± 0.14); (2) no single model was flawless in reconstructing the TWS_GRACE over all grids or basins, so a combination of models is necessary; (3) basin-scale models outperform grid-scale models; (4) the DNN model outperforms both GLM and GBM at the basin scale, whereas the GBM outperforms at the grid scale; (5) among other inputs, the Global Land Data Assimilation System (GLDAS)-derived TWS controls the model performance on both basin and grid scales; and (6) the reconstructed TWS_GRACE data captured extreme climatic events over the investigated basins and grid cells. The developed approach is robust, effective, and could be used to accurately reconstruct TWS_GRACE for any hydrologic system across the globe.

Land subsidence and sea level rise are well-known, ongoing problems that are negatively impacting the entire Texas coast. Although ground-based monitoring techniques using long-term global navigation satellite systems (GNSS) records provide accurate subsidence rates, they are labor intensive, expensive, time-consuming, and spatially limited. In this study, interferometric synthetic aperture radar (InSAR) data and techniques were used to map the locations and quantify rates of land subsidence in the Texas Coastal Bend region during the period from October 2016 to July 2019. InSAR-derived land subsidence rates were then validated and calibrated against GNSS-derived rates. The factors controlling the observed land subsidence rates and locations were investigated. The consequences of spatial variability in land subsidence rates in Coastal Bend were also examined. The results indicated that: (1) land subsidence rates in the Texas Coastal Bend exhibited spatial variability, (2) InSAR-derived land subsidence rates were consistent with GNSS-derived deformation rates, (3) land subsidence in the Texas Coastal Bend could be attributed mainly to hydrocarbon and groundwater extraction as well as vertical movements along growth faults, and (4) land subsidence increased both flood frequency and severity in the Texas Coastal Bend. Our results provide valuable information regarding not only land deformation rates in the Texas Coastal Bend region, but also the effectiveness of interferometric techniques for other coastal rural areas around the globe.