Search NASA⌕ Search

SEARCH · Search NASA

Results for “Meteorologic data”

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

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

At least 307 records · Page 17

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL National Security Test Range

This assessment uses screening level models to calculate potential environmental impacts from proposed tests at two locations at the Idaho National Laboratory (INL) National Security Test Range (NSTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K 2 O, LaBr 3 , KBr, Cu, Zr, F, Ga, Ga 2 O 3 , NaNO 2 , Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface soil impacts were calculated with the computer program Mixing-Cell Model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated based on the assumption that 12 tests are conducted annually using all 11 material types for 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer code input, output, executable files, batch files, scripts, and spreadsheet files, are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS).

99 GENERAL AND MISCELLANEOUS↗

Super Resolution for Renewable Energy Resource Data With Wind From Reanalysis Data (Sup3rWind) and Application to Ukraine [Slides]

In this work we present a novel deep learning-based downscaling method, using generative adversarial networks (GANs), for generating high-resolution wind resource data from ECMWF Reanalysis v5 data (ERA5). We show that by training a GAN model on ERA5, as opposed to coarsened high-resolution data, we achieve results that are competitive with conventional dynamical downscaling. This GAN-based downscaling method additionally reduces computational costs over dynamical downscaling by two orders of magnitude. All GANs are trained on data sampled from CONUS, selected to provide a diverse sampling of terrain conditions, and validated on observational data along with data held out from training. This cross-validation shows low error and high correlations with observations and excellent agreement with hold out data across physical distributions. Our approach is finally used to downscale 30km hourly ERA5 to 2-km 5-minute wind data, for January 2000 through December 2023, at multiple hub heights, over Ukraine, Moldova, and part of Romania. Comparisons against observational data from Meteorological Assimilation Data Ingest System (MADIS) and multiple wind farms show the same level of performance as for CONUS validation. This 24 year data record is the first member of the "super resolution for renewable energy resource data with wind from reanalysis data" dataset (Sup3rWind).

17 WIND ENERGY↗

Building research capacity in Southern Appalachian mountain wetlands

This project increased research capacity in Southern Appalachian wetlands by facilitating new collaborations between the PI (UNC Asheville) and DOE and academic research partners. Mountain wetlands are among the rarest natural communities in the Southern Appalachians; most of the sites contain rare or endangered species and have been incorporated into or are being considered for the new Mountain Bogs National Wildlife Refuge (established in 2015). PI Wilcox’s research since 2010 has focused on wetland hydrology at over a dozen priority field sites in western North Carolina, eastern Tennessee, and northern Georgia. This project provided time and funding to learn about DOE programs and resources, including meteorological data, watershed modeling programs, analytical capabilities, and computing facilities. PI Wilcox was able to meet potential research partners at professional conferences and the DOE-BER-ESS PI meeting and invite several of those potential partners to visit UNC Asheville and nearby field sites in Western North Carolina. These connections have resulted in three new collaborative research proposals, and a new cooperative agreement with the US Fish and Wildlife Service that includes two new field sites.

54 ENVIRONMENTAL SCIENCES↗

The U.S. Agrivoltaic Shading Tool: A National-Scale Interface for Modeling Light and Shade Patterns in Ten Common Agrivoltaic Configurations

Agrivoltaic systems are dual-use configurations that co-locate agriculture and photovoltaic (PV) infrastructure and require careful design to balance crop performance and energy generation. A critical element of agrivoltaic design is the spatial and temporal distribution of irradiance and shade within and around PV arrays. To support research, planning, and stakeholder decision-making, we introduce the U.S. Agrivoltaic Shading Tool, a novel web-based application that delivers high-resolution irradiance and photosynthetically active radiation (PAR) modeling for ten standardized PV configurations across the conterminous United States. The tool leverages the National Laboratory of the Rockies (NLR) System Advisor Model (SAM) to perform detailed irradiance simulations, using meteorological data from the National Solar Radiation Database (NSRDB). Outputs include seasonal, monthly, weekly, and diurnal patterns of available sunlight, amount of shade, irradiance, and PAR at ground level within agrivoltaic system footprints. For a user's selected location, these results are visualized through interactive visualizations, heatmaps, and time-series plots, designed to be accessible to both technical and non-technical users. In addition to facilitating rapid spatial exploration of agrivoltaic light environments, the tool will offer seamless integration with the InSPIRE Agrivoltaics Design and Analysis Model (ADAM). This optional workflow will allow users to port selected site and configuration parameters into a more advanced modeling environment for further customization of structural layouts, crop-system compatibility, power generation, and technoeconomic performance. Finally, to promote open science, the entire dataset will be hosted and available for open access through the OpenEI platform. By standardizing and disseminating high-quality irradiance data and design tools, the U.S. Agrivoltaic Shading Tool supports a wide range of users, including researchers, landowners, energy developers, and policymakers, in evaluating the agronomic and energetic feasibility of agrivoltaic systems across the United States.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

National Park Air Quality Index Dataset

The National Park Air Quality Index dataset (NPS-AQI) consists of webcam images taken from the National Park Service's publicly available air quality web cameras and associated measurements for air pollutants, AQI, and meteorological data obtained via the publicly available NPS Gaseous Pollutant Monitoring Program. The full dataset is a collection of 146,822 images paired with air quality measurements. The specific measurements reported are: ozone ppm, 8-hour running average ozone ppm, so2 ppm, AQI (derived from ozone), temperature, and humidity. The images are 1500X1000 pixel PNG files arranged into folders by NPS site and named according to the time and date the image was taken. There are three CSV files (representing "training", "validation", and "testing" images splits) containing image names and associated NPS site names, air pollutant measurements, and meteorlogical data.

Svinth, Christian N↗

Similarity of the turbulent kinetic energy dissipation rate distribution in the upper mixed layer of the tropical Indian Ocean

Turbulence within the upper ocean mixed layer plays a key role in various physical, biological, and chemical processes. Between September and November 2011, a dataset of 570 vertical profiles of the turbulent kinetic energy (TKE) dissipation rate, as well as conventional hydrological and meteorological data, were collected in the upper layer of the tropical Indian Ocean. These data were used to statistically analyze the vertical distribution of the TKE dissipation rate in the mixed layer. The arithmetic-mean method made the statistical TKE dissipation rate profile more scattered than the median and geometric-mean methods. The statistical TKE dissipation rate were respectively scaled by the surface buoyancy flux and the TKE dissipation rate at the mixed-layer base. It was found that the TKE dissipation rate scaled by that at the mixed-layer base exhibited better similarity characteristics than that scaled by the surface buoyancy flux, whether the stability parameter D/|L MO | was greater or less than 10, indicating that the TKE dissipation rate at the mixed-layer base is a better characteristic scaling parameter for reflecting the intrinsic structure of the TKE dissipation rate in the mixed layer, where D and L MO are respectively the mixed-layer thickness and the Monin-Obukhov length scale. The parameterization of the TKE dissipation rate at the mixed-layer base on the shear-driven dissipation rate and the surface buoyancy flux was further explored. It was found that the TKE dissipation rate at the mixed-layer base could be well fitted by a linear combination of three terms: the wind-shear-driven dissipation rate, the surface buoyancy flux, and a simple nonlinear coupling term of these two.

54 ENVIRONMENTAL SCIENCES↗

Assessing shellfish water exposure to fecal bacteria pollution in Salish Sea: three-dimensional modeling and implications for monitoring

Fecal bacteria (FB) contamination poses significant risks to shellfish safety and management in coastal and estuarine waters. Despite extensive pollution identification and correction efforts, FB contamination in shellfish-growing areas persists in the Salish Sea, highlighting the need to identify overlooked sources and better understand FB transport from riverine and shoreline inputs to shellfish beds. To address this, a high-resolution three-dimensional hydrodynamic model coupled with FB kinetics was developed and applied to a case study site in Salish Sea—Portage Bay—to simulate freshwater plume circulation, flushing dynamics, and bacterial transport. Daily FB loading from the major freshwater inflow—Nooksack River was generated by both linear interpolation and integrating a machine learning approach (XGBoost), trained on historical hydrological and meteorological data. The model successfully reproduced both the magnitude and seasonal variation of FB concentrations in Portage Bay for the year of 2021, demonstrating that simplified FB kinetics with first-order decay due to mortality was effective in this dynamic coastal environment with short flushing time. Model results identified the Nooksack River as the dominant far-field FB source, while scenario simulations showed that near-field coastal stormwater outfalls elevated local FB levels following rainfall, particularly under low-flow conditions. The XGBoost prediction provided comparable or superior accuracy to linear interpolation, particularly during periods of missing observational data, by capturing short-term variability and event-driven loading more effectively. Integrating data-driven riverine FB inputs with mechanistic coastal numerical modeling provides a robust framework for operational forecasting of shellfish bed exposure risk and supports adaptive monitoring and management of shellfish growing areas in the Salish Sea and similar coastal systems.

Salish Sea↗

Impacts of a Prescribed Fire on Air Quality in Central New Mexico

A short-duration but high-impact air quality event occurred on 28 November 2018 along the Rio Grande Valley of New Mexico. This fire occurred outside the typical wildfire season, and greatly impacted the air quality in Socorro, NM, and the surroundings. Measurements were taken during the event using an aerosol light scattering technique (integrating nephelometer) and a particulate mass concentration monitor (DustTrak PM optical monitor). The instruments sampled the ambient air during the event on the campus of the New Mexico Institute of Mining and Technology in Socorro, New Mexico. The peak values on a 5-min basis of light scattering and the PM mass concentration reached 470 Mm -1 and 270 µg/m 3 , respectively. We examined the meteorological context of the event using local meteorological data and back trajectories using the NOAA HYSPLIT model to determine atmospheric transport and possible sources. Several fires, both prescribed and wildfires, occurred in the region including a prescribed burn at Bosque del Apache National Wildlife Refuge (17 km south-southeast of the receptor site). The data suggest that the prescribed burn at Bosque del Apache was the dominant contributor due to transport evidence and the event’s narrow spatiotemporal extent. The increasing importance of restoring ecosystem function using prescribed fire in wildland fire management will likely lead to more frequent air quality impacts and sets up policy tradeoffs that require a balance between these public goals. This study examines the evidence of the effects of a prescribed fire in a protected wildland area impacting the air quality in a nearby populated area.

54 ENVIRONMENTAL SCIENCES↗

Computational Fluid Dynamics-Based Modeling of Methane Flows Around Oil and Gas Equipment

Recent studies estimate that emissions from oil and gas production facilities contribute between 20 and 50% of the total methane ( CH 4 ) emitted in the US; therefore, quantifying and reducing these emissions are crucial for achieving climate goals. Methane quantification depends on both measuring methane concentrations and converting them to emissions through a modeling framework. Currently, simple atmospheric dispersion models are primarily used to quantify emissions and concentrations, but these estimates are highly uncertain when quantifying emissions from complex aerodynamic sources, such as oil and gas facilities. This investigation used a CFD modeling approach, which can account for aerodynamic complexity but has hitherto not been used to model methane concentrations downwind of a methane release of a known rate, and compared it against in situ measurements. High-time-resolution (1 Hz) methane concentration and meteorological data were measured during experiments conducted at the METEC on 21 March and 11 July 2024. The METEC site configuration, measured wind data, and controlled emission rates were used as input for the CONVERGE CFD model to model downwind CH 4 concentration. The modeling was carried out between 20 and 70 m, from two different points of release in two separate controlled-release experiments, one from a separator and another from a wellhead. In these experiments, we found that the CFD model could predict the CH 4 concentrations downwind of the release to a good degree. The model was evaluated on multiple metrics to assess its performance in estimating methane concentrations at typical fence line distances (∼30 m). These results help us to understand external flows and the ability of CFD models to predict downwind concentrations in aerodynamically complex environments.

03 NATURAL GAS↗

Comparative Performance of Gaussian Plume and Backward Lagrangian Stochastic Models for Near-Field Methane Emission Estimation Using a Single Controlled Release Experiment

Methane (CH 4 ) is a major component of natural gas and a potent greenhouse gas. Increasing atmospheric methane concentrations are attributed to emissive anthropogenic activities by an average of 13 ppb per yr since 2020 and are linked to a changing global climate. Mitigating CH 4 emissions from oil and gas production sites has recently become a target to reduce overall greenhouse gas emissions; however, monitoring the efficacy of mitigation strategies depends on accurate quantification of CH 4 emissions at the facility-level. Near-field quantification of methane (CH 4 ) emissions from oil and gas (O&G) facilities remains challenging due to the effects of atmospheric variability and sensor configuration on atmospheric dispersion models. This study evaluates the performance of two atmospheric dispersion models, the Gaussian plume (GP) and backward Lagrangian stochastic (bLS), by comparing calculated CH 4 emissions to controlled single-point emissions between 0.4 and 5.2 kg CH 4 h −1 . Emissions were calculated by both models using 121 individual sets of measurements comprising five-minute averaged downwind methane mixing ratios and matching meteorological data. The comparison shows that the bLS approach achieved a higher proportion of emission estimates within a factor of two (FAC2) of the known emission rates compared to the GP approach. The emissions calculated by the bLS model also had a lower multiplicative error and reduced bias relative to GP. Other error-based metrics further confirmed the bLS model performed better, as it yielded lower RMSE and MAE than GP. Statistical analysis of the emission data shows that the lateral and vertical alignment of the source and the sensor plays a critical role in emission estimations, as measurements made closer to the plume centerline and at a distance between 40 and 80 m downwind yielded the best FAC2 agreement. High wind meander degraded the ability of both approaches to generate representative emissions, particularly with the GP approach, as it violates the modeling approach’s assumption of steady-state emissions. Data suggest emissions calculated by the bLS model are comprehensively in better agreement, but the computational demands of the modeling approach and integration into fenceline systems limit real-time applicability. While these results provide insight into model performance under controlled near-field conditions, their applicability to more complex or heterogeneous oil and gas production environments (e.g., the regions Marcellus or Unita Basins) remains limited and uncertain.

gaussian plume↗

Synchronicity from synchronized chaos

The synchronization of loosely-coupled chaotic oscillators, a phenomenon investigated intensively for the last two decades, may realize the philosophical concept of “synchronicity”—the commonplace notion that related events mysteriously occur at the same time. When extended to continuous media and/or large discrete arrays, and when general (non-identical) correspondences are considered between states, intermittent synchronous relationships indeed become ubiquitous. Meaningful synchronicity follows naturally if meaningful events are identified with coherent structures, defined by internal synchronization between remote degrees of freedom; a condition that has been posited as necessary for synchronizability with an external system. The important case of synchronization between mind and matter is realized if mind is analogized to a computer model, synchronizing with a sporadically observed system, as in meteorological data assimilation. Evidence for the ubiquity of synchronization is reviewed along with recent proposals that: (1) synchronization of different models of the same objective process may be an expeditious route to improved computational modeling and may also describe the functioning of conscious brains; and (2) the nonlocality in quantum phenomena implied by Bell’s theorem may be explained in a variety of deterministic (hidden variable) interpretations if the quantum world resides on a generalized synchronization “manifold”.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Solar-Tracker Diffuse-Response Algorithm for Balancing Energy Gain and Mechanical Wear in Photovoltaic Systems

Single-axis solar tracking maximizes photovoltaic energy production under clear-sky conditions; however, its effectiveness decreases under cloudy and overcast skies, where diffuse irradiance dominates and the optimal module orientation changes. Conventional tracking algorithms either neglect sky conditions or rely on simplified diffuse-response strategies that may trigger frequent tracker repositioning under variable cloud cover, leading to increased mechanical wear with marginal energy gains. This work proposes an enhanced diffuse-response tracking algorithm that explicitly accounts for both the intensity and temporal persistence of cloudiness. By requiring overcast conditions to persist for a minimum duration before reorienting the tracker to a diffuse-stow position, the proposed approach reduces unnecessary movements while preserving the benefits of diffuse-response operation. The algorithm is evaluated through numerical simulations based on historical meteorological data and validated using field measurements on monofacial and bifacial photovoltaic strings. The results show that the proposed strategy reduces excess tracker movement from 114% to 0.16% while maintaining nearly the same energy yield. Compared to a conventional diffuse-response algorithm, the associated energy reduction is minimal (≈0.17%) relative to the ≈0.37% yield gain observed at the studied location. These findings demonstrate that incorporating cloudiness duration enables a practical compromise between energy performance and tracker durability, particularly for monofacial photovoltaic systems.

14 SOLAR ENERGY↗

Upscaling Gross Primary Production in Corn-Soybean Rotation Systems in the Midwest

The Midwestern US is dominated by corn (Zea mays L.) and soybean (Glycine max [L.] Merr.) production, and the carbon dynamics of this region are dominated by these production systems. An accurate regional estimate of gross primary production (GPP) is imperative and requires upscaling approaches. The aim of this study was to upscale corn and soybean GPP (referred to as GPPcalc) in four counties in Central Iowa in the 2016 growing season (DOY 145–269). Eight eddy-covariance (EC) stations recorded carbon dioxide fluxes of corn (n = 4) and soybean (n = 4), and net ecosystem production (NEP) was partitioned into GPP and ecosystem respiration (RE). Additional field-measured NDVI was used to calculate radiation use efficiency (RUEmax). GPPcalc was calculated using 16 MODIS satellite images, ground-based RUEmax and meteorological data, and improved land use maps. Seasonal NEP, GPP, and RE ( x ¯ ± SE) were 678 ± 63, 1483 ± 100, and −805 ± 40 g C m−2 for corn, and 263 ± 40, 811 ± 53, and −548 ± 14 g C m−2 for soybean, respectively. Field-measured NDVI aligned well with MODIS fPAR (R2 = 0.99), and the calculated RUEmax was 3.24 and 1.90 g C MJ−1 for corn and soybean, respectively. The GPPcalc vs. EC-derived GPP had a RMSE of 2.24 and 2.81 g C m−2 d−1, for corn and soybean, respectively, which is an improvement to the GPPMODIS product (2.44 and 3.30 g C m−2 d−1, respectively). Corn yield, calculated from GPPcalc (12.82 ± 0.65 Mg ha−1), corresponded well to official yield data (13.09 ± 0.09 Mg ha−1), while soybean yield was overestimated (6.73 ± 0.27 vs. 4.03 ± 0.04 Mg ha−1). The approach presented has the potential to increase the accuracy of regional corn and soybean GPP and grain yield estimates by integrating field-based flux estimates with remote sensing reflectance observations and high-resolution land use maps.

Dold, Christian (ORCID:0000000260355597)↗

Determining Off-Cycle Fuel Economy Benefits of 2-Layer HVAC Technology

This work presents a methodology to determine the off-cycle fuel economy benefit of a 2-Layer HVAC system which reduces ventilation and heat rejection losses of the heater core versus a vehicle using a standard system. Experimental dynamometer tests using EPA drive cycles over a broad range of ambient temperatures were conducted on a highly instrumented 2016 Lexus RX350 (3.5L, 8 speed automatic). These tests were conducted to measure differences in engine efficiency caused by changes in engine warmup due to the 2-Layer HVAC technology in use versus the technology being disabled (disabled equals fresh air-considered as the standard technology baseline). These experimental datasets were used to develop simplified response surface and lumped capacitance vehicle thermal models predictive of vehicle efficiency as a function of thermal state. These vehicle models were integrated into a database of measured on road testing and coupled with U.S. typical meteorological data to simulate vehicle efficiency across seasonal thermal and operational conditions for hundreds of thousands of drive cycles. Fuel economy benefits utilizing the 2-Layer HVAC technology are presented in addition to goodness of fit statistics of the modeling approach relative to the experimental test data.

2-layer technology↗

Evaluation of column-averaged methane in models and TCCON with a focus on the stratosphere

The distribution of methane (CH 4 ) in the stratosphere can be a major driver of spatial variability in the dry-air column-averaged CH 4 mixing ratio (XCH 4 ), which is being measured increasingly for the assessment of CH 4 surface emissions. Chemistry-transport models (CTMs) therefore need to simulate the tropospheric and stratospheric fractional columns of XCH 4 accurately for estimating surface emissions from XCH 4 . Simulations from three CTMs are tested against XCH 4 observations from the Total Carbon Column Network (TCCON). We analyze how the model–TCCON agreement in XCH 4 depends on the model representation of stratospheric CH 4 distributions. Model equivalents of TCCON XCH 4 are computed with stratospheric CH 4 fields from both the model simulations and from satellite-based CH 4 distributions from MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) and MIPAS CH 4 fields adjusted to ACE-FTS (Atmospheric Chemistry Experiment Fourier Transform Spectrometer) observations. Using MIPAS-based stratospheric CH 4 fields in place of model simulations improves the model–TCCON XCH 4 agreement for all models. For the Atmospheric Chemistry Transport Model (ACTM) the average XCH 4 bias is significantly reduced from 38.1 to 13.7 ppb, whereas small improvements are found for the models TM5 (Transport Model, version 5; from 8.7 to 4.3 ppb) and LMDz (Laboratoire de Météorologie Dynamique model with zooming capability; from 6.8 to 4.3 ppb). Replacing model simulations with MIPAS stratospheric CH 4 fields adjusted to ACE-FTS reduces the average XCH 4 bias for ACTM (3.3 ppb), but increases the average XCH 4 bias for TM5 (10.8 ppb) and LMDz (20.0 ppb). These findings imply that model errors in simulating stratospheric CH 4 contribute to model biases. Current satellite instruments cannot definitively measure stratospheric CH 4 to sufficient accuracy to eliminate these biases. Applying transport diagnostics to the models indicates that model-to-model differences in the simulation of stratospheric transport, notably the age of stratospheric air, can largely explain the inter-model spread in stratospheric CH 4 and, hence, its contribution to XCH 4 . Furthermore, it would be worthwhile to analyze how individual model components (e.g., physical parameterization, meteorological data sets, model horizontal/vertical resolution) impact the simulation of stratospheric CH 4 and XCH 4 .

54 ENVIRONMENTAL SCIENCES↗

FLUXNET-CH 4 : a global, multi-ecosystem dataset and analysis of methane seasonality from freshwater wetlands

Methane (CH 4 ) emissions from natural landscapes constitute roughly half of global CH4 contributions to the atmosphere, yet large uncertainties remain in the absolute magnitude and the seasonality of emission quantities and drivers. Eddy covariance (EC) measurements of CH4 flux are ideal for constraining ecosystem-scale CH 4 emissions due to quasi-continuous and high-temporal-resolution CH 4 flux measurements, coincident carbon dioxide, water, and energy flux measurements, lack of ecosystem disturbance, and increased availability of datasets over the last decade. Here, we (1) describe the newly published dataset, FLUXNET-CH 4 Version 1.0, the first open-source global dataset of CH 4 EC measurements (available at https://fluxnet.org/data/fluxnet-ch4-community-product/, last access: 7 April 2021). FLUXNET-CH 4 includes half-hourly and daily gap-filled and non-gap-filled aggregated CH4 fluxes and meteorological data from 79 sites globally: 42 freshwater wetlands, 6 brackish and saline wetlands, 7 formerly drained ecosystems, 7 rice paddy sites, 2 lakes, and 15 uplands. Then, we (2) evaluate FLUXNET-CH 4 representativeness for freshwater wetland coverage globally because the majority of sites in FLUXNET-CH4 Version 1.0 are freshwater wetlands which are a substantial source of total atmospheric CH 4 emissions; and (3) we provide the first global estimates of the seasonal variability and seasonality predictors of freshwater wetland CH4 fluxes. Our representativeness analysis suggests that the freshwater wetland sites in the dataset cover global wetland bioclimatic attributes (encompassing energy, moisture, and vegetation-related parameters) in arctic, boreal, and temperate regions but only sparsely cover humid tropical regions. Seasonality metrics of wetland CH 4 emissions vary considerably across latitudinal bands. In freshwater wetlands (except those between 20°S to 20°N) the spring onset of elevated CH 4 emissions starts 3 d earlier, and the CH 4 emission season lasts 4 d longer, for each degree Celsius increase in mean annual air temperature. On average, the spring onset of increasing CH 4 emissions lags behind soil warming by 1 month, with very few sites experiencing increased CH 4 emissions prior to the onset of soil warming. In contrast, roughly half of these sites experience the spring onset of rising CH 4 emissions prior to the spring increase in gross primary productivity (GPP). The timing of peak summer CH 4 emissions does not correlate with the timing for either peak summer temperature or peak GPP. Our results provide seasonality parameters for CH 4 modeling and highlight seasonality metrics that cannot be predicted by temperature or GPP (i.e., seasonality of CH 4 peak). FLUXNET-CH 4 is a powerful new resource for diagnosing and understanding the role of terrestrial ecosystems and climate drivers in the global CH 4 cycle, and future additions of sites in tropical ecosystems and site years of data collection will provide added value to this database. All seasonality parameters are available at https://doi.org/10.5281/zenodo.4672601 (Delwiche et al., 2021). Additionally, raw FLUXNET-CH 4 data used to extract seasonality parameters can be downloaded from https://fluxnet.org/data/fluxnet-ch4-community-product/ (last access: 7 April 2021), and a complete list of the 79 individual site data DOIs is provided in Table 2 of this paper.

54 ENVIRONMENTAL SCIENCES↗

Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra site

Numerical simulations are essential tools for understanding the complex hydrologic response of Arctic regions to a warming climate. However, strong coupling among thermal and hydrological processes on the surface and in the subsurface and the significant role that subtle variations in surface topography have in regulating flow direction and surface storage lead to significant uncertainties. Careful model evaluation against field observations is thus important to build confidence.We evaluate the integrated surface/subsurface permafrost thermal hydrology models in the Advanced Terrestrial Simulator (ATS) against field observations from polygonal tundra at the Barrow Environmental Observatory. ATS couples a multiphase, 3D representation of subsurface thermal hydrology with representations of overland nonisothermal flows, snow processes, and surface energy balance. We simulated thermal hydrology of a 3D ice-wedge polygon with geometry that is abstracted but broadly consistent with the surface microtopography at our study site. The simulations were forced by meteorological data and observed water table elevations in ice-wedge polygon troughs. With limited calibration of parameters appearing in the soil evaporation model, the 3-year simulations agreed reasonably well with snow depth, summer water table elevations in the polygon center, and high-frequency soil temperature measurements at several depths in the trough, rim, and center of the polygon. Upscaled evaporation is in good agreement with flux tower observations. The simulations were found to be sensitive to parameters in the bare soil evaporation model, snowpack, and the lateral saturated hydraulic conductivity.Timing of fall freeze-up was found to be sensitive to initial snow density, illustrating the importance of including snow aging effects. The study provides new support for an emerging class of integrated surface/subsurface permafrost simulators.

54 ENVIRONMENTAL SCIENCES↗

Thermal modeling of three lakes within the continuous permafrost zone in Alaska using the LAKE 2.0 model

Abstract. Lakes in the Arctic are important reservoirs of heat with much lower albedo in summer and greater absorption of solar radiation than surrounding tundra vegetation. In the winter, lakes that do not freeze to their bed have a mean annual bed temperature >0 ∘C in an otherwise frozen landscape. Under climate warming scenarios, we expect Arctic lakes to accelerate thawing of underlying permafrost due to warming water temperatures in the summer and winter. Previous studies of Arctic lakes have focused on ice cover and thickness, the ice decay process, catchment hydrology, lake water balance, and eddy covariance measurements, but little work has been done in the Arctic to model lake heat balance. We applied the LAKE 2.0 model to simulate water temperatures in three Arctic lakes in northern Alaska over several years and tested the sensitivity of the model to several perturbations of input meteorological variables (precipitation, shortwave radiation, and air temperature) and several model parameters (water vertical resolution, sediment vertical resolution, depth of soil column, and temporal resolution). The LAKE 2.0 model is a one-dimensional model that explicitly solves vertical profiles of water state variables on a grid. We used a combination of meteorological data from local and remote weather stations, as well as data derived from remote sensing, to drive the model. We validated modeled water temperatures with data of observed lake water temperatures at several depths over several years for each lake. Our validation of the LAKE 2.0 model is a necessary step toward modeling changes in Arctic lake ice regimes, lake heat balance, and thermal interactions with permafrost. The sensitivity analysis shows us that lake water temperature is not highly sensitive to small changes in air temperature or precipitation, while changes in shortwave radiation and large changes in precipitation produced larger effects. Snow depth and lake ice strongly affect water temperatures during the frozen season, which dominates the annual thermal regime of Arctic lakes. These findings suggest that reductions in lake ice thickness and duration could lead to more heat storage by lakes and enhanced permafrost degradation.

Clark, Jason A. (ORCID:0000000248105261)↗