GRL

Knowledge regarding the abundance and distribution of solar wind (SW)-sourced water (OH/H2O) on the Moon in the shallow subsurface remains limited. Here, we report the NanoSIMS measurements of H abundances and D/H ratios on soil grains from three deepest sections of the Chang'E-5 drill core sampled at depths of 0.45–0.8 m. High water contents of 0.13–1.3 wt.% are present on approximately half of the grain surfaces (topmost ∼100 nm), comparable to the values of Chang'E-5 scooped soils. The extremely low δD values (as low as −995‰) and negative correlations between δD and water contents indicate that SW implantation is an important source of water beneath the lunar surface. The results are indicative of homogeneous distribution of SW-derived water in the vertical direction, providing compelling evidence for the well-mixed nature of the lunar regolith. Moreover, the findings demonstrate that the shallow subsurface regolith of the Moon contains a considerable amount of water.

The sub-seasonal CO2 flux (FCO2) variability across the Southern Ocean is poorly understood due to sparse observations at the required temporal and spatial scales. Twinned surface and profiling gliders experiments were used to investigate how storms influence FCO2 through the air-sea gradient in partial pressure of CO2 (ΔpCO2) in the sub-Antarctic zone. Winter-spring storms caused ΔpCO2 to weaken (by 22–37 μatm) due to mixing/entrainment and weaker stratification. This weakening in ΔpCO2 was in phase with the increase in wind stress resulting in a reduction of the storm-driven CO2 uptake by 6%–27%. During summer, stronger stratification explained the weaker sensitivity of ΔpCO2 to storms, instead temperature changes dominated the ΔpCO2 variability. These results highlight the importance of observing synoptic-scale variability in ΔpCO2, the absence of which may propagate significant biases to the mean annual FCO2 estimates from large-scale observing programmes and reconstructions.

The Tibetan Plateau (TP) is an area highly sensitive to climate change and is warming faster than the global average. The TP temperature change has a significant impact on the local ecological environment and the downstream weather and climate. The TP will undoubtedly warm in the future, but the warming extent is uncertain. Using the Coupled Model Inter-comparison Project Phase 6 multi-model ensemble, we found that models simulating smaller TP temperature increases in recent decades tend to project weaker warming in the future. This relationship is driven by the simulation of snowmelt response to greenhouse gas increases, as snow-related albedo feedback dominates the TP temperature changes in both historical and future periods. Based on a two-step emergent constraint approach, the rectified TP warming magnitude increases by about 0.3°C compared to the unconstrained result under both the medium and high emission scenarios, and the inter-model uncertainty is reduced by about 60%.

No abstract is available for this article.

Using the high-time-resolution data from the Magnetospheric Multiscale mission, precursor waves upstream of foreshock transient (FT) shocks are statistically investigated using the four-spacecraft timing method. The wave frequencies and wave vectors determined in the plasma rest frame (PRF) are shown to follow the cold plasma dispersion relation for whistler waves. Combining with the feature of the right-hand polarization in the PRF, the precursors are identified as whistler-mode waves around the lower hybrid frequency. The occurrence of whistler precursors is independent of the Alfvén Mach number and the FT shock normal angle. More importantly, all the whistler precursors have group velocities pointing upstream in the shock frame, suggesting the dispersive radiation to be a possible generation mechanism. The study improves the understanding of not only the whistler precursors but also the overall FT shock dynamics.

The Surface Ocean CO2 Atlas (SOCAT) of CO2 fugacity (fCO2) observations is a key resource supporting annual assessments of CO2 uptake by the ocean and its side effects on the marine ecosystem. SOCAT data are usually released with a lag of up to 1.5 years which hampers timely quantification of recent variations of carbon fluxes between the Earth System components, not only with the ocean. This study uses a statistical ensemble approach to analyze fCO2 with a latency of one month only based on the previous SOCAT release and a series of predictors. Results indicate a modest degradation in a retrospective prediction test for 2021–2022. The generated fCO2 and fluxes for January–August 2023 show a progressive reduction in the Equatorial Pacific source following the La Niña retreat. A breaking-record decrease in the northeastern Atlantic CO2 sink has been diagnosed on account of the marine heatwave event in June 2023.

Electromagnetic ion cyclotron waves in the Earth's outer radiation belt drive rapid electron losses through wave-particle interactions. The precipitating electron flux can be high in the hundreds of keV energy range, well below the typical minimum resonance energy. One of the proposed explanations relies on nonresonant scattering, which causes pitch-angle diffusion away from the fundamental cyclotron resonance. Here we propose the fractional sub-cyclotron resonance, a second-order nonlinear effect that scatters particles at resonance order n = 1/2, as an alternate explanation. Using test-particle simulations, we evaluate the precipitation ratios of sub-MeV electrons for wave packets with various shapes, amplitudes, and wave normal angles. We show that the nonlinear sub-cyclotron scattering produces larger ratios than the nonresonant scattering when the wave amplitude reaches sufficiently large values. The ELFIN CubeSats detected several events with precipitation ratio patterns matching our simulation, demonstrating the importance of sub-cyclotron resonances during intense precipitation events.

A 2-D GCPIC simulation in a dipole field system has been conducted to explore the excitation of oblique whistler mode chorus waves driven by energetic electrons with temperature anisotropy. The rising tone chorus waves are initially generated near the magnetic equator, consisting of a series of subpackets, and become oblique during their propagation. It is found that electron holes in the wave phase space, which are formed due to the nonlinear cyclotron resonance, oscillate in size with time during subpacket formation. The associated inhomogeneity factor varies accordingly, giving rise to various frequency chirping in different phases of subpackets. Distinct nongyrotropic electron distributions are detected in both wave gyrophase and stationary gyrophase. Landau resonance is found to coexist with cyclotron resonance. This study provides multidimensional electron distributions involved in subpacket formation, enabling us to comprehensively understand the nonlinear physics in chorus wave evolution.

Arctic near-surface wind speed (NWS) plays an increasingly crucial role in influencing the local air-sea interactions and the safety of trans-Arctic shipping, but its potential changes in a warming climate and underlying causes remain unclear. Using reanalysis and model simulation data sets, we reveal that the Arctic NWS has increased remarkably since the 1960s, with the strongest increase in the Arctic Ocean surface. We propose that the acceleration of Arctic NWS is primarily driven by reduced stability in the lower troposphere due to increased upward heat fluxes and decreased surface roughness owing to the losses of Arctic glaciers and sea ice in a warming climate. In addition, the coupled climate models project a robust increase in the Arctic NWS under various warming scenarios during the 21st century, especially in the vicinity of the Kara Sea and the Beaufort Sea.

Determination of the key vertical level for cloud condensation nuclei concentration (CCNC) explosions has been a long-term issue in CCN-cloud interaction studies. An idealized hailstorm is simulated with 37 sensitivity runs, including an initial CCNC grouping vertically from the ground to the cloud top, increasing from 100 to 3,000 mg−1. The results reveal a key zone from 750 to 800 hPa near the median boundary layer, where an explosion of CCNC plays a dominant role in the nonmonotonic response of the hail precipitation rate. The explosion of CCNC in this zone could initially result in the condensation of more water vapor into the clouds, which could be transported to a greater vertical extent to significantly affect the riming collection efficiency. However, the dominant zone for the total precipitation rate is wider at heights of 700–800 hPa due to the lower sensitivity of the riming collection efficiency.

Global surface temperature short-term trends fluctuate between cooling and fast-warming under the combined action of external forcing and internal variability, significantly influencing the detectability of near-term climate change. A key driver of these variations is anthropogenic aerosols (AAs), which have undergone a non-monotonic evolution with rapid reduction in recent decades. However, their reduction is projected to decelerate under a high carbon emission scenario, yet the impact on surface temperature trends remains unknown. Here, using initial-perturbation large ensembles, we find that future slowdown in AA reduction over Europe and North America expedites the subpolar North Atlantic surface warming by intensifying the Atlantic meridional overturning circulation. Further, it accelerates the South Indian Ocean and Southern Ocean surface warming through positive low-cloud feedback and oceanic dynamical adjustment, triggered by the poleward migration of westerlies under interhemispheric energy constraint. These AA-driven warmings exacerbate greenhouse warming, significantly enhancing the detectability of local decadal warming trends.

The ionospheric convection electric field is often assumed to be a potential field. This assumption is not always valid, especially when the ionosphere changes on short time scales T≲5 $T\lesssim 5$ min. We present a technique for estimating the induction electric field using ground magnetometer measurements. The technique is demonstrated on real and simulated data for sudden increases in solar wind dynamic pressure of ∼ ${\sim} $1 and 10 nPa, respectively. For the real data, the ionospheric induction electric field is 0.15 ± $\pm $ 0.015 mV/m, and the corresponding compressional flow is 2.5 ± $\pm $ 0.3 m/s. For the simulated data, the induction electric field and compressional flow reach 3 mV/m and 50 m/s, respectively. The induction electric field can locally constitute tens of percent of the total electric field. Inclusion of the induction electric field increased the total Joule heating by 2.4%. Locally the Joule heating changed by tens of percent. This corresponds to energy dissipation that is not accounted for in existing models.

While remote sensing has provided extensive insights into the global terrestrial water, carbon, and energy cycles, space-based retrievals remain limited in observing the belowground influence of the full soil moisture (SM) profile on ecosystem function. We show that this gap can be addressed when coupling 70 m resolution ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station retrievals of land surface temperature (LST) with in-situ SM profile measurements. These data sets together reveal that ecosystem water use decreases with depth with 93% of sites showing significant LST coupling with SM shallower than 20 cm while 34% of sites have interactions with SM deeper than 50 cm. Furthermore, the median depth of peak ecosystem water use is estimated to be 10 cm, though forests have more common peak interactions with deeper soil layers (50–100 cm) in 37% of cases. High spatial resolution remote sensing coupled with field-level data can thus elucidate the role of belowground processes on land surface behavior.

Sulfur-rich volcanic eruptions happen sporadically. If Stratospheric Aerosol Injection (SAI) were to be deployed, it is likely that explosive volcanic eruptions would happen during such a deployment. Here we use an ensemble of Earth System Model simulations to show how changing the injection strategy post-eruption could be used to reduce the climate risks of a large volcanic eruption; the risks are also modified even without any change to the strategy. For a medium-size eruption (10 Tg-SO2) comparable to the SAI injection rate, the volcanic-induced cooling would be reduced if it occurs under SAI, especially if artificial sulfur dioxide injections were immediately suspended. Alternatively, suspending injection only in the eruption hemisphere and continuing injection in the opposite would reduce shifts in precipitation in the tropical belt and thus mitigate eruption-induced drought. Finally, we show that for eruptions much larger than the SAI deployment, changes in SAI strategy would have minimal effect.

The lower limb of the Atlantic meridional overturning circulation (AMOC) is the equatorward flow of dense waters formed through the cooling and freshening of the poleward-flowing upper limb. In the subpolar North Atlantic (SPNA), upper limb variability is primarily set by the North Atlantic Current, whereas lower limb variability is less well understood. Using observations from a SPNA mooring array, we show that variability of the AMOC's lower limb is connected to poleward flow in the interior Irminger Sea. We identify this poleward flow as the northward branch of the Irminger Gyre (IG), accounting for 55% of the AMOC's lower limb variability. Over 2014–2018, wind stress curl fluctuations over the Labrador and Irminger Seas drive this IG and AMOC variability. On longer (>annual) timescales, however, an increasing trend in the thickness of intermediate water, from 2014 to 2020, within the Irminger Sea coincides with a decreasing trend in IG strength.

The two phases of El-Niño-Southern Oscillation (ENSO) influence both regional and global terrestrial vegetation productivity on inter-annual scales. However, the major drivers for the regional vegetation productivity and their controlling strengths during different phases of ENSO remain unclear. We herein disentangled the impacts of two phases of ENSO on regional carbon cycle using multiple data sets. We found that soil moisture predominantly accounts for ∼40% of the variability in regional vegetation productivity during ENSO events. Our results showed that the satellite-derived vegetation productivity proxies, gross primary productivity from data-driven models (FLUXCOM) and observation-constrained ecosystem model (Carbon Cycle Data Assimilation System) generally agree in depicting the contribution of soil moisture and air temperature in modulating regional vegetation productivity. However, the ensemble of weakly constrained ecosystem models exhibits non-negligible discrepancies in the roles of vapor pressure deficit and radiation over extra-tropics. This study highlights the significance of water in regulating regional vegetation productivity during ENSO.

Despite its importance for the global oxidative capacity, spatially resolved trends and variability of the hydroxyl radical (OH) are poorly constrained. We demonstrate the utility of a tropospheric column OH (TCOH) product, created from machine learning and satellite proxy data, in determining the spatial variability in trends of tropical OH over the oceans during September through November. While OH increases domain-wide by 2.1%/decade from 2005–2019, we find significant spatial heterogeneity in regional trends, with decreases in some areas of 2.5%/decade. Our analysis of the trends in the proxy data indicate anthropogenic-driven changes in emissions of OH drivers as well as increasing temperatures cause these trends. This OH product is potentially a significant advance in constraining OH spatial variability and serves as a useful complement to existing tools in understanding the atmospheric oxidative capacity. Comprehensive observations of TCOH are required to assess the fidelity of this method.

Drainage basins are fundamental units of Earth's surface, describing how flows accumulate across landscapes. They are direct expressions of how tectonics and climatic forces alter Earth's surface morphology. Here, we measure the width-to-length ratios (WLRs) of 386,931 drainage basins (average area ∼157 km2), covering all continents except Antarctica and Greenland. Global variations in WLRs are correlated with climatic aridity, whole-basin slope, and local topographic roughness. Basins in arid landscapes tend to be narrower, potentially reflecting a higher prevalence of surface runoff and therefore a stronger slope-parallel component of the transporting flow. Local topographic roughness is associated with wider basins, potentially reflecting greater dispersion of flow directions. Conversely, whole-basin topographic gradients, potentially reflecting gradients in uplift, are associated with narrower basins. However, steeper basins are also often rougher, so revealing the effects of whole-basin slope requires correcting for the confounding effects of roughness variations.

The influence of Atlantic Niño on the following El Niño–Southern Oscillation becomes significant since mid-1970s. However, exact mechanisms for this inter-decadal change are still unclear. Here, we perform a set of model pacemaker experiments to probe the relative contributions of the changes in the Atlantic Niño itself and the mean-state under global warming. The results suggest that the warmer background of the tropical Atlantic plays an essential role in enhancing local mean precipitation, inducing stronger divergence and low-level easterlies in the Pacific. Under a favorable condition in the Pacific, even a weak Atlantic Niño-related warming could promote the development of La Niña through cross-basin Walker circulation and the Indian Ocean-relayed Kelvin wave response. In contrast, the Atlantic Niño pattern change itself induces feeble convection anomalies in the western Atlantic, which cannot induce significant atmospheric response in the Pacific. These results imply an important modulation of global warming on the inter-basin connection.

This study examines marine boundary layer cloud regime transition during a cold air outbreak (CAO) over the Norwegian Sea, simulated by a global storm-resolving model (GSRM) known as the Simple Cloud-Resolving Energy Exascale Earth System Model Atmosphere Model (SCREAM). By selecting observational references based on a combination of large-scale conditions rather than strict time-matched comparisons, this study finds that SCREAM qualitatively captures the CAO cloud transition, including boundary layer growth, cloud mesoscale structure, and phase partitioning. SCREAM also accurately locates the greatest ice and liquid in the mesoscale updrafts, however, underestimates supercooled liquid water in cumulus clouds. The model evaluation approach adopted by this study takes advantages of the existing computational-expensive global simulations of GSRM and the available observations to understand model performance and can be applied to assessments of other cloud regimes in different regions. Such practice provides valuable guidance on the future effort to correct and improve biased model behaviors.