Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reanalysis”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Intense windstorms in the northeastern United States

Abstract. Windstorms are a major natural hazard in many countries. The objective of this study is to identify and characterize intense windstorms during the last 4 decades in the US Northeast and determine both the sources of cyclones responsible for these events and the manner in which those cyclones differ from the cyclone climatology. The windstorm detection is based on the spatial extent of locally extreme wind speeds at 100 m height from the ERA5 reanalysis database. During the top 10 windstorms, wind speeds exceed their local 99.9th percentile over at least one-third of land-based ERA5 grid cells in this high-population-density region of the USA. Maximum sustained wind speeds at 100 m during these windstorms range from 26 to over 43 ms−1, with wind speed return periods exceeding 6.5 to 106 years (considering the top 5 % of grid cells during each storm). Property damage associated with these storms, with inflation adjusted to January 2020, ranges from USD 24 million to over USD 29 billion. Two of these windstorms are linked to decaying tropical cyclones, three are Alberta clippers, and the remaining storms are Colorado lows. Two of the 10 re-intensified off the east coast, leading to development of nor'easters. These windstorms followed frequently observed cyclone tracks but exhibit maximum intensities as measured using 700 hPa relative vorticity and mean sea level pressure that is 5–10 times the mean values for cyclones that followed similar tracks over this 40-year period. The time evolution of wind speeds and concurrent precipitation for those windstorms that occurred after the year 2000 exhibit good agreement with in situ ground-based and remote sensing observations, plus storm damage reports, indicating that the ERA5 reanalysis data have a high degree of fidelity for large, damaging windstorms such as these. A larger pool of the top 50 largest windstorms exhibit evidence of only weak serial clustering, which is in contrast to the relatively strong serial clustering of windstorms in Europe.

58 GEOSCIENCES↗
DOE OSTI2026

Evaluation of a high-resolution regional climate simulation for surface and hub-height wind climatology over North America

Assessing the availability of key wind resources requires augmenting observations to support the implementation of wind energy infrastructure. However, observations are limited, necessitating the development of high-resolution, long-term gridded datasets. This study presents a robust, dynamically downscaled climatological dataset, offering 20 years of hourly wind data at a 4 km spatial resolution across North America, and evaluates its performance against observations, including meteorological towers and automated surface-observing system (ASOS) stations, as well as coarse-resolution reanalysis data (the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis version 5 (ERA5)). Results demonstrate that the downscaled high-resolution wind data outperform ERA5 in regions of complex terrain and coastal areas, with improved overlap coefficients for wind data distributions and reduced root mean square errors (RMSEs) for hub-height and near-surface diurnal wind patterns. The downscaled simulation also captures the synoptic drivers of seasonal wind direction patterns reasonably well, indicated by high wind rose similarity indices. This study also provides an analysis of interannual variability, utilizing the dataset's full 20-year period, and model uncertainty, generated by varying model initial conditions and physics parameterizations across 1-year ensemble members, which are key considerations for wind resource assessment in wind farm development.

17 WIND ENERGY↗
DOE OSTI2022

Evaluation of obstacle modelling approaches for resource assessment and small wind turbine siting: case study in the northern Netherlands

Abstract. Growth in adoption of distributed wind turbines for energy generation is significantly impacted by challenges associated with siting and accurate estimation of the wind resource. Small turbines, at hub heights of 40 m or less, are greatly impacted by terrestrial obstacles such as built structures and vegetation that can cause complex wake effects. While some progress in high-fidelity complex fluid dynamics (CFD) models has increased the potential accuracy for modelling the impacts of obstacles on turbulent wind flow, these models are too computationally expensive for practical siting and resource assessment applications. To understand the efficacy of available models in situ, this study evaluates classic and commonly used methods alongside new state-of-the-art lower-order models derived from CFD simulations and machine learning approaches. This evaluation is conducted using a subset of an extensive original dataset of measurements from more than 300 operational wind turbines in the northern Netherlands. The results show that data-driven methods (e.g. machine learning and statistical modelling) are most effective at predicting production at real sites with an average error in annual energy production of 2.5 %. When sufficient data may not be available de novo to support these data-driven approaches, models derived from high-fidelity simulations show promise and reliably outperform classic methods. On average these models have 6.3 %–11.5 % error compared with 26 % for classic methods and 27 % baseline error for reanalysis data without obstacle correction. While more performant on average, these methods are also sensitive to the quality of obstacle descriptions and reanalysis inputs.

17 WIND ENERGY↗
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↗
DOE OSTI2020

How Wrong Can the Operational AEP Uncertainty Estimate Be When We Ignore the Correlations Between the Uncertainty Components?

Calculations of wind farm annual energy production (AEP) on operational data are essential for a variety of financial transactions during the life of wind plants. The AEP estimate is associated with an uncertainty value, which is calculated by combining contributions connected to on-site measurements, long-term reference measurements, losses, regression, windiness adjustment, and wind resource interannual variability. Although very limited documentation on the topic exists, the conventional approach currently used by the wind energy community to assess the uncertainty connected to the operational AEP estimate assume that the different uncertainty components are uncorrelated and therefore calculates the overall uncertainty with a sum of squares approach. In this analysis, we contrast the traditional technique to estimate the overall AEP uncertainty by ignoring the correlation between its different components with a novel Monte Carlo based approach, which can instead take into account the correlation between different uncertainty categories. We consider monthly operational data from 472 wind farms from the Energy Information Administration (EIA) 923 database to assess the difference between the two approaches. Long-term wind speed data needed for the AEP assessment are taken from three reanalysis products: the Modern-Era Retrospective analysis for Research and Applications v2 (MERRA-2), the European Reanalysis Interim (ERA-interim), and the National Centers for Environmental Prediction v2 (NCEP-2). The results of the Monte Carlo approach show that three pairs of AEP uncertainty components do show a statistically significant correlation: the uncertainty connected with wind resource inter-annual variability is positively correlated with the one related to the windiness correction and negatively correlated with the one due to the regression, and the wind measurement uncertainty is positively correlated with the regression uncertainty. All these correlations, which are found between uncertainty components which are not only part of an operational analysis, but also of a wind resource assessment, are currently ignored in the conventional techniques used as industry standard. We further investigate the causes of these correlations, in terms of common dependencies of different uncertainty components on wind resource variability, number of data points, and quality of the regression between wind speed and energy production data. Next, we quantify the error in the current industry standard technique, in terms of the percentage difference in total uncertainty calculated with the two considered approaches, for all the analyzed wind farms. We find a mean absolute percentage difference of about 6%, with the largest differences being greater than 20%. The data clearly confirm that ignoring the actual correlation between the uncertainty components can lead to large errors in the assessment of the operational AEP uncertainty, and the proposed Monte Carlo approach should be preferred.

Monte Carlo↗
DOE OSTI2020

How many offshore wind turbines does New England need?

The proliferation of countries and regions with 100% clean or renewable energy targets necessitates an analysis to determine the number of generating units and storage needed to meet real-time electricity demand on the electric grid. The coastal areas of New England have the capacity to produce a large percentage of the region's energy needs with offshore wind turbines. Here we model offshore wind turbine power production data using MERRA-2 reanalysis and lidar wind speed data sets. We compare this power production to the New England hourly grid demand over the course of one year. 2,000 10 MW offshore wind turbines could satisfy New England's grid demand for about 37% of the year. When combined with 55 GWh of storage, 2,000 turbines could satisfy grid demand for about 72% of the year.

17 WIND ENERGY↗
DOE OSTI2022

A regional (land–ocean) comparison of the seasonal to decadal variability of the Northern Hemisphere jet stream 1871–2011

Seasonal to decadal variations in Northern Hemisphere jet stream latitude and speed over land (Eurasia, North America) and oceanic (North Atlantic, North Pacific) regions are presented for the period 1871–2011 from the Twentieth Century Reanalysis dataset. Significant regional differences are seen on seasonal to decadal timescales. Seasonally, the jet latitude range is lower over the oceans compared to land, reduced from 20° over Eurasia to 10° over the North Atlantic where the ocean meridional heat transport is greatest. The mean jet latitude range is at a minimum in winter (DJF), particularly along the western boundary of the North Pacific and North Atlantic, where the land-sea contrast and SST gradients are strongest. The 141-year trends in jet latitude and speed show differences on a regional basis. The North Atlantic has significant increasing jet latitude trends in all seasons, up to 3° in winter. Eurasia has significant increasing trends in winter and summer, however, no increase is seen across the North Pacific or North America. Jet speed shows significant increases evident in winter (up to 4.7 ms –1 ), spring and autumn over the North Atlantic, Eurasia and North America however, over the North Pacific no increase is observed. Long term trends are generally overlaid by multidecadal variability, particularly evident in the North Pacific, where 20-year variability in jet latitude and jet speed are seen, associated with the Pacific Decadal Oscillation which explains 50% of the winter variance in jet latitude since 1940. The results highlight that northern hemisphere jet variability and trends differ on a regional basis (North Atlantic, North Pacific, Eurasia and North America) on seasonal to decadal timescales, suggesting that different mechanisms are influencing the jet latitude and speed. This is important from a climate modelling perspective and for climate predictions in the near and longer term.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Disentangling error structures of precipitation datasets using decision trees

Characterizing error structures in precipitation products not only facilitates their proper applications for scientific and practical purposes but also helps improve their retrieval algorithms and processing methods. Despite the fact that multiple precipitation products have been assessed in the literature, factors that affect their error structures remain inadequately addressed. By interpreting 60 binary decision trees, this study disentangles the error characteristics of precipitation products in terms of their spatiotemporal patterns and geographical factors. Three independent precipitation products - two satellite-based and one reanalysis datasets: the Integrated Multi-satellitE Retrievals for GPM (Global Precipitation Measurement) late run (IMERG-L), Soil Moisture to Rain-Advanced SCATterometer (SM2RAIN-ASCAT), and the Modern-Era Retrospective analysis for Research and Applications, Version 2 uncorrected precipitation output (MERRA2-UC), are evaluated across the contiguous United States from 2010 to 2019. Here, the ground-based Stage IV precipitation dataset is used as the ground truth. Results indicate that the MERRA2-UC outperforms the IMERG-L and SM2RAIN-ASCAT with higher accuracy and more stable interannual patterns for the analysis period. Decision trees cross-assess three spatiotemporal factors and find that the underestimation of MERRA2-UC occurs in the east of the Rocky Mountains, and SM2RAIN-ASCAT underestimates precipitation over high latitudes, especially in winter. Additionally, the decision tree method ascribes system errors to nine different geographical characteristics, of which the distance to the coast, soil type, and DEM are the three dominant features. On the other hand, the land cover type, topography position index, and aspect are three relatively weak factors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Can ERA5 Be Used to Study Mesoscale Convective System Climatological Characteristics?

Mesoscale convective systems (MCSs) produce more than half of tropical rainfall and are central to the global hydrologic cycle. As the climate warms, environments favorable for MCSs may become more common; however, limited observational records hamper understanding of how MCSs respond to variations and changes in their environments. Here, we evaluate how well MCSs are represented in ERA5, a widely used global high‐resolution reanalysis product. Using PyFLEXTRKR, which jointly tracks top‐of‐atmosphere infrared brightness temperature and surface precipitation, we identified MCSs in ERA5 and compared them with those identified in satellite observations using the same detection algorithm. This comparison analysis spans 2007–2020 using hourly data at 0.25° horizontal resolution focusing over the tropics. ERA5 reproduces observed brightness‐temperature statistics and captures the geographic distribution and seasonal and diurnal cycles of MCS cold cloud shields. However, ERA5 precipitation exhibits an intensity bias—too much light rain and too little heavy rain—which shifts the rain‐rate distribution and reduces the frequency of MCSs relative to observations. Within MCSs, ERA5 precipitation exhibits the same pattern of bias, yielding a systematic underestimation of MCS precipitation intensity. Consistent with these biases, ERA5 underestimates the contribution of MCS to tropical rainfall by 25%–34% in key regions. Overall, ERA5 is suitable for studying MCS cold cloud‐shield climatology and evolution, but precipitation‐based MCS characteristics (including event‐level precipitation features and the geospatial distribution of MCS precipitation) should be interpreted with caution. These findings clarify which aspects of MCS behavior are robustly represented in ERA5 for climatological applications.

mesoscale convection↗
DOE OSTI2022

Summertime Marine Boundary Layer Cloud, Thermodynamic, and Drizzle Morphology over the Eastern North Atlantic: A Four-Year Study

Abstract Summertime remote sensor and in situ data from 2016 to 2019 collected at the ARM Eastern North Atlantic (ENA) Observatory are combined with aircraft measurements from the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign to quantify marine boundary layer (MBL) cloud, thermodynamic, and drizzle morphology in the region. A radar reflectivity–rainfall rate relationship ( Z – R ) is developed from aircraft data and 6-h cloud morphological regimes are identified from ENA data using a k -means algorithm driven by three independent inputs quantifying cloud thickness, drizzle intensity, and cloud field geometric complexity. Four separate MBL structural regimes representing non- or weakly drizzling single-layer stratocumulus, drizzling stratocumulus and cumulus-coupled stratocumulus, deep convection, and broken clouds embedded in northerly flow are identified. Single-layer stratocumulus is indicated when weak subtropical anticyclones are significantly west of the ENA site, and the MBL is cooler and drier than when drizzling and cumulus-coupled stratocumulus and broken clouds are observed. Drizzling and cumulus-coupled stratocumulus clouds are observed on the eastern flank of strong subtropical anticyclones in deep warm moist air masses with wind speeds exceeding 7 m s −1 and strong near-surface wind shear. Broken clouds exhibit strong wind shear near the inversion, while single-layer stratocumulus clouds have lower wind speeds and minimal shear. Net latent heat fluxes in the subcloud layer resulting from a combination of the ocean surface heat flux and evaporating drizzle average near zero over long periods in drizzling and cumulus-coupled stratocumulus. The ECMWF reanalysis version 5 (ERA5) is found to accurately represent single-layer stratocumulus properties, while producing significant discrepancies when drizzling stratocumulus and cumulus-coupled stratocumulus are observed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Do Machine Learning Approaches Offer Skill Improvement for Short-Term Forecasting of Wind Gust Occurrence and Magnitude?

Abstract Wind gusts, and in particular intense gusts, are societally relevant but extremely challenging to forecast. This study systematically assesses the skill enhancement that can be achieved using artificial neural networks (ANNs) for forecasting of wind gust occurrence and magnitude. Geophysical predictors from the ERA5 reanalysis are used in conjunction with an autoregressive term in regression and ANN models with different predictors, and varying model complexity. Models are derived and assessed for the warm (April–September) and cold (October–March) seasons for three high passenger volume airports in the United States. Model uncertainty is assessed by deriving models for 1000 different randomly selected training (70%) and testing (30%) subsets. Gust prediction fidelity in independent test samples is critically dependent on inclusion of an autoregressive term. Gust occurrence probabilities derived using five-layer ANNs exhibit consistently higher fidelity than those from regression models and shallower ANNs. Inclusion of the autoregressive term and increasing the number of hidden layers in ANNs from 1 to 5 also improve the model performance for gust magnitudes (lower RMSE, increased correlation, and model standard deviations that more closely approximate observed values). Deeper ANNs (e.g., 20 hidden layers) exhibit higher skill in forecasting strong (17–25.7 m s −1 ) and damaging (≥25.7 m s −1 ) wind gusts. However, such deep networks exhibit evidence of overfitting and still substantially underestimate (by 50%) the frequency of strong and damaging wind gusts at the three airports considered herein. Significance Statement Improved short-term forecasting of wind gusts will enhance aviation safety and logistics and may offer other societal benefits. Here we present a rigorous investigation of the relative skill of models of wind gust occurrence and magnitude that employ different statistical methods. It is shown that artificial neural networks (ANNs) offer considerable skill enhancement over regression methods, particularly for strong and damaging wind gusts. For wind gust magnitudes in particular, application of deeper learning networks (e.g., five or more hidden layers) offers tangible improvements in forecast accuracy. However, deeper networks are vulnerable to overfitting and exhibit substantial variability with the specific training and testing data subset used. Also, even deep ANNs reproduce only half of strong and damaging wind gusts. These results indicate the need for future work to elucidate the dynamical mechanisms of intense wind gusts and advance solutions to their prediction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Multivariate Space‐Time Dynamic Model for Characterizing the Atmospheric Impacts Following the Mt. Pinatubo Eruption

The June 1991 Mt. Pinatubo eruption resulted in a massive increase of sulfate aerosols in the atmosphere, absorbing radiation and leading to global changes in surface and stratospheric temperatures. A volcanic eruption of this magnitude serves as a natural analog for stratospheric aerosol injection, a proposed solar radiation modification method to combat a warming climate. The impacts of such an event are multifaceted and region-specific. Our goal is to characterize the multivariate and dynamic nature of the atmospheric impacts following the Mt. Pinatubo eruption. We developed a multivariate space-time dynamic linear model to understand the full extent of the spatially- and temporally-varying impacts. Specifically, spatial variation is modeled using a flexible set of basis functions for which the basis coefficients are allowed to vary in time through a vector autoregressive (VAR) structure. This novel model is cast in a Dynamic Linear Model (DLM) framework and estimated via a customized MCMC approach. We demonstrate how the model quantifies the relationships between key atmospheric parameters prior to and following the Mt. Pinatubo eruption with reanalysis data from MERRA-2 and highlight when such a model is advantageous over univariate models.

Dynamic Linear Model↗
DOE OSTI2023

A convection-permitting dynamically downscaled dataset over the Midwestern United States

Climate change is expected to have far-reaching effects at both the global and regional scale, but local effects are difficult to determine from coarse-resolution climate studies. Dynamical downscaling can provide insight into future climate projections on local scales. Here, we present a new dynamically downscaled dataset for Indiana and the surrounding regions. Output from the Community Earth System Model (CESM) version 1 is downscaled using the Weather Research and Forecasting model (WRF). Simulations are run with a 24-hr reinitialization strategy and a 12-hr spin-up window. WRF output is bias corrected to the National Centers for Environmental Protection/National Center for Atmospheric Research 40-year Reanalysis project (NCEP) using a modified quantile mapping method. Bias-corrected 2-m air temperature and accumulated precipitation are the initial focus, with additional variables planned for future releases. Regional climate change signals agree well with larger global studies, and local fine-scaled features are visible in the resulting dataset, such as urban heat islands, frontal passages, and orographic temperature gradients. This high-resolution climate dataset could be used for down-stream applications focused on impacts across the domain, such as urban planning, energy usage, water resources, agriculture and public health.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Projecting climate change in South America using variable‐resolution Community Earth System Model: An application to Chile

Abstract We introduce variable‐resolution enabled Community Earth System Model (VR‐CESM) results simulating historical and future climate conditions at 28 km over South America and 14 km over the Andes. Three 30‐year simulations are performed: a historic (1985–2014), a near future (2030–2059), and an end‐century (2070–2099) simulation under the RCP8.5 scenario. Historic results compare favourably to several temperature and precipitation reanalysis products, though local biases are present, particularly during austral summer. Future simulations highlight broad warming patterns (+3–6°C by end‐century) and heterogeneous precipitation responses across South America that qualitatively agree with prior modelling efforts. Our results reveal that the interaction between temperature and precipitation changes produce shifts in several Köppen–Geiger climates. Notable changes include the near‐elimination of the Andean Tundra or Alpine climates, a 15% decrease in Tropical Rainforests and a Tropical Savannah expansion of 20%. To provide a regionally focused analysis of projected climate change and to illustrate the benefits of variable resolution modelling, we analyse changes in the magnitude and trend in seasonal and daily temperature and precipitation in Chile. We also examined several metrics [e.g., snow water equivalent (SWE), temperatures on wet days, and days below 0°C] to evaluate potential impacts of climate change on the Chilean cryosphere between the end‐of‐century and historic periods, finding wide‐ranging indications of cryospheric decline. These changes are interpreted through reductions in the timing (1–2.5 months earlier peak SWE) and magnitude (200–1,000 mm SWE decreases) of water stored as snow in the Andes, a 10–30% decrease in number of cool season wet days with temperatures below 1°C, and 50–200 fewer days (annually) with minimum temperatures below 0°C. Our aim in producing a high‐resolution dataset of climate projections from VR‐CESM is to support analyses of climate change throughout South America but especially in vulnerable montane regions and to provide additional results for comparison with previous, ongoing, and upcoming modelling efforts.

variable-resolution global climate models↗
DOE OSTI2022

Deficient precipitation sensitivity to Sahel land surface forcings among CMIP5 models

Abstract The overall performance of the simulated seasonal precipitation response to local terrestrial forcings, namely vegetation abundance and soil moisture, in the Sahel among the Coupled Model Intercomparison Project Phase Five (CMIP5) Earth System Models (ESMs) is systematically investigated and compared with its observational counterpart using a multivariate statistical method. The observed seasonal precipitation response is evaluated against a large ensemble of observational, reanalysis, and satellite data sets to provide quantification of uncertainties. The behaviour of models with and without a Dynamic Global Vegetation Model (DGVM) component is also explored, along with the mechanisms responsible for terrestrial feedback on rainfall. In general, the CMIP5 models can reasonably capture the seasonal evolution of Sahel precipitation and soil moisture, albeit with wet biases during the pre‐monsoon period and dry biases during the peak monsoon period. The non‐DGVM ESMs simulate comparable leaf area indices (LAIs) with observations, while DGVM‐enabled ESMs simulate too much year‐round LAI. The variance of precipitation that is attributed to oceanic forcings in CMIP5 is comparable with observations; however, the variance of precipitation that is attributed to terrestrial forcings is smaller in CMIP5 models than observed, especially for non‐DGVM ESMs. CMIP5 models, especially those without DGVMs, undervalue precipitation's observed response strength to soil moisture anomalies. In both observations and CMIP5 models, none of the atmospheric variables show significant responses to direct vegetation forcing, except for the response in transpiration. Although vegetation has minimal direct effect on the atmospheric state, it can affect the atmosphere by modifying soil moisture and transpiration rate indirectly, which helps explain the more realistic simulation of rainfall in DGVM‐enabled ESMs than non‐DGVM ESMs. Coupling of an ESM to a DGVM is critical in generating reasonable land–atmosphere feedback and examining future ecological and climatic changes over the Sahel.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A new climatology of South American extratropical cyclogenesis with an intercomparison among ERA5 , JRA55 and the Brazilian Navy

Abstract In South America, the most destructive extratropical cyclones occur over the southeastern quadrant of the continent, a region that includes some of the world's largest population centres. However, there are few studies of cyclogenesis in this region, and little is known about how the varying origins of these storms impact subsequent behaviour. By supplementing an observational record with current reanalysis data, this study reveals the general characteristics of cyclones that develop in the lee regions of the Andes and those that develop nearer the Atlantic coast. The comparative climatologies demonstrate that cyclone development in these separate contexts is driven by topography, local atmospheric circulations and basin‐scale climate oscillations. The analysis further reveals a low‐frequency seesaw of lee and coastal cyclogenesis, with cross‐basin ocean temperature differences playing a teleconnectional role.

Meteorology & Atmospheric Sciences↗
DOE OSTI2020

Evaluation of the Quasi‐Biennial Oscillation in global climate models for the SPARC QBO‐initiative

Abstract Quasi‐biennial oscillations (QBOs) in thirteen atmospheric general circulation models forced with both observed and annually repeating sea surface temperatures (SSTs) are evaluated. In most models the QBO period is close to, but shorter than, the observed period of 28 months. Amplitudes are within ±20 % of the observed QBO amplitude at 10 hPa, but typically about half of that observed at lower altitudes (50 and 70 hPa). For almost all models, the oscillation's amplitude profile shows an overall upward shift compared to reanalysis and its meridional extent is too narrow. Asymmetry in the duration of eastward and westward phases is reasonably well captured, though not all models replicate the observed slowing of the descending westward shear. Westward phases are generally too weak, and most models have an eastward time mean wind bias throughout the depth of the QBO. The intercycle period variability is realistic and in some models is enhanced in the experiment with observed SSTs compared to the experiment with repeated annual cycle SSTs. Mean periods are also sensitive to this difference between SSTs, but only when parametrized non‐orographic gravity wave (NOGW) sources are coupled to tropospheric parameters and not prescribed with a fixed value. Overall, however, modelled QBOs are very similar whether or not the prescribed SSTs vary interannually. A portrait of the overall ensemble performance is provided by a normalized grading of QBO metrics. To simulate a QBO, all but one model used parametrized NOGWs, which provided the majority of the total wave forcing at altitudes above 70 hPa in most models. Hence the representation of NOGWs either explicitly or through parametrization is still a major uncertainty underlying QBO simulation in these present‐day experiments.

Bushell, A. C.↗
DOE OSTI2024

Characterizing climate pathways using feature importance on echo state networks

The 2022 National Defense Strategy of the United States listed climate change as a serious threat to national security. Climate intervention methods, such as stratospheric aerosol injection, have been proposed as mitigation strategies, but the downstream effects of such actions on a complex climate system are not well understood. The development of algorithmic techniques for quantifying relationships between source and impact variables related to a climate event (i.e., a climate pathway) would help inform policy decisions. Data-driven deep learning models have become powerful tools for modeling highly nonlinear relationships and may provide a route to characterize climate variable relationships. In this paper, we explore the use of an echo state network (ESN) for characterizing climate pathways. ESNs are a computationally efficient neural network variation designed for temporal data, and recent work proposes ESNs as a useful tool for forecasting spatiotemporal climate data. However, ESNs are noninterpretable black-box models along with other neural networks. The lack of model transparency poses a hurdle for understanding variable relationships. We address this issue by developing feature importance methods for ESNs in the context of spatiotemporal data to quantify variable relationships captured by the model. We conduct a simulation study to assess and compare the feature importance techniques, and we demonstrate the approach on reanalysis climate data. In the climate application, we consider a time period that includes the 1991 volcanic eruption of Mount Pinatubo. This event was a significant stratospheric aerosol injection, which acts as a proxy for an anthropogenic stratospheric aerosol injection. Furthermore, we are able to use the proposed approach to characterize relationships between pathway variables associated with this event that agree with relationships previously identified by climate scientists.

black-box models↗