Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied meteorology”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Detection of atmospheric pressure loading using very long baseline interferometry measurements

Loading of the Earth by the temporal redistribution of global atmospheric mass is likely to displace the positions of geodetic monuments by tens of millimeters both vertically and horizontally. Estimates of these displacements are determined by convolving National Meteorological Center (NMC) global values of atmospheric surface pressure with Farrell's elastic Green's functions. An analysis of the distances between radio telescopes determined by very long baseline interferometry (VLBI) between 1984 and 1992 reveals that in many of the cases studied there is a significant contribution to baseline length change due to atmospheric pressure loading. Our analysis covers intersite distances of between 1000 and 10,000 km and is restricted to those baselines measured more than 100 times. Accounting for the load effects (after first removing a best fit slope) reduces the weighted root-mean-square (WRMS) scatter of the baseline length residuals on 11 of the 22 baselines investigated. The slight degradation observed in the WRMS scatter on the remaining baselines is largely consistent with the expected statistical fluctuations when a small correction is applied to a data set having a much larger random noise. The results from all baselines are consistent with approximately 60% of the computed pressure contribution being present in the VLBI length determinations. Site dependent coefficients determined by fitting local pressure to the theoretical radial displacement are found to reproduce the deformation caused by the regional pressure to within 25% for most inland sites. The coefficients are less reliable at near coastal and island stations.

Vandam, T. M.↗
DOE OSTI2023

Multi-scale planning model for robust urban drought response

Increasingly severe droughts are straining municipal water resources and jeopardizing urban water security, but uncertainty in their duration, frequency, and intensity challenges drought planning and response. We develop the Drought Resilient Interscale Portfolio Planning model (DRIPP) to generate optimal planning responses to urban drought. DRIPP is a generalizable multi-scale framework for optimizing dynamic planning strategies of long-term infrastructure deployment and short-term drought response. It integrates climate and hydrological variability with high-fidelity representations of urban water distribution, available technology options, and demand reduction measures to yield robust and cost-effective water supply portfolios that are location-specific. We apply DRIPP in Santa Barbara, California to assess how least cost water supply portfolios vary under different drought scenarios and identify portfolios that are robust across drought scenarios. In Santa Barbara, we find that drought intensity, not duration or frequency, drives cost increases, reliability risk, and regret of overbuilding infrastructure. Under uncertain drought conditions, a diversified technology portfolio that includes both rapidly deployable, decentralized technologies alongside larger centralized technologies minimizes water supply cost while maintaining high robustness to climate uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Estimated impacts of forest restoration scenarios on smoke exposures among outdoor agricultural workers in California

Abstract As wildfires continue to worsen across western United States, forest managers are increasingly employing prescribed burns as a way to reduce excess fuels and future wildfire risk. While the ecological benefits of these fuel treatments are clear, little is known about the smoke exposure tradeoffs of using prescribed burns to mitigate wildfires, particularly among at-risk populations. Outdoor agricultural workers are a population at increased risk of smoke exposure because of their time spent outside and the physical demands of their work. Here, we assess the smoke exposure impacts among outdoor agricultural workers resulting from the implementation of six forest management scenarios proposed for a landscape in the Central Sierra, California. We leverage emissions estimates from LANDIS-II to model daily PM 2.5 concentrations with the Hybrid Single Particle Lagrangian Integrated Trajectory Model (HYSPLIT) and link those to agricultural employment data from the Bureau of Labor Statistics. We find a u-shaped result, in that moderate amounts of prescribed burning result in the greatest reduction in total smoke exposure among outdoor agricultural workers, particularly during months of peak agricultural activity due to wildfire-specific smoke reductions. The reduction in total smoke exposure, relative to scenarios with the least amount of management, decreases as more prescribed burning is applied to the landscape due to the contributions of the fuel treatments themselves to overall smoke burden. The results of this analysis may contribute to preparedness efforts aimed at reducing smoke exposures among outdoor agricultural workers, while also informing forest management planning for this specific landscape.

Environmental Sciences & Ecology↗
DOE OSTI2021

Estimating and Forecasting Time-Varying Groundwater Recharge in Fractured Rock: A State-Space Formulation with Preferential and Diffuse Flow to the Water Table

Rapid infiltration following precipitation may result in groundwater contamination from surface contaminants or pathogens. In fractured rock, contaminants can migrate rapidly to points of groundwater withdrawals. In contrast to the temporal availability of groundwater quality chemical indicators, meteorological and groundwater level observations are available in real-time to estimate time-varying recharge, which can act as a surrogate to identify periods of rapid infiltration that may indicate contamination susceptibility. Estimating recharge using methods, such as base-flow recession, unsaturated infiltration models, or Water-Table Fluctuations (WTF), cannot capitalize on currently available technologies and telecommunication infrastructure to conduct real-time recharge estimation at scales relevant to characterizing rapid infiltration. We present a linear, physics-based State-Space (SS) model of one-dimensional infiltration to estimate recharge, which includes preferential and diffuse-flow to the water table. The model can take advantage of real-time data for water-table altitude, precipitation, and evapotranspiration. Model parameters are calibrated over an observation period, and the Kalman Filter (KF) is subsequently applied to continuously update the observed (water-table altitude) and unobserved (groundwater recharge) system states and predict future states as new data become available. The SS/KF algorithm is demonstrated at the Masser Groundwater Recharge Site in Pennsylvania, USA and comparisons are made with recharge estimates from WTF methods. Model results indicate that the frequency of observations (daily versus sub-daily) dictates the allocation between preferential and diffuse flow. Additionally, because infiltration processes encompass many nonlinearities, model parameters estimated from observation periods need to be updated at least seasonally to account for changing recharge conditions.

groundwater, Recharge, infiltration, recursive est↗
DOE OSTI2024

Daily water stable isotopes, transpiration, and matrix potential data for an aspen and engelmann stand in the East River Watershed (version 2)

We provide daily stable isotope (2H & 18O) ratios in soil water and xylem (plant stem) water, as well as the sap flow (transpiration) and the soil's matric potential at a forested site near Gothic, Colorado, in the East River catchment. We measured the stable isotopic composition of the transpiration and the daily transpiration flux sum of three aspen and three engelmann spruce. In both forest stands, we installed a soil profile and measured the soil matric potential at 15, 30, and 60 cm depth as well as the stable isotopes of soil pore water at 5, 10, 30, 60, and 90 cm depths. All isotope measurements were done in situ via vapor probes connected to a cavity ring down spectrometer (Picarro L1240i).We further report the daily meteorological data observed at billy barr near our study site. We also provide for each tree the relative share of root water uptake derived from the isotope measurements via a Bayesian mixing model (MixSIAR).The daily data is provided as a time series in "Iso_MP_Sap_DataDaily_ESSDiveUpload.csv" and the units are provided in "dd.csv"; the location of the instrumented trees and soil profiles are given as latitude and longitude coordinates saved as CSV and KMZ files; and a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata.The data was gathered to investigate the short-term changes of the water sources (i.e., variation of root water uptake from different soil depths) of the studied subalpine trees.Update 07/16/2025: The relative and absolute plant water uptake depths were grouped to ensure that the MixSIAR model was applied with endmembers that differed in their d2H value by at least 3 permill and at least 1 permill for d18O. Whenever the difference between observed d2H values for two or more probes at neighboring depths was less than 3 permill, we used the average value for the source water endmember. For days at which probes that were not next to each other measurements did not differ at least 3 permill, the average of all probes between these two depths was used as the water source endmember representing the depth range between these two probes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2007

The METOP-A Orbit Acquisition Strategy and its LEOP Operational Experience

Europe's first polar-orbiting weather satellite, METOPA, was launched by a Soyuz launcher from Baikonur Cosmodrome on the 19th of October of 2006. The routine operations of METOP-A are conducted by EUMETSAT (European Organization for Exploitation of Meteorological Satellites) in the frame of the European Polar System mission (EPS). The METOP-A Launch and Early Orbit Phase (LEOP) operations have been performed by ESA/ESOC. The Flight Dynamics Orbit Determination and Control team (OD&C) at ESOC was in charge of correcting the S/C orbit as delivered by the launcher in such a way that EUMETSAT would be able to acquire the reference orbit with a drift-stop manoeuvre approximately two weeks after a LEOP of 3 days and Hand-Over to the EUMETSAT Control Centre (EUMETSAT-CC) in Darmstadt, Germany. The various strict constraints and the short amount of time available for ESOC operations made this task challenging. Several strategies were prepared before launch and analysed during LEOP based on the achieved injection orbit. This paper presents the different manoeuvre strategies investigated and finally applied to acquire the operational orbit, reporting as well the details of its execution and final achieved state.

Merz, K.↗
DOE OSTI2025

Examining future changes in coastal low-level jet properties offshore California through dynamical downscaling

The coastal low-level jet, or coastal low-level jet (CLLJ), is a synoptically-forced meteorological feature frequently present offshore the western United States (U.S.). Characterized by a wind speed maximum that resides at the top of the marine boundary layer, the CLLJ is largely controlled by the location and strength of the North Pacific High (NPH) as well as the coastal geometry. Considering the rich wind resource available in this offshore region, the Bureau of Ocean Energy Management identified wind energy lease areas offshore California and supported the deployment of two U.S. Department of Energy wind lidar buoys near Morro Bay and Humboldt. Despite our relatively good understanding of the fundamental mechanisms responsible for large-scale CLLJ properties offshore the western U.S., future changes in CLLJ characteristics are less clear. To address this research challenge, and ultimately to better inform future wind turbine deployments, we use simulations driven by three global climate models (GCMs). We apply self-organizing maps to the model outputs for a historical and two future climate periods to show the range of NPH regimes that support CLLJ conditions during the warm seasons, as well as the subtle contribution from land-falling cyclones approaching the mainland during the cold seasons. Compared to the historical period, the three GCM-driven simulations agree that CLLJ conditions will become more (less) prevalent from central California northward (southward). They agree less with respect to future changes in maximum CLLJ wind speeds and CLLJ heights. However, after considering model biases present during the historical period, wind speeds between the models are actually more similar during the 2070–2095 period than during the historical period. The potential combination of more frequent CLLJ conditions characterized by relatively consistent wind speeds occurring at lower heights across northern California suggests that the Humboldt lease area may be ideal for a long-term wind turbine deployment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Stratospheric aerosol forcing for CMIP7 – Part 1: optical properties for pre-industrial, historical, and scenario simulations

Stratospheric aerosols, most of which originate from explosive volcanic sulfur emissions into the stratosphere, are a key natural driver of climate variability. They are thus a forcing provided by the Coupled Model Intercomparison Project (CMIP) Climate Forcings Task Team to climate modelling groups participating in phase 7 of CMIP. For the historical period, we provide two datasets covering 1750–2023: (i) a volcanic upper tropospheric-stratospheric sulfur emission dataset, documented in a companion paper; and (ii) a stratospheric sulfate aerosol optical property dataset, documented in this paper. For the satellite era (from 1979 onwards), stratospheric aerosol optical properties are derived from the Global Space-based Stratospheric Aerosol Climatology (GloSSAC) dataset. For the pre-satellite era (1750–1978), optical properties are derived from our volcanic SO 2 emission dataset using a new version of the reduced-complexity volcanic aerosol model Easy Volcanic Aerosol (Height) (EVA_H v2). A background, non-volcanic stratospheric aerosol climatology is derived from the 1998–2001 period with a trend over 1850–1978 accounting for increasing anthropogenic aerosols. A monthly stratospheric aerosol climatology is derived from the 1850–2021 average for both pre-industrial and scenario (future) simulations, with a 9-year ramp over 2022–2030 for scenario simulations to ensure a smooth transition from the historical period. CMIP7 simulations should use version 2.2.1 of our pre-industrial and historical datasets, and version 2.2.2 of our scenario dataset. Our methodology to produce historical aerosol optical properties significantly differs from CMIP6 for the pre-satellite era, and the resulting forcings in turn largely differ. In particular, the CMIP6 dataset was mostly based on the sparse and uncertain pyrheliometer record, which resulted in strongly underrepresented emissions from small-to-moderate magnitude eruptions. The resulting bias is addressed in CMIP7, which is entirely emission-derived in the pre-satellite era and uses more recent ice-core-based volcanic sulfur emission inventories than CMIP6. Our approach results in an overall larger volcanic aerosol forcing for CMIP7, with the 1850–2014 mean mid-visible global mean stratospheric aerosol optical depth (SAOD) in CMIP7 (0.0138) being 29 % higher than in CMIP6 (0.0107). The pre-industrial mean of the same variable is 26 % higher in CMIP7 (0.0135, derived from the historical 1850–2021) than CMIP6 (0.0107, derived from the historical 1850–2014 mean). Using a reduced-complexity climate model, we simulate a global mean surface temperature that is 0.07 °C colder for 1850–1900 when using the CMIP7 dataset instead of CMIP6, whereas 2000–2014 is 0.03 °C warmer in CMIP7. Our dataset also exhibits lower forcing for 1960–1980, resulting in temperatures 0.06 °C warmer when averaged over 1960–1990, a period for which CMIP6 climate models exhibit a cold bias. Given the large uncertainties characterizing the dataset, in particular for the pre-satellite era, we advise against treating the CMIP7 or CMIP6 dataset as uniquely superior for any specific year and highlight the need for further evaluation. We conclude the study by discussing sources of uncertainty for the dataset, future research avenues to improve it, as well as requirements to operationalize the production of the dataset, i.e. extend it and update it on an annual basis instead of every 5–7 years following CMIP cycles.

Aubry, Thomas J. [Univ. of Exeter, Devon (United K↗
NASA NTRSDate not supplied

Overview of GSICS visible calibration methods.

Global Space-based Inter-Calibration System (GSICS) is an international collaborative effort initiated in 2005 by the World Meteorological Organization (WMO) and the Coordination Group for Meteorological Satellites (CGMS) to 1-monitor, 2-improve and 3-harmonize the quality of observations from operational weather and environmental satellites of the Global Observing System (GOS). GSICS aims at ensuring consistent accuracy among space-based observations worldwide for climate monitoring, weather forecasting, and environmental applications by 1-monitoring instrument performance, 2-operational inter-calibration of satellite instruments, 3-radiometrically scaling observations to absolute reference standards, 4-recalibration of archived datasets. Although many of the GSIC visible inter-calibration methods are based on medium pixel-level imagers such as MODIS, VIIRS, AVHRR, and geostationary sensors, the methods are applicable to high spatial resolution imagery such as Landsat. For sensor inter-calibration efforts, coincident ray-matched reflectance pairs are used to transfer the calibration from the reference to the target instrument. Spectral band adjustment factors (SBAF) are applied to account for any spectral band differences. For sensor stability monitoring, PICS, deep convective cloud (DCC) and polar ice Earth viewed targets are utilized. GSICS uses the moon for stability monitoring by using the GSICS Implementation of the USGS Robotic Lunar Observatory (ROLO) model (GIRO). GSICS is looking forward to the launch of CLARREO on the international space station (ISS) in 2023. The CLARREO hyperspectral reflective solar band (RSB) sensor will establish an absolute calibration reference in space, which will characterize invariant targets and will inter-calibrate concurrent sensors. Examples of these calibration methods will be presented at the workshop.

Visible calibration methods↗
DOE OSTI2025

Best estimate of the planetary boundary layer height from multiple remote sensing measurements

Remote sensing measurements have been widely used to estimate the planetary boundary layer height (PBLHT). Each remote sensing approach offers unique strengths and faces different limitations. In this study, we use machine learning (ML) methods to produce a best-estimate PBLHT (PBLHT-BE-ML) by integrating four PBLHT estimates derived from remote sensing measurements at the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) observatory. Three ML models – random forest (RF) classifier, RF regressor, and light gradient-boosting machine (LightGBM) – were trained on a dataset from 2017 to 2023 that included radiosonde, various remote sensing PBLHT estimates, and atmospheric meteorological conditions. Evaluations indicated that PBLHT-BE-ML from all three models improved alignment with the PBLHT derived from radiosonde data (PBLHT-SONDE), with LightGBM demonstrating the highest accuracy under both stable and unstable boundary layer conditions. Feature analysis revealed that the most influential input features at the SGP site were the PBLHT estimates derived from (a) potential temperature profiles retrieved using Raman lidar (RL) and atmospheric emitted radiance interferometer (AERI) measurements (PBLHT-THERMO), (b) vertical velocity variance profiles from Doppler lidar (PBLHT-DL), and (c) aerosol backscatter profiles from micropulse lidar (PBLHT-MPL). The trained models were then used to predict PBLHT-BE-ML at a temporal resolution of 10 min, effectively capturing the diurnal evolution of PBLHT and its significant seasonal variations, with the largest diurnal variation observed over summer at the SGP site. We applied these trained models to data from the ARM Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) field campaign (EPC), where the PBLHT-BE-ML, particularly with the LightGBM model, demonstrated improved accuracy against PBLHT-SONDE. Analyses of model performance at both the SGP and EPC sites suggest that expanding the training dataset to include various surface types, such as ocean and ice-covered areas, could further enhance ML model performance for PBLHT estimation across varied geographic regions.

Zhang, Damao [Pacific Northwest National Laborator↗
NASA NTRSDate not supplied

Performance of Heterogeneous Sulfate Aerosol Formation Mechanisms Against Aircraft and Ground-Based Observations in Seoul, South Korea

Haze events are a common occurrence in East Asia, where local meteorology combines with anthropogenic emissions under certain conditions to create rapid increases in secondary aerosol. Often these events are driven by rapid production of secondary sulfate-nitrate-ammonium (SIA) aerosol. Models have typically included simple parameterizations of the production of sulfate during these events which vary widely in their magnitude. There are a range of different mechanistic descriptions of the chemical processes occurring during these events in large part due to poor constraints on aerosol pH and transition metal content, both of which are important to determining the key pathway(s) driving aerosol production. Here, we use the GEOS-Chem chemical transport model to interpret detailed ground and aircraft observations from May and June 2016 during the Korea-United States Air Quality (KORUS-AQ) field study to evaluate the mechanistic pathways governing sulfate production during a haze event. We apply our model of sulfate production mechanisms to ground-based speciated aerosol measurements in winter, summer, and fall in Seoul.

Katherine R Travis↗
NASA NTRS1993

Modeling and data assimilation for Mars Observer

The Mars Observer mission will return a substantial set of meteorological data for the atmosphere of Mars, principally in the form of thermal infrared soundings from the Pressure Modulator Infrared Radiometer (PMIRR). The dataset will enable global studies of a range of dynamical and physical phenomena in the martian atmosphere. The solar orbital configuration of Mars Observer, however, means that the coverage of the martian atmosphere and surface by remote sounding instruments will be asynchronous. This leads to difficulties of interpretation when studying synoptic-scale phenomena with timescales of a few days. In an attempt to overcome such problems we propose to use data assimilation techniques, as currently employed for operational weather forecasting on the Earth, in conjunction with a Martian General Circulation Model (MGCM) that is under development at Oxford and Reading Universities. An overview of the model and assimilation scheme will be given, with reference to problems encountered so far in applying such methods to remote-sensing observations of the martian atmosphere.

Lewis, Stephen R.↗
NASA NTRS2016

Evaluation of Different MODIS AOD Retrieval Algorithms for PM (sub 2.5) Estimation in the Western, Midwestern and Southeastern United States with Implications for Public Health

This study was part of the research activities of the Center for Applied Atmospheric Research and Education (CAARE) funded by the NASA MUREP (Minority University Research and Education Project) Institutional Research Opportunity (MIRO) Program. Satellite measurements of Aerosol Optical Depth (AOD) have been shown to be correlated with ground measurements of fine particulate matter less than 2.5 microns PM (sub 2.5), which in turn has been linked to respiratory and heart diseases. The strength of the correlation between AOD and PM (sub 2.5) varies for different AOD retrieval algorithms and geographic regions. We evaluated several Moderate Resolution Imaging Spectrometer (MODIS) AOD products from different satellites (Aqua vs. Terra), retrieval algorithms (Dark Target versus Deep Blue), Collections (5.1 versus 6) and spatial resolutions (10-kilometers versus 3-kilometers) for cities in the Western, Midwestern and Southeastern U.S. We developed and validated PM (sub 2.5) prediction models using remotely-sensed AOD data, which were improved by incorporating meteorological variables (temperature, relative humidity, precipitation, wind speed, and wind direction) from the North American Land Data Assimilation System Phase 2 (NLDAS-2). Adding these meteorological data significantly improved the predictive power of all the PM (sub 2.5) models, especially in the Western U.S. Temperature, relative humidity and wind speed were the most significant meteorological variables throughout the year in the Western U.S. Wind speed was the most significant meteorological variable for the cold season while temperature was the most significant variable for the warm season in the Midwestern and Southeastern U.S. Our study re-establishes the connection between PM (sub 2.5) and public health concerns including respiratory and cardiovascular diseases (asthma, high blood pressure, coronary heart disease, heart attack, and stroke). Using PM (sub 2.5) data and health data from the Centers for Disease Control and Prevention (CDC)'s Behavioral Risk Factor Surveillance System (BRFSS), our statistical analysis showed that heart attack and stroke occurrences had the strongest correlations with PM (sub 2.5).

Environmental Public Health↗
NASA NTRS1994

Ice surface temperature retrieval from AVHRR, ATSR, and passive microwave satellite data: Algorithm development and application

One essential parameter used in the estimation of radiative and turbulent heat fluxes from satellite data is surface temperature. Sea and land surface temperature (SST and LST) retrieval algorithms that utilize the thermal infrared portion of the spectrum have been developed, with the degree of success dependent primarily upon the variability of the surface and atmospheric characteristics. However, little effort has been directed to the retrieval of the sea ice surface temperature (IST) in the Arctic and Antarctic pack ice or the ice sheet surface temperature over Antarctica and Greenland. The reason is not one of methodology, but rather our limited knowledge of atmospheric temperature, humidity, and aerosol vertical, spatial and temporal distributions, the microphysical properties of polar clouds, and the spectral characteristics of snow, ice, and water surfaces. Over the open ocean the surface is warm, dark, and relatively homogeneous. This makes SST retrieval, including cloud clearing, a fairly straightforward task. Over the ice, however, the surface within a single satellite pixel is likely to be highly heterogeneous, a mixture of ice of various thicknesses, open water, and snow cover in the case of sea ice. Additionally, the Arctic is cloudy - very cloudy - with typical cloud cover amounts ranging from 60-90 percent. There are few observations of cloud cover amounts over Antarctica. The goal of this research is to increase our knowledge of surface temperature patterns and magnitudes in both polar regions, by examining existing data and improving our ability to use satellite data as a monitoring tool. Four instruments are of interest in this study: the AVHRR, ATSR, SMMR, and SSM/I. Our objectives are as follows. Refine the existing AVHRR retrieval algorithm defined in Key and Haefliger (1992; hereafter KH92) and applied elsewhere. Develop a method for IST retrieval from ATSR data similar to the one used for SST. Further investigate the possibility of estimating surface temperature from passive microwave data (in conjunction with AVHRR clear sky samples) through the use of 'effective emissivities' and physical relationships between skin temperature and subsurface temperature. Use the general method outlined in MK93 to calculate a 12-year record of clear sky equivalent surface temperatures, or possibly all-sky snow-ice interface physical temperatures, from SMMR and SSM/I, compare these temperatures to climatologies, ECMWF modeled surface temperatures, and surface temperatures predicted by a 2-D ice model. And intercompare several ice surface retrieval methods and validate them against ground measurements from the Swiss Camp on the Greenland ice sheet. Additionally, we intend to develop a surface temperature product based on AVHRR data and possibly blended with drifting buoy and meteorological station temperatures.

Key, Jeff↗
DOE OSTI2020

EBSD seed LDRD project: Does Corona Virus – 2019 (COVID-19) and Seasonal Flu have similar meteorology and air quality controls driving their spread?

Seasonal influenza and Influenza like Illnesses (ILI) pose a serious public health risk and in turn affect the economy. Various factors affect ILI cases and mortality, including the pathogen and its interaction with the host, as well as environmental and socioeconomic factors such as meteorology, household structure, air pollution, urbanization, and population. Despite the growing number of studies on influenza and ILI, challenges remain in forecasting the timing of seasonal onset, outbreak patterns, and key factors affecting transmission. In particular, the impacts of meteorology and air quality on ILI have been challenging to understand, with linear regression studies focused on different geographic regions producing contradictory results. For example, influenza seasonality has been associated with cold-dry conditions in temperate mid-latitudes, but with humid-rainy conditions in tropical climates. These apparently contradictory results imply that the relationships between influenza cases and atmospheric variables may be too complex to be captured by linear regression models. In this seed project, we analyze meteorology and air quality variables from numerical models to determine which atmospheric variables are most helpful in predicting weekly changes in recorded flu cases. In contrast to most previous studies that relied on linear regression analysis to predict the timing of the flu onset or peak, we employ a robust machine learning algorithm to evaluate the contribution of atmospheric variables to weekly changes in recorded ILI cases. These results may also be relevant to the spread of other respiratory illnesses such as Corona Virus Disease – 2019 (COVID-19).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Effects of random forest modeling decisions on biogeochemical time series predictions

Abstract Random forests (RF) are an increasingly popular machine learning approach used to model biogeochemical processes in the Earth system. While RF models are robust to many assumptions that complicate deterministic models, there are several important parameterization decisions for appropriate use and optimal model fit. We explored the role that parameter decisions, including training/testing data splitting strategies, variable selection, and hyperparameters play on RF goodness‐of‐fit by constructing models using 1296 unique parameter combinations to predict concentrations of nitrate, a key nutrient for biogeochemical cycling in aquatic ecosystems. Models were built on long‐term, publicly available water quality and meteorology time series collected by the National Estuarine Research Reserve monitoring network for two contrasting ecosystems representing freshwater and brackish estuaries. We found that accounting for temporal dependence when splitting data into training and testing subsets was key for avoiding over‐estimation of model predictive power. In addition, variable selection, the ratio of training to testing data, and to a lesser degree, variables per split and number of trees, were significant parameters for optimizing RF goodness‐of‐fit. We also explored how model parameter decisions influenced interpretation of the relative importance of predictors to the model, and model predictor‐dependent variable relationships, with results suggesting that both data structure and model parameterization influence these factors. Because much of the current RF literature is written for the computational and statistical science communities, the primary goal of this study is to provide guidelines for aquatic scientists new to machine learning to apply RF techniques appropriately to aquatic biogeochemical datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Assessing the large-scale drivers of precipitation in the northeastern United States via linear orthogonal decomposition

This study examines the linear orthogonal modes associated with monthly precipitation in the northeastern United States, from CESM1 LENS (35 ensemble members, 1979–2005) and two reanalysis datasets (ERA5, 1950–2018 and NOAA-CIRES-DOE 20CRv3, 1950–2015). Calendar months are aggregated together, and any linear trends in data are removed. Using region-averaged precipitation anomaly time series and monthly anomalies for several global 2D atmospheric fields, a linear orthogonal decomposition method is implemented to iteratively extract time series (based on field and geographic location) of absolute maximum correlation. Linear modes associated with this method are then projected onto the full set of 2D fields to provide physical insight into the mechanisms involved in generating precipitation. In this region, the first mode is associated with vapor transport from the Atlantic seaboard, the second mode is characterized by westward vapor transport associated with extratropical cyclones, and the third mode captures vapor transport from the Gulf of Mexico during the fall and winter. However, the third mode is less robust in the spring and summer. Results are generally consistent across the datasets, and applying multiple linear regression with the linear modes to predict the precipitation anomalies produces R-squared values of around 0.54–0.65 for CESM1 LENS, and around 0.58–0.88 for reanalysis, with the lowest values generally in the spring and late summer. The influence of low-frequency climate variability on the modes is considered for CESM1 LENS, and the modes in late winter can be predicted with some success via a combination of several, prominent large-scale teleconnection patterns.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Validation of Atmosphere/Ionosphere Signals Associated with Major Earthquakes by Multi-Instrument Space-Borne and Ground Observations

The latest catastrophic earthquake in Japan (March 2011) has renewed interest in the important question of the existence of pre-earthquake anomalous signals related to strong earthquakes. Recent studies have shown that there were precursory atmospheric/ionospheric signals observed in space associated with major earthquakes. The critical question, still widely debated in the scientific community, is whether such ionospheric/atmospheric signals systematically precede large earthquakes. To address this problem we have started to investigate anomalous ionospheric / atmospheric signals occurring prior to large earthquakes. We are studying the Earth's atmospheric electromagnetic environment by developing a multisensor model for monitoring the signals related to active tectonic faulting and earthquake processes. The integrated satellite and terrestrial framework (ISTF) is our method for validation and is based on a joint analysis of several physical and environmental parameters (thermal infrared radiation, electron concentration in the ionosphere, lineament analysis, radon/ion activities, air temperature and seismicity) that were found to be associated with earthquakes. A physical link between these parameters and earthquake processes has been provided by the recent version of Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model. Our experimental measurements have supported the new theoretical estimates of LAIC hypothesis for an increase in the surface latent heat flux, integrated variability of outgoing long wave radiation (OLR) and anomalous variations of the total electron content (TEC) registered over the epicenters. Some of the major earthquakes are accompanied by an intensification of gas migration to the surface, thermodynamic and hydrodynamic processes of transformation of latent heat into thermal energy and with vertical transport of charged aerosols in the lower atmosphere. These processes lead to the generation of external electric currents in specific regions of the atmosphere and the modifications, by dc electric fields, in the ionosphere-atmosphere electric circuit. We retrospectively analyzed temporal and spatial variations of four different physical parameters (gas/radon counting rate, lineaments change, long-wave radiation transitions and ionospheric electron density/plasma variations) characterizing the state of the lithosphere/atmosphere coupling several days before the onset of the earthquakes. Validation processes consist in two phases: A. Case studies for seven recent major earthquakes: Japan (M9.0, 2011), China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and, Chile (M8.8, 2010) and B. A continuous retrospective analysis was preformed over two different regions with high seismicity- Taiwan and Japan for 2003-2009. Satellite, ground surface, and troposphere data were obtained from Terra/ASTER, Aqua/AIRS, POES and ionospheric variations from DEMETER and COSMIC-I data. Radon and GPS/TEC were obtaining from monitoring sites in Taiwan, Japan and Italy and from global ionosphere maps (GIM) respectively. Our analysis of ground and satellite data during the occurrence of 7 global earthquakes has shown the presence of anomalies in the atmosphere. Our results for Tohoku M9.0 earthquake show that on March 7th, 2011 (4 days before the main shock and 1 day before the M7.2 foreshock of March 8, 2011) a rapid increase of emitted infrared radiation was observed by the satellite data and an anomaly was developed near the epicenter. The GPS/TEC data indicate an increase and variation in electron density reaching a maximum value on March 8. From March 3 to 11 a large increase in electron concentration was recorded at all four Japanese ground-based ionosondes, which returned to normal after the main earthquake. Similar approach for analyzing atmospheric and ionospheric parameters has been applied for China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and Chile (M8.8, 2010) eahquakes. Results have revealed the presence of related variations of these parameters implying their connection with the earthquake process. The second phase (B) of this validation included 102 major earthquakes (M>5.9) in Taiwan and Japan. We have found anomalous behavior before all of these events with no false negatives. False alarm ratio for false positives is less then 10% and has been calculated for the same month of the earthquake occurrence for the entire period of analysis (2003-2009). The commonalities for detecting atmospheric/ionospheric anomalies are: i.) Regularly appearance over regions of maximum stress (i.e., along plate boundaries); ii.) Anomaly existence over land and sea; and iii) association with M>5.9 earthquakes not deeper than 100km. Due to their long duration over the same region these anomalies are not consistent with a meteorological origin. Our initial results from the ISTF validation of multi-instrument space-borne and ground observations show a systematic appearance of atmospheric anomalies near the epicentral area, one to seven (average) days prior to the largest earthquakes, and suggest that it could be explained by a coupling process between the observed physical parameters and the pre-earthquake preparation processes.

Ouzounov, Dimitar↗