Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hindcast”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Evaluation of the Impacts of Biomass Burning Aerosols on the NASA GEOS Sub-Seasonal Climate Forecasts

Emission sources of trace gases and aerosol particles in the South American (SA)and African (Af) continents have a strong seasonal and space variability associated with the extensive vegetation fires activities. In both continents, during the austral winter, the fires affect mainly tropical forest and savannah-type biomes and are mostly associated with deforestation and agricultural/pasture land management. Smoke aerosol particles, on average, contribute to at least 90% of the total aerosol optical depth (AOD) in the visible spectrum in the case of the South America regional smoke. Smoke aerosols also act as cloud condensation nuclei affecting cloud microphysics properties and therefore, changing the radiation budget, hydrological cycle and global circulation patterns over disturbed areas (Kaufman, 1995; Rosenfeld, 1999; Andreae,et al., 2004; Koren et al., 2004, Zhang, 2008; Ott et al., 2010; Randles et al., 2013). This study aims to evaluate and quantify the impact of including a comprehensive emission field of biomass burning aerosol on the performance of a seasonal climate forecast system, not only regarding the AOD itself but mainly on the meteorological state variable (e.g., precipitation and temperature). To address the questions put above, we designed two numerical experiments: 1- named"AERO_CTL" which applies the Quick Fire Emissions Dataset (QFED) emissions estimated with intra-diurnal variation (hereafter, BBE), and 2- named "AERO_CLM" where the sourcee mission is based on a climatology of the QFED emissions, with only monthly variation(hereafter, BBCLIM). Hindcast simulations were produced using the Goddard Earth ObservingSystem global circulation model, version 5, sub-seasonal to seasonal (GEOS5-S2S) system with a nominal spatial resolution of 56km (Rienecker et al., 2008). In both experiments, the aerosol feedbacks from cloud developments and radiation interactions were accounted. The two experiments consisted of 4 members each and ran from June to November spanning over the years 2000 to 2015. Model performance was evaluated by calculating statistical metrics on the mean area of SA and Af. Our results demonstrated that the skill model in predicting AOD is significantly improve when BBE source emission is applied over SA, but not over the Afcontinent. Over SA, the correlation between the AERO_CTL model configuration and MERRA-2 is 0.93 (R2= 0.86, RMS=0.02, BIAS=0.01), while the AERO_CLM model presents a value of0.81 (R2= 0.65, RMS=0.04, BIAS=0.06). However, the AERO_CTL experiment better represents the inter-annual variability of the AOS in both regions. The gain of the skill in predicting the AOD by the AERO_CTL experiment is also seen in some meteorological variables. We observed an increase in the model skill in predicting the 2-meter temperature and precipitation of up to 0.3 for the AERO_CTL experiment in comparison to the AERO_CLM. AERO_CLM. According to the analyzed hindcast, we inferred that representing the BBE more realistically implies in a significant gain of skills in the seasonal climate forecasting over SA and Af continents.

Aerosol↗
NASA NTRSDate not supplied

Investigating seasonal ENSO forecast amplitude calibration for GEOS-S2S-2

The GEOS-S2S-2 is a global coupled model and assimilation system, encompassing many aspects of the Earth climate system. This project looked at seasonal (nine-month) ensemble forecasts, over the period from 1982 through the present. We explored and validated a forecast amplitude correction technique used by the North American Multi-Model Ensemble (NMME), of which GEOS is a member, for Niño3.4 sea-surface temperature (SST) anomaly predictions. The method relies on deriving a set of correction factors based on the standard deviations of hindcast and observed SST anomalies. This algorithm was implemented at the Global Modeling and Assimilation Office (GMAO) and applied to the GEOS-S2S-2 Niño 3.4 hindcasts. In general, the correction reduced the forecast amplitude. This improved the quality of the forecasts when measured by the root mean square error (RMSE) of the ensemble mean forecast. RMSE decreased for most initialization months and forecast leads. The correction had a larger impact on seasons with strong El Niño-Southern Oscillation (ENSO) events than neutral seasons. The effects of the correction on the ensemble characteristics were also examined, and a variation of the algorithm was proposed that preserves the ensemble spread. The outcome of this project is a tool that can be used to calibrate GEOS-S2S-2 ENSO forecasts and potentially improve the precision of multi-model El Niño outlooks.

ENSO↗
DOE OSTI2023

The June 2012 North American Derecho: A Testbed for Evaluating Regional and Global Climate Modeling Systems at Cloud‐Resolving Scales

Abstract In this paper, we introduce a testbed for evaluating and comparing climate modeling systems at cloud resolving scales using hindcasts of the June 2012 North American derecho. To demonstrate its utility for model intercomparison, the testbed is applied to two models: the regionally‐refined Simple Cloud‐Resolving E3SM Atmosphere Model (SCREAM) at 6.5, 3.25 and 1.625 km grid spacing and the Weather Research and Forecasting (WRF) model with 3.2 and 1.6 km grid spacing. We find the simulation results to be highly sensitive to the initial conditions (ICs), initialization time, and model configurations, with ICs from the Rapid Refresh producing the best simulation. Significant improvement is identified in both models as horizontal grid spacing is refined. While a propagation delay of approximately 2 hr is found in both models, SCREAM at 1.625 km simulates the observed bow echo structure of the derecho well and predicts strong surface gusts that exceed 30 m/s. In comparison, WRF has difficulty producing surface wind over 25 m/s, with wind gusts in WRF 42%–46% lower than in SCREAM. However, WRF has a lower bias in simulating cloud top temperature and extent, but overestimates precipitation intensity. Both models reproduce the observed outgoing longwave radiation spatial patterns well (Pearson correlation >0.88), but, compared with NEXRAD observations, simulate generally larger areas of composite radar reflectivity >40 dBZ and underestimate the precipitating area by ∼47%.

2012 North American derecho↗
DOE OSTI2023

Diurnal cycle of precipitation over the tropics and central United States: intercomparison of general circulation models

Diurnal precipitation is a fundamental mode of variability that climate models have difficulty in accurately simulating. Here, in this work, the diurnal cycle of precipitation (DCP) in participating climate models from the Global Energy and Water Exchanges' DCP project is evaluated over the tropics and central United States. Common model biases such as excessive precipitation over the tropics, too frequent light-to-moderate rain, and the failure to capture propagating convection in the central United States still exist. Over the central United States, the issues of too weak rainfall intensity in climate runs is well improved in their hindcast runs with initial conditions from numerical weather prediction analyses. But the improvement is minimal over the central Amazon. Incorporating the role of the large-scale environment in convective triggering processes helps resolve the phase-locking issue in many models where precipitation often incorrectly peaks near noon due to maximum insolation over land. Allowing air parcels to be lifted above the boundary layer improves the simulation of nocturnal precipitation which is often associated with the propagation of mesoscale systems. Including convective memory in cumulus parameterizations acts to suppress light-to-moderate rain and promote intense rainfall; however, it also weakens the diurnal variability. Simply increasing model resolution (with cumulus parameterizations still used) cannot fully resolve the biases of low-resolution climate models in DCP. The hierarchy modeling framework from this study is useful for identifying the missing physics in models and testing new development of model convective processes over different convective regimes.

58 GEOSCIENCES↗
DOE OSTI2021

Anthropogenic influences on the African easterly jet–African easterly wave system

The African easterly jet (AEJ) and African easterly waves (AEWs) can have both local and far-reaching impacts on weather. It is therefore crucial to understand how the AEJ and AEWs will respond to future climate change. In this study, we examine anthropogenic influences on the AEJ–AEW system using the Weather Research and Forecasting (WRF) model configured as a tropical channel model (TCM). Hindcast simulations for the years 2001–2010 were performed using the WRF TCM, and ten additional years of simulations were performed using the pseudo-global warming method with the initial and boundary conditions of the model modified as if it were the late twenty-first century. A comparison of the simulations from the two climate scenarios indicates robust changes to both the AEJ and AEWs. For the AEJ, the jet is weaker and shifted northwards and upwards in the future climate, in association with an increase in precipitation over the Sahel and a strengthening of the meridional temperature gradient. For the AEWs, there is an increase in the number and strength of the waves in the future climate, in association with an increase in the baroclinic and barotropic energy conversions. Here, the barotropic energy conversion in particular has a larger contribution in the future climate, which manifests in the southern AEW track experiencing greater future strengthening than the northern track.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Shallow-Water Habitat in the Lower Columbia River Estuary: A Highly Altered System

Abstract Decreases in shallow-water habitat area (SWHA) in the Lower Columbia River and Estuary (LCRE) have adversely affected salmonid populations. We investigate the causes by hindcasting SWHA from 1928 to 2004, system-wide, based on daily higher high water (HHW) and system hypsometry. Physics-based regression models are used to represent HHW along the system as a function of river inflow, tides, and coastal processes, and hypsometry is used to estimate the associated SWHA. Scenario modeling is employed to attribute SWHA losses to levees, flow regulation, diversion, navigational development, and climate-induced hydrologic change, for subsidence scenarios of up to 2 m, and for 0.5 m fill. For zero subsidence, the system-wide annual-average loss of SWHA is 55 ± 5%, or 51 × 10 5 ha/year; levees have caused the largest decrease ( $${54}_{-14}^{+5}$$ 54 - 14 + 5 %, or ~ 50 × 10 5 ha/year). The loss in SWHA due to operation of the hydropower system is small, but spatially and seasonally variable. During the spring freshet critical to juvenile salmonids, the total SWHA loss was $${63}_{-3}^{+2}$$ 63 - 3 + 2 %, with the hydropower system causing losses of 5–16% (depending on subsidence). Climate change and navigation have caused SWHA losses of $${5}_{-5}^{+16}$$ 5 - 5 + 16 % and $${4}_{-6}^{+14}$$ 4 - 6 + 14 %, respectively, but with high spatial variability; irrigation impacts have been small. Uncertain subsidence causes most of the uncertainty in estimates; the sum of the individual factors exceeds the total loss, because factors interact. Any factor that reduces mean or peak flows (reservoirs, diversion, and climate change) or alters tides and along-channel slope (navigation) becomes more impactful as assumed historical elevations are increased to account for subsidence, while levees matter less.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Unraveling the 2021 Central Tennessee flood event using a hierarchical multi-model inundation modeling framework

Flood prediction systems need hierarchical atmospheric, hydrologic, and hydraulic models to predict rainfall, runoff, streamflow, and floodplain inundation. The accuracy of such systems depends on the error propagation through the modeling chain, sensitivity to input data, and choice of models. In this study, we used multiple precipitation forcings (hindcast and forecast) to drive hydrologic and hydrodynamic models to analyze the impacts of various drivers on the estimates of flood inundation depth and extent. We implement this framework to unravel the August 2021 extreme flooding event that occurred in Central Tennessee, USA. We used two radar-based quantitative precipitation estimates (STAGE4 and MRMS) as well as quantitative precipitation forecasts (QPF) from the National Weather Service Weather Prediction Center (WPC) to drive a series of models in the hierarchical framework, including the Variable Infiltration Capacity (VIC) land surface model, the Routing Application for Parallel Computation of Discharge (RAPID) river routing model, and the AutoRoute and TRITON inundation models. An evaluation with observed high-water marks demonstrates that the framework can reasonably simulate flood inundation. Despite the complex error propagation mechanism of the modeling chain, we show that inundation estimates are most sensitive to rainfall estimates. Most notably, QPF significantly underestimates flood magnitudes and inundations leading to unanticipated severe flooding for all stakeholders involved in the event. Finally, we discuss the implications of the hydrodynamic modeling framework for real-time flood forecasting.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Machine-learning predictions of the shale wells’ performance

The ultra-low permeability nature of shale reservoirs leads to an extended linear flow and necessitates horizontal wells with multi-stage engineered fractures to efficiently extract hydrocarbons resources. These artificially-generated and naturally-occurring fractures form complex networks that create complex flow regimes which control oil production. These fractures are neither identical nor equally-spaced, which leads to a production profile with a masked onset of the boundary-dominated flow. The combination of the extended linear flow with the indeterminate onset of the boundary-dominated flow challenges the current deterministic analytic approaches to forecast the estimated ultimate recovery (EUR). In this work, we propose a novel machine-learning approach which overcomes these challenges and provides reliable EUR estimates based on field-wide analyses. We implement a novel unsupervised machine learning (ML) methodology, which allows for automatic identification of the optimal number of features (signals) present in the data based on non-negative matrix/tensor factorization coupled with k-means clustering incorporating regularization and physics constraints. In the presented analyses, the input data to the ML algorithm is the available (public) production history from the field collected at existing unconventional reservoirs. We validate our approach through hindcasting of the production data, where we achieved an excellent agreement. In addition, our approach is able to identify the poorly-performing wells, which could benefit from early refracing. Our approach provides fast and accurate estimations of the well performance without presumptions about the state of the well or the flow regime.

03 NATURAL GAS↗
DOE OSTI2020

Non-stationary historical trends in wave energy climate for coastal waters of the United States

The thirty-year non-stationary historical trends in the wave energy climate for United States coastal waters between 1980 and 2009 are investigated using spectral partitioned wave data generated from a WaveWatch III® (version 5.05) hindcast. In addition to historical trends in the omni-directional wave power, frequency and directionally resolved wave power, frequency and directional spreading, and seasonal variability, are examined for the first time, including their geographical distribution. These historical wave energy climate trends are linked to changes to the dominant wave systems and commensurate trends in the historical wind climate. Total wave power trends are consistent with other studies, but the present study identifies regions, and specific frequency and direction bands, where significant wave energy climate changes have occurred. For some regions, changes to wave climate parameters exceeded one-percent annually, more than thirty-percent over the study period. In conclusion, non-stationary trends of this magnitude have significant implications for ocean and coastal engineering projects designed assuming wave climates are stationary and warran

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Evaluating Multimodel Ensemble Seasonal Climate Forecasts on Rangeland Plant Production in the California Annual Grassland

Seasonal precipitation and temperature directly affect total plant production in the California Annual Grassland (CAG). Technological advances have resulted in skillful seasonal climate forecasts (i.e., significant correlations between actual and forecasted climate), which could be input into plant production models to inform stocking and other rangeland management decisions. This study presents a procedure for forecasting plant production in the CAG ecosystem to predict annual plant production for grazing, restoration, or other rangeland management practices using a combination of historical gridMET climate data and seasonal hindcasts (i.e., retrospective forecasts, from the North American Multi-Model Ensemble program). Further, the results of this study first confirmed high forecast skill, throughout the growing season at all sites. We also identified skillful plant production forecasts across most of the growing season at two sites and in three of the seven forecasting months at one study site. Forecasting climate and end-of-year plant production across the growing season at three CAG sites allowed us to identify the places and times in the growing season when forecasting might be most helpful in informing management decisions. Integrating plant production forecasting into rangeland management practices could significantly improve rangeland management outcomes. These procedures provide a user guide for creating plant production forecasts for any given area of interest and may be applicable across a wide range of other agricultural and rangeland management systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Characterization of Gulf Stream energy resources: A 30-year high-resolution simulation study

The strength and proximity of the Gulf Stream (GS) along the U.S. coast offers opportunities for marine energy extraction. Leveraging high-resolution (up to 400 m) long-term (30-year) hindcasts based on an unstructured grid model, this study characterizes the energy potential of the GS in two viable locations: the Florida Straits (FS) and Cape Hatteras (CH). Results reveal distinct characteristics in the kinetic energy flux (KEF) and kinetic energy density (KED) for the two regions. The 30-year mean KEF at CH is 27.50 GW, surpassing the 19.74 GW estimated for the FS, which is attributed to the wider cross-section of the GS in CH. Additionally, CH exhibits higher KEF variability due to GS’s larger path shifts and meanders. Estimated KEDs at ADCP site locations reach up to 2,908 W/m 2 in the FS and 1,512 W/m 2 in the CH at a depth of 20 m. High-energy zones are closer to shore in the FS, whereas energy-rich areas at CH are located further offshore. A depth-based KED analysis shows significant decreases across typical energy extraction depths (20 to 100 m), with reductions up to 44%. Finally, this study provides high-resolution resource characterization, forming a foundation for optimized site selection and marine hydrokinetic system design.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Evaluating the effect of meso/submesoscale current–wave interactions on wave energy resource characterization at northeast U.S. coast

Wave energy is a promising renewable resource, but accurate assessment is difficult in regions with strong currents due to wave–current interactions (WCI). Here, this study develops a two-way coupled WCI model within the Coupled Ocean Atmosphere Wave Sediment Transport (COAWST) framework at 2 km resolution to improve wave energy characterization along the northeastern U.S. coast, including the Mid-Atlantic Bight and Gulf of Maine. The model integrates WaveWatchIII (WWIII) and the Regional Ocean Modeling System (ROMS) to enhance wave hindcasting by accounting for Doppler-shift, refraction, and nonlinear energy exchanges. Validation against buoy and satellite observations confirms model accuracy. Analysis shows that Doppler-shifting can alter wave power density by over 20%, while strong current gradients and shear distort wave crests via focusing/defocusing and stretching/squeezing, modifying wave direction and frequency. These processes together can induce wave power fluctuations of up to 40% on synoptic scales. Applying a 2.5 MW Ocean Energy Converter power matrix shows that WCI may change harvested energy by up to 100% in shallow-waters and 60% in deep-waters. These results underscore the importance of incorporating nonlinear WCI for reliable wave climate predictions and resource assessments in energetic coastal regions.

doppler-shift↗
DOE OSTI2020

Toward Understanding the Simulated Phase Partitioning of Arctic Single-Layer Mixed-Phase Clouds in E3SM

21 Mixed-phase clouds are widespread in the Arctic and have substantial impacts on the 22 surface energy budget, sea ice extent, and regional climate. Significant changes are found in the 23 modeled phase partitioning of Arctic mixed-phase clouds in the newly developed U.S. 24 Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model 25 version 1 (EAMv1) compared to its predecessor (CAM5/EAMv0). The mixed-phase clouds 26 simulated by EAMv1 is overly dominated by super-cool liquid which is often underestimated by 27 CAM5. Compared to CAM5, EAMv1 has made several important changes in its cloud physical 28 parameterizations, including using CLUBB to unify shallow convection, PBL turbulence, and 29 cloud macrophysics, as well as updating the MG two-moment cloud microphysics scheme and 30 using the CNT heterogeneous ice nucleation scheme to replace the Meyers scheme for mixed- 31 phase clouds. To understand which of the newly-implemented physical parameterizations is 32 responsible for the change of modeled mixed-phase cloud phase partitioning in EAMv1, 33 sensitivity tests are conducted by utilizing the short-term hindcast approach and results are 34 compared to the observations from the ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) 35 field campaign. 36 It is found that although each scheme contributes to the decrease of cloud ice water 37 content (IWC) and increase of cloud liquid water content (LWC), only the combined effect of 38 CLUBB and CNT can explain the significant underestimation of IWC. By examining the budgets 39 of cloud microphysical processes, CLUBB is found to play a dominant role due to the fact that 40 CLUBB does not diagnoses cloud ice phase water in its assumed joint PDF method and thus may 41 not provide sufficient cloud ice to initialize ice phase microphysical processes such as the growth 42 of ice crystals through vapor deposition. Although a more effective heterogeneous ice nucleation 43 may alleviate the underestimation of IWC, the ice phase process in CLUBB requires urgent 44 attention in the future model development.

Zhang, Meng↗
DOE OSTI2022

Seasonality of Interbasin SST Contributions to Atlantic Tropical Cyclone Activity

Abstract Recent studies have demonstrated that the difference in sea surface temperature anomalies (SSTAs) between the tropical Atlantic main development region (MDR) and the tropical Pacific (Niño 3) modulates Atlantic tropical cyclone activity. This study further explores the seasonality of Pacific and Atlantic contributions to Atlantic hurricane activity. Our analysis shows that while MDR and Niño 3 SSTAs are equally important for late‐season (September–November) activity, early‐season (June‐August) activity is largely modulated by MDR SSTAs. This reflects the increased (reduced) variance of MDR (Niño 3) SSTAs in the early‐season due to their phase locking to the seasonal cycle. Further analysis yields skillful forecasts using an MDR‐Niño 3 interbasin index derived from hindcasts of the North American Multi‐Model Ensemble with May initial conditions. However, the prediction skill for MDR SSTAs is lower than that of Niño 3 SSTAs, suggesting that increasing the prediction skill for MDR SSTAs is key to improving seasonal outlooks.

West, Robert↗
DOE OSTI2022

Anthropogenic Contributions to the 2021 Pacific Northwest Heatwave

Abstract Daily maximum temperatures during the 2021 heatwave in the Pacific Northwest United States and Canada shattered century old records. Multiple causal factors, including anthropogenic climate change, contributed to these high temperatures, challenging traditional methods of attributing human influence. We demonstrate that the observed 2021 daily maximum temperatures are far above the bounds of Generalized Extreme Value distributions fitted from historical data. Hence, confidence in Granger causal inference statements about the human influence on this heatwave is low. Alternatively, we present a more conditional hindcast attribution study using two regional models. We performed ensembles of simulations of the heatwave to investigate how the event would have changed if it had occurred without anthropogenic climate change and with future warming. We found that global warming caused a ∼0.8°C–1°C increase in heatwave temperatures. Future warming would lead to a ∼5°C increase in heatwave temperature by the end of the 21st century.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Captured QBO‐MJO Connection in a Subseasonal Prediction System

Abstract The quasi‐biennial oscillation (QBO) impacts the Madden‐Julian Oscillation (MJO) activity with a stronger MJO in QBO easterly (QBOE) than QBO westerly (QBOW) winters. However, this relationship is poorly represented in the current generation climate models. For the first time, this paper applies a stratospheric zonal‐mean nudging in a subseasonal prediction system to capture it. Two strong MJO cases in a QBO‐neutral winter are investigated. The QBO temperature and zonal wind anomalies are added separately as well as together to the stratosphere using nudging in MJO case hindcast. Only by nudging the QBO temperature anomalies while leaving the zonal wind free, can the prediction system capture the observed QBO‐MJO connection. The tropopause instability is found positively correlated to the MJO amplitude, but it cannot fully explain the captured connection. The free‐evolving zonal wind anomalies in the stratosphere due to the nudged QBO temperature are crucial for the captured connection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Future Projections of EURO‐CORDEX Raw and Bias‐Corrected Daily Maximum Wind Speeds Over Scandinavia

Twenty historical and future Coordinated Regional Climate Downscaling Experiment ensemble for Europe simulations are bias corrected to investigate the future changes in the daily maximum wind speed over Scandinavia. We use quantile mapping to adjust the wind for the historical period (1985–2014) and quantile delta mapping for two future periods (2041–2070, 2071–2100, RCP8.5) with the 3-km spatial resolution NORA3 hindcast as the reference data set. Decomposing the variance, we find that most of the inter-model spread in the bias and the response to climate change is due to the regional climate model over land and mainly to the general climate model over sea. On average, the mean daily maximum wind speed is projected to increase everywhere except along the western coast of Norway and Denmark. In summer, we see an opposite sign over Sweden and Finland. The Norwegian Sea experiences weaker mean and high wind whereas the Baltic Sea experiences a strengthening. Bias correction influences the amplitude of the response, not the response pattern. Wind speed distributions can have different shapes in the future, with, for example, a flattening of the distribution off the coast of Norway with more frequent weak and very strong winds. Apart from a few locations in Norway, we find an increase in the number of days in the local highest historical wind category. Overall, summer exhibits opposite signals compared to the three other seasons. At country scale, the sign of the change in the mean and 98th percentile varies greatly depending on the region and among the simulations, especially for Norway.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Revealing the Statistics of Extreme Events Hidden in Short Weather Forecast Data

Extreme weather events have significant consequences, dominating the impact of climate on society. While high-resolution weather models can forecast many types of extreme events on synoptic timescales, long-term climatological risk assessment is an altogether different problem. A once-in-a-century event takes, on average, 100 years of simulation time to appear just once, far beyond the typical integration length of a weather forecast model. Therefore, this task is left to cheaper, but less accurate, low-resolution or statistical models. But there is untapped potential in weather model output: despite being short in duration, weather forecast ensembles are produced multiple times a week. Integrations are launched with independent perturbations, causing them to spread apart over time and broadly sample phase space. Collectively, these integrations add up to thousands of years of data. We establish methods to extract climatological information from these short weather simulations. Using ensemble hindcasts by the European Center for Medium-range Weather Forecasting archived in the subseasonal-to-seasonal (S2S) database, we characterize sudden stratospheric warming (SSW) events with multi-centennial return times. Consistent results are found between alternative methods, including basic counting strategies and Markov state modeling. By carefully combining trajectories together, we obtain estimates of SSW frequencies and their seasonal distributions that are consistent with reanalysis-derived estimates for moderately rare events, but with much tighter uncertainty bounds, and which can be extended to events of unprecedented severity that have not yet been observed historically. These methods hold potential for assessing extreme events throughout the climate system, beyond this example of stratospheric extremes.

58 GEOSCIENCES↗