SEARCH · Engineering Papers

Results for “Reanalysis data”

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.

SourceResearch connectionFrom yearThrough year

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

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Subseasonal Clustering of Atmospheric Rivers Over the Western United States

Abstract The serial occurrence of atmospheric rivers (ARs) along the US West Coast can lead to prolonged and exacerbated hydrologic impacts, threatening flood‐control and water‐supply infrastructure due to soil saturation and diminished recovery time between storms. Here a statistical approach for quantifying subseasonal temporal clustering among extreme events is applied to a 41‐year (1979–2019) wintertime AR catalog across the western United States (US). Observed AR occurrence, compared against a randomly distributed AR timeseries with the same average event density, reveals temporal clustering at a greater‐than‐random rate across the western US with a distinct geographical pattern. Compared to the Pacific Northwest, significant AR clusters over the northern Coastal Range of California and Sierra Nevada are more frequent and occur over longer time periods. Clusters along the California Coastal Range typically persist for 2 weeks, are composed of 4–5 ARs per cluster, and account for over 85% of total AR occurrence. Across the northwest Coast‐Cascade Ranges, clusters account for ∼50% of total AR occurrence, typically last 8–10 days, and contain 3–4 individual AR events. Based on precipitation data from a high‐resolution dynamical downscaling of reanalysis, the fractions of total and extreme hourly precipitation attributable to AR clusters are largest along the northern California coast and in the Sierra Nevada. Interannual variability among clusters highlights their importance for determining whether a particular water year is anomalously wet or dry. The mechanisms behind this unusual clustering are unclear and require further research.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

Evaluating downscaled products with expected hydroclimatic co-variances

Abstract. There has been widespread adoption of downscaled products amongst practitioners and stakeholders to ascertain risk from climate hazards at the local scale (e.g., ∼ 5 km resolution). Such products must nevertheless be consistent with physical laws to be credible and of value to users. Here we evaluate statistically and dynamically downscaled products by examining local co-evolution of downscaled temperature and precipitation during convective and frontal precipitation events (two mechanisms testable with just temperature and precipitation). We find that two widely used statistical downscaling techniques (Localized Constructed Analogs version 2, LOCA2, and Seasonal Trends and Analysis of Residuals Empirical Statistical Downscaling Model, STAR-ESDM) generally preserve expected co-variances during convective precipitation events over the historical and future projected intervals as compared to European Centre for Medium-Range Weather Forecasts Reanalysis v5 (ERA5) and two observation-based data products (Livneh and nClimGrid-Daily). However, both techniques dampen future intensification of frontal precipitation that is otherwise robustly captured in global climate models (i.e., prior to downscaling) and with process-based dynamical downscaling across five different regional climate models. In the case of LOCA2, this leads to appreciable underestimation of future frontal precipitation event intensity. This study is one of the first to quantify a likely ramification of the stationarity assumption underlying statistical downscaling methods and identify a phenomenon where projections of future change diverge depending on data production method employed. Finally, our work proposes expected co-variances during convective and frontal precipitation as useful evaluation diagnostics that can be universally applied to a wide range of statistically downscaled products.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Assessment of Model Hub Height Wind Speed Performance Using DOE Lidar Buoy Data

As the offshore wind industry continues to gain momentum both in the U.S. and globally, reanalysis products remain essential tools for the lifecycle of a wind farm. From project siting and production estimates to determining construction and maintenance windows, reanalysis products play an important role. These regional and global long-term models also serve as inputs to higher resolution models. With such wide-ranging impacts on the offshore wind industry, the need for validation of reanalysis products is strong, particularly at typical turbine hub heights. Pacific Northwest National Laboratory (PNNL) operates two AXYS WindSentinelTM buoys for the U.S. Department of Energy (DOE) in order to collect meteorological and oceanographic data, including hub height wind speed, in areas of interest for offshore wind development. This work uses DOE buoy observations off the coasts of New Jersey and Virginia to evaluate the performance of Modern Era Retrospective Analysis for Research and Applications-2 (MERRA-2), the North American Regional Reanalysis (NARR), the European Center for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5th generation (ERA5), and the National Oceanic and Atmospheric Administration (NOAA) assimilation system Rapid Refresh (RAP). Biases and degrees of correspondence are determined for each reanalysis product in order to provide insights on the performance and uncertainty for long-term wind resource characterization. With the baseline metrics determined, this work proceeds with an investigation into the sources of large deviations between the modeled and observed hub height wind speeds during the East Coast lidar buoy deployments. Consistent sources of reanalysis model error include stably stratified flow conditions, high wind shear, flow parallel to the coastline, and tropical and winter storms. For the near shore New Jersey location, offshore winds (i.e. blowing from land to water) tend to be associated with model overestimation of observed hub height wind speed, while winds parallel to the coastline are correlated with model underestimation of wind speed. Similarly, near the Virginia buoy, during southerly winds (parallel to the coastline) models underestimate the observed wind speed due to coastal upwelling. Reanalysis model wind speeds underestimate the hub height winds with increasing stable atmospheric stratification. Reanalysis model biases are seen when the observed wind speeds approach or exceed typical turbine cut-out speeds (around 20 m s-1). Large model errors during winter and tropical storm events, high wind speed events, and ramp events are often the result of wind speed timing and magnitude offsets. Overall, ERA5 provides the most successful representation of observed offshore hub height wind speeds at the U.S. East Coast buoy locations and is therefore best suited as an input boundary condition for model case studies. ERA5 performs well for a variety of atmospheric phenomena, including storms and sea breezes, however RAP, the next most successful reanalysis model, performs best in capturing the frequency of ramp events.

17 WIND ENERGY↗

DOE OSTI2025

A unified ensemble soil moisture dataset across the continental United States

Abstract A unified ensemble soil moisture (SM) package has been developed over the Continental United States (CONUS). The data package includes 19 products from land surface models, remote sensing, reanalysis, and machine learning models. All datasets are unified to a 0.25-degree and monthly spatiotemporal resolution, providing a comprehensive view of surface SM dynamics. The statistical analysis of the datasets leverages the Koppen-Geiger Climate Classification to explore surface SM’s spatiotemporal variabilities. The extracted SM characteristics highlight distinct patterns, with the western CONUS showing larger coefficient of variation values and the eastern CONUS exhibiting higher SM values. Remote sensing datasets tend to be drier, while reanalysis products present wetter conditions. In-situ SM observations serve as the basis for wavelet power spectrum analyses to explain discrepancies in temporal scales across datasets facilitating daily SM records. This study provides a comprehensive soil moisture data package and an analysis framework that can be used for Earth system model evaluations and uncertainty quantification, quantifying drought impacts and land–atmosphere interactions and making recommendations for drought response planning.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Tropical Cyclone Wind Shear-Relative Asymmetry in Reanalyses

Abstract While tropical cyclones (TCs) are axisymmetric vortices to the first order, they often exhibit noteworthy structural asymmetries. These often result from environmental vertical wind shear, which tilts the vortex and induces a wavenumber 1 pattern in the circulation and precipitation fields. Reanalyses and climate models have improved in representing the TC structure and climatology, but their relatively coarse resolution and dependence on parameterized physics cast doubt on their ability to capture the asymmetric TC structure. We perform the most comprehensive process-oriented assessment of TC asymmetry to date in reanalyses. Specifically, we analyze the composite shear-relative TC structure in ERA5 and Climate Forecast System Reanalysis (CFSR), which vary in their resolutions, physical parameterization suites, and data assimilation techniques. These structures are compared with aircraft reconnaissance radar observations. In agreement with the observations, the strongest tangential winds are usually found left-of-shear, while inner core rainfall, ascent, vortex tilt, and low-level inflow are favored directly downshear or in the downshear-left quadrant. Outer rainband convection generally peaks in the downshear-right quadrant. Thermodynamic asymmetries are also apparent, with anomalous low-level moisture right-of-shear, midlevel warmth in the upshear-right quadrant (uptilt), and cloud properties suggestive of a realistic precipitation life cycle from growth to fallout. We also decompose rainfall contributions from the convective parameterization and large-scale cloud schemes and highlight the roles of vorticity advection, buoyancy advection, and diabatic processes in driving asymmetric vertical motions in the inner core and outer rainband regions. Our results suggest that process-level studies of TC asymmetry and TC–wind shear interaction under future warming are viable using climate models. Significance Statement Asymmetries are common in tropical cyclones (TCs), influencing their intensity, track, and hazards. Vertical wind shear often plays a leading-order role in causing these asymmetries. It is uncertain how well asymmetric structures and processes are captured in reanalyses and global climate models (GCMs) with grid spacings of 0.25° and coarser. In this study, we first evaluate TC asymmetry in reanalyses, which have the benefit of being forced by observations. This helps to assess whether the resolutions associated with GCMs sufficiently capture asymmetric structures and processes and motivates upcoming work with free-running GCMs to study how TC asymmetry may change in a warming climate.

Carstens, Jacob D.↗

DOE OSTI2024

Cloud Properties and Boundary Layer Stability Above Southern Ocean Sea Ice and Coastal Antarctica

Significant variability in climate predictions originates from the simulated cloud cover over the Southern Ocean. Historically, Southern Ocean cloud and aerosol properties have been less studied than their northern hemisphere counterparts, and cloud-sea-ice interactions over the Southern Ocean also remain largely unexamined. We used data from combined radar, lidar, radiometer, radiosonde, and ERA5 reanalysis profiles to investigate cloud property relationships to cloud temperature, sea-ice concentration, and boundary layer stability. Our findings show correlations between both cloud macrophysical properties and radiative effects and sea-ice concentration, and that the marine atmospheric boundary layer is more stable over higher sea-ice concentrations. Mixed-phase cloud frequency of occurrence was highest over the sea-ice zone at 15%, three times higher than over cold water south of the Antarctic Polar Front. For temperatures greater than –15°C, low-level, single-layer clouds were more likely to precipitate ice if they were coupled to cold-water or sea-ice surfaces than if they were decoupled from these surfaces, with the highest percentage of clouds precipitating ice observed over sea ice. These findings suggest a surface source of ice-nucleating particles at high southern latitudes that increases cloud glaciation probability. We discuss the implications of our results for future studies into the relationship between cloud properties, aerosols, sea ice, and boundary layer stability at high latitudes over the Southern Ocean.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

An Evaluation of Dynamical Downscaling Methods Used to Project Regional Climate Change

In the past decade, dynamical downscaling using “pseudo-global-warming” (PGW) techniques has been applied frequently to project regional climate change. Such techniques generate signals by adding mean global climate model (GCM)-simulated climate change signals in temperature, moisture, and circulation to lateral and surface boundary conditions derived from reanalysis. An alternative to PGW is to downscale GCM data directly. This technique should be advantageous, especially for simulation of extremes, since it incorporates the GCM's full spectrum of changing synoptic-scale dynamics in the regional solution. Here, we test this assumption, by comparing simulations in Europe and Western North America. We find that for warming and changes in temperature extremes, PGW often produces similar results to direct downscaling in both regions. For mean and extreme precipitation changes, PGW generally also performs surprisingly well in many cases. Moisture budget analysis in the Western North America domain reveals why. Large fractions of the downscaled hydroclimate changes arise from mean changes in large-scale thermodynamics and circulation, that is, increases in temperature, moisture, and winds, included in PGW by design. The one component PGW may have difficulty with is the contribution from changes in synoptic-scale variability. When this component is large, PGW performance could be degraded. Global analysis of GCM data shows there are regions where it is large or dominant. Hence, our results provide a road map to identify, through GCM analyses, the circumstances when PGW would not be expected to accurately regionalize GCM climate signals.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Arctic extreme precipitation changes from 1980 to 2022 in response to sea ice decline and enhanced atmospheric rivers

Arctic extreme precipitation (EP) broadly impacts permafrost degradation, glacier and snow cover changes, and ice sheet mass balance as well as ecosystems. However, investigation of EP spatiotemporal variations over the Arctic remains challenging, and their primary drivers are still poorly understood. Performance estimation of three state-of-the-art reanalysis products (Climate Forecast System Reanalysis (CFSR), European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA-5), and Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2)) against gauge-based precipitation observations reveals that MERRA-2 outperforms other reanalysis for Arctic EP changes. Based on MERRA-2 data, both annual EP amount and occurrence days averaged over the Arctic show statistically significant positive trends during 1980–2022 (3.37 ± 1.03 mm dec −1 and 0.42 ± 0.17 d dec −1 , respectively), with the most pronounced increase in the autumn. Spatial heterogeneity in annual and seasonal EP trends is found across the Arctic, with the largest positive annual trends of 30 mm dec −1 over the Bering Sea and the Denmark Strait. The significant EP increase is closely associated with intensified atmospheric rivers (ARs) and widespread decline in sea ice concentration (SIC). Specially, SIC and ARs are responsible for 12% and 50% of Arctic EP inter-annual variance, respectively, while ARs directly contribute 28.3% of the total annual EP amounts. These findings explain the mechanistic controls on Arctic EP, providing critical insights for projecting polar weather and climate extremes and their impacts on the Arctic environment.

atmospheric rivers↗

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 OSTI2023

Deep Learning Parameterization of Vertical Wind Velocity Variability via Constrained Adversarial Training

Atmospheric models with typical resolution in the tenths of kilometers cannot resolve the dynamics of air parcel ascent, which varies on scales ranging from tens to hundreds of meters. Small-scale wind fluctuations are thus characterized by a subgrid distribution of vertical wind velocity W with standard deviation σ W . The parameterization of σ W is fundamental to the representation of aerosol–cloud interactions, yet it is poorly constrained. Using a novel deep learning technique, this work develops a new parameterization for σ W merging data from global storm-resolving model simulations, high-frequency retrievals of W , and climate reanalysis products. The parameterization reproduces the observed statistics of σ W and leverages learned physical relations from the model simulations to guide extrapolation beyond the observed domain. Incorporating observational data during the training phase was found to be critical for its performance. The parameterization can be applied online within large-scale atmospheric models, or offline using output from weather forecasting and reanalysis products.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Evaluating the potential of short-term instrument deployment to improve distributed wind resource assessment

Distributed wind projects, which are connected at the distribution level of an electricity system or in off-grid applications to serve specific or local energy needs, often rely solely on wind resource models to establish wind speed and energy generation expectations. Historically, anemometer loan programs have provided an affordable avenue for more accurate onsite wind resource assessment, and the lowering cost of lidar systems has shown similar advantages for more recent assessments. While a full 12 months of onsite wind measurement is the standard for correcting model-based long-term wind speed estimates for utility-scale wind farms, the time and capital investment involved in gathering onsite measurements must be reconciled with the energy needs and funding opportunities that drive expedient deployment of distributed wind projects. Much literature exists to quantify the performance of correcting long-term wind speed estimates with 1 or more years of observational data, but few studies explore the impacts of correcting with months-long observational periods. This study aims to answer the question of how short you can go in terms of the observational time period needed to make impactful improvements to model-based long-term wind speed estimates. Three algorithms, multivariable linear regression, adaptive regression splines, and regression trees, are evaluated for their skill at correcting long-term wind resource estimates from the European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5) using months-long periods of observational data from 66 locations across the US. On average, correction with even 1 month of observations provides significant improvement over the baseline ERA5 wind speed estimates and produces median bias magnitudes and relative errors within 0.22 m s −1 and 4 percentage points of the median bias magnitudes and relative errors achieved using the standard 12 months of data for correction. However, in cases when the shortest observational periods (1 to 2 months) used for correction are not well correlated with the overlapping ERA5 reference, the resultant long-term wind speed errors are worse than those produced using ERA5 without correction. Summer months, which are characterized by weaker relative wind speeds and standard deviations for most of the evaluation sites, tend to produce the worst results for long-term correction using months-long observations. The three tested algorithms perform similarly for long-term wind speed bias; however, regression trees perform notably worse than multivariable linear regression and adaptive regression splines in terms of correlation when using 6 months or less of observational data for correction. Translating the analysis to wind energy, median relative errors in the capacity factor are on average within 10 % using 1 month of training. If the observation period used for correction is not well correlated with the reference data, however, misrepresentation of the observed capacity factor can be substantial. The risk associated with poor correlation between the observed and reference datasets decreases with increasing training period length. In the worst-correlation scenarios, the median capacity factor relative errors from using 1, 3, and 6 months are within 47 %, 26 %, and 16 %, respectively.

17 WIND ENERGY↗

DOE OSTI2021

A database of global storm surge reconstructions

Storm surges are among the deadliest coastal hazards and understanding how they have been affected by climate change and variability in the past is crucial to prepare for the future. However, tide gauge records are often too short to assess trends and perform robust statistical analyses. Here we use a data-driven modeling framework to simulate daily maximum surge values at 882 tide gauge locations across the globe. We use five different atmospheric reanalysis products for the storm surge reconstruction, the longest one going as far back as 1836. The data that we generate can be used, for example, for long-term trend analyses of the storm surge climate and identification of regions where changes in the intensity and/or frequency of storms surges have occurred in the past. It also provides a better basis for robust extreme value analysis, especially for tide gauges where observational records are short. The data are made available for public use through an interactive web-map as well as a public data repository.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

An Experimental and Theoretical Study of the High Temperature Reactions of the Four Butyl Radical Isomers

The high temperature gas phase chemistry of the four butyl radical isomers (n-butyl,sec-butyl, iso-butyl, andtert-butyl) was investigated in a combined experimental and theoretical study. Organic nitrites were used as convenient and clean sources of each of the butyl radical isomers. Rate coefficients for dissociation of each nitrite were obtained experimentally and are at, or close to, the high pressure limit. Low pressure experiments were performed in a diaphragmless shock tube with laser schlieren densitometry at post-shock pressures of 65, 130, and 260 Torr and post-shock temperatures of 700-1000 K. Additional experiments were conducted with iso-butyl radicals at 805 K and 8.7 bar to elucidate changes in mechanism at higher pressures. These experiments were performed in a miniature shock tube with synchrotron-based photoionization mass spectrometry. The mass spectra confirmed that scission of the O-NO bond is the primary channel by which the precursors dissociate, but they also provided evidence of a minor channel (<7.7%) through HNO loss and formation of an aldehyde. These high pressure experiments were also used to determine the disproportionation/recombination ratio for iso-butyl radicals as 0.3. Reanalysis of the lower-temperature literature and the present data yielded rate constants for the disproportionation reaction, iso-butyl + iso-butyl = iso-butene + iso-butane. A chemical kinetics model was developed for the reactions of the butyl isomers that included new paths for highly energized adducts. These adducts are formed by the addition of H, CH 3 or C 2 H 5 to the butyl radicals. Accompanying theoretical investigations show that chemically activated pathways are competitive with stabilization of the adduct by collision under the conditions of the laser schlieren experiments. These calculations also show that at 10 bar and T< 1000 K stabilization is the only important reaction, but at higher temperatures, even at 10 bar, chemically activated product channels should also be considered. Branching fractions and rate coefficients are presented for these reactions. Furthermore, this study also highlights the importance of the radical structure for determining branching ratios for disproportionation and recombination of alkyl radicals, and these were facilitated by theoretical calculations of recombination rate coefficients for the four butyl radical isomers. Overall, the results reveal previously unknown features of butyl radical chemistry under conditions that are relevant to a wide range of applications and reaction mechanisms are presented that incorporate pressure dependent rate coefficients for the key steps.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Wintertime extreme warming events in the high Arctic: characteristics, drivers, trends, and the role of atmospheric rivers

Abstract. An extreme warming event near the North Pole, with 2 m temperature rising above 0 °C, was observed in late December 2015. This specific event has been attributed to cyclones and their associated moisture intrusions. However, little is known about the characteristics and drivers of similar events in the historical record. Here, using data from European Centre for Medium-Range Weather Forecasts Reanalysis, version 5 (ERA5), we study these winter extreme warming events with 2 m temperature over a grid point above 0 °C over the high Arctic (poleward of 80° N) that occurred during 1980–2021. In ERA5, such wintertime extreme warming events can only be found over the Atlantic sector. They occur rarely over many grid points, with a total absence during some winters. Furthermore, even when occurring, they tend to be short-lived, with the majority of the events lasting for less than a day. By examining their surface energy budget, we found that these events transition with increasing latitude from a regime dominated by turbulent heat flux into the one dominated by downward longwave radiation. Positive sea level pressure anomalies which resemble blocking over northern Eurasia are identified as a key ingredient in driving these events, as they can effectively deflect the eastward propagating cyclones poleward, leading to intense moisture and heat intrusions into the high Arctic. Using an atmospheric river (AR) detection algorithm, the roles of ARs in contributing to the occurrence of these extreme warming events defined at the grid-point scale are explicitly quantified. The importance of ARs in inducing these events increases with latitude. Poleward of about 83° N, 100 % of these events occurred under AR conditions, corroborating that ARs were essential in contributing to the occurrence of these events. Over the past 4 decades, both the frequency, duration, and magnitude of these events have been increasing significantly. As the Arctic continues to warm, these events are likely to increase in both frequency, duration, and magnitude, with great implications for the local sea ice, hydrological cycle, and ecosystem.

54 ENVIRONMENTAL SCIENCES↗

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 OSTI2022

Dynamically Downscaled (WRF) 1km, Hourly Meteorological Conditions 1987-2020. East/Taylor Watersheds

This dataset contains meteorological output from the Weather Research and Forecasting (WRF) version 3.8.1. This dataset has been created to 1) investigate hydrometeorological processes impacting the East River and water-delivery to the Critical Zone, and 2) provide meteorological forcing data for distributed Earth-science modeling applications in the East River watershed. Variables have a 1 kilometer spatial resolution and hourly temporal resolution and encompass a rectangular region encompassing the East and Taylor River watersheds, Colorado, near the town of Crested Butte. WRF was forced using Climate Forecast System Reanalysis (CFSR) lateral boundary conditions. Each .zip file contains one "water year" of data (October 1 -- September 30; i.e. water year 2017 starts October 1, 2016 and ends September 30, 2017). Each zip folder contains 12 netcdf (.nc) files containing one month of hourly data each and are approximately 250mb. Model timestamps are in UTC time.The files contain the following data variables:EAST_MASK: binary mask of the watershed regionTAYLOR_MASK: binary mask of the watershed regionGLW downwelling longwave radiation (w/m2)HR_PRCP: Hourly Precipitation Rate (mm/hr). Includes all hydrometeors (solid+liquid). HFX NoahMP LSM total grid-cell modelled sensible heat flux (w/m2) [positive towards atmosphere]LH NoahMP LSM total grid-cell modelled latent heat flux (w/m2) [positive towards atmosphere; can be converted to ET]PSFC Surface Barometric Pressure (hPa)Q2 Two-meter specific humidity (kg/kg)SWDOWN Downwelling shortwave solar radiation (w/m2)SWNORM Terrain-normal downwelling shortwave radiation (w/m2)T2 Two-meter air temperature (deg K)U10 10-m U-component of wind velocity (m/s)V10 10m V-component of wind velocity (m/s)XLAT Latitude of grid-center point XLONG Longitude of grid-center point XTIME Model timestamp, **in UTC**

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Surface temperature comparison of the Arctic winter MOSAiC observations, ERA5 reanalysis, and MODIS satellite retrieval

Atmospheric model systems, such as those used for weather forecast and reanalysis production, often have significant and systematic errors in their representation of the Arctic surface energy budget and its components. The newly available observation data of the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition (2019/2020) enable a range of model analyses and validation in order to advance our understanding of potential model deficiencies. In the present study, we analyze deficiencies in the surface radiative energy budget over Arctic sea ice in the ERA5 global atmospheric reanalysis by comparing against the winter MOSAiC campaign data, as well as, a pan-Arctic level-2 MODIS ice surface temperature remote sensing product. We find that ERA5 can simulate the timing of radiatively clear periods, though it is not able to distinguish the two observed radiative Arctic winter states, radiatively clear and opaquely cloudy, in the distribution of the net surface radiative budget. The ERA5 surface temperature over Arctic sea ice has a conditional error with a positive bias in radiatively clear conditions and a negative bias in opaquely cloudy conditions. The mean surface temperature error is 4°C for radiatively clear situations at MOSAiC and up to 15°C in some parts of the Arctic. The spatial variability of the surface temperature, given by 4 observation sites at MOSAiC, is not captured by ERA5 due to its spatial resolution but represented in the level-2 satellite product. The sensitivity analysis of possible error sources, using satellite products of snow depth and sea ice thickness, shows that the positive surface temperature errors during radiatively clear events are, to a large extent, caused by insufficient sea ice thickness and snow depth representation in the reanalysis system. A positive bias characterizes regions with ice thickness greater than 1.5 m, while the negative bias for thinner ice is partly compensated by the effect of snow.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Evaluation of precipitation indices in suites of dynamically and statistically downscaled regional climate models over Florida

Abstract The present work evaluates historical precipitation and its indices defined by the Expert Team on Climate Change Detection and Indices (ETCCDI) in suites of dynamically and statistically downscaled regional climate models (RCMs) against NOAA’s Global Historical Climatology Network Daily (GHCN-Daily) dataset over Florida. The models examined here are: (1) nested RCMs involved in the North American CORDEX (NA-CORDEX) program, (2) variable resolution Community Earth System Models (VR-CESM), (3) Coupled Model Intercomparison Project phase 5 (CMIP5) models statistically downscaled using localized constructed analogs (LOCA) technique. To quantify observational uncertainty, three in situ-based (PRISM, Livneh, CPC) and three reanalysis (ERA5, MERRA2, NARR) datasets are also evaluated against the station data. The reanalyses and dynamically downscaled RCMs generally underestimate the magnitude of the monthly precipitation and the frequency of the extreme rainfall in summer. The models forced with CanESM2 miss the phase of the seasonality of extreme precipitation. All models and reanalyses severely underestimate both the mean and interannual variability of mean wet-day precipitation (SDII), consecutive dry days (CDD), and overestimate consecutive wet days (CWD). Metric analysis suggests large uncertainty across NA-CORDEX models. Both the LOCA and VR-CESM models perform better than the majority of models. Overall, RegCM4 and WRF models perform poorer than the median model performance. The performance uncertainty across models is comparable to that in the reanalyses. Specifically, NARR performs poorer than the median model performance in simulating the mean indices and MERRA2 performs worse than the majority of models in capturing the interannual variability of the indices.

54 ENVIRONMENTAL SCIENCES↗