Search NASA⌕ Search

SEARCH · Search NASA

Results for “Soil Temperature”

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 163 records · Page 9

Improved ELMv1-ECA simulations of zero-curtain periods and cold-season CH 4 and CO 2 emissions at Alaskan Arctic tundra sites

Field measurements have shown that cold-season methane (CH 4 ) and carbon dioxide (CO 2 ) emissions contribute a substantial portion to the annual net carbon emissions in permafrost regions. However, most earth system land models do not accurately reproduce cold-season CH 4 and CO 2 emissions, especially over the shoulder (i.e., thawing and freezing) seasons. Here we use the Energy Exascale Earth System Model (E3SM) land model version 1 (ELMv1-ECA) to tackle this challenge and fill the knowledge gap of how cold-season CH 4 and CO 2 emissions contribute to the annual totals at Alaska Arctic tundra sites. Specifically, we improved the ELMv1-ECA soil water phase-change scheme, environmental controls on microbial activity, and the methane module. Results demonstrate that both soil temperature and the duration of zero-curtain periods (i.e., the fall period when soil temperatures linger around 0°C) simulated by the updated ELMv1-ECA were greatly improved; e.g., the mean absolute error (MAE) in zero-curtain durations at 12 cm depth was reduced by 62 % on average. Furthermore, the MAEs of simulated cold-season carbon emissions at three tundra sites were improved by 72 % and 70 % on average for CH 4 and CO 2 , respectively. Overall, CH 4 emitted during the early cold season (September and October), which often includes most of the zero-curtain period in Arctic tundra, accounted for more than 50 % of the total emissions throughout the entire cold season (September to May) in the model, compared with around 49.4 % (43 %–58 %) in observations. From 1950 to 2017, both CO 2 emissions during the zero-curtain period and during the entire cold season showed increasing trends, for example, of 0.17 and 0.36 gC m -2 yr -1 at Atqasuk. This study highlights the importance of zero-curtain periods in facilitating cold-season CH 4 and CO 2 emissions from tundra ecosystems.

58 GEOSCIENCES↗

Investigation of remote sensing techniques of measuring soil moisture

Major activities described include development and evaluation of theoretical models that describe both active and passive microwave sensing of soil moisture, the evaluation of these models for their applicability, the execution of a controlled field experiment during which passive microwave measurements were acquired to validate these models, and evaluation of previously acquired aircraft microwave measurements. The development of a root zone soil water and soil temperature profile model and the calibration and evaluation of gamma ray attenuation probes for measuring soil moisture profiles are considered. The analysis of spatial variability of soil information as related to remote sensing is discussed as well as the implementation of an instrumented field site for acquisition of soil moisture and meteorologic information for use in validating the soil water profile and soil temperature profile models.

Newton, R. W.↗

Ecosystem Carbon Emissions from 2015 Forest Fires in Interior Alaska

In the summer of 2015, hundreds of wildfires burned across the state of Alaska, and consumed more than 1.6 million ha of boreal forest and wetlands in the Yukon-Koyukuk region. Mapping of 113 large wildfires using Landsat satellite images from before and after 2015 indicated that nearly 60% of this area was burned at moderate-to-high severity levels. Field measurements near the town of Tanana on the Yukon River were carried out in July of 2017 in both unburned and 2015 burned forested areas (nearly adjacent to one-another) to visually verify locations of different Landsat burn severity classes (low, moderate, or high). Results: Field measurements indicated that the loss of surface organic layers in boreal ecosystem fires is a major factor determining post-fire soil temperature changes, depth of thawing, and carbon losses from the mineral topsoil layer. Measurements in forest sites showed that soil temperature profiles to 30 cm depth at burned forest sites increased by an average of 8o - 10o C compared to unburned forest sites. Sampling and laboratory analysis indicated a 65% reduction in soil carbon content and a 58% reduction in soil nitrogen content in severely burned sample sites compared to soil mineral samples from nearby unburned spruce forests. Conclusions: Combined with nearly unprecedented forest areas severely burned in the Interior region of Alaska in 2015, total ecosystem fire emission of carbon to the atmosphere exceeded most previous estimates for the state.

Ecosystem↗

Quantifying the influence of different biocrust community states and their responses to warming temperatures on soil biogeochemistry in field and mesocosm studies

Biocrusts influence soil biogeochemistry by fixing carbon (C) and nitrogen (N) and through leachate inputs to soils. Functional rates can vary among biocrust community states and in response to edaphic properties, heterotrophic microbial activity, and global change. Using soils and biocrusts from the Colorado Plateau, Utah, USA, we aimed to quantify the influence of early-successional (ES) and late-successional (LS) biocrust community states on soil biogeochemistry. In a field setting, we found soil was less “fertile” under ES than LS biocrust, but ES biocrust had a relative influence 1.3 times greater than LS biocrust on soil fertility. Leachate collected from LS biocrust had, on average, 6 times more organic C and 1.7 times more dissolved N than ES biocrust, but concentrations of phosphorus (P) and inorganic N did not differ among the two biocrust types. Here, to disentangle influences of biocrusts and soil properties on biogeochemical pools, we constructed mesocosms from homogenized soil and left the surface bare or covered with ES or LS biocrust, before assignment to ambient or warmed (+5 °C) temperature treatments for 3 months. Multivariate biogeochemical properties differed among cover types, yet all exhibited losses of P, N, and organic C and nearly half of the biogeochemical variables considered did not differ among cover types. Mesocosms with LS biocrust retained more dissolved N, supported 8 additional, significant correlations among biogeochemical pools of the biocrust and mineral soil layer on average, and lost fewer of these correlations under warming. Overall, while soils under LS biocrusts were more fertile (i.e., had higher nutrient concentrations) than under ES, we did not find evidence implicating leachate as the primary driver of this difference. Biocrust influences on soil fertility were greater when mineral soil nutrients were in lower concentrations, highlighting the value of even incipient biocrusts for dryland functioning.

54 ENVIRONMENTAL SCIENCES↗

Modeling microbial carbon fluxes and stocks in global soils from 1901 to 2016

Abstract. Soil microbes play a crucial role in the carbon (C) cycle; however, they have been overlooked in predicting the terrestrial C cycle. We applied a microbial-explicit Earth system model – the Community Land Model-Microbe (CLM-Microbe) – to investigate the dynamics of soil microbes during 1901 to 2016. The CLM-Microbe model was able to reproduce the variations of gross (GPP) and net (NPP) primary productivity, heterotrophic (HR) and soil (SR) respiration, microbial (MBC) biomass C in fungi (FBC) and bacteria (BBC) in the top 30 cm and 1 m, and dissolved (DOC) and soil organic C (SOC) in the top 30 cm and 1 m during 1901–2016. During the study period, simulated C variables increased by approximately 12 PgC yr−1 for HR, 25 PgC yr−1 for SR, 1.0 PgC for FBC and 0.4 PgC for BBC in 0–30 cm, and 1.2 PgC for FBC and 0.7 PgC for BBC in 0–1 m. Increases in microbial C fluxes and pools were widely found, particularly at high latitudes and in equatorial regions, but we also observed their decreases in some grids. Overall, the area-weighted averages of HR, SR, FBC, and BBC in the top 1 m were significantly correlated with those of soil moisture and soil temperature in the top 1 m. These results suggested that microbial C fluxes and pools were jointly governed by vegetation C input and soil temperature and moisture. Our simulations revealed the spatial and temporal patterns of microbial C fluxes and pools in response to environmental change, laying the foundation for an improved understanding of soil microbial roles in the global terrestrial C cycle.

He, Liyuan↗

Experimental Whole-Ecosystem Warming Enables Novel Estimation of Snow Cover and Depth Sensitivities to Temperature, and Quantification of the Snow-Albedo Feedback Effect

Climate change is reducing the amount, duration, and extent of snow across high-latitude ecosystems. But, in landscapes where persistent winter snow cover develops, experimental platforms to specifically investigate interactions between warming and changes in snowpack, and impacts on ecosystem processes, have been lacking. We leveraged a whole-ecosystem warming experiment in a boreal peatland forest to quantify how snow duration, depth, and fractional cover vary with warming of up to +9°C. We found that every snow-related quantity we examined declined precipitously as the amount of warming increased. The importance of deep, continuous snow cover for moderating shallow soil temperature is highlighted by an increase in soil temperature variance and the frequency of short-duration freeze-thaw cycles in the warmer plots. We used a paired-plot approach to estimate the magnitude of the snow-albedo feedback effect, and demonstrate that albedo-driven warming linked to reduced snow cover varies between December (+0.4°C increase in maximum air temperature) and March (+1.2°C increase) because of differences in insolation. Overall, results show that even modest future warming will have profound impacts on northern winters and cold-season ecosystem processes. Plot-level data from this warming experiment, and emergent relationships between warming and quantities related to snow cover and duration, could be of enormous value for testing and improving the representation of snow processes in simulation models, especially under future climate scenarios that are outside of the range of historically observed variability.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of simulated soil carbon dynamics in Arctic-Boreal ecosystems

Given the magnitude of soil carbon stocks in northern ecosystems, and the vulnerability of these stocks to climate warming, land surface models must accurately represent soil carbon dynamics in these regions. We evaluate soil carbon stocks and turnover rates, and the relationship between soil carbon loss with soil temperature and moisture, from an ensemble of eleven global land surface models. We focus on the region of NASA's Arctic-Boreal vulnerability experiment (ABoVE) in North America to inform data collection and model development efforts. Models exhibit an order of magnitude difference in estimates of current total soil carbon stocks, generally under- or overestimating the size of current soil carbon stocks by greater than 50 PgC. We find that a model's soil carbon stock at steady-state in 1901 is the prime driver of its soil carbon stock a hundred years later—overwhelming the effect of environmental forcing factors like climate. The greatest divergence between modeled and observed soil carbon stocks is in regions dominated by peat and permafrost soils, suggesting that models are failing to capture the frozen soil carbon dynamics of permafrost regions. Using a set of functional benchmarks to test the simulated relationship of soil respiration to both soil temperature and moisture, we find that although models capture the observed shape of the soil moisture response of respiration, almost half of the models examined show temperature sensitivities, or Q10 values, that are half of observed. Significantly, models that perform better against observational constraints of respiration or carbon stock size do not necessarily perform well in terms of their functional response to key climatic factors like changing temperature. This suggests that models may be arriving at the right result, but for the wrong reason. The results of this work can help to bridge the gap between data and models by both pointing to the need to constrain initial carbon pool sizes, as well as highlighting the importance of incorporating functional benchmarks into ongoing, mechanistic modeling activities such as those included in ABoVE.

54 ENVIRONMENTAL SCIENCES↗

Soil moisture from temperature measurements at the Earth's surface, update

Soil moisture budgets at the Earth's surface were investigated based on soil and atmospheric temperature variations. A number of data sets were plotted and statistically analyzed in order to accentuate the existence and the characteristics of mesoscale soil temperature extrema variations and their relations to other parameters. The correlations between diurnal temperature extrema for air and soil in drought and non-drought periods appear to follow different characteristic patterns, allowing an inference of soil moisture content from temperature data. The recovery of temperature extrema after a precipitation event also follows a characteristic power curve rise between two limiting values which is an indicator of evaporation rates. If these indicators are applied universally to regional temperature data, soil moisture content or drought conditions can be inferred directly from temperature measurements.

Welker, J. E.↗

Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series: Supporting Data

This dataset consists of soil temperature time series that were used to estimate soil thermal diffusivity and its uncertainty trough the probabilistic modelling approach developed and presented in the article named "Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series" and published in Earth Surface Dynamics. There are two compressed (.zip) files that contains synthetic (Synthetic_soiltemp_Teller.zip) and field (Field_soiltemp_Teller.zip) data. There is one MATLAB file that requires MATLAB to execute but any text editor can open it. The Synthetic_soiltemp_Teller.zip file includes 5 comma-delimited data files (.csv) each of which contains soil temperature time series generated through forward modeling (i.e., heat-conduction process in a heterogeneous medium using an explicit finite difference method) to mimic various types of temperature gradients, trend and fluctuations. These synthetic soil temperatures were used to investigate the impact of different environmental conditions on the uncertainty of thermal diffusivity estimates. The Field_soiltemp_Teller.zip file contains 28 comma-delimited data files (.csv) out of which (a) 27 files includes soil temperature time series recorded from 27 temperature probes located in a site along Teller Road about 40 km northwest of Nome, Alaska (64.72°N, 165.94°W), (b) one includes the name and coordinates of the 27 probes. These field soil temperatures were used to infer soil thermal diffusivity at numerous locations and depths in a discontinuous permafrost environment, and to evaluate the links between the estimated soil thermal diffusivity values and soil physical properties. The comma-delimited data files (.csv) of the synthetic and field soil temperature time series includes date and time (UTC) in the first column and soil temperature from 5 cm below the ground surface to 1.05 m depth (with 5 or 10 cm spacing between sensors) in the other columns. The measurements were acquired every 15 minutes. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Refinements to SSiB with an Emphasis on Snow-Physics: Evaluation and Validation Using GSWP and Valdai Data

Refinements to the snow-physics scheme of SSiB (Simplified Simple Biosphere Model) are described and evaluated. The upgrades include a partial redesign of the conceptual architecture to better simulate the diurnal temperature of the snow surface. For a deep snowpack, there are two separate prognostic temperature snow layers - the top layer responds to diurnal fluctuations in the surface forcing, while the deep layer exhibits a slowly varying response. In addition, the use of a very deep soil temperature and a treatment of snow aging with its influence on snow density is parameterized and evaluated. The upgraded snow scheme produces better timing of snow melt in GSWP-style simulations using ISLSCP Initiative I data for 1987-1988 in the Russian Wheat Belt region. To simulate more realistic runoff in regions with high orographic variability, additional improvements are made to SSiB's soil hydrology. These improvements include an orography-based surface runoff scheme as well as interaction with a water table below SSiB's three soil layers. The addition of these parameterizations further help to simulate more realistic runoff and accompanying prognostic soil moisture fields in the GSWP-style simulations. In intercomparisons of the performance of the new snow-physics SSiB with its earlier versions using an 18-year single-site dataset from Valdai Russia, the version of SSiB described in this paper again produces the earliest onset of snow melt. Soil moisture and deep soil temperatures also compare favorably with observations.

Mocko, David M.↗

Data from: Lowland Tropical Forests Remain a Methane Sink Under Warming and Long-Term Hurricane Disturbance Recovery

The repository folder contains spreadsheets and script for soil greenhouse gas (GHG) fluxes, soil moisture, soil temperature, air temperature, and precipitation measurements collected from the Tropical Responses to Altered Climate Experiment (TRACE) at the Sabana Research Field Station, El Yunque National Forest (USDA Forest Service; 18°19′28.74″ N, 65°43′50.09″ W) — an open-air field warming experiment located in a lowland tropical forest in Puerto Rico within the Luquillo Experimental Forest (LEF) — six to seven years after Hurricanes Irma and Maria (2017). All spreadsheets for soil and air microclimate data, as well as soil greenhouse gas data, are included as csv files. Air temperature data are also included as Excel spreadsheets (.xlsx). The script is built in R Studio, which is the only software required to run data analysis. This dataset is associated with the manuscript “Larocca Conte G ; Zuvela L ; Cruz-Pérez R ; Barreto-Vélez T ; Becerra-Santillan N ; Campbell S ; Chu H ; Dam T ; Grullón-Penkova I ; Kleit M ; Ortiz-Iglesias D ; Rubio-Lebrón L ; Cavaleri M ; Reed S ; Sihi D ; Wood T ; O'Connell C., 2026. Lowland Tropical Forests Remain a Methane Sink Under Warming and Long-Term Hurricane Disturbance Recovery. Agricultural and Forest Meteorology. In review". The dataset was used to test the effect of warming on soil CH4 dynamics following long-term legacy effects of hurricane disturbance. The dataset includes: - An overall README file in word and pdf format describing methodology and spreadsheets’ structure. - Continuous measurements of soil temperature and moisture from January 2023 to July 2024 measured with Campbell CS655 probes (“TRACE_soil_temperature_and_moisture_2023_cleaned(in).csv” and “TRACE_soil_temperature_and_moisture_2024_cleaned. csv”). - Air temperature data measured with a HOBO MX23O1A data logger (“Hobo air temperature 2023 Sep 2024” and “Hobo air temperature 2023 Sep 2024” – “CSV FILES folders”). - Precipitation data from a nearby weather tower downloaded from González et al. (2025; “sabana_2020-2025.csv”). - Soil CH4 and CO2 effluxes measured intermittently in two summer campaigns (June – August 2023 and June – July 2024) with a LI-COR 8200-01S Portable Smart Chamber coupled with a LI-COR LI-7810 CH4/ CO2/H2O Trace Gas Analyzer (“23_24COMBO2.0.csv”). - R markdown script for data analysis (“Trace new_PLOTS.Rmd”).

54 ENVIRONMENTAL SCIENCES↗

Evaluation of HCMM data for assessing soil moisture and water table depth

Data were analyzed for variations in eastern South Dakota. Soil moisture in the 0-4 cm layer could be estimated with 1-mm soil temperatures throughout the growing season of a rainfed barley crop (% cover ranging from 30% to 90%) with an r squared = 0.81. Empirical equations were developed to reduce the effect of canopy cover when radiometrically estimating the 1-mm soil temperature, r squared = 0.88. The corrective equations were applied to an aircraft simulation of HCMM data for a diversity of crop types and land cover conditions to estimate the 0-4 cm soil moisture. The average difference between observed and measured soil moisture was 1.6% of field capacity. HCMM data were used to estimate the soil moisture for four dates with an r squared = 0.55 after correction for crop conditions. Location of shallow alluvial aquifers could be accomplished with HCMM predawn data. After correction of HCMM day data for vegetation differences, equations were developed for predicting water table depths within the aquifer (r=0.8).

Moore, D. G.↗

Modeling Modern Methane Emissions from Natural Wetlands: Model Description and Results - 1

Methane is an important greenhouse gas which contributes about 22 percent to the present greenhouse effect. Natural wetlands currently constitute the biggest methane source and were the major source in preindustrial times. Wetland emissions depend highly on the climate, i.e., on soil temperature and water table. To investigate the response of methane emissions from natural wetlands to climate variations, a process-based model that derives methane emissions from natural wetlands as a function of soil temperature, water table, and net primary productivity is used. For its application on the global scale, global data sets for all model parameters are generated. In addition, a simple hydrologic model is developed in order to simulate the position of the water table in wetlands. The hydrologic model is tested against data from different wetland sites, and the sensitivity of the hydrologic model to changes in precipitation is examined. The global methane­ hydrology model constitutes a tool to study temporal and spatial variations in methane emissions from natural wetlands. The model is applied using high-frequency atmospheric forcing fields from European Center for Medium-range Weather Forecasts (ECMWF) re-analyses of the period from 1982 to 1993. We calculate global annual methane emissions from wetlands to be 260 teragrams per year. Twenty-five percent of these methane emissions originate from wetlands north of 30 degrees North Latitude. Only 60 percent of the produced methane is emitted, while the rest is re-oxidized. A comparison of zonal integrals of simulated global wetland emissions and results obtained by an inverse modeling approach shows good agreement. In a test with data from two wetlands the seasonality of simulated and observed methane emissions agrees well.

Biogeochemical processes↗

Changes in precipitation and air temperature contribute comparably to permafrost degradation in a warmer climate

Abstract Surface energy budgets of high-latitude permafrost systems are poorly represented in Earth system models (ESMs), yet permafrost is rapidly degrading and these dynamics are critical to future carbon-climate feedback predictions. A potentially important factor in permafrost degradation neglected so far by ESMs is heat transfer from precipitation, although increases in soil temperature and thaw depth have been observed following increases in precipitation. Using observations and a mechanistic ecosystem model, we show here that increases in precipitation hasten active layer development beyond that caused by surface air warming across the North Slope of Alaska (NSA) under recent and 21st century climate (RCP8.5). Modeled active layer depth (ALD) in simulations that allow precipitation heat transfer agreed very well with observations from 28 Circumpolar Active Layer Monitoring sites (R 2 = 0.63; RMSE = 10 cm). Simulations that ignored precipitation heat transfer resulted in lower spatially-averaged soil temperatures and a 39 cm shallower ALD by 2100 across the NSA. The results from our sensitivity analysis show that projected increases in 21st century precipitation deepen the active layer by enhancing precipitation heat transfer and ground thermal conductivity, suggesting that precipitation is as important an environmental control on permafrost degradation as surface air temperature. We conclude that ESMs that do not account for precipitation heat transfer likely underestimate ALD rates of change, and thus likely predict biased ecosystem responses.

54 ENVIRONMENTAL SCIENCES↗

The Land Surface in Current and Planned MERRA Reanalysis Products

Current global atmospheric reanalysis products such as the European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5), the NASA Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2), and the Japanese Reanalysis for Three Quarters of a Century (JRA-3Q) provide estimates of land surface states and fluxes, including soil moisture, soil temperature, snow mass, latent and sensible heat fluxes, and runoff, that are widely used in research and applications. These land surface estimates are based on land surface process models and, depending on the reanalysis product, on precipitation observations or the assimilation of land surface observations of soil moisture, soil temperature, snow conditions, and screen-level air temperature and humidity from satellite observations and in situ measurements. In this presentation, we review the land surface modeling and data assimilation components of the suite of current and planned MERRA reanalysis products. In addition to MERRA-2, we will discuss the latest NASA reanalysis, MERRA for the 21st century (M21C), which is currently under production, as well as the development and planning of the next version of the MERRA reanalysis, tentatively labeled MERRA-3. In MERRA-2, observations-based precipitation data products are used to correct the precipitation falling on the land surface. Outside of the high-latitudes and Africa, the daily, 0.5-degree, gauge-based Climate Prediction Center (CPC) Unified (CPCU) product is used. In Africa, the pentad, 2.5-degree, satellite- and gauge-based CPC Merged Analysis of Precipitation (CMAP) product is used. Poleward of 62.5 degrees latitude, the land surface sees the precipitation generated by the atmospheric model in the cycling data assimilation system. This configuration provides improved soil moisture estimates compared to those of the original (version 1) MERRA estimates, which did not benefit from the use of precipitation observations. Moreover, the use of precipitation observations facilitates a seamless spin-up of the land surface initial conditions across the MERRA-2 production streams. The use of a gauge-only precipitation product in MERRA-2 across much of the globe, however, adversely impacts the quality of the MERRA-2 land surface estimates in regions with poor gauge coverage, including most of South America and Australia. Therefore, the forthcoming M21C reanalysis uses satellite- and gauge-based precipitation from the Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement Mission (IMERG). This change results in significant improvements in the quality of the M21C soil moisture estimates in the Southern Hemisphere compared to those from MERRA-2. Planning for MERRA-3 focuses on the assimilation of soil moisture observations from the Soil Moisture Active Passive (SMAP) mission and the Advanced Scatterometer (ASCAT), along with snow cover area fraction observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) to further improve the quality of the land surface estimates from the reanalysis. As a first step towards the assimilation of land surface observations in MERRA-3, the offline (land-only) M21C-Land reanalysis is currently under development as a supplemental M21C product that includes the assimilation of SMAP, ASCAT, and MODIS observations. Preliminary results from M21C and M21C-Land will be discussed in the context of MERRA-2 and plans for MERRA-3.

Rolf Reichle↗

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

58 GEOSCIENCES↗

A Process-based, Climate-Sensitive Model to Derive Methane Emissions from Natural Wetlands: Application to 5 Wetland Sites, Sensitivity to Model Parameters and Climate

Methane emissions from natural wetlands constitutes the largest methane source at present and depends highly on the climate. In order to investigate the response of methane emissions from natural wetlands to climate variations, a 1-dimensional process-based climate-sensitive model to derive methane emissions from natural wetlands is developed. In the model the processes leading to methane emission are simulated within a 1-dimensional soil column and the three different transport mechanisms diffusion, plant-mediated transport and ebullition are modeled explicitly. The model forcing consists of daily values of soil temperature, water table and Net Primary Productivity, and at permafrost sites the thaw depth is included. The methane model is tested using observational data obtained at 5 wetland sites located in North America, Europe and Central America, representing a large variety of environmental conditions. It can be shown that in most cases seasonal variations in methane emissions can be explained by the combined effect of changes in soil temperature and the position of the water table. Our results also show that a process-based approach is needed, because there is no simple relationship between these controlling factors and methane emissions that applies to a variety of wetland sites. The sensitivity of the model to the choice of key model parameters is tested and further sensitivity tests are performed to demonstrate how methane emissions from wetlands respond to climate variations.

Walter, Bernadette P.↗

Ground Heat Flux Reconstruction Using Bayesian Uncertainty Quantification Machinery and Surrogate Modeling

Ground heat flux (G 0 ) is a key component of the land-surface energy balance of high-latitude regions. Despite its crucial role in controlling permafrost degradation due to global warming, G 0 is sparsely measured and not well represented in the outputs of global scale model simulation. In this study, an analytical heat transfer model is tested to reconstruct G 0 across seasons using soil temperature series from field measurements, Global Climate Model, and climate reanalysis outputs. The probability density functions of ground heat flux and of model parameters are inferred using available G 0 data (measured or modeled) for snow-free period as a reference. When observed G 0 is not available, a numerical model is applied using estimates of surface heat flux (dependent on parameters) as the top boundary condition. These estimates (and thus the corresponding parameters) are verified by comparing the distributions of simulated and measured soil temperature at several depths. Aided by state-of-the-art uncertainty quantification methods, the developed G 0 reconstruction approach provides novel means for assessing the probabilistic structure of the ground heat flux for regional permafrost change studies.

54 ENVIRONMENTAL SCIENCES↗