Search

Climate change is expected to increase waterborne diseases especially in developing countries. However, we lack understanding of how different types of water sources (both improved and unimproved) are affected by climate change, and thus, where to prioritize future investments and improvements to maximize health outcomes. This is due to limited knowledge of the relationships between source water quality and the observed variability in climate conditions. To address this gap, a 20‐month observational study was conducted in Tanzania, aiming to understand how water quality changes at various types of sources due to short‐term climate variability. Nine rounds of microbiological water quality sampling were conducted for Escherichia coli and total coliforms, at three study sites within different climatic regions. Each round included approximately 233 samples from water sources and 632 samples from households. To identify relationships between water quality and short‐term climate variability, Bayesian hierarchical modeling was adopted, allowing these relationships to vary with source types and sampling regions to account for potentially different physical processes. Across water sources, increases in E. coli/total coliform levels were most closely related to increases in recent heavy rainfall. Our key recommendations to future longitudinal studies are (a) demonstrated value of high sampling frequency and temporal coverage (a minimum of 3 years) especially during wet seasons; (b) utility of the Bayesian hierarchical models to pool data from multiple sites while allowing for variations across space and water sources; and (c) importance of a multidisciplinary team approach with consistent commitment and sharing of knowledge., Key Points We present a longitudinal study in a developing country on water quality changes at a range of water sources under climate variabilityIncreases in E. coli and Total Coliform levels were most closely related to recent heavy rainfallRecommendations for future cross disciplinary studies on drinking water quality and relationships with climate variability are made

Three different methods of specifying ozone in an atmosphere‐only version of the HadGEM3‐A global circulation model are compared to the coupled chemistry configuration of this model. These methods include a specified zonal‐mean ozone climatology, a specified 3‐D ozone climatology, and a calculated‐asymmetry scheme in which a specified zonal‐mean ozone field is adapted online to be consistent with dynamically produced zonal asymmetries. These simulations all use identical boundary conditions and, by construction, have the same climatological zonal‐mean ozone, that of the coupled chemistry configuration of the model. Prescribing ozone, regardless of scheme, results in a simulation which is 3–4 times faster than the coupled chemistry‐climate model (CCM). Prescribing climatological zonal asymmetries leads to a vortex which is the correct intensity but which is systematically displaced over regions with lower prescribed ozone. When zonal asymmetries in ozone are free to evolve interactively with model dynamics, the modeled wintertime stratospheric vortex shape and mean sea level pressure patterns closely resemble that produced by the full CCM in both hemispheres, in terms of statistically significant differences. Further, we separate out the two distinct pathways by which zonal ozone asymmetries influence modeled dynamics. We present this interactive‐ozone zonal‐asymmetry scheme as an inexpensive tool for accurately modeling the impacts of dynamically consistent ozone fields as seen in a CCM which ultimately influence mean sea level pressure and tropospheric circulation (particularly during wintertime in the Northern Hemisphere, when ozone asymmetries are generally largest), without the computational burden of simulating interactive chemistry., Key Points Dynamical feedbacks resulting from ozone zonal asymmetries are important for accurately representing Northern Hemisphere circulation patterns during wintertimeThree simulations using different methods of specifying ozone mixing ratios are compared with a fully coupled CCM simulation in order to examine the pathways by which zonal ozone asymmetries influence modeled dynamicsDynamical feedbacks occurring in a coupled CCM can be reproduced, at a greatly reduced computational cost, by appropriately placing a prescribed ozone climatology onto a dynamically evolving coordinate

Vb cyclones typically emerge in the Western Mediterranean and propagate to the Northeast into Central Europe. This paper explores the temporal characteristics of Vb cyclone occurrence based on cyclone tracks identified at the atmospheric levels of Z700 and sea level pressure, using JRA‐55 reanalysis data for the period 1959–2015. The risk of Vb occurrence was significantly high in the 1960s and has remained at a lower level since then. Vb cyclones do not occur fully randomly according to a Poisson point process. Eleven well‐separated and distinct clusters as well as 11 hiatus periods are identified, with average occurrence rates of 21.5 and 5.2 yrea−1, respectively. During the event of Vb, the large‐scale atmospheric circulation is changed into a state favoring the development of successive Vb cyclones. Clustering is very prominent in the case of Genoan Vb cyclones in summer as well as those Vb cyclones developing over the Iberian Peninsula or the North African Coast in winter. Superposition of the polar and the subtropical jet stream over the Western Mediterranean is identified as a main feature at the onset of Vb cyclones. Vb cyclone occurrence appears to be synchronized with the Northern Atlantic Oscillation (NAO; at Z500) and Arctic Oscillation (AO; at Z1000). Clusters have occurred when both NAO and AO were negative. This relation applies to Western Mediterranean cyclones not following a Vb track as well, however to a much weaker extent. In contrast, Vb cyclone frequency was particularly low from 1988 to 1997 during a sustained positive phase of both NAO and AO., Key Points In the past 50 years, the frequency of Vb cyclones was high in the 1960s and lower since thenVb cyclones were strongly synchronized with NAO and AO; high frequency clusters occurred when NAO and AO were negativeThe coupling of the polar and the subtropical jet stream over the Western Mediterranean is a driving mechanism of the onset of Vb cyclones

The NASA Cyclone Global Navigation Satellite System (CYGNSS) was launched in December 2016, providing an unprecedented opportunity to obtain ocean surface wind speeds including wind estimates over the hurricane inner‐core region. This study demonstrates the influence of assimilating an early version of CYGNSS observations of ocean surface wind speeds on numerical simulations of two notable landfalling hurricanes, Harvey and Irma (2017). A research version of the National Centers for Environmental Prediction operational Hurricane Weather Research and Forecasting model and the Gridpoint Statistical Interpolation‐based hybrid ensemble three‐dimensional variational data assimilation system are used. It is found that the assimilation of CYGNSS data results in improved track, intensity, and structure forecasts for both hurricane cases, especially for the weak phase of a hurricane, implying potential benefits of using such data for future research and operational applications., Key Points The NASA Cyclone Global Navigation Satellite System (CYGNSS) provides an unprecedented opportunity to obtain ocean surface wind data over a hurricane inner‐core regionThis study found that the assimilation of CYGNSS data results in improved track, intensity, and structure forecasts for two notable landfalling hurricanes, Harvey and Irma (2017)

The ability of FLake, WRF‐Lake, and CoLM‐Lake models in simulating the thermal features of Lake Nam Co in Central Tibetan Plateau has been evaluated in this study. All the three models with default settings exhibited distinct errors in the simulated vertical temperature profile. Then model calibration was conducted by adjusting three (four) key parameters within FLake and CoLM‐Lake (WRF‐Lake) in a series of sensitive experiments. Results showed that each model's performance is sensitive to the key parameters and becomes much better when adjusting all the key parameters relative to tuning single parameter. Overall, setting the temperature of maximum water density to 1.1 °C instead of 4 °C in the three models consistently leads to improved vertical thermal structure simulation during cold seasons; reducing the light extinction coefficient in FLake results in much deeper mixed layer and warmer thermocline during warm seasons in better agreement with the observation. The vertical thermal structure can be clearly improved by decreasing the light extinction coefficient and increasing the turbulent mixing in WRF‐Lake and CoLM‐Lake during warm seasons. Meanwhile, the modeled water temperature profile in warm seasons can be significantly improved by further replacing the constant surface roughness lengths by a parameterized scheme in WRF‐Lake. Further intercomparison indicates that among the three calibrated models, FLake (WRF‐Lake) performs the best to simulate the temporal evolution and intensity of temperature in the layers shallower (deeper) than 10 m, while WRF‐Lake is the best at simulating the amplitude and pattern of the temperature variability at all depths., Key Points The performance of three one‐dimensional lake models in simulating the thermal structure of Nam Co Lake is evaluated and improvedKey processes related to the simulated thermal regime of alpine lakes on the Tibetan Plateau are indicated and revealed

The development of a unified similarity scaling has so far failed over complex surfaces, as scaling studies show large deviations from the empirical formulations developed over flat and horizontally homogeneous terrain as well as large deviations between the different complex terrain data sets. However, a recent study of turbulence anisotropy for flat and horizontally homogeneous terrain has shown that separating the data according to the limiting states of anisotropy (isotropic, two‐component axisymmetric and one‐component turbulence) improves near‐surface scaling. In this paper we explore whether this finding can be extended to turbulence over inclined and horizontally heterogeneous surfaces by examining near‐surface scaling for 12 different data sets obtained over terrain ranging from flat to mountainous. Although these data sets show large deviations in scaling when all anisotropy types are examined together, the separation according to the limiting states of anisotropy significantly improves the collapse of data onto common scaling relations, indicating the possibility of a unified framework for turbulence scaling. A measure of turbulence complexity is developed, and the causes for the breakdown of scaling and the physical mechanisms behind the turbulence complexity encountered over complex terrain are identified and shown to be related to the distance to the isotropic state, prevalence of directional shear with height in mountainous terrain, and the deviations from isotropy in the inertial subrange., Key Points A pathway toward a unified theory of near‐surface atmospheric turbulence is providedThe separation according to the limiting states of anisotropy significantly improves the collapse of data onto common scaling relationsPhysical processes responsible for turbulence complexity are identified

Terrestrial gamma ray flashes (TGFs) are sub‐millisecond bursts of high energetic gamma radiation associated with intracloud flashes in thunderstorms. In this paper we use the simultaneity of lightning detections by World Wide Lightning Location Network to find TGFs in the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) data that are too faint to be identified by standard search algorithms. A similar approach has been used in an earlier paper, but here we expand the data set to include all years of RHESSI + World Wide Lightning Location Network data and show that there is a population of observationally weak TGFs all the way down to 0.22 of the RHESSI detection threshold (three counts in the detector). One should note that the majority of these are “normal” TGFs that are produced further away from the subsatellite point (and experience a 1/r 2 effect) or produced at higher latitudes with a lower tropoause and thus experience increased atmospheric attenuation. This supports the idea that the TGF production rate is higher than currently reported. We also show that compared to lightning flashes, TGFs are more partial to ocean and coastal regions than over land., Key Points We show that there is a fluence distribution of TGFs down to 0.22 of the detection threshold, and we do not see a cutoff in the RHESSI TGF distributionWe find that there are six times as many TGFs inside of RHESSI's field of view, than can be identified by current search algorithmsWe find that observationally weak TGFs largely originate at larger radial distances and higher latitudes from the subsatellite point

Global multiconstituent concentration and emission fields obtained from the assimilation of the satellite retrievals of ozone, CO, NO2, HNO3, and SO2 from the Ozone Monitoring Instrument (OMI), Global Ozone Monitoring Experiment 2, Measurements of Pollution in the Troposphere, Microwave Limb Sounder, and Atmospheric Infrared Sounder (AIRS)/OMI are used to understand the processes controlling air pollution during the Korea‐United States Air Quality (KORUS‐AQ) campaign. Estimated emissions in South Korea were 0.42 Tg N for NOx and 1.1 Tg CO for CO, which were 40% and 83% higher, respectively, than the a priori bottom‐up inventories, and increased mean ozone concentration by up to 7.5 ± 1.6 ppbv. The observed boundary layer ozone exceeded 90 ppbv over Seoul under stagnant phases, whereas it was approximately 60 ppbv during dynamical conditions given equivalent emissions. Chemical reanalysis showed that mean ozone concentration was persistently higher over Seoul (75.10 ± 7.6 ppbv) than the broader KORUS‐AQ domain (70.5 ± 9.2 ppbv) at 700 hPa. Large bias reductions (>75%) in the free tropospheric OH show that multiple‐species assimilation is critical for balanced tropospheric chemistry analysis and emissions. The assimilation performance was dependent on the particular phase. While the evaluation of data assimilation fields shows an improved agreement with aircraft measurements in ozone (to less than 5 ppbv biases), CO, NO2, SO2, PAN, and OH profiles, lower tropospheric ozone analysis error was largest at stagnant conditions, whereas the model errors were mostly removed by data assimilation under dynamic weather conditions. Assimilation of new AIRS/OMI ozone profiles allowed for additional error reductions, especially under dynamic weather conditions. Our results show the important balance of dynamics and emissions both on pollution and the chemical assimilation system performance., Key Points Multiconstituent data assimilation during KORUS‐AQ showed that emissions in South Korea were 0.42 Tg N for NOx and 1.1 Tg CO for COThese emissions were 40% and 83% higher, respectively, than the a priori bottom‐up inventories and increased ozone by up to 7.5 ± 1.6 ppbvMean ozone concentration was persistently higher over Seoul (75.1 ± 7.6 ppbv) than the broader KORUS‐AQ domain (70.5 ± 9.2 ppbv) at 700 hPa

Extreme temperature and pressure conditions on the surface of Venus present formidable technological challenges against performing ground‐based seismology. Efficient coupling between the Venusian atmosphere and the solid planet theoretically allows the study of seismically generated acoustic waves using balloons in the upper atmosphere, where conditions are far more clement. However, earthquake detection from a balloon has never been demonstrated. We present the first detection of an earthquake from a balloon‐borne microbarometer near Ridgecrest, CA in July 2019 and include a detailed analysis of the dependence of seismic infrasound, as measured from a balloon on earthquake source parameters, topography, and crustal and atmospheric structure. Our comprehensive analysis of seismo‐acoustic phenomenology demonstrates that seismic activity is detectable from a high‐altitude platform on Earth, and that Rayleigh wave‐induced infrasound can be used to constrain subsurface velocities, paving the way for the detection and characterization of such signals on Venus., Key Points First detection of a natural earthquake using balloon‐borne infrasound dataRayleigh wave‐induced infrasound dispersion characteristics provide constraints on subsurface velocitiesShallow waveguides, focal mechanism, and subwavelength topographic changes control infrasound amplitude and dispersion by weak earthquakes

The Tien Shan and Pamir mountains host over 28,000 glaciers providing essential water resources for increasing water demand in Central Asia. A disequilibrium between glaciers and climate affects meltwater release to Central Asian rivers, challenging the region's water availability. Previous research has neglected temporal variability. We present glacier mass balance estimates based on transient snowline and geodetic surveys with unprecedented spatiotemporal resolution from 1999/00 to 2017/18. Our results reveal spatiotemporal heterogeneity characterized by two mass balance clusters: (a) positive, low variability, and (b) negative, high variability. This translates into variable glacial meltwater release (≈1–16%) of annual river runoff for two watersheds. Our study reveals more complex climate forcing‐runoff responses and importance of glacial meltwater variability for the region than suggested previously., Key Points Annual glacier mass balance for Central Asia (1999/00–2017/18) is derived by combining transient snowlines, geodetic surveys, and modelingStrong spatiotemporal heterogeneity with contrasting patterns of mass gain and loss are foundHot spots of heterogeneous mass balance variability are associated with highly variable glacier melt water runoff

24 pages, 8 figures, supporting information https://doi.org/10.1029/2021JD034811.-- Data Availability Statement: The data set used in this study are available in (1) Antoine et al. (2019) available at https://doi.org/10.5281/zenodo.3406983; (2) Chen et al. (2019) available at https://doi.org/10.5281/zenodo.3559982; (3) Landwehr et al. (2020) available at https://doi.org/10.5281/zenodo.3836439; (4) Landwehr, Thurnherr et al. (2020) available at https://doi.org/10.5194/amt-13-3487-2020; (5) Schmale et al. (2019) available at https://doi.org/10.5281/zenodo.2636709; (6) Tatzelt et al. (2020) available at https://doi.org/10.5281/zenodo.3922147; (7) Thomalla et al., (2020) available at https://doi.org/10.5281/zenodo.3859515; (8) Thurnherr et al. (2020) available at https://doi.org/10.5281/zenodo.4031705; In addition, the fluorescent aerosol and gel–like POM measurements have been uploading as supplementary information supporting this article, In this study, we investigate the occurrence of primary biological aerosol particles (PBAP) over all sectors of the Southern Ocean (SO) based on a 90-day data set collected during the Antarctic Circumnavigation Expedition (ACE) in austral summer 2016–2017. Super-micrometer PBAP (1–16 µm diameter) were measured by a wide band integrated bioaerosol sensor (WIBS-4). Low (3σ) and high (9σ) fluorescence thresholds are used to obtain statistics on fluorescent and hyper-fluorescent PBAP, respectively. Our focus is on data obtained over the pristine ocean, that is, more than 200 km away from land. The results indicate that (hyper-)fluorescent PBAP are correlated to atmospheric variables associated with sea spray aerosol (SSA) particles (wind speed, total super-micrometer aerosol number concentration, chloride and sodium concentrations). This suggests that a main source of PBAP over the SO is SSA. The median percentage contribution of fluorescent and hyper-fluorescent PBAP to super-micrometer SSA was 1.6% and 0.13%, respectively. We demonstrate that the fraction of (hyper-)fluorescent PBAP to total super-micrometer particles positively correlates with concentrations of bacteria and several taxa of pythoplankton measured in seawater, indicating that marine biota concentrations modulate the PBAP source flux. We investigate the fluorescent properties of (hyper-)fluorescent PBAP for several events that occurred near land masses. We find that the fluorescence signal characteristics of particles near land is much more variable than over the pristine ocean. We conclude that the source and concentration of fluorescent PBAP over the open ocean is similar across all sampled sectors of the SO, The authors acknowledge funding for ACE-SPACE, ACE-SORPASSO by the Swiss Polar Institute, and Ferring Pharmaceuticals. ACE-SPACE was carried out with additional support from the European FP7 project BACCHUS (Grant Agreement 49603445). Gang Chen received funding from the cost action of Chemical On-Line cOmpoSition and Source Apportionment of fine aerosol (COLOSSAL, CA16109), a COST related project of the Swiss National Science Foundation, Source apportionment using long-term Aerosol Mass Spectrometry and Aethalometer Measurements (SAMSAM, IZCOZ0_177063). Sebastian Landwehr received funding from the Swiss Data Science Center project c17-02. Andrea Baccarin received funding from the Swiss National Science Foundation (Grant 200021_169090). Julia Schmale holds the Ingvar Kamprad Chair. Martin Gysel-Beer received funding from the ERC under Grant ERC-CoG-615922-BLACARAT. Charlotte Robinson received funding from the Australian Research Council (Grand ARC DP160103387), With the funding support of the ‘Severo Ochoa Centre of Excellence’ accreditation (CEX2019-000928-S), of the Spanish Research Agency (AEI)

Responses to COVID‐19 have resulted in unintended reductions of city‐scale carbon dioxide (CO2) emissions. Here, we detect and estimate decreases in CO2 emissions in Los Angeles and Washington DC/Baltimore during March and April 2020. We present three lines of evidence using methods that have increasing model dependency, including an inverse model to estimate relative emissions changes in 2020 compared to 2018 and 2019. The March decrease (25%) in Washington DC/Baltimore is largely supported by a drop in natural gas consumption associated with a warm spring whereas the decrease in April (33%) correlates with changes in gasoline fuel sales. In contrast, only a fraction of the March (17%) and April (34%) reduction in Los Angeles is explained by traffic declines. Methods and measurements used herein highlight the advantages of atmospheric CO2 observations for providing timely insights into rapidly changing emissions patterns that can empower cities to course‐correct CO2 reduction activities efficiently., Key Points Atmospheric CO2 observations can be used to detect the onset of the COVID‐19 response in Los Angeles and Washington DC/BaltimoreRelative reductions in April 2020 associated with COVID‐19 are ∼30% when compared to emissions in 2018 and 2019Decreases in vehicular traffic do not completely explain observed emissions reductions in both Los Angeles and Washington DC/Baltimore

Responding to the 2020 COVID‐19 outbreak, China imposed an unprecedented lockdown producing reductions in air pollutant emissions. However, the lockdown driven air pollution changes have not been fully quantified. We applied machine learning to quantify the effects of meteorology on surface air quality data in 31 major Chinese cities. The meteorologically normalized NO2, O3, and PM2.5 concentrations changed by −29.5%, +31.2%, and −7.0%, respectively, after the lockdown began. However, part of this effect was also associated with emission changes due to the Chinese Spring Festival, which led to ∼14.1% decrease in NO2, ∼6.6% increase in O3 and a mixed effect on PM2.5 in the studied cities that largely resulted from festival associated fireworks. After decoupling the weather and Spring Festival effects, changes in air quality attributable to the lockdown were much smaller: −15.4%, +24.6%, and −9.7% for NO2, O3, and PM2.5, respectively., Key Points The Chinese Spring Festival led to a 14.1% decrease in NO2, 6.6% increase in O3 and a mixed effect on PM2.5 across China in 2015–2019The 2020 lockdown resulted in −15.4%, −17.0%, −14.5%, −7.6%, −9.7%, and +24.6% changes for NO2, SO2, CO, PM10, PM2.5, and O3, respectivelyFireworks emissions contributed to haze in many cities such as an additional 29 μg m−3 PM2.5 in Beijing on the first lunar day of 2020

The Mesozoic era (∼252 to 66 million years ago) was a key interval in Earth's evolution toward its modern state, witnessing the breakup of the supercontinent Pangaea and significant biotic innovations like the early evolution of mammals. Plate tectonic dynamics drove a fundamental climatic transition from the early Mesozoic supercontinent toward the Late Cretaceous fragmented continental configuration. Here, key aspects of Mesozoic long‐term environmental changes are assessed in a climate model ensemble framework. We analyze so far the most extended ensemble of equilibrium climate states simulated for evolving Mesozoic boundary conditions covering the period from 255 to 60 Ma in 5 Myr timesteps. Global mean temperatures are generally found to be elevated above the present and exhibit a baseline warming trend driven by rising sea levels and increasing solar luminosity. Warm (Triassic and mid‐Cretaceous) and cool (Jurassic and end‐Cretaceous) anomalies result from pCO2 changes indicated by different reconstructions. Seasonal and zonal temperature contrasts as well as continental aridity show an overall decrease from the Late Triassic‐Early Jurassic to the Late Cretaceous. Meridional temperature gradients are reduced at higher global temperatures and less land area in the high latitudes. With systematic sensitivity experiments, the influence of paleogeography, sea level, vegetation patterns, pCO2, solar luminosity, and orbital configuration on these trends is investigated. For example, long‐term seasonality trends are driven by paleogeography, but orbital cycles could have had similar‐scale effects on shorter timescales. Global mean temperatures, continental humidity, and meridional temperature gradients are, however, also strongly affected by pCO2., Key Points We assess global long‐term climate trends through the Mesozoic era with an ensemble of climate model simulationsVarying carbon dioxide levels cause anomalies around an overall warming trend due to changing paleogeography and increasing insolationSeasonal and zonal temperature contrasts as well as aridity decrease with time, while meridional gradients vary with paleogeography

Understanding the complex interrelationships between wildfire and its environmental and anthropogenic controls is crucial for wildfire modeling and management. Although machine learning (ML) models have yielded significant improvements in wildfire predictions, their limited interpretability has been an obstacle for their use in advancing understanding of wildfires. This study builds an ML model incorporating predictors of local meteorology, land‐surface characteristics, and socioeconomic variables to predict monthly burned area at grid cells of 0.25° × 0.25° resolution over the contiguous United States. Besides these predictors, we construct and include predictors representing the large‐scale circulation patterns conducive to wildfires, which largely improves the temporal correlations in several regions by 14%–44%. The Shapley additive explanation is introduced to quantify the contributions of the predictors to burned area. Results show a key role of longitude and latitude in delineating fire regimes with different temporal patterns of burned area. The model captures the physical relationship between burned area and vapor pressure deficit, relative humidity (RH), and energy release component (ERC), in agreement with the prior findings. Aggregating the contribution of predictor variables of all the grids by region, analyses show that ERC is the major contributor accounting for 14%–27% to large burned areas in the western US. In contrast, there is no leading factor contributing to large burned areas in the eastern US, although large‐scale circulation patterns featuring less active upper‐level ridge‐trough and low RH two months earlier in winter contribute relatively more to large burned areas in spring in the southeastern US., Key Points US fire burned areas are well predicted by a machine learning model and SHAP improves interpretation of the contributing factorsIncluding large‐scale circulation patterns conducive to wildfires as predictors improve prediction of burned areas in several regionsFire‐season fuel dryness and dry winters are important contributors to large burned areas in western and southeastern US, respectively

The Madden‐Julian Oscillation (MJO) is the leading mode of intraseasonal climate variability, having profound impacts on a wide range of weather and climate phenomena. Here, we use a wavelet‐based spectral Principal Component Analysis (wsPCA) to evaluate the skill of 20 state‐of‐the‐art CMIP6 models in capturing the magnitude and dynamics of the MJO. By construction, wsPCA has the ability to focus on desired frequencies and capture each propagative physical mode with one principal component (PC). We show that the MJO contribution to the total intraseasonal climate variability is substantially underestimated in most CMIP6 models. The joint distribution of the modulus and angular frequency of the wavelet PC series associated with MJO is used to rank models relatively to the observations through the Wasserstein distance. Using Hovmöller phase‐longitude diagrams, we also show that precipitation variability associated with MJO is underestimated in most CMIP6 models for the Amazonia, Southwest Africa, and Maritime Continent., Key Points A wavelet‐based spectral principal component analysis is used to examine CMIP6 models in reproducing the Madden‐Julian Oscillation (MJO)CMIP6 models capture the average MJO propagation speed but significantly underestimate the MJO contribution to the total intraseasonal climate variabilityPrecipitation variability related to MJO over the Amazonia, Southwest Africa, and Maritime Continent is underestimated in CMIP6 models

A new automated method to retrieve charge layer polarity from flashes, named Chargepol, is presented in this paper. Using data from the NASA Lightning Mapping Array (LMA) deployed during the Remote sensing of Electrification, Lightning, And Mesoscale/microscale Processes with Adaptive Ground Observations (RELAMPAGO) field campaign in Cordoba, Argentina, from November 2018 to April 2019, this method estimates the polarity of vertical charge distributions and their altitudes and thicknesses (or vertical depth) using the very‐high frequency (VHF) source emissions detected by LMAs. When this method is applied to LMA data for extended periods of time, it is capable of inferring a storm's bulk electrical charge structure throughout its life cycle. This method reliably predicted the polarity of charge within which lightning flashes propagated and was validated in comparison to methods that require manual assignment of polarities via visual inspection of VHF lightning sources. Examples of normal and anomalous charge structures retrieved using Chargepol for storms in Central Argentina during RELAMPAGO are presented for the first time. Application of Chargepol to five months of LMA data in Central Argentina and several locations in the United States allowed for the characterization of the charge structure in these regions and for a reliable comparison using the same methodology. About 13.3% of Cordoba thunderstorms were defined by an anomalous charge structure, slightly higher than in Oklahoma (12.5%) and West Texas (11.1%), higher than Alabama (7.3%), and considerably lower than in Colorado (82.6%). Some of the Cordoba anomalous thunderstorms presented enhanced low‐level positive charge, a feature rarely if ever observed in Colorado thunderstorms., Key Points A new automated method to estimate the polarity, altitude, and vertical depth of charge layers from flashes is presentedThunderstorm charge structures in Central Argentina are characterized from months of Lightning Mapping Array (LMA) data, with 13.3% anomalous charge structureCordoba anomalously charged thunderstorms have a distinct charge altitude distribution when compared to those in Colorado

Coastal boulder deposits (CBD) provide what are sometimes the only remaining signatures of wave inundation on rocky coastlines; in recent decades, CBD combined with initiation of motion (IoM) analyses have repeatedly been used as primary evidence to infer the existence of ancient tsunamis. However, IoM storm wave heights inferred by these studies have been shown to be highly inaccurate, bringing some inferences into question. This work develops a dimensionless framework to relate CBD properties with storm‐wave hindcasts and measurements, producing data‐driven relations between wave climate and boulder properties. We present an elevation‐density‐size‐inland distance‐wave height analysis for individual storm‐transported boulders which delineates the dynamic space where storm‐wave CBD occur. Testing these new relations against presumed tsunami CBD demonstrates that some fall well within the capabilities of storm events, suggesting that some previous studies might be fruitfully reexamined within the context of this new framework., Key Points Size, elevation, and distance inland in storm‐wave‐generated coastal boulder deposits are all functions of climatological wave heightsClear data ranges exist where storm wave coastal boulder deposits do/do not existSome coastal boulder deposits previously identified as having tsunami origin may have been generated by storm waves

As of July 27, 2020, COVID‐19 has caused 640,000 deaths worldwide and has had a major impact on people's productivity and lives. Analyzing the spatial distribution characteristics of COVID‐19 cases and their relationships with meteorological and environmental factors might help enrich our knowledge of virus transmission and formulate reasonable epidemic prevention strategies. Taking the cumulative confirmed cases in Hubei province from January 23, 2020, to April 8, 2020, as an example, this study analyzed the spatial evolution characteristics of confirmed COVID‐19 cases in Hubei province using exploratory spatial data analysis and explored the spatial relationship between the main environmental and meteorological factors and confirmed COVID‐19 cases using a geographically weighted regression (GWR) model. Results show that there was no obvious spatial clustering of confirmed COVID‐19 cases in Hubei province, while the decline and end of the newly confirmed cases revealed relatively obvious negative spatial correlations. Due to the lockdown in Hubei province, the main air quality indexes (e.g., AQI and PM2.5) decreased significantly and environmental quality was better than historical contemporaneous levels. Meanwhile, the results of the GWR model suggest that the impacts of environmental and meteorological factors on the development of COVID‐19 were not significant. These findings indicate that measures such as social distancing and isolation played the primary role in controlling the development of the COVID‐19 epidemic., Key Points The spatiotemporal patterns of COVID‐19 cases in Hubei Province were analyzedAnomalies of environmental and meteorological factors contributed little to COVID‐19 developmentSocial distancing and isolation played a key role in controlling the epidemic

The energetic particle precipitation (EPP) indirect effect (IE) refers to the downward transport of reactive odd nitrogen (NOx = NO + NO2) produced by EPP (EPP‐NOx) from the polar winter mesosphere and lower thermosphere to the stratosphere where it can destroy ozone. Previous studies of the EPP IE examined NOx descent averaged over the polar region, but the work presented here considers longitudinal variations. We report that the January 2009 split Arctic vortex in the stratosphere left an imprint on the distribution of NO near the mesopause, and that the magnitude of EPP‐NOx descent in the upper mesosphere depends strongly on the planetary wave (PW) phase. We focus on an 11‐day case study in late January immediately following the 2009 sudden stratospheric warming during which regional‐scale Lagrangian coherent structures (LCSs) formed atop the strengthening mesospheric vortex. The LCSs emerged over the north Atlantic in the vicinity of the trough of a 10‐day westward traveling planetary wave. Over the next week, the LCSs acted to confine NO‐rich air to polar latitudes, effectively prolonging its lifetime as it descended into the top of the polar vortex. Both a whole atmosphere data assimilation model and satellite observations show that the PW trough remained coincident in space and time with the NO‐rich air as both migrated westward over the Canadian Arctic. Estimates of descent rates indicate five times stronger descent inside the PW trough compared to other longitudes. This case serves to set the stage for future climatological analysis of NO transport via LCSs., Key Points First demonstration of the impact of the split Arctic vortex on the geographic distribution of nitric oxide at the winter mesopauseFirst evidence that a Lagrangian coherent structure inhibits horizontal transport of nitric oxide at the polar winter mesopauseDescent of nitric oxide is five times stronger between 80 and 90 km in a westward traveling planetary wave trough compared to the ridge

The ongoing coronavirus disease 2019 (COVID‐19) pandemic has caused more than 150 million cases of infection to date and poses a serious threat to global public health. In this study, global COVID‐19 data were used to examine the dynamical variations from the perspectives of immunity and contact of 84 countries across the five climate regions: tropical, arid, temperate, and cold. A new approach named Yi Hua Jie Mu is proposed to obtain the transmission rates based on the COVID‐19 data between the countries with the same climate region over the Northern Hemisphere and Southern Hemisphere. Our results suggest that the COVID‐19 pandemic will persist over a long period of time or enter into regular circulation in multiple periods of 1–2 years. Moreover, based on the simulated results by the COVID‐19 data, it is found that the temperate and cold climate regions have higher infection rates than the tropical and arid climate regions, which indicates that climate may modulate the transmission of COVID‐19. The role of the climate on the COVID‐19 variations should be concluded with more data and more cautions. The non‐pharmaceutical interventions still play the key role in controlling and prevention this global pandemic., Key Points A new approached is proposed to predict the future COVID‐19 variations rather than relying on information on other corona virusesCOVID‐19 pandemic will persist in multiple periods of 1–2 yearsThe temperate and cold climate regions have higher infection rates than the tropical and arid climate regions

Earth observations offer potential pathways for accurately closing the water and energy balance of watersheds, a fundamental challenge in hydrology. However, previous attempts based on purely satellite‐based estimates have focused on closing the water and energy balances separately. They are hindered by the lack of estimates of key components, such as runoff. Here, we posit a novel approach based on Budyko’s water and energy balance constraints. The approach is applied to quantify the degree of long‐term closure at the watershed scale, as well as its associated uncertainties, using an ensemble of global satellite data sets. We find large spatial variability across aridity, elevation, and other environmental gradients. Specifically, we find a positive correlation between elevation and closure uncertainty, as derived from the Budyko approach. In mountainous watersheds the uncertainty in closure is 3.9 ± 0.7 (dimensionless). Our results show that uncertainties in terrestrial evaporation contribute twice as much as precipitation uncertainties to errors in the closure of water and energy balance. Moreover, our results highlight the need for improving satellite‐based precipitation and evaporation data in humid temperate forests, where the closure error in the Budyko space is as high as 1.1 ± 0.3, compared to only 0.2 ± 0.03 in tropical forests. Comparing the results with land surface model‐based data sets driven by in situ precipitation, we find that Earth observation‐based data sets perform better in regions where precipitation gauges are sparse. These findings have implications for improving the understanding of global hydrology and regional water management and can guide the development of satellite remote sensing‐based data sets and Earth system models., Key Points A Budyko‐based approach to water and energy balance closure mitigates the need for runoff dataErrors in water and energy balance closure are influenced more by uncertainties in evaporation rather than precipitationInability of Earth observations to close the water and energy balance of temperate forests

The ongoing Coronavirus Disease 2019 (COVID‐19) has posed a serious threat to human public health and global economy. Population mobility is an important factor that drives the spread of COVID‐19. This study aimed to quantitatively evaluate the impact of population flow on the spread of COVID‐19 from a spatiotemporal perspective. To this end, a case study was carried out in Hubei Province, which was once the most affected area of COVID‐19 outbreak in Mainland China. The geographically and temporally weighted regression (GTWR) model was applied to model the spatiotemporal association between COVID‐19 epidemic and population mobility. Two patterns of population flows, including the population inflow from Wuhan and intra‐city population movement, were considered to construct explanatory variables. Results indicate that the GTWR model can reveal the spatial–temporal‐varying relationships between COVID‐19 and population mobility. Moreover, the association between COVID‐19 case counts and population movements presented three stages of temporal variation characteristics due to the virus incubation period and implementation of strict lockdown measures. In the spatial dimension, evident geographical disparities were observed across Hubei Province. These findings can provide policymakers useful knowledge about the impact of population movement on the spatio‐temporal transmission of COVID‐19. Thus, targeted interventions, if necessary in certain time periods, can be implemented to restrict population flow in cities with high transmission risk., Key Points To quantitatively evaluate how the effect of population flows on COVID‐19 spread varies over time and space in Hubei Province, ChinaThe association showed three stages of temporal variation characteristics and strong geographical disparities across Hubei ProvinceThe implementation of control measures by restricting population movements can play a key role in mitigating the spread of this epidemic

Because cross‐polarized radar returns are highly associated with volume scatter, radar polarimetry returns tend to show strong evidence of wildfire scars and recovery in forest and chaparral. We focus on the polarimetry images from UAVSAR (PolSAR) line SanAnd_08525, which covers a roughly 20 km wide swath over the Transverse Range including parts of the Santa Monica, San Gabriel and San Bernardino Mountains. We select images from four acquisition dates from October 2009 to September 2020, very roughly 4 years apart. These are compared to fire perimeters from the national Geospatial Multi‐Agency Coordination and NIFC databases for years 2003–2020, which shows the areas affected by the major fires (west to east) Springs2013, Woolsey2018, Topanga2005, LaTuna2017, Station2009, BlueCut2016, Pilot2016, Slide2007, Butler2007, and many smaller fires. PolSAR images are shown to be helpful in identifying types and boundaries of fire, 50‐meter scale details of vegetation loss, and variability of vegetation recovery in post‐fire years., Key Points The urban‐facing part of the Transverse Range is imaged by 27 Uninhabited Aerial Vehicle Synthetic Aperture Radar visits from 2009 to 2020, including locations of over 20 wildfiresPolSAR images identify types and boundaries of fire, 50 m scale details of vegetation loss, and variability of post‐fire recoveryObserved signs of biomass recovery in gap areas between recent fires provide clues to the locations and extents of near‐future fires

Six conceptually different transport models were applied to the macrodispersion experiment (MADE)-1 field tracer experiment as a first major attempt for model comparison. The objective was to show that complex mass distributions in heterogeneous aquifers can be predicted without calibration of transport parameters, solely making use of structural and flow data. The models differ in their conceptualization of the heterogeneous aquifer structure, computational complexity, and use of conductivity data obtained from various observation methods (direct push injection logging, DPIL, grain size analysis, pumping tests and flowmeter). They share the same underlying physical transport process of advection by the velocity field solely. Predictive capability is assessed by comparing results to observed longitudinal mass distributions of the MADE-1 experiment. The decreasing mass recovery of the observed plume is attributed to sampling and no physical process like mass transfer is invoked by the models. Measures like peak location and strength are used in comparing the modeled and measured plume mass distribution. Comparison of models reveals that the predictions of the solute plume agree reasonably well with observations, if the models are underlain by a few parameters of close values: mean velocity, a parameter reflecting log-conductivity variability, and a horizontal length scale related to conductivity spatial correlation. The robustness of the results implies that conservative transport models with appropriate conductivity upscaling strategies of various observation data provide reasonable predictions of plumes longitudinal mass distribution, as long as key features are taken into account., The authors thank Marco Bianchi and Boris Baeumer for their support during the development of the study. A. Fiori and A. Bellin acknowledge funding from the Italian Ministry of Education, University and Research (MIUR) in the frame of the Departments of Excellence Initiative 2018–2022 granted to Dept. of Engineering of Roma Tre University, and to the Dept. of Civil, Environmental and Mechanical Engineering of the University of Trento, respectively. M. Dentz acknowledges funding of the European Research Council (ERC) through the project MHetScale (contract number 617511), and the Spanish Research Agency (AEI) through the project HydroPore (contract number PID2019-106887GB-C31).

Cloud‐top heights (CTH) from the Multiangle Imaging Spectroradiometer (MISR) and the Moderate Resolution Imaging Spectroradiometer (MODIS) on Terra constitute our longest‐running single‐platform CTH record from a stable orbit. Here, we provide the first evaluation of the Terra Level 2 CTH record against collocated International Space Station Cloud‐Aerosol Transport System (CATS) lidar observations between 50ºN and 50ºS. Bias and precision of Terra CTH relative to CATS is shown to be strongly tied to cloud horizontal and vertical heterogeneity and altitude. For single‐layered, unbroken, optically thick clouds observed over all altitudes, the uncertainties in MODIS and MISR CTH are −540 ± 690 m and −280 ± 370 m, respectively. The uncertainties are generally smaller for lower altitude clouds and larger for optically thin clouds. For multi‐layered clouds, errors are summarized herein using both absolute CTH and CATS‐layer‐altitude proximity to Terra CTH. We show that MISR detects the lower cloud in a two‐layered system, provided top‐layer optical depth, Key Points We present the first semi‐global (50°N–50°S) comparison of Terra cloud‐top heights with coincident samples from a space‐based lidarUsing lidar as truth, Terra cloud‐top height bias and precision are summarized as a function of cloud geometrical and optical propertiesWith the first measurement of stereo‐opacity bias (−60 to −260 m, depending on cloud type), MISR cloud height error‐budget is closed

Improved understanding of how our coasts will evolve over a range of time scales (years‐decades) is critical for effective and sustainable management of coastal infrastructure. A robust knowledge of the spatial, directional and temporal variability of the inshore wave climate is required to predict future coastal evolution and hence vulnerability. However, the variability of the inshore directional wave climate has received little attention, and an improved understanding could drive development of skillful seasonal or decadal forecasts of coastal response. We examine inshore wave climate at 63 locations throughout the United Kingdom and Ireland (1980–2017) and show that 73% are directionally bimodal. We find that winter‐averaged expressions of six leading atmospheric indices are strongly correlated (r = 0.60–0.87) with both total and directional winter wave power (peak spectral wave direction) at all studied sites. Regional inshore wave climate classification through hierarchical cluster analysis and stepwise multi‐linear regression of directional wave correlations with atmospheric indices defined four spatially coherent regions. We show that combinations of indices have significant skill in predicting directional wave climates (R 2 = 0.45–0.8; p, Key Points Over 70% of inshore wave climates analyzed throughout the United Kingdom and Ireland were directionally bimodalCombinations of winter atmospheric indices NAO, WEPA, SCAND, and EA are significantly correlated with directional wave climate in all regionsRegression models using multiple winter atmospheric indices enable skillful reconstructions of directional wave climate in all regions

Many nations responded to the corona virus disease‐2019 (COVID‐19) pandemic by restricting travel and other activities during 2020, resulting in temporarily reduced emissions of CO2, other greenhouse gases and ozone and aerosol precursors. We present the initial results from a coordinated Intercomparison, CovidMIP, of Earth system model simulations which assess the impact on climate of these emissions reductions. 12 models performed multiple initial‐condition ensembles to produce over 300 simulations spanning both initial condition and model structural uncertainty. We find model consensus on reduced aerosol amounts (particularly over southern and eastern Asia) and associated increases in surface shortwave radiation levels. However, any impact on near‐surface temperature or rainfall during 2020–2024 is extremely small and is not detectable in this initial analysis. Regional analyses on a finer scale, and closer attention to extremes (especially linked to changes in atmospheric composition and air quality) are required to test the impact of COVID‐19‐related emission reductions on near‐term climate., Key Points Lockdown restrictions during COVID‐19 have reduced emissions of aerosols and greenhouse gases12 CMIP6 Earth system models have performed coordinated experiments to assess the impact of this on climateAerosol amounts are reduced over southern and eastern Asia but there is no detectable change in annually averaged temperature or precipitation

Geoenergy and geoengineering applications usually involve fluid injection into and production from fractured media. Accounting for fractures is important because of the strong poromechanical coupling that ties pore pressure changes and deformation. A possible approach to the problem uses equivalent porous media to reduce the computational cost and model complexity instead of explicitly including fractures in the models. We investigate the validity of this simplification by comparing these two approaches. Simulation results show that pore pressure distribution significantly differs between the two approaches even when both are calibrated to predict identical values at the injection and production wells. Additionally, changes in fracture stability are not well captured with the equivalent porous medium. We conclude that explicitly accounting for fractures in numerical models may be necessary under some circumstances to perform reliable coupled thermohydromechanical simulations, which could be used in conjunction with other tools for induced seismicity forecasting., A.Z. acknowledges the financial support received from the “Iran's Ministry of Science, Research and Technology” (PhD students' sabbatical grants) for visiting IDAEA‐CSIC. The contribution of F.P. is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—project number PA 3451/1‐1. R.M. is thankful for the support from US DOE through Carbon SAFE Macon County Project DE‐FE0029381. V.V. acknowledges funding from the Spanish National Research Council (CSIC) through the Intramural project 201730I100 and from the European Research Council (ERC) under the European Union's Horizon 2020 Research and Innovation Program through the Starting Grant GEoREST (www.georest.eu) (Grant agreement No. 801809). IDAEA‐CSIC is a Center of Excellence Severo Ochoa (Spanish Ministry of Science and Innovation, Project CEX2018‐000794‐S). The authors declare no conflict of interest.

Due to the specific hydrothermal conditions of dry‐hot valleys, temperature changes caused by the development of large‐scale hydropower projects may be more extreme than they are in other regions. In this study, we analyzed these temperature changes at four hydropower stations in both dry‐hot and non‐dry‐hot valleys. Based on the calculated relative temperatures of the downstream river and the areas surrounding the reservoirs, we employed two indices to quantify the influence of the reservoirs on the temperatures of these two regions: the downstream river temperature change and the reservoir effect change intensity. Our results are as follows: (a) In the downstream rivers, the temperature regulation effect was more pronounced in the wet season; in the regions surrounding the reservoirs, the temperature regulation effect was more pronounced in the dry season. (b) The downstream river temperature in both the dry‐hot and wet‐hot valleys exhibited noticeable warming in both the wet and dry seasons, while the cold‐dry valley was characterized by cooling in the dry season and warming in the wet season. With the exception of the Liyuan station (where the influence of the reservoir on the downstream temperatures only extended to a distance of 9 km from the dam) during the dry season, the existence of the hydropower stations affected the temperatures of the entire downstream region. (c) For the areas surrounding the reservoir, the presence of a hydropower station mainly caused the temperatures in the dry‐hot valleys to rise and the temperatures in the non‐dry‐hot valleys to decrease., Key Points The difference of surface temperature affected by reservoirs was compared in dry‐hot and non‐dry‐hot valleysTo facilitate comparison of the different downstream river temperature changes, downstream river temperature change index was establishedThe reservoir effect change intensity was established to explore how reservoirs regulate the surface temperature in the surrounding areas

Micromagnetic tomography aims at reconstructing large numbers of individual magnetizations of magnetic particles from combining high‐resolution magnetic scanning techniques with micro X‐ray computed tomography (microCT). Previous work demonstrated that dipole moments can be robustly inferred, and mathematical analysis showed that the potential field of each particle is uniquely determined. Here, we describe a mathematical procedure to recover higher orders of the magnetic potential of the individual magnetic particles in terms of their spherical harmonic expansions (SHE). We test this approach on data from scanning superconducting quantum interference device microscopy and microCT of a reference sample. For particles with high signal‐to‐noise ratio of the magnetic scan we demonstrate that SHE up to order n = 3 can be robustly recovered. This additional level of detail restricts the possible internal magnetization structures of the particles and provides valuable rock magnetic information with respect to their stability and reliability as paleomagnetic remanence carriers. Micromagnetic tomography therefore enables a new approach for detailed rock magnetic studies on large ensembles of individual particles., Key Points Micromagnetic Tomography uniquely recovers higher‐order multipole terms for several individual grains in a sampleHigher order multipole moments are an expression of the internal domain structure of magnetic grainsUltimately, this enables to select individual grains for rock‐ and paleomagnetic studies based on domain configuration

Intensive observations and WRF‐Chem simulations are applied in this study to investigate the adverse impacts of regional transport on the PM2.5 (fine particulate matter; diameter ≤2.5 μm) changes in Shanghai during the Coronavirus Disease 2019 lockdown. As the local atmospheric oxidation capacity was observed to be generally weakened, strong regional transport carried by the frequent westerly winds is suggested to be the main driver of the unexpected pollution episodes, increasing the input of both primary and secondary aerosols. Contributing 40%–80% to the PM2.5, the transport contributed aerosols are simulated to exhibit less decreases (13.2%–21.8%) than the local particles (37.1%–64.8%) in urban Shanghai due to the lockdown, which largely results from the less decreased industrial and residential emissions in surrounding provinces. To reduce the influence of the transport, synergetic emission control, especially synergetic ammonia control, measures are proved to be effective strategies, which need to be considered in future regulations., Key Points PM2.5 pollution in Shanghai were mostly attributed to the regional transport carried by frequent westerly winds during the COVID‐lockdownThe transport contributed aerosols in Shanghai are less affected by the lockdown measures due to the less decreased surrounding emissionsSynergetic emission control measures, especially those on ammonia emissions, are effective in mitigating the adverse transport impacts

We aim at developing a simple model as an interpretative framework for the water vapor isotopic variations in the tropical troposphere over the ocean. We use large‐eddy simulations of disorganized convection in radiative‐convective equilibrium to justify the underlying assumptions of this simple model, to constrain its input parameters and to evaluate its results. We also aim at interpreting the depletion of the water vapor isotopic composition in the lower and midtroposphere as precipitation increases, which is a salient feature in tropical oceanic observations. This feature constitutes a stringent test on the relevance of our interpretative framework. Previous studies, based on observations or on models with parameterized convection, have highlighted the roles of deep convective and mesoscale downdrafts, rain evaporation, rain‐vapor diffusive exchanges, and mixing processes. The interpretative framework that we develop, valid in case of disorganized convection, is a two‐column model representing the net ascent in clouds and the net descent in the environment. We show that the mechanisms for depleting the troposphere as the precipitation rate increases all stem from the higher tropospheric relative humidity. First, when the relative humidity is larger, less snow sublimates before melting and a smaller fraction of rain evaporates. Both effects lead to more depleted rain evaporation and eventually more depleted water vapor. This mechanism dominates in regimes of large‐scale ascent. Second, the entrainment of dry air into clouds reduces the vertical isotopic gradient and limits the depletion of tropospheric water vapor. This mechanism dominates in regimes of large‐scale descent., Key Points Isotopic enrichment of tropospheric water vapor by rain evaporation is stronger when drier air enhances sublimation and evaporationEntrainment of dry air weakens the vertical isotopic gradient and limits the depletion of tropospheric water vaporThese mechanisms explain the increased depletion of tropospheric water vapor as tropospheric relative humidity increases

Although uncommon, nontuberculous mycobacterial (NTM) pulmonary infection in the Hawaiian Islands has a relatively high incidence and mortality compared to the mainland U.S. As a result, this study examines the possible geological and hydrological pathways by which NTM patients may become infected, including the environmental conditions that may favor growth and transport. Previously suggested infection routes include the inhalation of NTM attached to micro‐droplets from infected home plumbing systems and aerosolized dust from garden soil. In this study, we evaluate the possible routes NTM may take from riparian environments, into groundwater, into public water supplies and then into homes. Because NTM are notoriously hydrophobic and prone to attach to surfaces, mineralogy, and surface chemistry of suspended sediment in streams, soils, and rock scrapings suggest that NTM may especially attach to Fe‐oxides/hydroxides, and be transported as particles from losing streams to the aquifer on time‐scales of minutes to days. Within the aquifer, flow models indicate that water may be drawn into production wells on time scales (months) that permit NTM to survive and enter domestic water supplies. These processes depend on the presence of interconnected fracture networks with sufficient aperture to preclude complete autofiltration. The common occurrence of NTM in and around streams, in addition to wells, implies that the natural and built environments are capable of introducing a source of NTM into domestic water supplies via groundwater withdrawals. This may produce a persistent source of NTM infection to individuals through the presence of NTM‐laden biofilms in home plumbing., Key Points Nontuberculous mycobacterial (NTM) are found in soils and biofilms of riparian environments in Hawai'iNTM are likely transported from losing streams to aquifersPumped wells draw NTM into culinary water supplies and into homes. NTM are chlorine resistant and outcompete other taxa in home plumbing

In a rapidly changing world, what is today an unprecedented extreme may soon become the norm. As a result, extreme‐related disasters are expected to become more frequent and intense. This will have widespread socio‐economic consequences and affect the ability of different societal groups to recover from and adapt to rapidly changing environmental conditions. Therefore, there is the need to decipher the relation between genesis of unprecedented events, accumulation and distribution of risk, and recovery trajectories across different societal groups. Here, we develop an analytical approach to unravel the complexity of future extremes and multiscalar societal responses—from households to national governments and from immediate impacts to longer term recovery. This requires creating new forms of knowledge that integrate analyses of the past—that is, structural causes and political processes of risk accumulation and differentiated recovery trajectories—with plausible scenarios of future environmental extremes grounded in the event‐specific literature. We specifically seek to combine the physical characteristics of the extremes with examinations of how culture, politics, power, and policy visions shape societal responses to unprecedented events, and interpret the events as social‐environmental extremes. This new approach, at the nexus between social and natural sciences, has the concrete advantage of providing an impact‐focused vision of future social‐environmental risks, beyond what is achievable within conventional disciplinary boundaries. In this paper, we focus on extreme flooding events and the societal responses they elicit. However, our approach is flexible and applicable to a wide range of extreme events. We see it as the first building block of a new field of research, allowing for novel and integrated theoretical explanations and forecasting of social‐environmental extremes., Key Points We conceptualize unprecedented extremes as social‐environmental processes shaped by institutional, political, and economic changeAs social‐environmental extremes become more frequent, there is an urgency to unravel their genesis and the possible societal responsesThis approach is the first building block of a new field of research in social‐environmental extreme event forecasts

A growing constellation of satellites is providing near‐global coverage of column‐averaged CO2 observations. Launched in 2019, NASA’s OCO‐3 instrument is set to provide XCO2 observations at a high spatial and temporal resolution for regional domains (100 × 100 km). The atmospheric column version of the Stochastic Time‐Inverted Lagrangian Transport (X‐STILT) model is an established method of determining the influence of upwind sources on column measurements of the atmosphere, providing a means of analysis for current OCO‐3 observations and future space‐based column‐observing missions. However, OCO‐3 is expected to provide hundreds of soundings per targeted observation, straining this already computationally intensive technique. This work proposes a novel scheme to be used with the X‐STILT model to generate upwind influence footprints with less computational expense. The method uses X‐STILT generated influence footprints from a key subset of OCO‐3 soundings. A nonlinear weighted averaging is applied to these footprints to construct additional footprints for the remaining soundings. The effects of subset selection, meteorological data, and topography are investigated for two test sites: Los Angeles, California, and Salt Lake City, Utah. The computational time required to model the source sensitivities for OCO‐3 interpretation was reduced by 62% and 78% with errors smaller than other previously acknowledged uncertainties in the modeling system (OCO‐3 retrieval error, atmospheric transport error, prior emissions error, etc.). Limitations and future applications for future CO2 missions are also discussed., Key Points Determining sources of spatially dense XCO2 observations with LPDM techniques can become time intensive and strain computational resourcesPresented in this work is an interpolation scheme that eases the computational burden of spatially dense XCO2 source determination studiesEvaluating the efficiency of this interpolation scheme revealed reductions of >50% in computational time at two testing locations

The intensification of heat stress reduces the labor capacity and hence poses a threat to socio‐economic development. The reliable projection of the changing climate and the development of sound adaptation strategies are thus desired for adapting to the decreasing labor productivity under climate change. In this study, an optimization modeling approach coupled with dynamical downscaling is proposed to design the optimal adaptation strategies for improving labor productivity under heat stress in China. The future changes in heat stress represented by the wet‐bulb globe temperature (WBGT) are projected with a spatial resolution of 25 × 25 km by a regional climate model (RCM) through the dynamical downscaling of its driving global climate model (GCM). Uncertain information such as system costs, environmental costs, and subsidies are also incorporated into the optimization process to provide reliable decision alternatives for improving labor productivity. Results indicate that the intensification of WBGT is overestimated by the GCM compared to the RCM. Such an overestimation can lead to more losses in working hours derived from the GCM than those from the RCM regardless of climate scenarios. Nevertheless, the overestimated heat stress does not alter the regional measures taken to adapt to decreasing labor productivity. Compared to inland regions, the monsoon‐affected regions tend to improve labor productivity by applying air conditioning rather than working overtime due to the cost differences. Consequently, decision‐makers need to optimally make a balance between working overtime and air conditioning measures to meet sustainable development goals., Key Points An optimization model coupled with dynamical downscaling is developed to improve labor productivity under intensified heat stress in ChinaCivilian workers in inland regions will be more vulnerable to the intensified heat stress in a changing climateCompared to the regional climate model, the global climate model overestimates heat stress and thus exaggerates related system costs to recover from the working hour loss

Supraglacial debris affects glacier mass balance as a thin layer enhances surface melting, while a thick layer reduces it. While many glaciers are debris‐covered, global glacier models do not account for debris because its thickness is unknown. We provide the first globally distributed debris thickness estimates using a novel approach combining sub‐debris melt and surface temperature inversion methods. Results are evaluated against observations from 22 glaciers. We find the median global debris thickness is ∼0.15 ± 0.06 m. In all regions, the net effect of accounting for debris is a reduction in sub‐debris melt, on average, by 37%, which can impact regional mass balance by up to 0.40 m water equivalent (w.e.) yr‐1. We also find recent observations of similar thinning rates over debris‐covered and clean ice glacier tongues is primarily due to differences in ice dynamics. Our results demonstrate the importance of accounting for debris in glacier modeling efforts., Key Points We produce the first distributed global debris thickness estimatesAccounting for debris significantly reduces regional glacier mass lossThe similar thinning rates of debris‐covered and clean ice glaciers in High Mountain Asia is primarily caused by differences in ice dynamics

The COVID‐19 pandemic led to a 72% reduction of air traffic over Europe in March–August 2020 compared to 2019. Modeled contrail cover declined similarly, and computed mean instantaneous radiative contrail forcing dropped regionally by up to 0.7 W m−2. Here, model predictions of cirrus optical thickness and the top‐of‐atmosphere outgoing longwave and reflected shortwave irradiances are tested by comparison to Meteosat‐SEVIRI‐derived data. The agreement between observations and modeled data is slightly better when modeled contrail cirrus contributions are included. The spatial distributions and diurnal cycles of the differences in these data between 2019 and 2020 are partially caused by differences in atmospheric and surface conditions, particularly for solar radiation in the spring of 2020. Aviation signals become discernible in the observed differences of these data between 2019 and 2020 when subtracting numerical weather prediction model results that approximate the atmosphere and surface conditions without contrails., Key Points Air traffic reduction during the COVID‐19 pandemic is used for testing aviation contrail cirrus and related radiative forcing modelingCirrus optical depth and top‐of‐atmosphere irradiances from satellite observations and model results show reasonable correlationsAviation contrail signatures are identifiable in spite of strongly different atmosphere and surface conditions in 2019 and 2020

The decline in global emissions of carbon dioxide due to the COVID‐19 pandemic provides a unique opportunity to investigate the sensitivity of the global carbon cycle and climate system to emissions reductions. Recent efforts to study the response to these emissions declines has not addressed their impact on the ocean, yet ocean carbon absorption is particularly susceptible to changing atmospheric carbon concentrations. Here, we use ensembles of simulations conducted with an Earth system model to explore the potential detection of COVID‐related emissions reductions in the partial pressure difference in carbon dioxide between the surface ocean and overlying atmosphere (ΔpCO2), a quantity that is regularly measured. We find a unique fingerprint in global‐scale ΔpCO2 that is attributable to COVID, though the fingerprint is difficult to detect in individual model realizations unless we force the model with a scenario that has four times the observed emissions reduction., Key Points COVID‐related emissions reductions will be imperceptible in surface ocean pH observationsThe CanESM5 COVID ensemble predicts a unique fingerprint of COVID‐related emissions reductions in global mean ΔpCO2 (pCO2oc ‐ pCO2atm)The fingerprint is potentially detectable in global‐scale observations of ΔpCO2, but only with large emissions reductions

Humboldt Gletscher is a 100‐km wide, slow‐moving glacier in north Greenland which holds a 19‐cm global sea level equivalent. Humboldt has been the fourth largest contributor to sea level rise since 1972 but the cause of its mass loss has not been elucidated. Multi‐beam echo sounding data collected in 2019 indicate a seabed 200 m deeper than previously known. Conductivity temperature depth data reveal the presence of warm water of Atlantic origin at 0°C at the glacier front and a warming of the ocean waters by 0.9 ± 0.1°C since 1962. Using an ocean model, we reconstruct grounded ice undercutting by the ocean, combine it with calculated retreat caused by ice thinning to floatation, and are able to fully explain the observed retreat. Two thirds of the retreat are caused by undercutting of grounded ice, which is a physical process not included in most ice sheet models., Key Points The 100‐km wide Humboldt Gletscher holds a 19‐cm sea level rise equivalent, lost 161 billion tons of mass, and retreated 13 km since 1972Warm waters at 0°C flood a 350–400 m deep trough on its northern flank that remains below sea level more than 100 km inlandWe explain the glacier retreat as 70% from ocean‐induced undercutting and 30% from thinning‐induced retreat

Exposure to wildfire smoke increases the risk of respiratory and cardiovascular hospital admissions. Health impact assessments, used to inform decision‐making processes, characterize the health impacts of environmental exposures by combining preexisting epidemiological concentration–response functions (CRFs) with estimates of exposure. These two key inputs influence the magnitude and uncertainty of the health impacts estimated, but for wildfire‐related impact assessments the extent of their impact is largely unknown. We first estimated the number of respiratory, cardiovascular, and asthma hospital admissions attributable to fire‐originated PM2.5 exposure in central California during the October 2017 wildfires, using Monte Carlo simulations to quantify uncertainty with respect to the exposure and epidemiological inputs. We next conducted sensitivity analyses, comparing four estimates of fire‐originated PM2.5 and two CRFs, wildfire and nonwildfire specific, to understand their impact on the estimation of excess admissions and sources of uncertainty. We estimate the fires accounted for an excess 240 (95% CI: 114, 404) respiratory, 68 (95% CI: −10, 159) cardiovascular, and 45 (95% CI: 18, 81) asthma hospital admissions, with 56% of admissions occurring in the Bay Area. Although differences between impact assessment methods are not statistically significant, the admissions estimates' magnitude is particularly sensitive to the CRF specified while the uncertainty is most sensitive to estimates of fire‐originated PM2.5. Not accounting for the exposure surface's uncertainty leads to an underestimation of the uncertainty of the health impacts estimated. Employing context‐specific CRFs and using accurate exposure estimates that combine multiple data sets generates more certain estimates of the acute health impacts of wildfires., Key Points PM2.5 exposure during the 2017 California fires accounted for an estimated 308 excess respiratory and cardiovascular hospital admissionsHealth impact assessments, used to inform decision‐making processes, are sensitive to the exposure and epidemiologic inputs specifiedAccurate exposure estimates and context‐specific health risk functions estimate fire‐attributable health impacts with greater certainty

We assess the value of machine learning as an accelerator for the parameterization schemes of operational weather forecasting systems, specifically the parameterization of nonorographic gravity wave drag. Emulators of this scheme can be trained to produce stable and accurate results up to seasonal forecasting timescales. Generally, networks that are more complex produce emulators that are more accurate. By training on an increased complexity version of the existing parameterization scheme, we build emulators that produce more accurate forecasts. For medium range forecasting, we have found evidence that our emulators are more accurate than the version of the parametrization scheme that is used for operational predictions. Using the current operational CPU hardware, our emulators have a similar computational cost to the existing scheme, but are heavily limited by data movement. On GPU hardware, our emulators perform 10 times faster than the existing scheme on a CPU., Key Points Nonorographic gravity wave drag parametrization can be accurately emulated with a neural networkThese emulators produce accurate and stable forecasts over long timescalesNeural networks can reduce the cost of an increased complexity scheme

17 pags., 9 figs., 1 tab., Ocean-going ships supply products from one region to another and contribute to the world’s economy. Ship exhaust contains many air pollutants and results in significant changes in marine atmospheric composition. The role of reactive halogen species (RHS) in the troposphere has received increasing recognition and oceans are the largest contributors to their atmospheric burden. However, the impact of shipping emissions on RHS and that of RHS on ship-originated air pollutants have not been studied in detail. Here, an updated Weather Research Forecasting coupled with Chemistry model is utilized to explore the chemical interactions between ship emissions and oceanic RHS over the East Asia seas in summer. The emissions and resulting chemical transformations from shipping activities increase the level of NO and NO at the surface, increase O in the South China Sea, but decrease O in the East China Sea. Such changes in pollutants result in remarkable changes in the levels of RHS (>200% increase of chlorine; ∼30% and ∼5% decrease of bromine and iodine, respectively) as well as in their partitioning. The abundant RHS, in turn, reshape the loadings of air pollutants (∼20% decrease of NO and NO; ∼15% decrease of O) and those of the oxidants (>10% reduction of OH and HO; ∼40% decrease of NO) with marked patterns along the ship tracks. We, therefore, suggest that these important chemical interactions of ship-originated emissions with RHS should be considered in the environmental policy assessments of the role of shipping emissions in air quality and climate., This study received funding from the European Research Council Executive Agency under the European Union’s Horizon 2020 Research and Innovation Program (Project ERC-2016- COG 726349 CLIMAHAL) and was supported by the Consejo Superior de Investigaciones Científicas (CSIC) of Spain and The Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning and Shanghai Thousand Talents Program. The development and maintenance of the WRF-Chem model are conducted by NOAA/ESRL/GSD in active collaboration with other institutes. Computing resources, support, and data storage were provided by the Climate Simulation Laboratory at NCAR’s Computational and Information Systems Laboratory (CISL), sponsored by the NSF. International cooperation between CSIC (Spain) and CONICET (Argentina) was supported by the i-COOP program (B20331). R.P.F. would like to thank ANPCyT (PICT 2015-0714), UNCuyo (SeCTyP M032/3853), and UTN (PID 4920-194/2018) for financial support. IITM is funded by the Ministry of Earth Sciences, Government of India.

Irrigated rice requires intense water management under typical agronomic practices. Cost effective tools to improve the efficiency and assessment of water use is a key need for industry and resource managers to scale ecosystem services. In this research we advance model‐based decomposition and machine learning to map inundated rice using time‐series polarimetric, L‐band Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) observations. Simultaneous ground truth observations recorded water depth inundation during the 2019 crop season using instrumented fields across the study site in Arkansas, USA. A three‐component model‐based decomposition generated metrics representing surface‐, double bounce‐, and volume‐scattering along with a shape factor, randomness factor, and the Radar Vegetation Index (RVI). These physically meaningful metrics characterized crop inundation status independent of growth stage including under dense canopy cover. Machine learning (ML) comparisons employed Random Forest (RF) using the UAVSAR derived parameters to identify cropland inundation status across the region. Outcomes show that RVI, proportion of the double‐bounce within total scattering, and the relative comparison between the double‐bounce and the volume scattering have moderate to strong mechanistic ability to identify rice inundation status with Overall Accuracy (OA) achieving 75%. The use of relative ratios further helped mitigate the impacts of far range incidence angles. The RF approach, which requires training data, achieved a higher OA and Kappa of 88% and 71%, respectively, when leveraging multiple SAR parameters. Thus, the combination of physical characterization and ML provides a powerful approach to retrieving cropland inundation under the canopy. The growth of polarimetric L‐band availability should enhance cropland inundation metrics beyond open water that are required for tracking water quantity at field scale over large areas., Key Points Cropland inundation assessment has largely focused on open waterQuad polarized L‐band SAR can help detect under canopy inundationThe underlying physical mechanisms driving scattering responses and machine learning algorithms had similar outcomes

10 pags., 4 figs., Many Chemistry-Climate Models (CCMs) include a simplified treatment of brominated very short-lived (VSL) species by assuming CHBr as a surrogate for VSL. However, neglecting a comprehensive treatment of VSL in CCMs may yield an unrealistic representation of the associated impacts. Here, we use the Community Atmospheric Model with Chemistry (CAM-Chem) CCM to quantify the tropospheric and stratospheric changes between various VSL chemical approaches with increasing degrees of complexity (i.e., surrogate, explicit, and full). Our CAM-Chem results highlight the improved accuracy achieved by considering a detailed treatment of VSL photochemistry, including sea-salt aerosol dehalogenation and heterogeneous recycling on ice-crystals. Differences between the full and surrogate schemes maximize in the lowermost stratosphere and midlatitude free troposphere, resulting in a latitudinally dependent reduction of ∼1–7 DU in total ozone column and a ∼5%–15% decrease of the OH/HO ratio. We encourage all CCMs to include a complete chemical treatment of VSL in the troposphere and stratosphere., This study has been funded by the European Union's Horizon 2020 Re-search and Innovation program (Project ‘ERC-2016-COG 726349 CLIMAHAL’), and supported by the Consejo Superior de Investigaciones Científicas of Spain. Computing resources and support are provided by the Computational and Information System Laboratory (CISL,2017). R. P. Fernandez would like to thank financial support from CONICET, ANPCyT (PICT 2015-0714), UNCuyo (SeCTyP M032/3853) and UTN (PID 4920-194/2018). NCAR is sponsored by NSF under grant number 1852977. R. J. Salawitch appreciates support from the NASA (grant ACMP 80NSSC19K0983). The authors thank two anonymous reviewers for their constructive comments

We present a new solar irradiance reference spectrum representative of solar minimum conditions between solar cycles 24 and 25. The Total and Spectral Solar Irradiance Sensor‐1 (TSIS‐1) Hybrid Solar Reference Spectrum (HSRS) is developed by applying a modified spectral ratio method to normalize very high spectral resolution solar line data to the absolute irradiance scale of the TSIS‐1 Spectral Irradiance Monitor (SIM) and the CubeSat Compact SIM (CSIM). The high spectral resolution solar line data are the Air Force Geophysical Laboratory ultraviolet solar irradiance balloon observations, the ground‐based Quality Assurance of Spectral Ultraviolet Measurements In Europe Fourier transform spectrometer solar irradiance observations, the Kitt Peak National Observatory solar transmittance atlas, and the semi‐empirical Solar Pseudo‐Transmittance Spectrum atlas. The TSIS‐1 HSRS spans 202–2730 nm at 0.01 to ∼0.001 nm spectral resolution with uncertainties of 0.3% between 460 and 2365 nm and 1.3% at wavelengths outside that range., Key Points The TSIS‐1 Spectral Irradiance Monitor and Compact SIM instruments observe the Sun's irradiance spectrum at high accuracyThe TSIS‐1 Hybrid Solar Reference Spectrum consists of high resolution solar line data normalized to the TSIS‐1 SIM irradiance spectrumThe TSIS‐1 Hybrid Solar Reference Spectrum has at least 0.01 nm spectral resolution, spans 202–2730 nm, and is accurate to 0.3%–1.3%

Low‐cloud based emergent constraints have the potential to substantially reduce uncertainty in Earth’s equilibrium climate sensitivity, but recent work has shown that previously developed constraints fail in the latest generation of climate models, suggesting that new approaches are needed. Here, we investigate the potential for emergent constraints to reduce uncertainty in regional cloud feedbacks, rather than the global‐mean cloud feedback. Strong relationships are found between the monthly and interannual variability of tropical clouds, and the tropical net cloud feedback. These relationships are combined with observations to substantially narrow the uncertainty in the tropical cloud feedback and demonstrate that the tropical cloud feedback is likely >0Wm−2K−1. Promising relationships are also found in the 90°–60°S and 30°–60°N regions, though these relationships are not robust across model generations and we have not identified the associated physical mechanisms., Key Points Global low‐cloud‐based emergent constraints on equilibrium climate sensitivity fail in CMIP6Strong relationships are found between unforced cloud variability and long‐term cloud feedbacks in several regionsRegional emergent constraints suggest the tropical cloud feedback is likely positive