Search NASA⌕ Search

SEARCH · Search NASA

Results for “Arctic, Earth System Modeling”

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

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

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Uncertain Spatial Pattern of Future Land Use and Land Cover Change and Its Impacts on Terrestrial Carbon Cycle Over the Arctic–Boreal Region of North America

Land use and land cover change (LULCC) represents a key process of human-Earth system interaction and has profound impacts on terrestrial ecosystem carbon cycling. As a key input for ecosystem models, future gridded LULCC data is typically spatially downscaled from regional LULCC projections by integrated assessment models, such as the Global Change Analysis Model (GCAM). The uncertainty associated with the different spatial downscaling methods and its impacts on the subsequent model projections have been historically ignored and rarely examined. This study investigated this problem using two representative spatial downscaling methods and focused on their impacts on the carbon cycle over the Arctic-Boreal Vulnerability Experiment (ABoVE) domain, where extensive LULCC is expected. Specifically, we used the Future Land Use Simulation model (FLUS) and the Demeter model to generate 0.25° gridded LULCC data (i.e., LULCC FLUS and LULCC Demeter , respectively) with the same input of regional LULCC projections from GCAM, under both the low (i.e., SSP126) and high (i.e., SSP585) greenhouse gas emission scenarios. The two sets of downscaled LULCC were used to drive the Community Land Model version 5 and prognostically simulate the terrestrial carbon cycle dynamics over the 21st century. The results suggest large spatial-temporal differences between LULCC FLUS and LULCC Demeter , and the spatial distributions of the needleleaf evergreen boreal tree, broadleaf deciduous boreal tree, broadleaf deciduous boreal shrub, and C3 arctic grass are particularly different under both SSP126 and SSP585. Additionally, the spatiotemporal differences are larger under SSP126 than SSP585, due to more intensive LULCC under SSP126 than SSP585 from GCAM projection. The differences in LULCC further lead to large discrepancies in the spatial patterns of projected gross primary productivity, ecosystem respiration, and net ecosystem exchange, which represent more than 79% of the contributions of future LULCC in 2100. Additionally, the difference in carbon flux under SSP126 is generally larger than those under SSP585. This study highlights the importance of considering the uncertainties induced by the spatial downscaling process in future LULCC projections and carbon cycle simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Nitrogen fixing shrubs advance the pace of tall-shrub expansion in low-Arctic tundra

Abstract Tall deciduous shrubs are critically important to carbon and nutrient cycling in high-latitude ecosystems. As Arctic regions warm, shrubs expand heterogeneously across their ranges, including within unburned terrain experiencing isometric gradients of warming. To constrain the effects of widespread shrub expansion in terrestrial and Earth System Models, improved knowledge of local-to-regional scale patterns, rates, and controls on decadal shrub expansion is required. Using fine-scale remote sensing, we modeled the drivers of patch-scale tall-shrub expansion over 68 years across the central Seward Peninsula of Alaska. Models show the heterogeneous patterns of tall-shrub expansion are not only predictable but have an upper limit defined by permafrost, climate, and edaphic gradients, two-thirds of which have yet to be colonized. These observations suggest that increased nitrogen inputs from nitrogen-fixing alders contributed to a positive feedback that advanced overall tall-shrub expansion. These findings will be useful for constraining and projecting vegetation-climate feedbacks in the Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Primary and secondary ice production: interactions and their relative importance

A discrepancy of up to 5 orders of magnitude between ice crystal and ice nucleating particle (INP) number concentrations was found in the measurements, indicating the potentially important role of secondary ice production (SIP) in the clouds. However, the interactions between primary and SIP processes and their relative importance remain unexplored. In this study, we implemented five different ice nucleation schemes as well as physical representations of SIP processes (i.e., droplet shattering during rain freezing, ice-ice collisional break-up, and rime splintering) in the Community Earth System Model version 2 (CESM2). We ran CESM2 in the single column mode for model comparisons with the DOE Atmospheric Radiation Measurement (ARM) Mixed-Phase Arctic Cloud Experiment (M-PACE) observations. We found that the model experiments with aerosol-aware ice nucleation schemes and SIP processes yield the best simulation results for the M-PACE single-layer mixed-phase clouds. We further investigated the relative importance of ice nucleation and SIP to ice number and cloud phase as well as interactions between ice nucleation and SIP in the M-PACE single-layer mixed-phase clouds. Our results show that SIP contributes 80 % to the total ice formation and transforms ~30 % of pure liquid-phase clouds simulated in the model experiments without considering SIP into mixed-phase clouds. The SIP is not only a result of ice crystals produced from ice nucleation, but also competes with the ice nucleation by reducing the number concentrations of cloud droplets and cloud-borne dust INPs. Conversely, strong ice nucleation also suppresses SIP by glaciating mixed-phase clouds and thereby reducing the amount of precipitation particles (rain and graupel).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Causes and evolution of winter polynyas north of Greenland

During the 42-year period (1979–2020) of satellite measurements, four major winter (December–March) polynyas have been observed north of Greenland: one in December 1986 and three in the last decade, i.e., February of 2011, 2017, and 2018. The 2018 polynya was unparalleled in its magnitude and duration compared to the three previous events. Given the apparent recent increase in the occurrence of these extreme events, this study aims to examine their evolution and causality, in terms of forced versus natural variability. The limited weather station and remotely sensed sea ice data are analyzed combining with output from the fully coupled Regional Arctic System Model (RASM), including one hindcast and two ensemble simulations. We found that neither the accompanying anomalous warm surface air intrusion nor the ocean below had an impact (i.e., no significant ice melting) on the evolution of the observed winter open-water episodes in the region. Instead, the extreme atmospheric wind forcing resulted in greater sea ice deformation and transport offshore, accounting for the majority of sea ice loss in all four polynyas. Our analysis suggests that strong southerly winds (i.e., northward wind with speeds greater than 10 m s –1 ) blowing persistently over the study region for at least 2 d or more were required over the study region to mechanically redistribute some of the thickest Arctic sea ice out of the region and thus to create open-water areas (i.e., a latent heat polynya). To assess the role of internal variability versus external forcing of such events, we carried out and examined results from the two RASM ensembles dynamically downscaled with output from the Community Earth System Model (CESM) Decadal Prediction Large Ensemble (DPLE) simulations. Out of 100 winters in each of the two ensembles (initialized 30 years apart: one in December 1985 and another in December 2015), 17 and 16 winter polynyas were produced north of Greenland, respectively. The frequency of polynya occurrence had no apparent sensitivity to the initial sea ice thickness in the study area pointing to internal variability of atmospheric forcing as a dominant cause of winter polynyas north of Greenland. We assert that dynamical downscaling using a high-resolution regional climate model offers a robust tool for process-level examination in space and time, synthesis with limited observations, and probabilistic forecasts of Arctic events, such as the ones being investigated here and elsewhere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effects of marine organic aerosols as sources of immersion-mode ice-nucleating particles on high-latitude mixed-phase clouds

Mixed-phase clouds are frequently observed in high-latitude regions and have important impacts on the surface energy budget and regional climate. Marine organic aerosol (MOA), a natural source of aerosol emitted over ~70 % of Earth's surface, may significantly modify the properties and radiative forcing of mixed-phase clouds. However, the relative importance of MOA as a source of ice-nucleating particles (INPs) in comparison to mineral dust, and MOA's effects as cloud condensation nuclei (CCN) and INPs on mixed-phase clouds are still open questions. In this study, we implement MOA as a new aerosol species into the Community Atmosphere Model version 6 (CAM6), the atmosphere component of the Community Earth System Model version 2 (CESM2), and allow the treatment of aerosol–cloud interactions of MOA via droplet activation and ice nucleation. CAM6 reproduces observed seasonal cycles of marine organic matter at Mace Head and Amsterdam Island when the MOA fraction of sea spray aerosol in the model is assumed to depend on sea spray biology but fails when this fraction is assumed to be constant. Model results indicate that marine INPs dominate primary ice nucleation below 400 hPa over the Southern Ocean and Arctic boundary layer, while dust INPs are more abundant elsewhere. By acting as CCN, MOA exerts a shortwave cloud forcing change of -2.78 W m-2 over the Southern Ocean in the austral summer. By acting as INPs, MOA enhances the longwave cloud forcing by 0.35 W m-2 over the Southern Ocean in the austral winter. The annual global mean net cloud forcing changes due to CCN and INPs of MOA are -0.35 and 0.016 W m-2, respectively. These findings highlight the vital importance for Earth system models to consider MOA as an important aerosol species for the interactions of biogeochemistry, hydrological cycle, and climate change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Vertical profiles of light absorption and scattering associated with black carbon particle fractions in the springtime Arctic above 79° N

Despite the potential importance of black carbon (BC) for radiative forcing of the Arctic atmosphere, vertically resolved measurements of the particle light scattering coefficient (σ sp ) and light absorption coefficient (σ ap ) in the springtime Arctic atmosphere are infrequent, especially measurements at latitudes at or above 80° N. Here, relationships among vertically distributed aerosol optical properties (σap, σsp and single scattering albedo or SSA), particle microphysics and particle chemistry are examined for a region of the Canadian archipelago between 79.9 and 83.4° N from near the surface to 50 hPa. Airborne data collected during April 2015 are combined with ground-based observations from the observatory at Alert, Nunavut and simulations from the Goddard Earth Observing System (GEOS) model, GEOS-Chem, coupled with the TwO-Moment Aerosol Sectional (TOMAS) model (collectively GEOS-Chem–TOMAS; Kodros et al., 2018) to further our knowledge of the effects of BC on light absorption in the Arctic troposphere. The results are constrained for σsp less than 15 Mm -1 , which represent 98% of the observed σsp, because the single scattering albedo (SSA) has a tendency to be lower at lower σsp, resulting in a larger relative contribution to Arctic warming. At 18.4 m 2 g -1 , the average BC mass absorption coefficient (MAC) from the combined airborne and Alert observations is substantially higher than the two averaged modelled MAC values (13.6 and 9.1 m2 g -1 ) for two different internal mixing assumptions, the latter of which is based on previous observations. The higher observed MAC value may be explained by an underestimation of BC, the presence of small amounts of dust and/or possible differences in BC microphysics and morphologies between the observations and model. In comparing the observations and simulations, we present σap and SSA, as measured, and σ ap 2 and the corresponding SSA to encompass the lower modelled MAC that is more consistent with accepted MAC values. Median values of the measured σap, rBC and the organic component of particles all increase by a factor of 1.8±0.1, going from near-surface to 750 hPa, and values higher than the surface persist to 600 hPa. Modelled BC, organics and σap agree with the near-surface measurements but do not reproduce the higher values observed between 900 and 600 hPa. The differences between modelled and observed optical properties follow the same trend as the differences between the modelled and observed concentrations of the carbonaceous components (black and organic). Model-observation discrepancies may be mostly due to the modelled ejection of biomass burning particles only into the boundary layer at the sources. For the assumption of the observed MAC value, the SSA range between 0.88 and 0.94, which is significantly lower than other recent estimates for the Arctic, in part reflecting the constraint of σsp<15 Mm -1 . The large uncertainties in measuring optical properties and BC, and the large differences between measured and modelled values here and in the literature, argue for improved measurements of BC and light absorption by BC and more vertical profiles of aerosol chemistry, microphysics and other optical properties in the Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

PAVC: The foundation for a Pan-Arctic Vegetation Cover database

Field-measured Arctic vegetation cover data is essential for creating accurate, high-quality vegetation structure and composition maps. Extrapolating field data into high-resolution cover maps provides detailed, function-specific information for use in Earth System Models, vegetation classifications, and monitoring vegetation change over time and space. However, field campaigns that collect plant cover vary substantially in scope, method, and purpose, which makes them difficult to unify across data stores, and they are often not designed to meet remote sensing needs. In this work, we synthesized and harmonized field-based fractional cover data from various data stores to create a high-quality, consistent repository schema for remote sensing-based vegetation cover mapping applications. We developed a reproducible workflow for synthesizing visual estimate and point-intercept fractional cover data. The resultant Pan-Arctic Vegetation Cover (PAVC) database contains synthesized fractional cover at both the species and plant functional type levels. The latter includes absolute foliar cover for deciduous shrubs and trees, evergreen shrubs and trees, forbs, graminoids, lichen, bryophytes, and “other” vegetation, as well as absolute cover for litter and top cover for water and bare ground.

Steckler, Morgan R. [Oak Ridge National Laboratory↗
DOE OSTI2022

Benthic Kinetic Energy in an Earth System Model and Observations

The benthic layer of the ocean is closely correlated with every upper layer above it. Close analysis of benthic layer stress can yield an accurate depiction of the state of the ocean, as well as the atmosphere above it. Through the use of the Department of Energy’s Energy Exascale Earth System Model (E3SM), we took a look at how these features changed over time in a simulation spanning the historical period from 1958 to 2018 in the Arctic Ocean north of 60° N. Making use of E3SM’s regional refinement capability to place higher resolution in regions of interest, the E3SM-Arctic configuration has ~10 km resolution in the ocean and sea ice components in this region, coarsening up to ~60 km outside of the Arctic (Veneziani et al. 2022). In particular, we narrowly focused on the benthic kinetic energy levels, and how trends changed over time. To validate the fidelity of the E3SM model, we compared model current speed outputs over a 369-day period in the Bering Strait to actual recorded observations during this time period. Our analysis of the benthic kinetic energy mean, standard deviation, and trend throughout this time period can be used to characterize the effects of environmental impacts in the Arctic Ocean.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Data Files for Runoff Evaluation in an Earth System Land Model for Permafrost Regions

Modeling of hydrological runoff is essential for accurately capturing spatiotemporal feedbacks within the land–atmosphere system, particularly in sensitive regions such as permafrost landscapes. However, substantial uncertainties persist in the terrestrial runoff parameterization schemes used in Earth system and land surface models. This is particularly true in permafrost regions, where landscape heterogeneity is high and reliable observational data are scarce.This data set includes all files that were produced and applied in the paper Runoff Evaluation in an Earth System Land Model for Permafrost Regions [Xiang et al. in review]. The paper is in review as of July 1 2025 in Geoscientific Model Development (GMD). In this study, we evaluate the performance of runoff parameterization schemes in the Energy Exascale Earth System Model (E3SM) land model (ELM). Our proposed framework leverages simulation results from the Advanced Terrestrial Simulator (ATS), which is a physics-rich integrated surface/subsurface hydrologic model that has been successfully evaluated previously in Arctic tundra regions. We used ATS to simulate runoff from 22 representative hillslopes in the Sagavanirktok River basin, located on the North Slope of Alaska, then compared the output with ELM’s parameterized representation of total runoff. This dataset contains 2 figure image files (*.png, *jpg) that describe the study site and methods, as well as folders (Figure*.zip) that contain the associated data files (*.csv, *.dat) and python code notebooks (*.ipynb) for figures 3-7 in the paper. Jupyter notebook (*.ipynb) files that produce the figure files using the associated data files will run within a python environment configured with Jupyter Lab or Notebook packages.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The MOSAiC Distributed Network: Observing the coupled Arctic system with multidisciplinary, coordinated platforms

Central Arctic properties and processes are important to the regional and global coupled climate system. The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) Distributed Network (DN) of autonomous ice-tethered systems aimed to bridge gaps in our understanding of temporal and spatial scales, in particular with respect to the resolution of Earth system models. By characterizing variability around local measurements made at a Central Observatory, the DN covers both the coupled system interactions involving the ocean-ice-atmosphere interfaces as well as three-dimensional processes in the ocean, sea ice, and atmosphere. The more than 200 autonomous instruments (“buoys”) were of varying complexity and set up at different sites mostly within 50 km of the Central Observatory. During an exemplary midwinter month, the DN observations captured the spatial variability of atmospheric processes on sub-monthly time scales, but less so for monthly means. They show significant variability in snow depth and ice thickness, and provide a temporally and spatially resolved characterization of ice motion and deformation, showing coherency at the DN scale but less at smaller spatial scales. Ocean data show the background gradient across the DN as well as spatially dependent time variability due to local mixed layer sub-mesoscale and mesoscale processes, influenced by a variable ice cover. The second case (May–June 2020) illustrates the utility of the DN during the absence of manually obtained data by providing continuity of physical and biological observations during this key transitional period. We show examples of synergies between the extensive MOSAiC remote sensing observations and numerical modeling, such as estimating the skill of ice drift forecasts and evaluating coupled system modeling. The MOSAiC DN has been proven to enable analysis of local to mesoscale processes in the coupled atmosphere-ice-ocean system and has the potential to improve model parameterizations of important, unresolved processes in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Simulating physical and biogeochemical dynamics of the Arctic Ocean [Slides]

Ocean physics (currents, mixing, air-sea transfer of heat and gases) and biogeochemistry (carbon and oxygen cycles, ecosystem dynamics, nutrient distributions) are tightly coupled. We used the DOE-developed Energy Exascale Earth System Model (E3SM) to simulate physical and biogeochemical dynamics of the global ocean, with an emphasis on the western North Pacific. In particular, we used the grid refinement capability of E3SM’s ocean component (LANL-developed MPASOcean) to quantify the effect of model resolution on the carbon and oxygen cycles using a pair of century-long simulations that differ only in the western North Pacific grid. Regional grid refinement has a noticeable impact on biogeochemistry. Distribution of important nutrients (such as nitrate) changes the amount of photosynthesis by phytoplankton, which in turn affects air-sea transfer of CO2 and the local carbon cycle.

58 GEOSCIENCES↗
DOE OSTI2025

Discrete global grid system-based flow routing datasets in the Amazon and Yukon basins

Abstract. Discrete global grid systems (DGGS) are emerging spatial data structures widely used to organize geospatial datasets across scales. While DGGS have found applications in various scientific disciplines, including atmospheric science and ecology, their integration into physically based hydrological models and Earth system models (ESMs) has been hindered by the lack of flow routing datasets based on DGGS. In response to this gap, this study pioneers the development of new flow routing datasets using icosahedral Snyder equal-area (ISEA) DGGS and a novel mesh-independent flow direction model. We present flow routing datasets for two large basins, the tropical Amazon River basin and the Arctic Yukon River basin. These datasets (1) facilitate the adoption of DGGS for hydrological models and (2) provide flow routing inputs for evaluation of DGGS-based flow routing in the Amazon and Yukon river basins. The data are available at https://doi.org/10.5281/zenodo.8377765 (Liao, 2023).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Runoff evaluation in an Earth System Land Model for permafrost regions in Alaska

Modeling of hydrological runoff is essential for accurately capturing spatiotemporal feedbacks within the land–atmosphere system, particularly in sensitive regions such as permafrost landscapes. However, substantial uncertainties persist in the terrestrial runoff parameterization schemes used in Earth system and land surface models. This is particularly true in permafrost regions, where landscape heterogeneity is high and reliable observational data are scarce. In this study, we evaluate the performance of runoff parameterization schemes in the Energy Exascale Earth System Model (E3SM) land model (ELM). Our proposed framework leverages simulation results from the Advanced Terrestrial Simulator (ATS), which is a physics-based integrated surface/subsurface hydrologic model that has been successfully evaluated previously in Arctic tundra regions. We used ATS to simulate runoff from 22 representative hillslopes in the Sagavanirktok River basin, located on the North Slope of Alaska, then compared the output with ELM's parameterized representation of total runoff. Results show that (1) ELM's total runoff was the same order of magnitude as the ATS simulations, and both models were similarly variable over time; (2) minor adjustments to coefficients in ELM's runoff parameterization improved the match between the ATS simulation and ELM's parameterized representation of annual and seasonal total runoff; (3) overall, runoff responses in ATS and ELM are more similar in flat hillslope environments compared to steep hillslopes; and (4) shallower active layer thicknesses and higher precipitation simulations resulted in lower correlations between the two models due to greater total runoff. By incorporating the optimized runoff coefficients from the Sagavanirktok River basin into ELM, the simulated total runoff better matched the streamflow observations at a small watershed located on the Seward Peninsula of Alaska. Our findings revealed important insights into the effectiveness of runoff parameterizations in land surface models and pathways for improving runoff coefficients in typical Arctic regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

How does humidity data impact land surface modeling of hydrothermal regimes at a permafrost site in Utqia$\mathrm{\dot{g}}$vik, Alaska?

Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. Here, we found that higher values of RH (up to 40%) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30% maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physically-based predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Climate Change Projection in the Twenty-First Century Simulated by NIMS-KMA CMIP6 Model Based on New GHGs Concentration Pathways

Abstract The National Institute of Meteorological Sciences-Korea Meteorological Administration (NIMS-KMA) has participated in the Coupled Model Inter-comparison Project (CMIP) and provided long-term simulations using the coupled climate model. The NIMS-KMA produces new future projections using the ensemble mean of KMA Advanced Community Earth system model (K-ACE) and UK Earth System Model version1 (UKESM1) simulations to provide scientific information of future climate changes. In this study, we analyze four experiments those conducted following the new shared socioeconomic pathway (SSP) based scenarios to examine projected climate change in the twenty-first century. Present day (PD) simulations show high performance skill in both climate mean and variability, which provide a reliability of the climate models and reduces the uncertainty in response to future forcing. In future projections, global temperature increases from 1.92 °C to 5.20 °C relative to the PD level (1995–2014). Global mean precipitation increases from 5.1% to 10.1% and sea ice extent decreases from 19% to 62% in the Arctic and from 18% to 54% in the Antarctic. In addition, climate changes are accelerating toward the late twenty-first century. Our CMIP6 simulations are released to the public through the Earth System Grid Federation (ESGF) international data sharing portal and are used to support the establishment of the national adaptation plan for climate change in South Korea.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Improved Understanding of Multicentury Greenland Ice Sheet Response to Strong Warming in the Coupled CESM2‐CISM2 With Regional Grid Refinement

The simulation of ice sheet‐climate interactions, such as surface mass balance fluxes, is sensitive to model grid resolution. Here we simulate the multi‐century evolution of the Greenland Ice Sheet (GrIS) and its interaction with the climate using the Community Earth System Model version 2.2 (CESM2.2) including an interactive GrIS component (the Community Ice Sheet Model v2.1 [CISM2.1]) under an idealized warming scenario (atmospheric CO 2 increases by 1% yr -1 until quadrupling the pre‐industrial level and then is held fixed). A variable‐resolution (VR) grid with 1/4° regional refinement over the broader Arctic and 1° resolution elsewhere is applied to the atmosphere and land components, and the results are compared with conventional 1° lat‐lon grid simulations to investigate the impact of grid refinement. Compared with the 1° runs, the VR run features a slower rate of surface melt, especially over the western and northern GrIS, where the ice surface slopes gently toward the periphery. This difference pattern originates primarily from higher snow albedo and, thus, weaker albedo feedback in the VR run. The VR grid better captures the CISM ice sheet topography by reducing elevation discrepancies between CAM and CISM and is, therefore, less reliant on the downscaling algorithm, which is known to underestimate albedo gradients. The sea level rise contribution from the GrIS in the VR run is 53 mm by year 150 and 831 mm by year 350, approximately 40% and 20% less than that of the 1° runs, respectively.

Earth System Model↗
DOE OSTI2025

A new biogeochemical modelling framework (FLaMe-v1.0) for lake methane emissions on the regional scale: development and application to the European domain

This study presents a new physical-biogeochemical modelling framework for simulating lake methane (CH 4 ) emissions at regional scales. The new model, FLaMe-v1.0 (Fluxes of Lake Methane), rests on an innovative, computationally efficient lake clustering approach that enables the simulation of CH 4 emissions across a large number of lakes. Building on the Canadian Small Lake Model (CSLM) that simulates the lake physics, we develop a suite of biogeochemical modules to simulate transient dynamics of organic Carbon (C), Oxygen (O 2 ), and CH 4 . We first test the performance of FLaMe-v1.0 by analyzing physical and biogeochemical processes in two theoretical lakes with characteristics that can be considered representative for many lakes (an oligotrophic, deep lake driven by cold climate versus a eutrophic, shallow lake driven by warm climate). Next, we evaluate the model by comparing simulated and observed timeseries of CH 4 emissions in four well-surveyed lakes. We then apply FLaMe-v1.0 at the European scale to evaluate simulated diffusive and ebullitive lake CH 4 fluxes against in-situ measurements in both boreal and central European regions. Finally, we provide a first assessment of the spatio-temporal variability in CH 4 emissions from European lakes with a surface area comprised between 0.1–1000 km 2 (n= 108 407, total area = 1.33 × 105 km 2 ), indicating a total emission of 0.97 ± 0.23 Tg CH 4 yr −1 , with the uncertainty constrained by combining FLaMe-v1.0 and machine learning techniques. Moreover, 30 % and 70 % of these CH 4 emissions are through diffusive and ebullitive pathways, respectively. Annually averaged CH 4 emission rates per unit lake area during 2010–2016 have a South-to-North decreasing gradient, resulting in a mean over the European domain as 7.39 g CH 4 m −2 yr −1 . Our simulations reveal a strong seasonality (with ice-blocking effects accounted for) in European lake CH 4 emissions, with nearly ten times higher emissions during late summer than during winter. This pronounced seasonal variation highlights the importance of accounting for the sub-annual variability in CH 4 emissions to accurately constrain regional CH 4 budgets. In the future, FLaMe-v1.0 could be embedded into Earth System Models to investigate the feedback between climate warming and global lake CH 4 emissions.

Maisonnier, Manon [Free Univ. of Brussels (Belgium↗
DOE OSTI2024

Opportunities for Earth Observation to Inform Risk Management for Ocean Tipping Points

Abstract As climate change continues, the likelihood of passing critical thresholds or tipping points increases. Hence, there is a need to advance the science for detecting such thresholds. In this paper, we assess the needs and opportunities for Earth Observation (EO, here understood to refer to satellite observations) to inform society in responding to the risks associated with ten potential large-scale ocean tipping elements: Atlantic Meridional Overturning Circulation; Atlantic Subpolar Gyre; Beaufort Gyre; Arctic halocline; Kuroshio Large Meander; deoxygenation; phytoplankton; zooplankton; higher level ecosystems (including fisheries); and marine biodiversity. We review current scientific understanding and identify specific EO and related modelling needs for each of these tipping elements. We draw out some generic points that apply across several of the elements. These common points include the importance of maintaining long-term, consistent time series; the need to combine EO data consistently with in situ data types (including subsurface), for example through data assimilation; and the need to reduce or work with current mismatches in resolution (in both directions) between climate models and EO datasets. Our analysis shows that developing EO, modelling and prediction systems together, with understanding of the strengths and limitations of each, provides many promising paths towards monitoring and early warning systems for tipping, and towards the development of the next generation of climate models.

Wood, Richard A. (ORCID:0000000239609513)↗