Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “temperature forecasting”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Diagnosing the Nature of Land-Atmosphere Coupling During the 2006-7 Dry/Wet Extremes in the U.S. Southern Great Plains

Land-atmosphere interactions play a critical role in determining the. diurnal evolution of both planetary boundary layer (PBL) and land surface temperature and moisture states. The degree of coupling between the land surface and PBL in numerical weather prediction and climate models remains largely unexplored and undiagnosed due to the complex interactions and feedbacks present across a range of scales. Further, uncoupled systems or experiments (e.g., the Project for Intercomparison of Land Parameterization Schemes, PILPS) may lead to inaccurate water and energy cycle process understanding by neglecting feedback processes such as PBL-top entrainment. In this study, a framework for diagnosing local land-atmosphere coupling (LoCo) is presented using a coupled mesoscale model with a suite of PBL and land surface model (LSM) options along with observations during the summers of 200617 in the U.S. Southern Great Plains. Specifically, the Weather Research and Forecasting (WRF) model has been coupled to NASA's Land Information System (LIS), which provides a flexible and high-resolution representation and initialization of land surface physics and states. A range of diagnostics exploring the links and feedbacks between soil moisture and precipitation are examined for the dry/wet extremes of this region, along with the sensitivity of PBL-LSM coupling to perturbations in soil moisture. As such, this methodology provides a potential pathway to study factors controlling local land-atmosphere coupling (LoCo) using the LIS-WRF system, which is serving as a testbed for LoCo experiments to evaluate coupling diagnostics within the community.

Santanello, Joseph A., Jr.↗
DOE OSTI2021

Exploring Sources of Surface Bias in HRRR Using New York State Mesonet

In recent years, there has been increasing demand for applications of short-term forecasting of renewable energy potential and assessments of the likelihood of extreme weather events using the High-Resolution Rapid Refresh (HRRR) model. Examining the biases in the newest version of HRRR is necessary to promote further model development. Using data from one of the most comprehensive and dense monitoring networks, New York State Mesonet (NYSM), we evaluate the HRRR version 3 meteorological fields for an entire year. In this work, the land-atmosphere-cloud coupling system is evaluated as an integrated whole. We investigate the physical processes influencing the soil hydrological balance and the thermodynamic interactions, from surface fluxes up to the level of boundary layer convection from both temporal (seasonal and diurnal) and spatial perspectives. Results show that the model 2 m temperature and humidity biases are seasonally dependent, with warm and dry bias present during the warm season, and an extreme nocturnal cold bias in winter. The summer warm bias includes both a land-surface-induced bias and a cloud-induced bias. Inaccurate representation of energy partition and soil hydrological process across different land use types as well as a hydrological bias in describing spring snowmelt are identified as the main source of the land-surface-induced bias. A feedback loop linking cloud presence, flux changes, and temperature contributes to the cloud-induced bias. The positive solar radiation bias increases from clear sky to overcast sky conditions. The most significant bias occurs during overcast and thick cloud conditions associated with frontal passage and thunderstorms.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1983

Low level water vapor fields from the VISSR Atmospheric Sounder (VAS) 'split window' channels

The use of the 11 and 12-micron IR-radiometer channels of the VISSR Atmosphere Sounder (VAS) on the GOES to detect water vapor in the lowest 300-400 mb of the troposphere is reported. An algorithm is developed to eliminate the background temperature, allowing the calculation of precipitable water (PW) over both land and water via a single-layer radiative model. This 'split-window' method is demonstrated in a case study, covering the US on July 13, 1981. PW values were calculated from five VAS images and compared with those from radiosonde data and surface measurements. It is shown that the VAS PW images have good resolution (15 km), reveal vivid, continuously evolving details, and differentiate deep, convection-supporting layers from shallow ones. PW ranged from 1.7 to 5.5 g/sq cm (+ or - 1.0 g/sq cm), in good agreement with other measurements. Since this method can detect mesoscale water-vapor fields in relatively clear air, it is considered of great potential value for numerical forecasting.

Robinson, W. D.↗
DOE OSTI2020

Visual Orbits of Spectroscopic Binaries with the CHARA Array. III. HD 8374 and HD 24546

We present the visual orbits of two long-period spectroscopic binary stars, HD 8374 and HD 24546, using interferometric observations acquired with the CHARA Array and the Palomar Testbed Interferometer. We also obtained new radial velocities from echelle spectra using the APO 3.5 m and Fairborn 2.0 m telescopes. By combining the visual and spectroscopic observations, we solve for the full, three-dimensional orbits and determine the stellar masses and distances to within 3% uncertainty. We then estimate the effective temperature and radius of each component star through Doppler tomography and spectral energy distribution analyses, in order to compare the observed stellar parameters to the predictions of stellar evolution models. For HD 8374, we find masses of M {sub 1} = 1.636 ± 0.050M {sub ⊙} and M {sub 2} = 1.587 ± 0.049M {sub ⊙}, radii of R {sub 1} = 1.84 ± 0.05R {sub ⊙} and R {sub 2} = 1.66 ± 0.12R {sub ⊙}, temperatures of T{sub eff1}=7280±110 K and T{sub eff2}=7280±120 K, and an estimated age of 1.0 Gyr. For HD 24546, we find masses of M {sub 1} = 1.434 ± 0.014M {sub ⊙} and M {sub 2} = 1.409 ± 0.014M {sub ⊙}, radii of R {sub 1} = 1.67 ± 0.06R {sub ⊙} and R {sub 2} = 1.60 ± 0.10R {sub ⊙}, temperatures of T{sub eff1}=6790±120 K and T{sub eff2}=6770±90 K, and an estimated age of 1.4 Gyr. HD 24546 is therefore too old to be a member of the Hyades cluster, despite its physical proximity to the group.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1998

A common-sense climate index: is climate changing noticeably?

We propose an index of climate change based on practical climate indicators such as heating degree days and the frequency of intense precipitation. We find that in most regions the index is positive, the sense predicted to accompany global warming. In a few regions, especially in Asia and western North America, the index indicates that climate change should be apparent already, but in most places climate trends are too small to stand out above year-to-year variability. The climate index is strongly correlated with global surface temperature, which has increased as rapidly as projected by climate models in the 1980s. We argue that the global area with obvious climate change will increase notably in the next few years. But we show that the growth rate of greenhouse gas climate forcing has declined in recent years, and thus there is an opportunity to keep climate change in the 21st century less than "business-as-usual" scenarios.

Time↗
DOE OSTI2021

Assessing dynamic vegetation model parameter uncertainty across Alaskan arctic tundra plant communities

Abstract As the Arctic region moves into uncharted territory under a warming climate, it is important to refine the terrestrial biosphere models (TBMs) that help us understand and predict change. One fundamental uncertainty in TBMs relates to model parameters, configuration variables internal to the model whose value can be estimated from data. We incorporate a version of the Terrestrial Ecosystem Model (TEM) developed for arctic ecosystems into the Predictive Ecosystem Analyzer (PEcAn) framework. PEcAn treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We examined how variation in 21 parameters in the equation for gross primary production influenced model sensitivity and uncertainty in terms of two carbon fluxes (net primary productivity and heterotrophic respiration) and two carbon (C) pools (vegetation C and soil C). We set up different parameterizations of TEM across a range of tundra types (tussock tundra, heath tundra, wet sedge tundra, and shrub tundra) in northern Alaska, along a latitudinal transect extending from the coastal plain near Utqiaġvik to the southern foothills of the Brooks Range, to the Seward Peninsula. TEM was most sensitive to parameters related to the temperature regulation of photosynthesis. Model uncertainty was mostly due to parameters related to leaf area, temperature regulation of photosynthesis, and the stomatal responses to ambient light conditions. Our analysis also showed that sensitivity and uncertainty to a given parameter varied spatially. At some sites, model sensitivity and uncertainty tended to be connected to a wider range of parameters, underlining the importance of assessing tundra community processes across environmental gradients or geographic locations. Generally, across sites, the flux of net primary productivity (NPP) and pool of vegetation C had about equal uncertainty, while heterotrophic respiration had higher uncertainty than the pool of soil C. Our study illustrates the complexity inherent in evaluating parameter uncertainty across highly heterogeneous arctic tundra plant communities. It also provides a framework for iteratively testing how newly collected field data related to key parameters may result in more effective forecasting of Arctic change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling Volatility-Based Aerosol Phase State Predictions in the Amazon Rainforest

Organic aerosol (OA) is a complex matrix of various constituents—fresh (primary organic aerosols—POA) and aged via oxidation (secondary organic aerosols—SOA), generated from biogenic, anthropogenic, and biomass burning sources. The viscosity of OA can be critical in influencing new particle formation, reactive uptake processes that impact evaporation-growth kinetics, and the lifetime of particles in the atmosphere. We report this work utilizes a well-defined relationship between volatility and viscosity for pure compounds, which we incorporated within the Weather Research and Forecasting Model coupled to chemistry (WRF-Chem) to simulate the phase state and viscosity of bulk OA during the dry-to-wet transition season (September–October) in the Amazon rainforest during 2014. Our simulations indicate spatial and temporal heterogeneity in aerosol phase state often not captured by global-scale models. We show the strong role of water associated with organic aerosol ($w_s$) as the dominant factor that can be used to quantitatively estimate OA viscosity. Analysis of WRF-Chem simulations across the entire atmospheric column indicates a strong inverse log-linear relationship between $w_s$ and OA viscosity with a correlation coefficient approaching 1, in the background and biomass burning-influenced conditions. At high altitudes where relative humidity (RH) and temperatures are low, our simulations indicate that OA exists in a semisolid-/solid-like phase state, consistent with previous studies. OA hygroscopicity is strongly correlated (ca. –0.8) with OA viscosity at RH ca. 30–50%, but this RH range is found mostly at low OA concentrations and the middle troposphere (ca. 6–10 km altitudes) in our simulated domain. OA hygroscopicity is uncorrelated with viscosity at higher-RH (near surface) and lower-RH (upper troposphere) regimes. At the urban site near surface, where day–night differences in RH are significant, RH is found to drive the phase state. At the background forested site near surface, where day–night RH differences are small, biomass burning-influenced OA is semisolid and a significant OA associated with background conditions is liquid-like. Simulations indicate a long tail of OA viscosity frequency distributions extending in the semisolid/solid regimes over background biogenic-influenced conditions due to the role of low-volatility OA components such as monoterpene oxidation products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

ERA5-Land Data for LASSO-CACTI Overview Paper

The European Centre for Medium-Range Weather Forecasts (ECMWF) generated a soil reanalysis dataset for the land component of the fifth generation of European ReAnalysis (ERA5), referred to as ERA5-Land. This is a model-generated dataset, with the original version available for the period 1950 to present. The version archived in this DOE ARM product is a subset of the data is for the period of the CACTI field campaign plus several preceding months, specifically from August 1, 2018 through March 22, 2019 with hourly intervals. The ARM copy is also a sub-region of the original global product; the ARM copy is for -60 to -5 °N by -105 to -30 °W. Only variables necessary to drive the WRF-Hydro model are included, which are the 2-m temperature and specific humidity, 10-m wind components, surface pressure, rain rate, and downward surface short and longwave radiation. These data have been obtained from the Copernicus Data Store.

10m wind u-component↗
NASA NTRS2006

Observing Ice in Clouds from Space

There are many satellite observations of cloud top properties and the liquid and rain content of clouds, however, we do not yet quantitatively understand the processes that control the water budget of the upper troposphere where ice is the predominant phase, and how these processes are linked to precipitation processes and the radiative energy budget. The ice in clouds in the upper troposphere either melts into rain or is detrained, and persists, as cirrus clouds affecting the hydrological and energy cycle, respectively. Fully modeling the Earth's climate and improving weather and climate forecasts requires accurate satellite measurements of various cloud properties at the temporal and spatial scales of cloud processes. These properties include cloud horizontal and vertical structure, cloud water content and some measure of particle sizes and shapes. The uncertainty in knowledge of these ice characteristics is reflected in the large discrepancies in model simulations of the upper tropospheric water budget. Model simulations are sensitive to the partition of ice between precipitation and outflow processes, i.e., to the parameterization of ice clouds and ice processes. One barrier to achieving accurate global ice cloud properties is the lack of adequate observations at millimeter and submillimeter wavelengths (183-874 GHz). Recent advances in instrumentation have allowed for the development and implementation of an airborne submillimeter-wave radiometer. The brightness temperatures at these frequencies are especially sensitive to cirrus ice particle sizes (because they are comparable to the wavelength). This allows for more accurate ice water path estimates when multiple channels are used to probe into the cloud layers. Further, submillimeter wavelengths offer simplicity in the retrieval algorithms because they do not probe into the liquid and near surface portions of clouds, thus requiring only one term of the radiative transfer equation (ice scattering) to relate brightness temperatures to ice. The next step is a satellite mission designed to acquire global Earth radiance measurements in the submillimeter-wave region, thus bridging the measurement gap between microwave sounders and shorter-wavelength infrared and visible sensors. This presentation provides scientific justification and an approach to measuring ice water path and particle size from a satellite platform that spans a range encompassing both the hydrologically active and radiatively active components of cloud systems.

Ackerman, S.↗
DOE OSTI2026

Hydrogenolysis Versus Hydrocracking for Polyolefin Upcycling

Global plastic production has reached 413.8 million metric tons in 2024 and is forecasted to surpass 1.2 billion metric tons by 2050. Polyolefins, mainly polyethylene (PE) and polypropylene (PP), dominate single-use packaging and account for approximately 55% of global plastic waste. The chemical inertness that makes these materials desirable for commercial applications also renders them persistent in the environment. Current recycling technologies have proven to be insufficient to divert plastic waste from landfills or environmental loss due to technical limitations and poor economic incentives. Conventional mechanical recycling is a form of downcycling, in which the polymer remelting process results in products with inferior material properties and reduced market value. Pyrolysis, a thermochemical route used for the chemical recycling of plastic wastes into refinery feedstock, requires severe conditions, typically 400–600 °C in an oxygen-free environment. The high temperature drives up energy costs and produces a wide range of poorly defined products, including undesirable light gases and heavy tars that pose challenges for downstream processing. As an alternative, catalytic chemical recycling offers a promising route for converting waste polyolefins back into value-added hydrocarbons—such as fuels, lubricant base oils, and other chemical feedstocks—at temperatures typically below 300 °C. This opinion article focuses on two dominant pathways in heterogeneous catalysis that are used to cleave the C–C bonds of polyolefins: hydrogenolysis and hydrocracking. Both catalytic pathways have the potential to selectively convert polyolefin waste into valuable fuels and chemical feedstocks under mild conditions. Notably, while hydrocracking primarily yields branched hydrocarbon products, hydrogenolysis predominantly produces linear hydrocarbons. Here, we compare the mechanisms and catalyst designs for hydrocracking and hydrogenolysis, analyze critical technical challenges from catalyst stability to process engineering, and provide an outlook on how these complementary pathways can be used to repurpose plastic waste into valuable products.

Zhang, Ruoxi [Iowa State Univ., Ames, IA (United S↗
DOE OSTI2025

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Wave

The U.S. Department of Energy and the National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence (AI) data centers and other variable loads. This dataset entry describes hydrogen production using a single, simulated wave energy conversion device. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel Hydrogen. While the unit supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack. For the wave energy, NLR used a wave energy converter model from PacWave. These devices can be equipped with accumulators and pressure relief values to smooth the power output by storing and releasing hydraulic energy. Using a peak power output of 10 MW, the model created two 25-minute profiles: one with and one without the accumulators and pressure relief valves. To down select the profile data from the native resolution of 20 Hz to 1 Hz, NLR took the mean of every 20 data points. NLR experimented with two simulated wave energy power plants: one that peaks at 10 MW, and one that peaks at 5 MW. These profiles were scaled for the physical 1.25 MW electrolyzer by multiplying the original profiles by one eighth and one quarter, respectively. The first profile matches the capacity rating of eight of the 1.25 MW electrolyzers, while the second matches four electrolyzers. Finally, NLR experimented with two settings for the electrolyzer power supply minimum and maximum current ramp rates (gain and slew): 200 and 400 amperes per second. The simulated profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single wave electrolysis experiment and is formatted as follows: {technology}-{accumulator?}_{number of 1.25 MW electrolyzers connected}-{electrolyzer ramp rate in amperes/second} For instance, “wavePacWave-Noacc_4-400.zip” represents the 25 minute-long experiment using the PacWave’s wave energy converter model, equipped with no accumulator, connected to four 1.25-MW electrolyzers with their power supplies set to a maximum current ramp rate (gain and slew) of 400 A/s. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. An experiment, labeled “characterization_200.zip”, demonstrates the MC250 electrolyzer steady-state response with 30 minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all wave profiles combined into one dataset labeled "combined_wave_experiments.csv". NLR also built an AI/machine-learning predictive model based on these datasets. The model ingests the electrolyzer current command in amperes, as well as various pressures and temperatures across the system, and predicts hydrogen output in kilograms per hour. The complete model can be found at https://huggingface.co/NatLabRockies/ptmelt-hydrogen-electrolysis.

08 HYDROGEN↗
NASA NTRS2019

Probability Index of Low Stratus and Fog at Dawn using Dual Geostationary Satellite Observations from COMS and FY-2D near the Korean Peninsula

We developed a new remote sensing method for detecting low stratus and fog (LSF) at dawn in terms of probability index (PI) of LSF from simultaneous stereo observations of two geostationary-orbit satellites; the Korean Communication, Ocean, and Meteorological Satellite (COMS; 128.2°E); and the Chinese FengYun satellite (FY-2D; 86.5°E). The algorithm was validated near the Korean Peninsula between the months of April and August from April 2012 to June 2015, by using surface observations at 45 meteorological stations in South Korea. The optical features of LSF were estimated by using satellite retrievals and simulated data from the radiative transfer model. The PI was calculated using the combination of three satellite-observed variables: 1) the reflectance at 0.67 μm (R(0.67)) from COMS, and 2) the FY-2D R(0.67) minus the COMS R(0.67) (△R(0.67)) and (3) the FY-2D-COMS difference in the brightness temperature difference between 3.7 and 11.0 μm (ΔBTD(3.7-11)). The three variables, adopted from the top three probability of detection (POD) scores for their fog detection thresholds: △R(0.67) (0.82) > ΔBTD(3.7-11) (0.73) > R(0.67) (0.70) > BTD(3.7-11) (0.51). The LSF PI for this algorithm was significantly better in the two case studies compared to that using COMS only (i.e., R(0.67) or BTD(3.7-11)), so that this improvement was due to △R(0.67) and ΔBTD(3.7-11). Overall, PI in the LSF spatial distribution has the merits of a high detection rate, a specific probability display, and a low rate of seasonality and variability in detection accuracy. Therefore, PI would be useful for monitoring LSF in near-real-time, and to further its forecast ability, using next-generation satellites.

fog↗
DOE OSTI2024

Sensitivity of Fine‐Resolution Urban Heat Island Simulations to Soil Moisture Parameterization

ABSTRACT Urban areas experience the impact of natural disasters, such as heatwaves and flash floods, disparately in different neighbourhoods across a city. The demand for precise urban hydrometeorological and hydroclimatological modelling to examine this disparity, and the interacting challenges posed by climate change and urbanisation, has thus surged. The Weather Research and Forecasting (WRF) model has served such operational and research purposes for decades. Recent advancements in WRF, including enhanced numerical schemes and sophisticated urban atmospheric‐hydrological parameterizations, have empowered the simulation of urban geophysical processes at high resolution (~1 km), but even this resolution misses significant urban microclimate variability. This study applies the large‐eddy simulations (LES) mode within WRF, coupled with single‐layer urban canopy models (SLUCM), to enable even finer‐scale modelling (150 m) of the Urban Heat Island (UHI) effect in the Baltimore metropolitan area. We run nine scenarios to evaluate various methods of initializing soil moisture and various spinup lead times, and to assess the impact of WRF's Mosaic approach in depicting subgrid‐scale processes. We evaluate the scenarios by comparing the WRF simulated land surface temperature (LST) against Landsat LST and the WRF simulated hourly 2‐m air temperatures (AT) with observations from eight weather stations across the domain. Results underscore the paramount influence of the lead spinup time on the spatiotemporal distribution of simulated soil moisture, consequently shaping WRF's efficacy in predicting the UHI. Furthermore, interpolating soil moisture‐related parameters from the parent for child domain initialization yields a notable reduction in mean and root‐mean‐squared errors. This improvement was particularly evident in simulations with the longest spinup time, affirming the importance of carefully designing the initialization of soil moisture for improved urban temperature predictions.

Talebpour, Mahdad↗
DOE OSTI2023

Global seamless tidal simulation using a 3D unstructured-grid model (SCHISM v5.10.0)

We present a new 3D unstructured-grid global ocean model to study both tidal and nontidal processes, with a focus on the total water elevation. Unlike existing global ocean models, the new model resolves estuaries and rivers down to ~8 m without the need for grid nesting. The model is validated with both satellite and in situ observations for elevation, temperature, and salinity. Tidal elevation solutions have a mean complex root-mean-square error (RMSE) of 4.2 cm for M2 and 5.4 cm for all five major constituents in the deep ocean. The RMSEs for the other four constituents, S2, N2, K1, and O1, are, respectively, 2.05, 0.93, 2.08, and 1.34 cm). The nontidal residual assessed by a tide gauge dataset (GESLA) has a mean RMSE of 7 cm. For the first time ever, we demonstrate the potential for seamless simulation on a single mesh from the global ocean into several estuaries along the US West Coast. The model is able to accurately capture the total elevation, even at some upstream stations. The model can therefore potentially serve as the backbone of a global tide surge and compound flooding forecasting framework.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2019

Impact of Assimilating Adaptively Thinned AIRS Cloud-Cleared Radiances in the GEOS

This work examines the assimilation of AIRS (Atmospheric Infrared Sounder) radiances from two points of view: the thinning strategy and the use of cloud-cleared radiances as opposed to clear-sky. Previous published work by this team, based on a very large set of Observing System Experiments performed with a 2014 3DVAR version of the GEOS (Goddard Earth Observing System), has shown that the assimilation of adaptively thinned AIRS cloud-cleared radiances (CCRs) improves the representation of tropical cyclones (TCs) without damaging the global forecast skill. The simple adaptive methodology is based on denser AIRS coverage in moving domains centered on TCs, and sparser coverage everywhere else. Subsequent experiments showed that the adaptive methodology produces good results also when applied to clear-sky CrIS (Cross-track Infrared Sounder) and IASI (Infrared Atmospheric Sounding Interferometer) radiances. In addition, the results indicate that the density of all hyperspectral data assimilated over meteorologically inactive areas is excessive, probably because of horizontal error correlation, suggesting that the global thinning should be more aggressive. More recent work focused on the polar regions has shown another positive impact of assimilating cloud-cleared AIRS radiances instead of clear-sky. The results show that high latitude atmospheric dynamics is very sensitive to the representation of the lower tropospheric temperature structure over the Arctic region. Specifically, assimilation of CCRs over areas that are data poor and also affected by broken stratus clouds, and as such minimally observed by AIRS clear-sky radiances, changes substantially the temperature structure over the Arctic low troposphere. Ingestion of CCRs over the region propagates, through hydrostatic adjustments, to mid-tropospheric geopotential height, allowing for better prediction of mid-latitude waves. In addition, adaptively thinned CCRs also improve the representation of mesoscale convective cyclones at high latitudes. An example of an Antarctic low is provided. Finally, recent ongoing work with the hybrid 4DenVAR GEOS, investigating the 2017 boreal TC season, has confirmed the previous results: namely that aggressively thinned cloud-cleared radiances improve TC structure with no loss of global skill.

CCR↗
NASA NTRS2019

Accelerated Simulation of Air Pollution Using NVIDIA RAPIDS

Atmospheric chemistry models are a central tool to study and forecast the impact of air pollution on the environment, vegetation, and human health. However, the numerical simulation of chemical kinetics is computationally expensive due to the stiffness of the system of ordinary differential equations that describes atmospheric chemistry. Here we present an alternative approach to the computation of atmospheric chemistry based on machine learning. Our training data set is produced using the NASA Goddard Earth Observing System (GEOS) model with GEOS-Chem chemistry, run on the NASA Center for Climate Simulation (NCCS) Discover supercomputing cluster on 384 Intel Xeon Haswell cores. This model spends more than 50% of total run time on solving atmospheric chemistry. The data set contains as input features the air pollution concentrations before solving the differential equations, together with some key physical parameters such as temperature and sun intensity. As target variables we define the air pollution concentrations after solving the differential equations. Using Dask-cuDF and Dask-XGBoost on the NVIDIA RAPIDS platform on 8 Tesla V100 GPUs, we generate from this training set gradient boosted decision tree models that can reproduce the simulation of chemical kinetics. We do this on the NCCS Advanced Data Analytics Platform (ADAPT) science cloud environment. Our application takes full advantage of recent advances in Dask-XGBoost, such as multi-node and multi-GPU scaling for distributed training with large data sets. The increase in training data size enabled by this is critical to capture the full range of chemical environments encountered across the globe and all annual seasons.The boosted tree models offer good predictability and show many of the features of the full chemistry reference simulation. Further improvements can be achieved through mass balance considerations and by accounting for error correlations. We incorporate the boosted tree models into the GEOS reference model using XGBoost's C API. This enables a seamless integration of the GPU trained models into GEOS-Chem, which is written in Fortran and optimized for use in a massively parallel CPU environment. We show the benefits of this approach and discuss the potential speedup of this machine learning accelerated atmospheric chemistry model.

Keller, Christoph A.↗
NASA NTRS2020

Climate change reshapes the drivers of false spring risk across European trees

(1) Temperate forests are shaped by late spring freezes after budburst—false springs—which may shift with climate change. Research to date has generated conflicting results, potentially because few studies focus on the multiple underlying drivers of false spring risk. (2) Here, we assessed the effects of mean spring temperature, distance from the coast, elevation and the North Atlantic Oscillation (NAO) using PEP725 leafout data for six tree species across 11648 sites in Europe, to determine which were the strongest predictors of false spring risk and how these predictors shifted with climate change. (3) All predictors influenced false spring risk before recent warming, but their effects have shifted in both magnitude and direction with warming. These shifts have potentially magnified the variation in false spring risk among species with an increase in risk for early‐leafout species (i.e., Aesculus hippocastanum , Alnus glutinosa , Betula pendula ) versus a decline or no change in risk among late‐leafout species (i.e., Fagus sylvatica , Fraxinus excelsior , Quercus robur ). (4) Our results show how climate change has reshaped the drivers of false spring risk, complicating forecasts of future false springs, and potentially reshaping plant community dynamics given uneven shifts in risk across species.

false spring↗
DOE OSTI2022

A multi-scale time-series dataset with benchmark for machine learning in decarbonized energy grids

The electric grid is a key enabling infrastructure for the ambitious transition towards carbon neutrality as we grapple with climate change. With deepening penetration of renewable resources, the reliable operation of the electric grid becomes increasingly challenging. In this paper, we present PSML, a first-of-its-kind open-access multi-scale time-series dataset, to aid in the development of data-driven machine learning (ML)-based approaches towards reliable operation of future electric grids. The dataset is synthesized from a joint transmission and distribution electric grid to capture the increasingly important interactions and uncertainties of the grid dynamics, containing power, voltage and current measurements over multiple spatio-temporal scales. Using PSML, we provide state-of-the-art ML benchmarks on three challenging use cases of critical importance to achieve: (i) early detection, accurate classification and localization of dynamic disturbances; (ii) robust hierarchical forecasting of load and renewable energy; and (iii) realistic synthetic generation of physical-law-constrained measurements. We envision that this dataset will provide use-inspired ML research in safety-critical systems, while simultaneously enabling ML researchers to contribute towards decarbonization of energy sectors.

54 ENVIRONMENTAL SCIENCES↗