Explore the vast range of titles available.

Saltwater intrusion is a critical concern for coastal communities due to its impacts on fresh ecosystems and civil infrastructure. Declining recharge and rising sea level are the two dominant drivers of saltwater intrusion along the land‐ocean continuum, but there are currently no global estimates of future saltwater intrusion that synthesize these two spatially variable processes. Here, for the first time, we provide a novel assessment of global saltwater intrusion risk by integrating future recharge and sea level rise while considering the unique geology and topography of coastal regions. We show that nearly 77% of global coastal areas below 60° north will undergo saltwater intrusion by 2100, with different dominant drivers. Climate‐driven changes in subsurface water replenishment (recharge) is responsible for the high‐magnitude cases of saltwater intrusion, whereas sea level rise and coastline migration are responsible for the global pervasiveness of saltwater intrusion and have a greater effect on low‐lying areas. Plain Language Summary Coastal watersheds around the globe are facing perilous changes to their freshwater systems. Driven by climatic changes in recharge and sea level working in tandem, sea water encroaches into coastal groundwater aquifers and consequently salinizes fresh groundwater, in a process called saltwater intrusion. To assess the vulnerability of coastal watersheds to future saltwater intrusion, we applied projections of sea level and groundwater recharge to a global analytical modeling framework. Nearly 77% of the global coast is expected to undergo measurable salinization by the year 2100. Changes in recharge have a greater effect on the magnitude of salinization, whereas sea level rise drives the widespread extensiveness of salinization around the global coast. Our results highlight the variable pressures of climate change on coastal regions and have implications for prioritizing management solutions. Key Points First global analysis of future saltwater intrusion vulnerability responding to spatially variable recharge and sea level rise is provided Recharge drives the extreme cases of saltwater intrusion, while sea level rise is responsible for its global pervasiveness Nearly 77% of global coastal areas below 60° north will undergo saltwater intrusion by 2100

The subduction of seamounts and ridge features at convergent plate boundaries plays an important role in the deformation of the overriding plate and influences geochemical cycling and associated biological processes. Active serpentinization of forearc mantle and serpentinite mud volcanism on the Mariana forearc (between the trench and active volcanic arc) provides windows on subduction processes.  Here, we present (1) the first observation of an extensive exposure of an undeformed Cretaceous seamount currently being subducted at the Mariana Trench inner slope; (2) vertical deformation of the forearc region related to subduction of Pacific Plate seamounts and thickened crust; (3) recovered Ocean Drilling Program and International Ocean Discovery Program cores of serpentinite mudflows that confirm exhumation of various Pacific Plate lithologies, including subducted reef limestone; (4) petrologic, geochemical and paleontological data from the cores that show that Pacific Plate seamount exhumation covers greater spatial and temporal extents; (5) the inference that microbial communities associated with serpentinite mud volcanism may also be exhumed from the subducted plate seafloor and/or seamounts; and (6) the implications for effects of these processes with regard to evolution of life. This article is part of a discussion meeting issue ‘Serpentine in the Earth system’.

Severe convective storms are among the most dangerous weather phenomena and accurate forecasts mitigate their impacts. The recently released suite of AI‐based weather models produces medium‐range forecasts within seconds, with a skill similar to state‐of‐the‐art operational forecasts for variables on single levels. However, predicting severe thunderstorm environments requires accurate combinations of dynamic and thermodynamic variables and the vertical structure of the atmosphere. Advancing the assessment of AI‐models toward process‐based evaluations lays the foundation for hazard‐driven applications. We assess the forecast skill of the top‐performing AI‐models GraphCast, Pangu‐Weather and FourCastNet for convective parameters at lead‐times up to 10 days against reanalysis and ECMWF's operational numerical weather prediction model IFS. In a case study and seasonal analyses, we see the best performance by GraphCast and Pangu‐Weather: these models match or even exceed the performance of IFS for instability and shear. This opens opportunities for fast and inexpensive predictions of severe weather environments. Plain Language Summary Over the past year, several global AI‐based weather models were released and produce a similar quality of forecasts as traditional weather models. AI‐models are very fast and computationally cheap to produce forecasts. The evaluation of AI‐models has largely focused on single atmospheric variables at certain heights. To forecast specific phenomena, such as thunderstorms, a combination of variables must be accurate at multiple heights. Here we use the output of AI‐models to derive thunderstorm‐related ingredients. We compare 10‐day‐forecasts between the AI‐models GraphCast, Pangu‐Weather and FourCastNet and a state‐of‐the‐art traditional weather model while using a reanalysis data set as the reference. The example of a tornado outbreak in the southern United States shows that all models are capable of forecasting thunderstorm ingredients multiple days in advance. To obtain a robust assessment, we evaluate the entire thunderstorm season in 2020 in North America, Europe, Argentina, and Australia, where severe thunderstorms occur frequently. Two of the three AI‐models achieve similar or even better results than the traditional weather model while being much cheaper to operate computationally. Forecasting thunderstorm parameters directly, instead of calculating them afterward, is likely to produce even better results. This opens opportunities for rapid and accessible forecasts for severe thunderstorm phenomena. Key Points AI‐based global weather models produce forecasts with sufficient accuracy to derive instability and shear metrics skillfully The best AI‐based weather models are capable of competing with state‐of‐the‐art numerical weather predictions of instability and shear This is a major step toward computationally inexpensive and fast convective outlooks

Volcanic activity has been shown to affect Earth's climate in a myriad of ways. One such example is that eruptions proximate to surface ice will promote ice melting. In turn, the crustal unloading associated with melting an ice sheet affects the internal dynamics of the underlying magma plumbing system. Geochronologic data from the Andes over the last two glacial cycles suggest that glaciation and volcanism may interact via a positive feedback loop. At present, accurate sea‐level predictions hinge on our ability to forecast the stability of the West Antarctic Ice Sheet, and thus require consideration of two‐way subglacial volcano‐deglaciation processes. The West Antarctic Ice Sheet is particularly vulnerable to collapse, yet its position atop an active volcanic rift is seldom considered. Ice unloading deepens the zone of melting and alters the crustal stress field, impacting conditions for dike initiation, propagation, and arrest. However, the consequences for internal magma chamber dynamics and long‐term eruption behavior remain elusive. Given that unloading‐triggered volcanism in West Antarctica may contribute to the uncertainty of ice loss projections, we adapt a previously published thermomechanical magma chamber model and simulate a shrinking ice load through a prescribed lithostatic pressure decrease. We investigate the impacts of varying unloading scenarios on magma volatile partitioning and eruptive trajectory. Considering the removal of km‐thick ice sheets, we demonstrate that the rate of unloading influences the cumulative mass erupted and consequently the heat released into the ice. These findings provide fundamental insights into the complex volcano‐ice interactions in West Antarctica and other subglacial volcanic settings. Plain Language Summary In regions like West Antarctica, volcanic eruptions occur underneath ice sheets. When hot magma comes in contact with ice, it can accelerate the melting of the ice cover. Beyond this, as climate change causes ice sheets to shrink, the decreasing weight on a volcano may affect its likelihood of erupting. The effects of ice loss above volcanoes on the underlying volcanic activity are not well understood. We conducted computer simulations to explore how gradual ice loss affects magma stored in the Earth's crust. We find that volcanoes beneath shrinking ice sheets are sensitive to the rate at which the ice sheet shrinks. As the ice melts away, the reduced weight on the volcano allows the magma to expand, applying pressure upon the surrounding rock that may facilitate eruptions. Additionally, the reduced weight from the melting ice above also allows dissolved water and carbon dioxide to form gas bubbles, which causes pressure to build up in the magma chamber and may eventually trigger an eruption. Under these conditions, we find that the removal of an ice sheet above a volcano results in more abundant and larger eruptions, which may potentially hasten the melting of overlying ice through complex feedback mechanisms. Key Points During deglaciation, the evolution of a crustal magma chamber beneath kilometers of ice is sensitive to the rate at which ice is removed A critical rate of unloading can trigger additional eruption events Ice unloading expedites the onset of volatile exsolution, with consequences for magma chamber pressurization and eruption size

A promising approach to improve cloud parameterizations within climate models and thus climate projections is to use deep learning in combination with training data from storm‐resolving model (SRM) simulations. The ICOsahedral Non‐hydrostatic (ICON) modeling framework permits simulations ranging from numerical weather prediction to climate projections, making it an ideal target to develop neural network (NN) based parameterizations for sub‐grid scale processes. Within the ICON framework, we train NN based cloud cover parameterizations with coarse‐grained data based on realistic regional and global ICON SRM simulations. We set up three different types of NNs that differ in the degree of vertical locality they assume for diagnosing cloud cover from coarse‐grained atmospheric state variables. The NNs accurately estimate sub‐grid scale cloud cover from coarse‐grained data that has similar geographical characteristics as their training data. Additionally, globally trained NNs can reproduce sub‐grid scale cloud cover of the regional SRM simulation. Using the game‐theory based interpretability library SHapley Additive exPlanations, we identify an overemphasis on specific humidity and cloud ice as the reason why our column‐based NN cannot perfectly generalize from the global to the regional coarse‐grained SRM data. The interpretability tool also helps visualize similarities and differences in feature importance between regionally and globally trained column‐based NNs, and reveals a local relationship between their cloud cover predictions and the thermodynamic environment. Our results show the potential of deep learning to derive accurate yet interpretable cloud cover parameterizations from global SRMs, and suggest that neighborhood‐based models may be a good compromise between accuracy and generalizability. Plain Language Summary Climate models, such as the ICOsahedral Non‐hydrostatic climate model, operate on low‐resolution grids, making it computationally feasible to use them for climate projections. However, physical processes –especially those associated with clouds– that happen on a sub‐grid scale (inside a grid box) cannot be resolved, yet they are critical for the climate. In this study, we train neural networks that return the cloudy fraction of a grid box knowing only low‐resolution grid‐box averaged variables (such as temperature, pressure, etc.) as the climate model sees them. We find that the neural networks can reproduce the sub‐grid scale cloud fraction on data sets similar to the one they were trained on. The networks trained on global data also prove to be applicable on regional data coming from a model simulation with an entirely different setup. Since neural networks are often described as black boxes that are therefore difficult to trust, we peek inside the black box to reveal what input features the neural networks have learned to focus on and in what respect the networks differ. Overall, the neural networks prove to be accurate methods of reproducing sub‐grid scale cloudiness and could improve climate model projections when implemented in a climate model. Key Points Neural networks can accurately learn sub‐grid scale cloud cover from realistic regional and global storm‐resolving simulations Three neural network types account for different degrees of vertical locality and differentiate between cloud volume and cloud area fraction Using a game theory based library we find that the neural networks tend to learn local mappings and are able to explain model errors

Serpentinization-fueled systems in the cool, hydrated forearc mantle of subduction zones may provide an environment that supports deep chemolithoautotrophic life. Here, we examine serpentinite clasts expelled from mud volcanoes above the Izu–Bonin–Mariana subduction zone forearc (Pacific Ocean) that contain complex organic matter and nanosized Ni–Fe alloys. Using time-of-flight secondary ion mass spectrometry and Raman spectroscopy, we determined that the organic matter consists of a mixture of aliphatic and aromatic compounds and functional groups such as amides. Although an abiotic or subduction slab-derived fluid origin cannot be excluded, the similarities between the molecular signatures identified in the clasts and those of bacteria-derived biopolymers from other serpentinizing systems hint at the possibility of deep microbial life within the forearc. To test this hypothesis, we coupled the currently known temperature limit for life, 122 °C, with a heat conduction model that predicts a potential depth limit for life within the forearc at ∼10,000 m below the seafloor. This is deeper than the 122 °C isotherm in known oceanic serpentinizing regions and an order of magnitude deeper than the downhole temperature at the serpentinized Atlantis Massif oceanic core complex, Mid-Atlantic Ridge. We suggest that the organic-rich serpentinites may be indicators for microbial life deep within or below the mud volcano. Thus, the hydrated forearc mantle may represent one of Earth’s largest hidden microbial ecosystems. These types of protected ecosystems may have allowed the deep biosphere to thrive, despite violent phases during Earth’s history such as the late heavy bombardment and global mass extinctions.

Fresh submarine groundwater discharge (FSGD) can deliver significant fluxes of water and solutes from land to sea. In the Arctic, which accounts for ∼34% of coastlines globally, direct observations and knowledge of FSGD are scarce. Through integration of observations and process‐based models, we found that regardless of ice‐bonded permafrost depth at the shore, summer SGD flow dynamics along portions of the Beaufort Sea coast of Alaska are similar to those in lower latitudes. Calculated summer FSGD fluxes in the Arctic are generally higher relative to low latitudes. The FSGD organic carbon and nitrogen fluxes are likely larger than summer riverine input. The FSGD also has very high CO2 making it a potentially significant source of inorganic carbon. Thus, the biogeochemistry of Arctic coastal waters is potentially influenced by groundwater inputs during summer. These water and solute fluxes will likely increase as coastal permafrost across the Arctic thaws. Plain Language Summary Groundwater flows from land to sea, transporting freshwater, organic matter, nutrients, and other solutes that impact coastal ecosystems. However, along coasts of the rapidly‐warming Arctic, there is limited knowledge regarding how much fresh groundwater enters the ocean. Using field observations and numerical models, we show that groundwater flowing from tundra in northern coastal Alaska carries large amounts of freshwater, organic matter, and carbon dioxide to the Arctic lagoons during summer. These inputs are likely significant to coastal biogeochemical cycling and marine food webs. Groundwater discharge and the associated transport of dissolved materials are expected to increase due to longer periods of above‐zero temperatures that thaw frozen soils below the tundra. Key Points Summer fresh submarine groundwater discharge (FSGD) to the Alaskan Beaufort Sea is only 3%–7% of rivers but carries as much organic matter Summer FSGD delivers a median of 116 (interquartile range: 35–405) and 6 (2–21) kg/d per km dissolved organic carbon and nitrogen Fresh groundwater at the beach of Simpson Lagoon (SL) has a median PCO2 of ∼33,000 μatm implying substantial CO2 flux

Global warming is shifting the thermal dynamics of lakes, with resulting climatic variability heavily affecting their mixing dynamics. We present a dual ensemble workflow coupling climate models with lake models. We used a large set of simulations across multiple domains, multi‐scenario, and multi GCM‐ RCM combinations from CORDEX data. We forced a set of multiple hydrodynamic lake models by these multiple climate simulations to explore climate change impacts on lakes. We also quantified the contributions from the different models to the overall uncertainty. We employed this workflow to investigate the effects of climate change on Lake Sevan (Armenia). We predicted for the end of the 21st century, under RCP 8.5, a sharp increase in surface temperature (4.3±0.7K)$(4.3\\pm 0.7\\,\\mathrm{K})$and substantial bottom warming (1.7±0.7K)$(1.7\\pm 0.7\\,\\mathrm{K})$ , longer stratification periods (+55 days) and disappearance of ice cover leading to a shift in mixing regime. Increased insufficient cooling during warmer winters points to the vulnerability of Lake Sevan to climate change. Our workflow leverages the strengths of multiple models at several levels of the model chain to provide a more robust projection and at the same time a better uncertainty estimate that accounts for the contributions of the different model levels to overall uncertainty. Although for specific variables, for example, summer bottom temperature, single lake models may perform better, the full ensemble provides a robust estimate of thermal dynamics that has a high transferability so that our workflow can be a blueprint for climate impact studies in other systems. Plain Language Summary Lakes are threatened by climate change because of effects related to the increasing temperature, long stratification, and ice disappearance. One of the best tools to predict these effects on lakes is numerical modeling of lakes that benefit from climate modeling. Climate modeling is normally done globally or in the so‐called general circulation model (GCM) or more detailed simulations on regional levels (RCM) like the CORDEX data set. In this study, we used the CORDEX data, which employed several climate models from several regions (domains) for several climatic scenarios (emissions scenarios) to force multiple lake models. This approach gave us an extensive prediction about various possible outputs. We applied this approach to Lake Sevan (Armenia), a large mountain lake. Our study predicted for the worst‐case scenario, an increase of the surface temperature by almost 4.3 K by the end of the 21st century, 1.75 K for bottom temperature, a total disappearance of ice cover, and about 55 extra days of stratification, showing its vulnerability for climate change. This optimized workflow uses the strength of a wide variety of models on the climate and lake levels to better understand the impact of climate change and quantify the sources of uncertainty in the workflow. Key Points Dual multi‐model ensemble of climate data and lake models is used for robust projections of climate change impacts Variance decomposition effectively identified the sources of uncertainty and contributions of different models to the overall uncertainty Significant warming, longer stratification periods, and loss of ice cover are predicted for Lake Sevan by the end of the 21st century

Anthropogenic aerosols (AER) and greenhouse gases (GHG)—the leading drivers of the forced historical change—produce different large‐scale climate response patterns, with correlations trending from negative to positive over the past century. To understand what caused the time‐evolving comparison between GHG and AER response patterns, we apply a low‐frequency component analysis to historical surface ocean changes from CESM1 single‐forcing large‐ensemble simulations. While GHG response is characterized by its first leading mode, AER response consists of two distinct modes. The first one, featuring long‐term global AER increase and global cooling, opposes GHG response patterns up to the mid‐twentieth century. The second one, featuring multidecadal variations in AER distributions and interhemispheric asymmetric surface ocean changes, appears to reinforce the GHG warming effect over recent decades. AER thus can have both competing and synergistic effects with GHG as their emissions change temporally and spatially. Plain Language Summary Anthropogenically forced climate change over the past century has been mainly caused by two types of emissions: greenhouse gases (GHG) and aerosols (AER). In general, sulfate aerosols from industrial sources can reflect shortwave radiation to yield a cooling effect opposite to the GHG warming effect. However, model simulations isolating GHG and AER forcings show that the large‐scale climate effect of AER does not always dampen the GHG effect. Instead, over recent decades, AER have produced surface ocean response patterns more like the GHG response. Using a novel low‐frequency statistical decomposion, we find that aerosols have driven two distinct modes of climate change patterns over the historical period. The first mode is associated with global aerosol increase, resulting in global‐wide cooling damping the GHG‐induced warming. The second mode is associated with the shift in aerosol emissions from north America/western Europe to southeast Asia, which drives regional changes enhancing the GHG effect. Our results highlight the importance of considering the temporal and spatial evolutions of AER emissions in assessing GHG and AER climate effects and attributing historical anthropogenic climate changes to GHG and AER forcings. Key Points Over the past century, GHG forced response is characterized by a single dominant mode while AER response consists of two distinct modes Monotonic global aerosol increases, mainly from Southeast Asia emissions, produce a global aerosol cooling mode opposing greenhouse warming Important in recent decades, geographic redistribution of AER emissions produces a second aerosol mode that reinforces greenhouse warming

Sea level rise (SLR) is a long-lasting consequence of climate change because global anthropogenic warming takes centuries to millennia to equilibrate for the deep ocean and ice sheets. SLR projections based on climate models support policy analysis, risk assessment and adaptation planning today, despite their large uncertainties. The central range of the SLR distribution is estimated by process-based models. However, risk-averse practitioners often require information about plausible future conditions that lie in the tails of the SLR distribution, which are poorly defined by existing models. Here, a community effort combining scientists and practitioners builds on a framework of discussing physical evidence to quantify high-end global SLR for practitioners. The approach is complementary to the IPCC AR6 report and provides further physically plausible high-end scenarios. High-end estimates for the different SLR components are developed for two climate scenarios at two timescales. For global warming of +2°C in 2100 (RCP2.6/SSP1-2.6) relative to pre-industrial values our high-end global SLR estimates are up to 0.9 m in 2100 and 2.5 m in 2300. Similarly, for a (RCP8.5/SSP5-8.5), we estimate up to 1.6 m in 2100 and up to 10.4 m in 2300. The large and growing differences between the scenarios beyond 2100 emphasize the long-term benefits of mitigation. However, even a modest 2°C warming may cause multi-meter SLR on centennial time scales with profound consequences for coastal areas. Earlier high-end assessments focused on instability mechanisms in Antarctica, while here we emphasize the importance of the timing of ice shelf collapse around Antarctica. This is highly uncertain due to low understanding of the driving processes. Hence both process understanding and emission scenario control high-end SLR.