Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Climate sensitivity”

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 235 records · Page 13

Critical habitat identification of peripheral Sage Thrashers under climate change

Abstract The Sage Thrasher ( Oreoscoptes montanus ) has been assessed as “Endangered” in Canada since 1992. Like other species with a geographic range that barely extends into Canada, Sage Thrashers are rare. Thirty‐one percent of Canada's bird species listed for recovery under the Canadian Species at Risk Act (SARA) are at the periphery of their range. A listing of “endangered” under SARA requires identification of critical habitat for the species. With anticipated climate change, recovery of species requires a more proactive intervention than relying on historical occurrence to locate suitable habitat. We synthesized 19 years of Sage Thrasher occurrence and related habitat data across the species' northern range in British Columbia (BC) and Washington (WA) to define critical habitat characteristics. We found Sage Thrashers selected less leaf litter and less grass cover in flat or low‐slope regions farther from anthropogenic or natural habitat breaks; habitat sensitive to the expected climate change impacts of fire, changes in precipitation, and invasive species establishment. By augmenting the BC data collected in the species' peripheral range with data from their core distribution in WA, we identified key habitat elements of an otherwise data‐poor species that do not breed in sufficient numbers in Canada to reliably characterize their habitat. These methods improve the identification of “critical habitat” for peripheral species like Sage Thrashers in preparation for climate‐induced range expansion northward. The framework developed demonstrates a useful template for conservation strategies for data‐limited peripheral populations in other regions. Focusing on the landscape‐level variables that indicate good habitat, and not the locations of habitat, can identify suitable future areas for conservation.

Millikin, Rhonda L.↗

How does uncertainty of soil organic carbon stock affect the calculation of carbon budgets and soil carbon credits for croplands in the U.S. Midwest?

Cropland carbon budget depicts the amount of carbon flowing in and out of agroecosystems and the changes in carbon stocks of soil and living biomass during the same period. Soil carbon credit is the additional change in soil carbon stock under certain farming practices compared with the business-as-usual practices. Accurately calculating cropland carbon budget and soil carbon credit is critical to assessing climate change mitigation potential in agroecosystems. The calculation of cropland carbon budget and soil carbon credit is sensitive to local soil and climatic conditions, especially initial soil organic carbon (SOC) stock, which is determined by both SOC concentration (SOC%) and bulk density (Bulk_Density). SOC stock data are either from soil sampling or gridded public survey data. In agroecosystem models, SOC stock data are a key model input for quantifying cropland carbon budget and soil carbon credit. However, various types and degrees of uncertainties exist in SOC stock datasets, which propagate to the quantification of SOC stock change. In particular, a large discrepancy is found in two widely used SOC stock datasets — Rapid Carbon Assessment dataset (RaCA) and Gridded Soil Survey Geographic Database (gSSURGO) — in the U.S. Midwest, with a relative difference (quantified using Normalized Root Mean Square Error, NRMSE) of 48.0% for 0–30 cm SOC stock between the two datasets. It remains largely unclear how uncertainty in SOC stocks affects the calculation of cropland carbon budget and soil carbon credit. To address this question, we used a well-validated process-based agroecosystem model, ecosys, to assess the impacts of SOC stock uncertainty on carbon budget and soil carbon credit calculation in the U.S. Midwestern corn-soybean rotation systems. Our results reveal the following findings: (1) A sizable discrepancy exists in simulated cropland carbon budget between using gSSURGO and using RaCA for their SOC% and Bulk_Density as model inputs, with a Pearson correlation coefficient (r) of only 0.4 for simulated change of SOC stock (ΔSOC) using these two different soil datasets. (2) Simulated cropland carbon budget components were more sensitive to initial SOC% than to Bulk_Density. For example, the upper and lower quartiles of multi-year averaged ΔSOC were –29.8 and 4.8 gC/m 2 /year for the selected counties respectively, with an uncertainty of 13.7 and 0.7 gC/m 2 /year induced by uncertainties in initial SOC% and Bulk_Density, respectively. (3) Both simulated ΔSOC and its uncertainty were negatively correlated with initial SOC%, whereas ΔSOC was negatively correlated with air temperature, and ΔSOC uncertainty was positively correlated with air temperature. (4) The uncertainty of calculated soil carbon credits was much smaller compared with the uncertainty of calculated absolute carbon budgets assuming the same SOC stock uncertainty level in the inputs. Specifically, in our assessment comparing planting cover crops vs no cover crop, the uncertainty of calculated soil carbon credits induced by initial SOC% uncertainty was less than 4% (relative to the quantified value of the soil carbon credits) for 90% of the cases. Our analysis highlights that high accuracy measurement of SOC% as inputs is needed for the calculation of cropland carbon budgets; however, soil carbon credit quantification is much less sensitive to the initial SOC% inputs, and the current publicly available soil datasets (e.g., gSSURGO) are largely suitable for the calculation of soil

54 ENVIRONMENTAL SCIENCES↗

The fully coupled regionally refined model of E3SM version 2: overview of the atmosphere, land, and river results

Abstract. This paper provides an overview of the United States (US) Department of Energy's (DOE's) Energy Exascale Earth System Model version 2 (E3SMv2) fully coupled regionally refined model (RRM) and documents the overall atmosphere, land, and river results from the Coupled Model Intercomparison Project 6 (CMIP6) DECK (Diagnosis, Evaluation, and Characterization of Klima) and historical simulations – a first-of-its-kind set of climate production simulations using RRM. The North American (NA) RRM (NARRM) is developed as the high-resolution configuration of E3SMv2 with the primary goal of more explicitly addressing DOE's mission needs regarding impacts to the US energy sector facing Earth system changes. The NARRM features finer horizontal resolution grids centered over NA, consisting of 25→100 km atmosphere and land, a 0.125∘ river-routing model, and 14→60 km ocean and sea ice. By design, the computational cost of NARRM is ∼3× of the uniform low-resolution (LR) model at 100 km but only ∼ 10 %–20 % of a globally uniform high-resolution model at 25 km. A novel hybrid time step strategy for the atmosphere is key for NARRM to achieve improved climate simulation fidelity within the high-resolution patch without sacrificing the overall global performance. The global climate, including climatology, time series, sensitivity, and feedback, is confirmed to be largely identical between NARRM and LR as quantified with typical climate metrics. Over the refined NA area, NARRM is generally superior to LR, including for precipitation and clouds over the contiguous US (CONUS), summertime marine stratocumulus clouds off the coast of California, liquid and ice phase clouds near the North Pole region, extratropical cyclones, and spatial variability in land hydrological processes. The improvements over land are related to the better-resolved topography in NARRM, whereas those over ocean are attributable to the improved air–sea interactions with finer grids for both atmosphere and ocean and sea ice. Some features appear insensitive to the resolution change analyzed here, for instance the diurnal propagation of organized mesoscale convective systems over CONUS and the warm-season land–atmosphere coupling at the southern Great Plains. In summary, our study presents a realistically efficient approach to leverage the fully coupled RRM framework for a standard Earth system model release and high-resolution climate production simulations.

54 ENVIRONMENTAL SCIENCES↗

Yesterday’s extremes, today’s new normal: flood risk in the Kathmandu Valley, Nepal

Unplanned urban growth has left many cities increasingly vulnerable to extreme rainfall events, particularly in regions with inadequate drainage infrastructures and development encroaching on natural floodplains. Here, in this perspective paper, we examine the September 2024 floods that struck Central Nepal, triggered by a persistent low-pressure system and enhanced by converging moisture flows from the Arabian Sea and the Bay of Bengal which led to widespread catastrophic damage. In the Kathmandu Valley, floodwaters expanded to more than 2.5 times the bankfull water extent, causing significant damage to housing, transportation network, and critical infrastructure, displacing thousands of residents, and severely disrupting urban services. This event highlights the urgent need for improved flood management strategies that integrate both structural and non-structural measures into the infrastructure development. While early warning systems provided critical lead time, challenges remain in reducing forecasting uncertainties and improving communication across government agencies and with local communities. A forward-looking approach is essential, including probabilistic flood forecasting systems, sustainable floodplain management, risk-sensitive land use planning, climate- and disaster- resilient infrastructure development, and the integration of nature-based solutions like urban green and blue spaces to mitigate flood impacts. By involving local communities in planning and preparedness efforts, particularly through citizen science initiatives, and engagement with underserved and disadvantaged communities, Nepal can better adapt to the growing risks posed by extreme rainfall and urban flooding and enhance long-term disaster resilience in rapidly urbanizing areas like Kathmandu Valley.

Kathmandu Valley↗

Residential energy demand, emissions, and expenditures at regional and income-decile level for alternative futures

Income and its distribution profile are important determinants of residential energy demand and carry direct implications for human well-being and climate. We explore the sensitivity of residential energy systems to income growth and distribution across shared socioeconomic pathway-representative concentration pathways scenarios using a global, integrated, multisector dynamics model, Global Change Analysis Model, which tracks national/regional household energy services and fuel choice by income decile. Nation/region energy use patterns across deciles tend to converge over time with aggregate income growth, as higher-income consumers approach satiation levels in floorspace and energy services. However, in some regions, existing within-region inequalities in energy consumption persist over time due to slow income growth in lower income groups. Due to continued differences in fuel types, lower income groups will have higher exposure to household air pollution, despite lower contributions to greenhouse gas emissions. We also find that the share of income dedicated to energy is higher for lower deciles, with strong regional differences.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Tropopause Characteristics Based on Long-Term ARM Radiosonde Data: A Fine-Scale Comparison at the Extratropical SGP Site and Arctic NSA Site

The variations in the characteristics of the tropopause are sensitive indicators for the climate system and climate change. By using Atmospheric Radiation Measurement (ARM) radiosonde data that were recorded at the extratropical Southern Great Plains (SGP) and Arctic North Slope of Alaska (NSA) sites over an 18-year period (January 2003 to December 2020), this study performs a fine-scale comparison of the climatological tropopause features between these two sites that are characterized by different climates. The static stability increases rapidly above the tropopause at both sites, indicating the widespread existence of a tropopause inversion layer. The structures of both the tropopause inversion layer and the stability transition layer are more obvious at NSA than at SGP, and the seasonal variation trends of the tropopause inversion layer and stability transition layer are distinctly different between the two sites. A fitting method was used to derive the fitted tropopause height and tropopause sharpness (λ). Although this fitting method may determine a secondary tropopause rather than the primary tropopause when multiple tropopause heights are identified on one radiosonde profile, the fitted tropopause heights generally agree well with the observed tropopause heights. Broad tropopause sharpness values (λ > 2 km) occur more frequently at SGP than at NSA, resulting in a greater average tropopause sharpness at SGP (1.0 km) than at NSA (0.6 km). Significant positive trends are exhibited by the tropopause heights over the two sites, with rates of increase of 23.7 ± 6.5 m yr ⁻1 at SGP and 28.0 ± 4.0 m yr ⁻1 at NSA during the study period.

54 ENVIRONMENTAL SCIENCES↗

Global Sensitivity Analysis Using the Ultra‐Low Resolution Energy Exascale Earth System Model

Abstract For decades, Arctic temperatures have increased twice as fast as average global temperatures. As a first step toward quantifying parametric uncertainty in Arctic climate, we performed a variance‐based global sensitivity analysis (GSA) using a fully coupled, ultra‐low resolution (ULR) configuration of version 1 of the U.S. Department of Energy's Energy Exascale Earth System Model (E3SMv1). Specifically, we quantified the sensitivity of six quantities of interests (QOIs), which characterize changes in Arctic climate over a 75 year period, to uncertainties in nine model parameters spanning the sea ice, atmosphere, and ocean components of E3SMv1. Sensitivity indices for each QOI were computed with a Gaussian process emulator using 139 random realizations of the random parameters and fixed preindustrial forcing. Uncertainties in the atmospheric parameters in the Cloud Layers Unified by Binormals (CLUBB) scheme were found to have the most impact on sea ice status and the larger Arctic climate. Our results demonstrate the importance of conducting sensitivity analyses with fully coupled climate models. The ULR configuration makes such studies computationally feasible today due to its low computational cost. When advances in computational power and modeling algorithms enable the tractable use of higher‐resolution models, our results will provide a baseline that can quantify the impact of model resolution on the accuracy of sensitivity indices. Moreover, the confidence intervals provided by our study, which we used to quantify the impact of the number of model evaluations on the accuracy of sensitivity estimates, have the potential to inform the computational resources needed for future sensitivity studies.

54 ENVIRONMENTAL SCIENCES↗

Local summer insolation and greenhouse gas forcing drove warming and glacier retreat in New Zealand during the Holocene

Geological climate archives from the Holocene Epoch provide baseline information concerning natural climate variability. Temperate mountain glacier extent is sensitive to summer air temperature, thus geological records of past glacier length changes are a useful proxy for this climatic variable. Here we present a new cosmogenic 10 Be chronology of glacier length changes at Dart Glacier in the Southern Alps, New Zealand. Prominent moraines deposited 321 ± 44 yr ago (n = 11) and 7.8 ± 0.3 ka (n = 5) show glaciers during the Little Ice Age were less extensive than during the early Holocene. This pattern of net Holocene glacier retreat is consistent with emerging data from other catchments in New Zealand and across the southern mid-latitudes. Using the physical framework of a transient global climate model simulation, we suggest that cool summers in the early Holocene were promoted by the local summer insolation minimum, together with low atmospheric greenhouse gas concentrations, causing an early Holocene austral glacial maximum. An insolation-driven reduction in seasonality at southern mid-latitudes may reconcile differences between early Holocene temperature reconstructions where climate proxies have different seasonal sensitivities. Furthermore, we suggest that rising greenhouse gas concentrations after 7 ka caused regional-scale glacier retreat and appear to be the dominant driver of multi-millennial summer temperature trends in the southern mid-latitudes during the present interglacial.

58 GEOSCIENCES↗

Increasing sensitivity of dryland vegetation greenness to precipitation due to rising atmospheric CO 2

Water availability plays a critical role in shaping terrestrial ecosystems, particularly in low- and mid-latitude regions. The sensitivity of vegetation growth to precipitation strongly regulates global vegetation dynamics and their responses to drought, yet sensitivity changes in response to climate change remain poorly understood. Here we use long-term satellite observations combined with a dynamic statistical learning approach to examine changes in the sensitivity of vegetation greenness to precipitation over the past four decades. We observe a robust increase in precipitation sensitivity (0.624% yr –1 ) for drylands, and a decrease (–0.618% yr –1 ) for wet regions. Using model simulations, we show that the contrasting trends between dry and wet regions are caused by elevated atmospheric CO 2 (eCO 2 ). eCO 2 universally decreases the precipitation sensitivity by reducing leaf-level transpiration, particularly in wet regions. However, in drylands, this leaf-level transpiration reduction is overridden at the canopy scale by a large proportional increase in leaf area. The increased sensitivity for global drylands implies a potential decrease in ecosystem stability and greater impacts of droughts in these vulnerable ecosystems under continued global change.

54 ENVIRONMENTAL SCIENCES↗

Numerically Relevant Timescales in the MG2 Microphysics Model

Climate models rely on parameterizations of a variety of processes in the atmospheric physics, but a common concern is that the temporal resolution is too coarse to consistently resolve the behavior that individual parameterizations are designed to capture. This study examines timescales numerically derived from the Morrison-Gettelman (MG2) microphysics as implemented within the Energy Exascale Earth System Model, version 1 (E3SMv1). Numerically relevant timescales in MG2 are derived by computing the eigenspectrum of its Jacobian. These timescales are found to often be smaller than the default 5 min time step used for MG2. The fast timescales are then heuristically connected to individual microphysics processes. By substepping a few particular rain processes within MG2, the time discretization error for those processes was considerably reduced with minimal additional expense to the overall microphysics. While this improvement has a substantial effect on the target processes and on the vertical distribution of stratiform-derived rain within E3SMv1, the overall model climate is found to not be sensitive to MG2 time step. Here, we hypothesize that this is because the surface climate does not depend strongly on certain process rates, especially MG2's rain evaporation rate.

58 GEOSCIENCES↗

Sensitivity of Arctic Clouds to Ice Microphysical Processes in the NorESM2 Climate Model

Abstract Ice formation remains one of the most poorly represented microphysical processes in climate models. While primary ice production (PIP) parameterizations are known to have a large influence on the modeled cloud properties, the representation of secondary ice production (SIP) is incomplete and its corresponding impact is therefore largely unquantified. Furthermore, ice aggregation is another important process for the total cloud ice budget, which also remains largely unconstrained. In this study, we examine the impact of PIP, SIP, and ice aggregation on Arctic clouds, using the Norwegian Earth System Model, version 2 (NorESM2). Simulations with both prognostic and diagnostic PIP show that heterogeneous freezing alone cannot reproduce the observed cloud ice content. The implementation of missing SIP mechanisms (collisional breakup, drop shattering, and sublimation breakup) in NorESM2 improves the modeled ice properties, while improvements in liquid content occur only in simulations with prognostic PIP. However, results are sensitive to the description of collisional breakup. This mechanism, which dominates SIP in the examined conditions, is very sensitive to the treatment of the sublimation correction factor, a poorly constrained parameter that is included in the utilized parameterization. Finally, variations in ice aggregation treatment can also significantly impact cloud properties, mainly through their impact on collisional breakup efficiency. Overall, enhancement in ice production through the addition of SIP mechanisms and the reduction in ice aggregation (in line with radar observations of shallow Arctic clouds) result in enhanced cloud cover and decreased TOA radiation biases, compared to satellite measurements, especially during the cold months. Significance Statement Arctic clouds remain a large source of uncertainty in projections of the future climate due to the poor representation of the microphysical processes that govern their life cycle. Ice formation is among the least understood processes. While it is widely recognized that better constraints on primary ice production (PIP) are needed to improve existing parameterizations, we show that secondary ice production (SIP) and ice aggregation can have also a significant impact on ice number concentrations. Constraining ice formation through the addition of missing SIP mechanisms and reducing ice aggregation can improve the representation of the cloud macrophysical properties and enhance total cloud cover in the Arctic region, which in turn contributes to decreased TOA radiation biases in the cold months.

Meteorology & Atmospheric Sciences↗

Modeled Results of Four Residential Energy Efficiency Measure Packages for Deriving Advanced Building Construction Research Targets

The Advanced Building Construction (ABC) Initiative from the U.S. Department of Energy Building Technologies Office is working to accelerate industrialized construction innovations for decarbonizing buildings. To inform performance and cost targets for research under the ABC Initiative, this analysis used the ResStock™ tool to evaluate the energy savings, utility bill impacts, and carbon emissions impacts of four simulated upgrade packages with specific target performance levels on a large sample of residential dwelling units (approximately 550,000) representative of the U.S. housing stock.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Stomatal closure as a driver of minimum leaf conductance declines at high temperature and vapor pressure deficit in Quercus

Abstract Rising global temperatures and vapor pressure deficits (VPDs) are increasing plant water demand and becoming major drivers of large-scale plant mortality. Controlling transient leaf water loss after stomatal closure (minimum stomatal conductance [gmin]) is recognized as a key trait determining how long plants survive during soil drought. Yet, substantial uncertainty remains regarding how gmin responds to elevated temperatures and VPD and the underlying mechanisms. We measured gmin in 24 Quercus species from temperate and Mediterranean climates to determine whether gmin was sensitive to a coupled temperature and VPD increase. We also explored mechanistic links to phenology, climate, evolutionary history, and leaf anatomy. We found that gmin in all species exhibited a nonlinear negative temperature and VPD dependence. At 25 °C (VPD = 2.2 kPa), gmin varied from 1.19 to 8.09 mmol m−2 s−1 across species but converged to 0.57 ± 0.06 mmol m−2 s−1 at 45 °C (VPD = 6.6 kPa). In a subset of species, the effect of temperature and VPD on gmin was reversible and linked to the degree of stomatal closure, which was greater at 45 °C than at 25 °C. Our results show that gmin is dependent on temperature and VPD, is highly conserved in Quercus species, and is linked to leaf anatomy and stomatal behavior.

Zailaa, Joseph (ORCID:000000019103190X)↗

Assessing Global Ocean Circulation Sensitivity to Changes in Antarctic Ice Discharge

The principal goal of this research is to assess how projected increases in meltwater discharge and iceberg calving from the Antarctic Ice Sheet will impact climate in the next 250 years. In short, the sensitivity of the global ocean circulation and climate, and in particular the Atlantic Meridional Overturning Circulation (AMOC), to large increases in iceberg calving and meltwater discharges to the Southern Ocean are poorly understood. The requirement to investigate this topic is heightened by growing evidence that the West Antarctic Ice Sheet (WAIS) is vulnerable to rapid retreat and collapse through marine-ice-sheet-instability (MISI) on multidecadal-to-centennial timescales. Future model simulations of the AIS show the potential for rapid ice sheet retreat in the next 50 – 300 years, with such a collapse associated with the discharge of enormous volumes of ice and meltwater to the Southern Ocean. This project uses a combination of an ice sheet/shelf model and a fully coupled (ocean-atmosphere-seaice) model (CESM 1.2) to explore the sensitive of global climate to future changes in runoff from the Antarctic continent under future greenhouse gas emissions scenarios IPCC RCP4.5 and 8.5. Overall, accounting for Antarctic meltwater reduces the projected global mean anthropogenic greenhouse warming by up to 2 deg C. during peak WAIS collapse and by ~1 deg. for most of the simulation. In addition, the increased meltwater contribution raises sub-surface ocean temperatures at the Antarctic ice sheet margin by more than 1ºC, such that the ice shelf basal melt rates increase beyond those currently projected. In contrast, the model predicts a strong cooling of both the surface air and ocean temperatures over the Southern Ocean that results in a massive expansion and thickening of the perennial sea-ice. The slow-down in anthropogenic warming also results in the onset of ice-free conditions in the Arctic being delayed in the model by several decades, and the projected slowdown in the AMOC to reduce.

58 GEOSCIENCES↗

Metabolic responses to increased temperatures in three semi-aquatic turtle species from the southeastern United States

Physiological maintenance of ectotherms is largely dependent on temperature. Abrupt changes in environmental conditions can cause shifts in energy budgets, with higher temperatures requiring more energy be devoted to maintenance functions, making less energy available for other activities. As a group, freshwater turtles are of conservation concern; thus, it is important to understand their physiological responses to temperature shifts, particularly given current climate change predictions. Here, we quantified metabolic rates (MR) for three species of semi-aquatic turtles that vary in their degree of terrestriality and propensity to aestivate: chicken turtles (Deirochelys reticularia), eastern mud turtles (Kinosternon subrubrum), and eastern musk turtles (Sternotherus odoratus). Metabolic trials were performed using a flow-through respirometer at three environmentally-relevant temperatures: 25°C, 30°C, and 35°C. As expected, MRs of turtles were significantly and positively associated with increases in mass and temperature. Pairwise comparisons revealed that D. reticularia exhibited significantly higher oxygen consumption rates (VO 2 ) relative to K. subrubrum, which had significantly higher VO 2 than S. odoratus. The combination of higher metabolic rates, the tendency to terrestrially aestivate, and declining quantity of high quality thermal habitats may place some semi-aquatic turtles at a disadvantage when considering future climate scenarios. Our findings, when combined with knowledge of other ecological traits, may be a useful tool for predicting relative sensitivities of turtle species to predicted climate change, particularly those species inhabiting seasonal wetlands.

Sternotherus odoratus↗

Nitrate loading projection is sensitive to freeze-thaw cycle representation

Climate change can have substantial impacts on nitrogen runoff, which is a major cause of eutrophication, harmful algal blooms, and hypoxia in freshwaters and coastal regions. We examined responses of nitrate loading to climate change in the Upper Mississippi River Basin (UMRB) with an enhanced Soil and Water Assessment Tool with physically based Freeze-Thaw cycle representation (SWAT-FT), as compared with the original SWAT model that employs an empirical equation. Driven by future climate projections from five General Circulation Models (GCMs) from 1960 to 2099 under the Representative Concentrations Pathways (RCP) 8.5 scenario, we analyzed changes in riverine nitrate loadings, as well as terrestrial surface and subsurface contributions of the UMRB in the 21st century relative to the baseline period of 1960–1999. By the end of the 21st century, the original SWAT model predicted about a 50% increase in riverine nitrate loadings which is nearly twice as much as that estimated by SWAT-FT (ca. 25%). Such a large difference in projected nitrate changes can potentially mislead mitigation strategies that aim to reduce nitrogen runoff from the UMRB. Further analysis shows that the difference between the original SWAT model and SWAT-FT led to substantial discrepancies in the spatial distribution of surface and subsurface nitrate loadings in the UMRB. In general, SWAT-FT predicted more nitrate leaching for northwestern parts of the UMRB which are more sensitive to freeze-thaw cycle, mainly because SWAT-FT simulated less frequent frozen soils. This study highlights the importance of using physically based freeze-thaw cycle representation in water quality modeling. Design of future nitrogen runoff reduction strategies should include careful assessment of effects that land management has on the freeze-thaw cycles to provide reliable projection of water quality under climate change.

Wang, Qianfeng↗

Projecting Permafrost Thaw of Sub-Arctic Tundra With a Thermodynamic Model Calibrated to Site Measurements

Northern circumpolar permafrost thaw affects global carbon cycling, as large amounts of stored soil carbon becomes accessible to microbial breakdown under a warming climate. The magnitude of carbon release is linked to the extent of permafrost thaw, which is locally variable and controlled by soil thermodynamics. Soil thermodynamic properties, such as thermal diffusivity, govern the reactivity of the soil-atmosphere thermal gradient, and are controlled by soil composition and drainage. In order to project permafrost thaw for an Alaskan tundra experimental site, we used seven years of site data to calibrate a soil thermodynamic model using a data assimilation technique. The model reproduced seasonal and interannual temperature dynamics for shallow (5–40 cm) and deep soil layers (2–4 m), and simulations of seasonal thaw depth closely matched observed data. The model was then used to project permafrost thaw at the site to the year 2100 using climate forcing data for three future climate scenarios (RCP 4.5, 6.0, and 8.5). Minimal permafrost thawing occurred until mean annual air temperatures rose above the freezing point, after which we measured over a 1 m increase in thaw depth for every 1 °C rise in mean annual air temperature. Under no projected warming scenario was permafrost remaining in the upper 3 m of soil by 2100. We demonstrated an effective data assimilation method that optimizes parameterization of a soil thermodynamic model. The sensitivity of local permafrost to climate warming illustrates the vulnerability of sub-Arctic tundra ecosystems to significant and rapid soil thawing.

54 ENVIRONMENTAL SCIENCES↗