Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Snow model”

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 937 records · Page 52

NCA-LDAS: Overview and Analysis of Hydrologic Trends for the National Climate Assessment

Terrestrial hydrologic trends over the conterminous United States are estimated for 1980-47 2015 using the National Climate Assessment-Land Data Assimilation System (NCA-LDAS) reanalysis. NCA-LDAS employs the uncoupled Noah Version 3.3 land surface model at 0.125 degreesx0.125 degrees forced with NLDAS-2 meteorology, rescaled Climate Prediction Center precipitation, and assimilated satellite-based soil moisture, snow depth and irrigation products. Mean annual trends are reported using the nonparametric Mann-Kendall test at p<0.1 significance. Results illustrate the interrelationship between regional gradients in forcing trends, and trends in other land energy and water stores and fluxes. Mean precipitation trends range from +3 to +9 mm/yr in the Upper Great Plains and Northeast to -1 to -9 in the West and South; net radiation flux trends range from +0.05 to +0.20 W/m(exp2)/yr in the East to -0.05 to -0.20 in the West; U.S.-wide temperature trends average about +0.03K/yr. Trends in soil moisture, snow cover, latent and sensible heat fluxes and runoff are consistent with forcings, contributing to increasing evaporative fraction trends from West to East. Evaluation of NCA-LDAS trends compared to independent data indicates mixed results. The RMSE of U.S.-wide trends in number of snow cover days improved from 3.13 to 2.89 days/yr while trend detection increased 11%; Trends in latent heat flux were hardly affected, RMSE decreasing only 0.17 to 0.16 W/m(exp 2)/yr, while trend detection increased 2%; NCA-LDAS runoff trends degraded significantly from 2.6 to 16.1 mm/yr while trend detection was unaffected. Analysis also indicated that NCA-LDAS exhibits relatively more skill in low precipitation station density areas, suggesting there are limits to the effectiveness of satellite data assimilation in densely gaged regions. Overall, NCA-LDAS demonstrates capability for quantifying physically consistent, U.S. hydrologic climate trends over the satellite era.

Snow depth↗

Polarimetric Radar Characteristics of Simulated and Observed Intense Convective Cores for a Midlatitude Continental and Tropical Maritime Environment

This study contrasts midlatitude continental and tropical maritime deep convective cores using polarimetric radar observables and retrievals from selected deep convection episodes during the MC3E and TWPICE field campaigns. The continental convective cores produce stronger radar reflectivities throughout the profiles, while maritime convective cores produce more positive differential reflectivity Z(dr) and larger specific differential phase K(dp) above the melting level. Hydrometeor identification retrievals revealed the presence of large fractions of rimed ice particles (snow aggregates) in the continental (maritime) convective cores, consistent with the Z(dr) and K(dp) observations. The regional cloud-resolving model simulations with bulk and size-resolved bin microphysics are conducted for the selected cases, and the simulation outputs are converted into polarimetric radar signals and retrievals identical to the observational composites. Both the bulk and the bin microphysics reproduce realistic land and ocean (L-O) contrasts in reflectivity, polarimetric variables of rain drops, and hydrometeor profiles, but there are still large uncertainties in describing Zdr and Kdp of ice crystals associated with the ice particle shapes/orientation assumptions. Sensitivity experiments are conducted by swapping background aerosols between the continental and maritime environments, revealing that background aerosols play a role in shaping the distinct L-O contrasts in radar reflectivity associated with raindrop sizes, in addition to the dominant role of background thermodynamics.

Toshihisa Matsui↗

NCA-LDAS: Overview and Analysis of Hydrologic Trends for the National Climate Assessment

Terrestrial hydrologic trends over the conterminous United States are estimated for 1980–2015 using the National Climate Assessment Land Data Assimilation System (NCA-LDAS) reanalysis. NCA-LDAS employs the uncoupled Noah version 3.3 land surface model at 0.125° 3 0.125° forced with NLDAS-2 meteorology, rescaled Climate Prediction Center precipitation, and assimilated satellite-based soil moisture, snow depth, and irrigation products. Mean annual trends are reported using the nonparametric Mann–Kendall test at p < 0.1 significance. Results illustrate the interrelationship between regional gradients in forcing trends and trends in other land energy and water stores and fluxes. Mean precipitation trends range from +3 to +9 mm/yr in the upper Great Plains and Northeast to -1 to -9 mm/yr in the West and South, net radiation flux trends range from +0.05 to +0.20 W/sq. m yr in the East to -0.05 to -0.20 W/sq. m yr in the West, and U.S.-wide temperature trends average about +0.03K/yr. Trends in soil moisture, snow cover, latent and sensible heat fluxes, and runoff are consistent with forcings, contributing to increasing evaporative fraction trends from west to east. Evaluation of NCA-LDAS trends compared to independent data indicates mixed results. The RMSE of U.S.-wide trends in number of snow cover days improved from 3.13 to 2.89 days/yr while trend detection increased 11%. Trends in latent heat flux were hardly affected, with RMSE decreasing only from 0.17 to 0.16 W/sq. m yr, while trend detection increased 2%. NCA-LDAS runoff trends degraded significantly from 2.6 to 16.1 mm/yr while trend detection was unaffected. Analysis also indicated that NCA-LDAS exhibits relatively more skill in low precipitation station density areas, suggesting there are limits to the effectiveness of satellite data assimilation in densely gauged regions. Overall, NCA-LDAS demonstrates capability for quantifying physically consistent, U.S. hydrologic climate trends over the satellite era.

Michael F. Jasinski↗

End-of-Winter Snow Depth, Temperature, Density, and SWE Measurements at Teller Road Site, Seward Peninsula, Alaska, 2023

Measurements of end-of-winter snow properties were collected at the NGEE Arctic Teller Road Site at mile marker 27 (TL_MM27) from March 31 to April 3, 2023. This dataset contains one *.pdf user guide and four *.csv data files of spatially distributed values of snow depth, snow water equivalent (SWE), snow temperature, and snow density. Snow temperature and snow density were collected throughout the snowpack at different depths. Data were collected toward the end of the winter season from late March to early April, when the snowpack is typically near its maximum. A Snow-Hydro MagnaProbe (http://www.snowhydro.com/products/column2.html) was used to improve collection efficiency and enhance spatial coverage. This dataset is a continuation of the previous end-of-winter snow surveys conducted at the Teller Road Site in 2016-2018 (Wilson et al., 2020; https://doi.org/10.5440/1592103), 2019 (Bennett et al., 2021; https://doi.org/10.5440/1798170), and 2022 (Bennett et al., 2022; https://doi.org/10.5440/1887250).The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) was a 15-year research effort (2012-2027) 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↗

Assessment of climate variability of the Greenland Ice Sheet: Integration of in situ and satellite data

The proposed research involves the application of multispectral satellite data in combination with ground truth measurements to monitor surface properties of the Greenland Ice Sheet which are essential for describing the energy and mass of the ice sheet. Several key components of the energy balance are parameterized using satellite data and in situ measurements. The analysis will be done for a ten year time period in order to get statistics on the seasonal and interannual variations of the surface processes and the climatology. Our goal is to investigate to what accuracy and over what geographic areas large scale snow properties and radiative fluxes can be derived based upon a combination of available remote sensing and meteorological data sets. Operational satellite sensors are calibrated based on ground measurements and atmospheric modeling prior to large scale analysis to ensure the quality of the satellite data. Further, several satellite sensors of different spatial and spectral resolution are intercompared to access the parameter accuracy. Proposed parameterization schemes to derive key component of the energy balance from satellite data are validated. For the understanding of the surface processes a field program was designed to collect information on spectral albedo, specular reflectance, soot content, grain size and the physical properties of different snow types. Further, the radiative and turbulent fluxes at the ice/snow surface are monitored for the parameterization and interpretation of the satellite data. The expected results include several baseline data sets of albedo, surface temperature, radiative fluxes, and different snow types of the entire Greenland Ice Sheet. These climatological data sets will be of potential use for climate sensitivity studies in the context of future climate change.

Steffen, K.↗

Improving the interface processes in the DOE/ACME model (Final Report)

This is the Final Report of our 4-year Energy Exascale Earth System Model (E3SM) project (3 years for the original project + 1-year no-cost extension). The overall objective of our project is to improve the interface processes in the E3SM. Two overarching questions have been addressed: 1) What are the major deficiencies of the (land-atmosphere, ocean-atmosphere, land-ocean, and snow-sea ice) interface processes in E3SM? 2) How can we improve the treatment of these deficiencies in E3SM? Four tasks have been carried out: to evaluate and improve the land–atmosphere coupling in E3SM; to evaluate and improve the ocean–atmosphere coupling in E3SM; to evaluate and improve the land–ocean coupling in E3SM; and to evaluate and improve the snow–sea ice coupling in E3SM.

54 ENVIRONMENTAL SCIENCES↗

Ultimate Bifacial Showdown: 75kW Field Results

This work presents one year results for a side-by-side comparison of 5 different bifacial technologies, deployed in a 75-kW single-axis tracked field at Golden, CO. Four PERC (multi and mono) and 1 Silicon Heterojunction manufacturers are represented. Reference strings of equivalent monofacial PV modules are also installed for a direct technology comparison and calculation of bifacial gain. Analysis indicates performance within expectation, with a cumulative model mean error within +/- 2% for both bifacial and monofacial models and cumulative bifacial energy gain between 6-9%. High accuracy string-level DC monitoring and module-level measurements are also presented, showing effects of snow and edge effects in the order of +/- 3% for non-uniform irradiance across the row. Data from a custom-module to measure shading loss from the torque-tube is installed, showing up to 6% irradiance non-uniformity for the month of January.

bifacial field validation↗

Drought Indicators Based on Model Assimilated GRACE Terrestrial Water Storage Observations

The Gravity Recovery and Climate Experiment (GRACE) twin satellites observe time variations in Earth's gravity field which yield valuable information about changes in terrestrial water storage (TWS). GRACE is characterized by low spatial (greater than 150,000 square kilometers) and temporal (greater than 10 day) resolution but has the unique ability to sense water stored at all levels (including groundwater) systematically and continuously. The GRACE Data Assimilation System (GRACE-DAS), based on the Catchment Land Surface Model (CLSM) enhances the value of the GRACE water storage data by enabling spatial and temporal downscaling and vertical decomposition into moisture 39 components (i.e. groundwater, soil moisture, snow), which individually are more useful for scientific applications. In this study, GRACE-DAS was applied to North America and GRACE-based drought indicators were developed as part of a larger effort that investigates the possibility of more comprehensive and objective identification of drought conditions by integrating spatially, temporally and vertically disaggregated GRACE data into the U.S. and North American Drought Monitors. Previously, the Drought Monitors lacked objective information on deep soil moisture and groundwater conditions, which are useful indicators of drought. Extensive datasets of groundwater storage from USGS monitoring wells and soil moisture from the Soil Climate Analysis Network (SCAN) were used to assess improvements in the hydrological modeling skill resulting from the assimilation of GRACE TWS data. The results point toward modest, but statistically significant, improvements in the hydrological modeling skill across major parts of the United States, highlighting the potential value of GRACE assimilated water storage field for improving drought detection.

Houborg, Rasmus↗

End-of-Winter Snow Depth, Temperature, Density and SWE Measurements at Teller Road Site, Seward Peninsula, Alaska, 2016-2018

Measurements of end-of-winter snow properties were collected at the NGEE Arctic Teller Road Site at mile marker 27 (TL_MM27) for three consecutive years - 2016, 2017, and 2018. This dataset contains one *.pdf user guide and four *.csv data files spatially distributed values of snow depth, snow water equivalent (SWE), temperature at the snow surface and at the bottom of the snowpack, and snow density as a function of depth from top of snowpack. Data was collected toward the end of the winter season, typically in late March or early April, when the snowpack would be near its maximum. In 2016, snow depth measurements were taken using either a thaw probe or avalanche probe. In 2017 and 2018 a Snow-Hydro MagnaProbe http://www.snowhydro.com/products/column2.html was used to improve collection efficiency and enhance spatial coverage. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), 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↗

End-of-Winter Snow Depth, Temperature, Density, and SWE Measurements at Kougarok Road Site, Seward Peninsula, Alaska, 2022

Measurements of end-of-winter snow properties were collected at the NGEE Arctic Kougarok Road Site at mile marker 64 (KG_MM64) from April 2 to April 5, 2022. This dataset contains three *.csv data files of spatially distributed values of snow depth, snow water equivalent (SWE), and snow temperature at the snow surface and at the bottom of the snowpack. Data were collected toward the end of the winter season during early April, when the snowpack would be near its maximum. A Snow-Hydro MagnaProbe (http://www.snowhydro.com/products/column2.html) was used to improve collection efficiency and enhance spatial coverage. A user guide is included as a *.pdf. This dataset is a continuation of the previous end-of-winter snow surveys conducted at the Kougarok Road Site in 2018 (Wilson et al., 2021; https://doi.org/10.5440/1593874).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↗

Empirical mass-loss rates for 25 O and early B stars, derived from Copernicus observations

Ultraviolet line profiles are fitted with theoretical line profiles in the cases of 25 stars covering a spectral type range from O4 to B1, including all luminosity classes. Ion column densities are compared for the determination of wind ionization, and it is found that the O VI/N V ratio is dependent on the mean density of the wind and not on effective temperature value, while the Si IV/N V ratio is temperature-dependent. The column densities are used to derive a mass-loss rate parameter that is empirically correlated against the mass-loss rate by means of standard stars with well-determined rates from IR or radio data. The empirical mass-loss rates obtained are compared with those derived by others and found to vary by as much as a factor of 10, which is shown to be due to uncertainties or errors in the ionization fractions of models used for wind ionization balance prediction.

Gathier, R.↗

Impacts of Microphysical Scheme on Convective and Stratiform Characteristics in Two High Precipitation Squall Line Events

This study investigates the impact of snow, graupel, and hail processes on simulated squall lines over the Southern Great Plains in the United States. The Weather Research and Forecasting (WRF) model is used to simulate two squall line events in Oklahoma during May 2007, and the simulations are validated against radar and surface observations. Several microphysics schemes are tested in this study, including the WRF 5-Class Microphysics (WSM5), WRF 6-Class Microphysics (WSM6), Goddard Cumulus Ensemble (GCE) Three Ice (3-ice) with graupel, Goddard Two Ice (2-ice), and Goddard 3-ice hail schemes. Simulated surface precipitation is sensitive to the microphysics scheme when the graupel or hail categories are included. All of the 3-ice schemes overestimate the total precipitation with WSM6 having the largest bias. The 2-ice schemes, without a graupel/hail category, produce less total precipitation than the 3-ice schemes. By applying a radar-based convective/stratiform partitioning algorithm, we find that including graupel/hail processes increases the convective areal coverage, precipitation intensity, updraft, and downdraft intensities, and reduces the stratiform areal coverage and precipitation intensity. For vertical structures, simulations have higher reflectivity values distributed aloft than the observed values in both the convective and stratiform regions. Three-ice schemes produce more high reflectivity values in convective regions, while 2-ice schemes produce more high reflectivity values in stratiform regions. In addition, this study has demonstrated that the radar-based convective/stratiform partitioning algorithm can reasonably identify WRF-simulated precipitation, wind, and microphysical fields in both convective and stratiform regions.

snow↗

The 20-22 January 2007 Snow Events over Canada: Microphysical Properties

One of the grand challenges of the Global Precipitation Measurement (GPM) mission is to improve precipitation measurements in mid- and high-latitudes during cold seasons through the use of high-frequency passive microwave radiometry. Toward this end, the Weather Research and Forecasting (WRF) model with the Goddard microphysics scheme is coupled with a Satellite Data Simulation Unit (WRF-SDSU) that has been developed to facilitate over-land snowfall retrieval algorithms by providing a virtual cloud library and microwave brightness temperature (Tb) measurements consistent with the GPM Microwave Imager (GMI). This study tested the Goddard cloud microphysics scheme in WRF for snowstorm events (January 20-22, 2007) that took place over the Canadian CloudSAT/CALIPSO Validation Project (C3VP) ground site (Centre for Atmospheric Research Experiments - CARE) in Ontario, Canada. In this paper, the performance of the Goddard cloud microphysics scheme both with 2ice (ice and snow) and 3ice (ice, snow and graupel) as well as other WRF microphysics schemes will be presented. The results are compared with data from the Environment Canada (EC) King Radar, an operational C-band radar located near the CARE site. In addition, the WRF model output is used to drive the Goddard SDSU to calculate radiances and backscattering signals consistent with direct satellite observations for evaluating the model results.

Tao. W.K.↗

Normalization Method and Application for MODIS TEB Assessments using Earth Scene Measurements

Selected Earth targets are commonly used for satellite sensor calibration assessments (e.g. sensor stability and inter-comparisons). Moreover, typical scenes used for the calibration assessment of the thermal emissive bands (TEBs) includeDome Concordia (Dome-C), ocean, desert, and deep convective clouds (DCC). Reference data used for these calibration assessments can come from another band, another instrument, or ground measurements. The Dome-C site, covered with uniformly-distributed permanent snow, is normally used for the assessment of the TEBsat cold temperatures. Furthermore, ocean and desert measurements prove useful for scenes with higher temperatures. The DCC, one of the most consistent and coldest targets, can be used for the TEBs calibration and product stability assessments. MODIS band 31 (~ 11 m) can be used as a reference for these scenes. However, measurements over these scenes have seasonal variations, and the DCC brightness temperatures (BTs) have asymmetrical distributions. These features can introduce additional uncertainty to the stability assessments. A normalization method is applied by using an empirical model to derive reference-dependent BTs. Using the developed empirical model, measurements can be normalized to a reference BT in order to enhance the calibration assessment’s accuracy. This method is evaluated using all four scene types (i.e. ocean, desert, snow (Dome-C), and DCC) and applied to all the Terra and Aqua MODIS TEBs. Stability assessments over the instruments’ entire data records are presented and discussed. The technique can be applied in future efforts to support MODIS TEBs calibration assessments.

MODIS TEB↗

Anticipating how rain-on-snow events will change through the 21st century: lessons from the 1997 new year’s flood event

The California-Nevada 1997 New Year’s flood was an atmospheric river (AR)-driven rain-on-snow (RoS) event and remains the costliest in their history. The joint occurrence of saturated soils, rainfall, and snowmelt generated inundation throughout northern California-Nevada. Although AR RoS events are projected to occur more frequently with climate change, the warming sensitivity of their flood drivers across scales remains understudied. We leverage the regionally refined mesh capabilities of the Energy Exascale Earth System Model (RRM-E3SM) to recreate the 1997 New Year’s flood with horizontal grid spacings of 3.5 km across California, with forecast lead times of up to 4 days, and across six warming levels ranging from pre-industrial conditions to +3.5° C. We describe the sensitivity of the flood drivers to warming including AR duration and intensity, precipitation phase, intensity and efficiency, snowpack mass and energy changes, and runoff efficiency. Our findings indicate current levels of climate change negligibly influence the flood drivers. At warming levels ≥ 1.7° C, AR hazard potential increases, snowpack nonlinearly decreases, antecedent soil moisture decreases (except where the snowline retreats), and runoff decreases (except in the southern Sierra Nevada where antecedent snowpack persists). Storm total precipitation increases, but at rates below warming-induced increases in saturation-specific humidity. Warming intensifies short-duration, high-intensity rainfall, particularly where snowfall-to-rainfall transitions occur. This study highlights the nonlinear tradeoffs in 21st-century RoS flood hazards with warming and provides water management and infrastructure investment adaptation considerations.

54 ENVIRONMENTAL SCIENCES↗

Quality Data Essential for Modeling Water Cycles Effectively

The experimental and observational data is the backbone of effective modeling of all aspects of water cycles. As the computer technology used for these data gathering efforts has improved dramatically in the recent decades, the environmental conditions these instruments are subjected to remain pretty much the same. These measuring instruments are exposed to the sun light, wind, rain, snow and ice; and they are immersed in dirt, water and mud. All of these are damaging to the electronics components used for data gathering, storage, and transmission, which creates many different kinds of data quality issues that affect the quality and reliability of modeling and prediction efforts.

54 ENVIRONMENTAL SCIENCES↗

Influence of Tundra Polygon Type and Climate Variability on CO 2 and CH 4 Fluxes Near Utqiagvik, Alaska

Arctic tundra has the potential to generate significant climate feedbacks, but spatial complexity makes it difficult to quantify the impacts of climate on ecosystem-atmosphere fluxes, particularly in polygonal tundra comprising wetter and drier polygon types on the scale of tens of meters. We measured CO 2 , CH 4 , and energy fluxes using eddy covariance for 7 yr (April to November, 2013–2019) in polygonal tundra near Utqiagvik, Alaska. This period saw the earliest snowmelt, latest snow accumulation, and hottest summer on record. To estimate fluxes by polygon type, we combined a polygon classification with a flux-footprint model. Methane fluxes were highest in the summer months but were also large during freeze-up and increased with the warming trend in August–November temperatures. While CO 2 respiration had a consistent, exponential relationship with temperature, net ecosystem exchange was more variable among years. CO 2 and CH 4 exchange (June–September) ranged between -0.83 (Standard error [SE] = 0.03) and -1.32 (SE = 0.04) μmol m -2 s-1 and 13.92 (SE = 0.26)—23.42 (SE = 0.45) nmol m -2 s -1 , respectively, and varied interannually (p ≤ 0.05). The maximum-influence method effectively attributed fluxes to polygon types. Areas dominated by low-centered polygons had higher CO 2 fluxes except in 2016–2017. Methane fluxes were highest in low-centered polygons 2013–2015 and in flat-centered polygons in subsequent years, possibly due to increasing temperature and precipitation. Sensible and latent heat fluxes also varied significantly among polygon types. Accurate characterization of Arctic fluxes and their climate dependencies requires spatial disaggregation and long term observations.

54 ENVIRONMENTAL SCIENCES↗

Near‐Surface Hydrology and Soil Properties Drive Heterogeneity in Permafrost Distribution, Vegetation Dynamics, and Carbon Cycling in a Sub‐Arctic Watershed

Abstract Discontinuous permafrost environments exhibit strong spatial heterogeneity at scales too small to be driven by weather forcing or captured by Earth System Models. Here we analyze effects of observed spatial heterogeneity in soil and vegetation properties, hydrology, and thermal dynamics on ecosystem carbon dynamics in a watershed on the Seward Peninsula in Alaska. We apply a Morris global sensitivity analysis to a process‐rich, successfully tested terrestrial ecosystem model (TEM), ecosys , varying soil properties, boundary conditions, and weather forcing. We show that landscape heterogeneity strongly impacts soil temperatures and vegetation composition. Snow depth, O‐horizon thickness, and near‐surface water content, which vary at scales of O(m), control the soil thermal regime more than an air temperature gradient corresponding to a 140 km north–south distance. High shrub productivity is simulated only in talik (perennially unfrozen) soils with high nitrogen availability. Through these effects on plant and permafrost dynamics, landscape heterogeneity impacts ecosystem productivity. Simulations with near‐surface taliks have higher microbial respiration (by 78.0 gC m −2 yr −1 ) and higher net primary productivity (by 104.9 gC m −2 yr −1 ) compared to runs with near‐surface permafrost, and simulations with high shrub productivity have outlying values of net carbon uptake. We explored the prediction uncertainty associated with ignoring observed landscape heterogeneity, and found that watershed net carbon uptake is 60% larger when heterogeneity is accounted for. Our results highlight the complexity inherent in discontinuous permafrost environments and demonstrate that missing representation of subgrid heterogeneity in TEMs could bias predictions of high‐latitude carbon budget.

54 ENVIRONMENTAL SCIENCES↗