Search NASA⌕ Search

SEARCH · Search NASA

Results for “meteorological modeling”

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

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

At least 91 records · Page 5

Quantifying mean, variability, and uncertainty in indoor radon exposure in Pennsylvania using random forest and quantile regression forest models

Radon is a naturally occurring radioactive gas that poses a serious health risk as the primary cause of lung cancer in non-smokers. Despite the well-known adverse association with health outcomes, current radon exposure assessments are limited to county-level or average-level estimates, which fail to capture regional variability. This study uses Machine Learning models, including Random Forest (RF) and Quantile Regression Forest (QRF), to estimate the indoor radon concentrations at the ZCTA (Zip code tabulation area)-level and characterize uncertainties in model estimates. Incorporating geological, meteorological, and building-specific data, the models aim to improve radon risk assessment by capturing mean exposure, variability, and extreme concentration levels. Processed radon test data (n = 718,111) were analyzed using average, variability, and quantile prediction methods. Models that estimate the average radon exposure at the ZCTA-level can yield promising model-fit results, but they do not capture the underlying variability of indoor radon exposure within a ZCTA. We utilize volatility analyses to identify characteristics indicative of high variability of indoor radon exposure. We also show that a QRF model can be used to estimate upper quantiles of residential radon exposure, thereby uncovering localized areas of elevated exposure that were not apparent in mean estimates. The results highlighted the need for a deep characterization of exposure risk and show that regions with moderate average exposure levels could still harbor extreme outliers with implications for evaluating health risks. Utilizing multiple radon exposure models allows for a deeper characterization of radon risk within a geographic area and can better identify high-risk areas. The results from this study provide a foundation for developing mitigation strategies and examining associations between radon exposure and health outcomes at fine scales. Future research should extend the geographic scope and incorporate additional environmental risk factors to establish a comprehensive framework for risk assessment.

Lee, Heechan [ORNL]↗

High-Resolution Meteorology with Climate Change Impacts from Global Climate Model Data Using Generative Machine Learning

As renewable energy generation increases, the impacts of weather and climate on energy generation and demand become critical to the reliability of the energy system. However, these impacts are often overlooked. Global climate models (GCMs) can be used to understand possible changes to our climate, but their coarse resolution makes them difficult to use in energy system modelling. Here we present open-source generative machine learning methods that produce meteorological data at a nominal spatial resolution of 4 km at an hourly frequency based on inputs from 100 km daily-average GCM data. These methods run 40 times faster than traditional downscaling methods and produce data that have high-resolution spatial and temporal attributes similar to historical datasets. We demonstrate that these methods can be used to downscale projected changes in wind, solar and temperature variables across multiple GCMs including projections for more frequent low-wind and high-temperature events in the Eastern United States.

climate change↗

Xanthos-Lake Dataset

The Xanthos-Lake v1.0 dataset provides the input data, trained machine-learning models, and simulation outputs needed to characterize lake water balance, snow and ice conditions, and mixing-layer temperature within the Xanthos global hydrological modeling framework. The dataset supports lake representation across a wide range of lake sizes and hydroclimatic conditions by combining xLSIM, a basin-specific machine-learning emulator of lake snow, ice, ice-cover fraction, and mixing-layer temperature, with the Xanthos-Lake water-balance model. The archive contains NetCDF datasets used to train and evaluate xLSIM, trained model weights, processed meteorological and lake-property inputs, and basin- and lake-category-specific simulation outputs. These materials are organized into four primary data groups, described below. Snowice_model_inputs: Contains the NetCDF input data used to train xLSIM. The xLSIM machine-learning framework uses three lake-based datasets. The meteorological forcing dataset provides monthly relative humidity, specific humidity, surface wind speed, maximum and minimum air temperature, downward longwave and shortwave radiation, snowfall, surface air pressure, and total precipitation. Lake surface area is included as an additional static predictor. The target-state dataset provides lake ice thickness, snow depth, snow cover, and lake mixing-layer temperature, while a companion lake-surface dataset provides the lake ice-cover fraction. Before training, ice thickness and snow depth are converted from meters to centimeters, mixing-layer temperature is converted from kelvin to degrees Celsius and constrained to nonnegative values, and ice-cover fraction is converted from a fraction to a percentage. The predictor variables are normalized using statistics calculated across the selected lakes and time steps. Snowice_model_outputs: Contains the NetCDF outputs generated by xLSIM. For each basin, xLSIM produces a file containing observed and predicted lake-state variables for the training, validation, and testing periods. The modeled variables include lake ice thickness, snow depth, snow cover, mixing-layer temperature, and lake ice-cover fraction. For basins without a sufficiently persistent snow-and-ice signal, the emulator predicts only mixing-layer temperature. The outputs also include training and validation loss histories, the selected model configuration, identifiers of the lakes used in training, and SHAP-based feature-importance information at the global, lake, and seasonal-regime levels. The trained machine-learning model weights are provided separately within the dataset archive. Together, these files support model evaluation and subsequent coupling with the Xanthos-Lake water-balance framework. XanthosLAKES: Contains the NetCDF input data used by the Xanthos-Lake framework. Monthly meteorological inputs include relative and specific humidity, downward shortwave and longwave radiation, mean, maximum, and minimum air temperature, wind speed, precipitation, snowfall, and surface air pressure. Static lake-property datasets provide lake identifiers, geographic locations, surface area, volume, mean depth, elevation, drainage area, fetch, outlet-routing information, and associated Xanthos grid-cell attributes. Separate bathymetric datasets provide the coefficients of the area–depth and volume–depth relationships for each aggregated lake unit. GLEV-based records provide observed lake surface area and evaporation data used to initialize lake states, define reference conditions, and calibrate and evaluate the model. Xanthos-Lake Outputs: Contains the basin- and lake-category-specific NetCDF outputs generated by Xanthos-Lake. Monthly variables include lake surface area, storage volume, outlet discharge, evaporation rate, evaporation volume, lake–groundwater exchange, lake inflow, ice thickness, snow depth, snow-cover fraction, ice-cover fraction, and mixing-layer temperature. The files also contain lake-specific calibration and validation statistics, including normalized root-mean-square error, mean absolute error, Nash–Sutcliffe efficiency, Kling–Gupta efficiency, and percent bias. Stored calibrated and derived parameters include the weir discharge coefficient, fractional freeboard, groundwater exchange coefficient, reference water level, corresponding reference surface area and storage volume, weir-width adjustment factor, and the fraction of routed inflow entering the lake. Basin identifiers, lake category, simulation period, calibration and validation periods, and parameter-schema information are retained as NetCDF metadata.

Abeshu, Guta [Pacific Northwest National Laborator↗

Monsoonal MCS Initiation, Rainfall, and Diurnal Gravity Waves over the Bay of Bengal: Observation and a Linear Model

Abstract Previous observational studies have indicated that mesoscale convective systems (MCSs) contribute the majority of precipitation over the Bay of Bengal (BoB) during the summer monsoon season, yet their initiation and propagation remain incompletely understood. To fill this knowledge gap, we conducted a comprehensive study using a combination of 20-yr satellite observations, MCS tracking, reanalysis data, and a theoretical linear model. Satellite observations reveal clear diurnal propagation signals of MCS initiation frequency and rainfall from the west coast of the BoB toward the central BoB, with the MCS rainfall propagating slightly slower than the MCS initiation frequency. Global reanalysis data indicate a strong association between the offshore-propagating MCS initiation frequency/rainfall and diurnal low-level wind perturbations, implying the potential role of gravity waves. To verify the hypothesis, we developed a 2D linear model that can be driven by realistic meteorological fields from reanalysis. The linear model realistically reproduces the characteristics of offshore-propagating diurnal wind perturbations. The wind perturbations, as well as the offshore propagation signals of MCS initiation frequency and rainfall, are associated with diurnal gravity waves emitted from the coastal regions, which in turn are caused by the diurnal land–sea thermal contrast. The ambient wind speed and vertical wind shear play crucial roles in modulating the timing, propagation, and amplitude of diurnal gravity waves. Using the linear model and satellite observations, we further show that the stronger monsoonal flows lead to faster offshore propagation of diurnal gravity waves, which subsequently control the offshore propagation signals of MCS initiation and rainfall. Significance Statement Rainfall over the Bay of Bengal (BoB) is primarily contributed by large and organized rainfall systems in the summer monsoon season. During this season, these systems are commonly observed over the east coast of India around midnight, the western BoB in the morning, and the central BoB in the afternoon. This eastward rainfall propagation is confirmed by observations, reanalysis data, and a theoretical model to have a strong association with atmospheric diurnal gravity waves. These waves are caused by land–sea thermal contrast and can trigger rainfall systems over the offshore regions. We also found that the diurnal gravity waves, as well as their triggered rainfall systems, can be greatly modulated by the large-scale monsoonal flows. All the above findings improve our understanding of diurnal rainfall cycle over the BoB and may contribute to the future improvement of rainfall forecast over the region.

Meteorology & Atmospheric Sciences↗

The vertical-velocity skewness in the inertial sublayer of turbulent wall flows

Empirical evidence is provided that within the inertial sublayer (i.e. logarithmic region) of adiabatic turbulent flows over smooth walls, the skewness of the vertical-velocity component 𝑆 𝑤 displays universal behaviour, being a positive constant and constrained within the range 𝑆 𝑤 ≈ 0.1–0.16, regardless of flow configuration and Reynolds number. A theoretical model is then proposed to explain this behaviour, including the observed range of variations of 𝑆 𝑤 . The proposed model clarifies why 𝑆 𝑤 cannot be predicted from down-gradient closure approximations routinely employed in large-scale meteorological and climate models. The proposed model also offers an alternative and implementable approach for such large-scale models.

turbulent boundary layers↗

Apparatus and methods for location and sizing of trace gas sources

A system for detecting gas leaks and determining their location and size. A data gathering portion of the system utilizes a chosen geometrical configuration to collect path-integrated spectroscopic data over multiple paths around an area. A processing portion of the system applies a transport model together with meteorological data of the area to generate an influence function of possible leak locations on gas detector measurement paths, and applies an inversion model to the influence function, prior data, and the spectroscopic data to generate gas source size and location.

03 NATURAL GAS↗

Processing Meteorological Data for the CAP-88 PC Model at Los Alamos National Laboratory

The Environmental Protection and Compliance-Compliance Programs (EPC-CP) group at Los Alamos National Laboratory (LANL) uses the Clean Air Act Assessment Package 1988 (CAP-88, Littleton 2020) PC model (Version 4.1) to estimate radiological doses for a set of areal sectors surrounding a release location, in order to satisfy the Environmental Protection Agency (EPA) National Emission Standards for Hazardous Air Pollutants (NESHAP) dose calculation requirement in 40 CFR 61 Subpart H. Among several types of data that must be prepared for CAP-88 input is a text file of meteorological data (“WIND” file), consisting of the joint frequency of wind direction, wind speed, and atmospheric stability categories. EPC-CP produces customized WIND files by running a CAP-88 utility program on a user generated text file of wind data in a different format, known as a STability ARray (STAR) file (Turner, 1964). At LANL, EPC-CP meteorologists prepare customized STAR files with data over desired time periods at selected meteorological towers. A custom program written in Precision Visuals -Workstation Analysis and Visualization Environment (PV-WAVE), a commercial Fortran-like language, is used to read LANL meteorological data and write a STAR file; the executable filename is “Star.out”. However, the outdated PV-WAVE utility program is being phased out by EPC-CP, due to the inefficient process to run it and an inability to modify the code. To preserve the ability to create customized meteorological data for CAP-88 in a way that will be easy to use and maintain, a new replacement utility program, written in the Python programming language, has been developed. The new, improved program reads a data file from any LANL meteorological tower, and at each desired observation time, determines the wind direction, wind speed, and stability categories defined in the CAP-88 documentation. The frequencies of all combinations of the three sets of categories are calculated and written to a file in the STAR format, which can later be converted to a WIND file for input into CAP-88.

54 ENVIRONMENTAL SCIENCES↗

Characterizing How Meteorological Forcing Selection and Parameter Uncertainty Influence Community Land Model Version 5 Hydrological Applications in the United States

Despite the increasing use of large-scale Land Surface Models (LSMs) in predicting hydrological responses in extreme conditions, there's a critical gap in understanding the uncertainties in these predictions. This study addresses this gap through a detailed diagnostic evaluation of the uncertainties arising from meteorological forcing selection and model parametrization in hydrological simulations of the Community Land Model version 5 (CLM5). CLM5 is configured at a spatial scale of about 12-km to simulate runoff processes for 464 headwater watersheds, selected from the Catchment Attributes for Large-Sample Studies (CAMELS) dataset to be representative of physiographic and climatic gradients across the conterminous United States. For each watershed, CLM5 is driven by five commonly used gridded forcing datasets in combination with a large ensemble (> 1200) of key CLM5 hydrologic parameters. Our results suggest that uncertainty in CLM5 runoff simulations resulting from both forcing and parametric sources is markedly higher in arid regions, e.g., Great Plains and Midwest regions. Uncertainty in low flow is dominated by parametric uncertainty, while the selection of meteorological forcing contributes more dominantly to high flow and seasonal flows during fall and spring. Our analysis also demonstrates that the selection of forcing datasets and the metrics used to calibrate CLM5 significantly impact the model’s predictive accuracy in extreme event severity for both floods and droughts. Overall, the results from this study highlight the need to understand and account for forcing and parametric uncertainties in CLM5 simulations, particularly for hazard and risk assessments addressing hydrologic extremes.

54 ENVIRONMENTAL SCIENCES↗

Estimation of Surface Sensible Heat Flux due to Precipitation over CONUS and Its Impact on Urban Extreme Precipitation Modeling

The surface sensible heat flux induced by precipitation (Q P ) is a consequence of the temperature difference between the surface and the rain droplets. Despite its seemingly negligible nature, Q P is frequently omitted from both meteorological and climatological models. Nevertheless, it is important to acknowledge the numerous occasions in which the instantaneous values of Q P can be significant, particularly during extreme precipitation events. This study undertakes a comprehensive assessment of Q P across the contiguous United States (CONUS) utilizing high-resolution reanalysis, observational data, and numerical modeling to examine the influence of Q P on precipitation and the surface energy budget. The findings indicate that the spatial distribution of Q P climatology is analogous to that of precipitation, with magnitudes ranging from 2 to 3 W m -2 predominantly over the Midwest and Southeast regions. A seasonal analysis of QP reveals that the highest values occurring during the June–August (JJA) period, averaging 3.18 W m -2 . Peak Q P values of approximately 4 W m -2 are observed during JJA over the Great Plains region. We hypothesize that the Q P during an extreme precipitation event would be nonnegligible and have a significant impact on the local weather. To test this conjecture, we perform high-resolution simulations with and without Q P during an extreme precipitation event over the Chicago Metropolitan Area (CMA). The results show that the Q P may be a dominant factor compared to other components of surface heat flux during the zenith of precipitation hours. Also, Q P has the potential to not only diminish precipitation but also alter and reconfigure the remaining surface energy budget components.

54 ENVIRONMENTAL SCIENCES↗

Meteorological characteristics of extreme ozone pollution events in China and their future predictions

Ozone (O 3 ) has become one of the most concerning air pollutants in China in recent decades. In this study, based on surface observations, reanalysis data, global atmospheric chemistry model simulations, and multi-model future predictions, meteorological characteristics conducive to extreme O 3 pollution in various regions of China are investigated, and their historical changes and future trends are analyzed. During the most severe O 3 polluted months, the chemical production of O 3 is enhanced under the hot and dry conditions over the North China Plain (NCP) in June 2018 and the Yangtze River Delta (YRD) in July 2017, while regional transport is the main reason for the severe O3 pollution over the Sichuan Basin (SCB) in July 2015 and the Pearl River Delta (PRD) in September 2019. Over the last 4 decades, the frequencies of high-temperature and low-relative-humidity conditions increased in 2000–2019 relative to 1980–1999, indicating that O 3 pollution in both the NCP and YRD has become more frequent under historical climate change. In the SCB and PRD, the occurrence of atmospheric circulation patterns similar to those during the most polluted months increased, together with the more frequent hot and dry conditions, contributing to the increases in severe O 3 pollution in the SCB and PRD during 1980–2019. In the future (by 2100), the frequencies of months with anomalous high temperature show stronger increasing trends in the high-forcing scenario (Shared Socioeconomic Pathway (SSP5-8.5)) compared to the sustainable scenario (SSP1-2.6) in China. It suggests that high anthropogenic forcing will not only lead to slow economic growth and climate warming but also likely result in environmental pollution issues.

54 ENVIRONMENTAL SCIENCES↗

Gap Flows Along the Columbia River Observed During the WFIP2 Field Campaign

Data collected in the Columbia River Gorge and Basin Area during the second Wind Forecasting Improvement Project (WFIP2) was used to study gap flow events in the region. In particular, the relationship between 200 m zonal wind speed and the zonal mesoscale pressure gradient at Troutdale (West of the gorge), Wasco (east of the gorge), and Boardman (farther east of the gorge) sites was examined. An objective criterion was used to identify 169 (Troutdale), 161 (Wasco), and 113 (Boardman) gap flow days that had weak synoptic forcing. The 200 m zonal wind speed exhibited a linear relationship with the zonal pressure gradient at Wasco during both easterly and westerly gap flow days at daily timescales. Such a relationship was only observed during easterly gap flow days at Troutdale and was not observed at Boardman. The zonal wind speed exhibited a diurnal cycle at Wasco and Boardman, with greater changes in the nighttime winds than the daytime winds due to changes in pressure gradient. The root mean square differences (RMSDs) between the observed pressure gradient at Wasco and that simulated by the High Resolution Rapid Refresh (HRRR) model in 3 km (HRRR-3km) and 750 m (HRRR-750m) horizontal resolution were 3.14 hPa and 0.57 hPa, respectively. However, the RMSDs for the 200 m zonal wind speed were 1.40 m s -1 for HRRR-3km and 1.20 m s -1 for HRRR- 750m. Improvement in simulating the zonal pressure gradient accompanied by negligible changes in simulating zonal wind speed due to higher horizontal resolution points to a complex interplay between different components of the model physics when simulating the meteorology. The HRRR-750m resolution model accurately mimicked the observed covariability (phase diagram) between the zonal pressure gradient and the 200 m zonal wind speed during the nighttime, but not during the daytime, suggesting model errors might be related to the representation of boundary layer processes.

54 ENVIRONMENTAL SCIENCES↗

Sensitivities of subgrid-scale physics schemes, meteorological forcing, and topographic radiation in atmosphere-through-bedrock integrated process models: a case study in the Upper Colorado River basin

Abstract. Mountain hydrology is controlled by interacting processes extending from the atmosphere through the bedrock. Integrated process models (IPMs), one of the main tools needed to interpret observations and refine conceptual models of the mountainous water cycle, require meteorological forcing that simulates the atmospheric process to predict hydroclimate then subsequently impacts surface–subsurface hydrology. Complex terrain and extreme spatial heterogeneity in mountainous environments drive uncertainty in several key considerations in IPM configurations and require further quantification and sensitivity analyses. Here, we present an IPM using the Weather Research and Forecasting (WRF) model which forces an integrated hydrologic model, ParFlow-CLM, implemented over a domain centered over the East River watershed (ERW), located in the Upper Colorado River basin (UCRB). The ERW is a heavily instrumented 300 km2 region in the headwaters of the UCRB near Crested Butte, CO, with a growing atmosphere-through-bedrock observation network. Through a series of experiments in the water year 2019 (WY19), we use four meteorological forcings derived from commonly used reanalysis datasets, three subgrid-scale physics scheme configurations in WRF, and two terrain shading options within WRF to test the relative importance of these experimental design choices for key hydrometeorological metrics including precipitation and snowpack, as well as evapotranspiration, groundwater storage, and discharge simulated by the ParFlow-CLM. Our hypothesis is that uncertainty from synoptic-scale forcings produces a much larger spread in surface–subsurface hydrologic fields than subgrid-scale physics scheme choice. Results reveal that the WRF subgrid-scale physics configuration leads to larger spatiotemporal variance in simulated hydrometeorological conditions, whereas variance across meteorological forcing with common subgrid-scale physics configurations is more spatiotemporally constrained. Despite reasonably simulating precipitation, a delay in simulated discharge peak is due to a systematic cold bias across WRF simulations, suggesting the need for bias correction. Discharge shows greater variance in response to the WRF simulations across subgrid-scale physics schemes (26 %) rather than meteorological forcing (6 %). The topographic radiation option has minor effects on the watershed-average hydrometeorological processes but adds profound spatial heterogeneity to local energy budgets (±30 W m−2 in shortwave radiation and 1 K air temperature differences in late summer). This is the first presentation of sensitivity analyses that provide support to help guide the scientific community to develop observational constraints on atmosphere-through-bedrock processes and their interactions.

54 ENVIRONMENTAL SCIENCES↗

Data-model files associated with the manuscript "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS)

This data package contains the model inputs and outputs used in "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS). The data.zip file contains the data used to drive the model simulations. The model.zip file contains the XML input file for ATS. The notebook.zip file contains the Jupyter notebooks for pre- and post- processing model results. The figures.zip file contains the raw figures associated with the manuscript.Meteorological forcing plays a critical role in accurately simulating the watershed hydrological cycle. With the advancement of high-performance computing and the development of integrated watershed models, simulating the watershed hydrological cycle at high temporal (hourly to daily) and spatial resolution (10s of meters) has become efficient and computationally affordable. These hyperresolution watershed models require high resolution of meteorological forcing as model input to ensure the fidelity and accuracy of simulated responses. In this study, we utilized the Advanced Terrestrial Simulator (ATS), an integrated watershed model, to simulate surface and subsurface flow and land surface processes using unstructured meshes at the Coal Creek Watershed near Crested Butte (Colorado). We compared simulated watershed hydrologic responses including streamflow, and distributed variables such as evapotranspiration, snow water equivalent (SWE), and groundwater table driven by three publicly available, gridded meteorological forcings (GMFs) -- Daily Surface Weather and Climatological Summaries (Daymet), Parameter-elevation Regressions on Independent Slopes Model (PRISM), and North American Land Data Assimilation System (NLDAS). By comparing various spatial resolutions (ranging from 400 m to 4 km) of PRISM, the simulated streamflow only becomes marginally worse when spatial resolution of meteorological forcing is coarsened to 4 km (or 30% of the watershed area). However, the 4 km resolution has much worse performance than finer resolution in spatially distributed variables such as SWE. Using temporally disaggregated PRISM, we compared models forced by different temporal resolutions (hourly to daily), sub-daily resolution preserves the dynamic watershed responses (e.g., diurnal fluctuation of streamflow) that are absent in results forced by daily resolution. Conversely, the simulated streamflow shows better performance using daily resolution compared to that using sub-daily resolution. Our findings suggest that the choice of GMF and its spatiotemporal resolution depends on the quantity of interest and its spatial and temporal scale, which may have important implications on model calibration and watershed management decisions.

54 ENVIRONMENTAL SCIENCES↗

How can biosphere models simulate enough vegetation biomass in the mountains of the western United States? Implications of meteorological forcing

Most carbon stocks and fluxes in the western United States are found in mountainous terrain, where observations and modeling are difficult. Terrestrial biosphere models generally underestimate above-ground biomass (AGB) over this region. In this report, we identify methods to reduce this underestimation by focusing upon 1) biases in meteorological datasets, 2) model representation of water stress, and 3) spatial resolution. We adopted the widely-used Community Land Model version 4.5 (CLM 4.5) with six different meteorological datasets and found a 6-fold variation in simulated AGB across Utah/Colorado. Simulations underestimated AGB because of warm and dry biases within the meteorological datasets that reduced water availability and restricted plant growth. To eliminate the AGB underestimation we adopted a meteorological dataset designed for complex terrain (gridMET), combined with a representation of plant hydraulic stress (CLM 5.0). Conversely, changes in spatial resolution (meteorological variables and land surface description) had negligible impact on simulated AGB.

54 ENVIRONMENTAL SCIENCES↗

A new parameterization of ice heterogeneous nucleation coupled to aerosol chemistry in WRF-Chem model version 3.5.1: evaluation through ISDAC measurements

In the Arctic, during polar night and early spring, ice clouds are separated into two leading types of ice clouds (TICs): (1) TIC1 clouds characterized by a large concentration of very small crystals and TIC2 clouds characterized by a low concentration of large ice crystals. Using a suitable parameterization of heterogeneous ice nucleation is essential for properly representing ice clouds in meteorological and climate models and subsequently understanding their interactions with aerosols and radiation. Here, we describe a new parameterization for ice crystal formation by heterogeneous nucleation in water-subsaturated conditions coupled to aerosol chemistry in the Weather Research and Forecasting model coupled with chemistry (WRF-Chem). The parameterization is implemented in the Milbrandt and Yau (2005a, b) two-moment cloud microphysics scheme, and we assess how the WRF-Chem model responds to the run-time interaction between chemistry and the new parameterization. Well-documented reference cases provided us with in situ data from the spring 2008 Indirect and Semi-Direct Aerosol Campaign (ISDAC) over Alaska. Our analysis reveals that the new parameterization clearly improves the representation of the ice water content (IWC) in polluted or unpolluted air masses and shows the poor performance of the reference parameterization in representing ice clouds with low IWC. The new parameterization is able to represent TIC1 and TIC2 microphysical characteristics at the top of the clouds, where heterogenous ice nucleation is most likely occurring, even with the known bias of simulated aerosols by WRF-Chem over the Arctic.

54 ENVIRONMENTAL SCIENCES↗

Machine Learning Analysis of Impact of Western US Fires on Central US Hailstorms

Fires, including wildfires, harm air quality and essential public services like transportation, communication, and utilities. These fires can also influence atmospheric conditions, including temperature and aerosols, potentially affecting severe convective storms. Here, we investigate the remote impacts of fires in the western United States (WUS) on the occurrence of large hail (size: $\geqslant$ 2.54 cm) in the central US (CUS) over the 20-year period of 2001–20 using the machine learning (ML), Random Forest (RF), and Extreme Gradient Boosting (XGB) methods. The developed RF and XGB models demonstrate high accuracy (> 90%) and F1 scores of up to 0.78 in predicting large hail occurrences when WUS fires and CUS hailstorms coincide, particularly in four states (Wyoming, South Dakota, Nebraska, and Kansas). The key contributing variables identified from both ML models include the meteorological variables in the fire region (temperature and moisture), the westerly wind over the plume transport path, and the fire features (i.e., the maximum fire power and burned area). Importantly, the results confirm a linkage between WUS fires and severe weather in the CUS, corroborating the findings of our previous modeling study conducted on case simulations with a detailed physics model.

54 ENVIRONMENTAL SCIENCES↗

Model estimates of dissolved organic carbon, runoff, snowmelt, and snow water equivalent from 1981-2010 across the western Arctic

This dataset contains estimates of dissolved organic carbon (DOC) yield (mg C m-2), runoff (mm day-1), snow water equivalent (SWE, mm day-1), and snowmelt (mm day-1) from a simulation of the Permafrost Water Balance Model (PWBM v4). The simulation effectively quantifies snowpack accumulation and melt, runoff, and DOC leaching, and broadly capture the seasonal cycle in DOC concentration and mass loadings. Daily air temperature, precipitation, and wind speed data from the Modern-Era Retrospective Analysis for Research and Applications (MERRA) were used for meteorological forcings. The model outputs are organized by grid and day of the month in gzipped ASCII text files contained within each archived ZIP file. No special software is required to work with these data. The basin and domain information files are also ASCII text. The gridded estimates are most useful for analyses of the spatial and temporal variations in terrestrial hydrology and leachate DOC concentrations and loadings across parts of the western Arctic. A manuscript describing the data and associated analysis has been accepted for publication in Journal of Geophysical Research - Biogeosciences (Rawlins et al., 2021).

54 ENVIRONMENTAL SCIENCES↗

Data Quality Objectives Supporting Radiological Air Emissions Monitoring for the PNNL-Richland Campus: North Campus Construction (Rev. 3)

Starting in fiscal year 2023, the north area of the PNNL-Richland Campus will undergo development. Initially, infrastructure (water, electrical) improvements will be installed. Later development includes the construction of several new office buildings. The area to be developed currently contains no buildings but does have two existing solar-powered ambient air sampling stations. Under the requirements of Washington State Department of Health Radioactive Air Emissions License -005, the PNNL-Richland Campus must operate and maintain a radiological air monitoring program. This revision documents and evaluates how the new North Campus construction impacts ambient air surveillance stations within the development region. This revision also considers current stack configurations and uses an updated environmental dispersion model and updated meteorological data. The DQO team concluded that one sampling station in the construction area will be relocated and the other remain in place. The evaluations conducted for this DQO also identified a third sampling station outside the development area that is recommended for relocation as a result of updated dispersion modeling. Recommended sites for the relocations are presented. Considerations for sampling in and around a construction zone are also addressed. Additionally, programmatic improvements to the ambient air sampling program were identified in the DQO process.

54 ENVIRONMENTAL SCIENCES↗