Engineering Papers⌕ Search

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.

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

At least 991 records · Page 55

Comparison of Observed Longwave, Shortwave Irradiance and Surface Temperature from “MOSAiC” to CERES Radiative Transfer Calculations and Inputs

As atmospheric temperatures rise due to increased anthropogenic forcing, the effect is expected to be larger over the arctic than midlatitude and tropics, known as polar amplification. A multi-national program, the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) was ran between Sept 2019 and Oct 2020. The project deployed scientific instruments on board the German research vessel Polarstern with the ship remaining across a year to observe all aspects of polar climate. The US DOE deployed the ARM AMF2 aboard the ship along with several off-ship sites for extended spatial observations. Of particular importance was the measurement of the energetics of the ice/ocean/atmosphere interactions through observations of surface irradiance and surface temperatures. NASA’s Clouds and the Earth’s Radiant Energy System (CERES) project produces the SYN1deg Edition 4 data product. This product provides high quality, hourly, globally gridded and temporally complete maps of top of atmosphere (TOA), in atmosphere, and surface irradiances. TOA fluxes are derived from CERES instruments and geostationary satellites. In addition, TOA, in atmosphere, and surface irradiances are computed using the Langley Fu and Liou radiative transfer model. The radiative transfer model is run hourly at 1 degree ~equal area spatial resolution. Meteorological profiles are provided by Global Modeling and Assimilation Office’s GEOS-541 reanalysis product. Cloud properties are derived solely from Terra and Aqua MODIS imagers northward of 60°. Here we compare SYN1deg hourly calculations of surface irradiance to observations from the AMF2 to validate the products estimates of surface irradiance in this challenging area. Along with the irradiance comparisons we take a close look at the surface temperature record in the GEOS-541 product and compare it over both time and space to several surface observations provided by MOSAiC. Along with the comparison to the re-analysis record we will compare these observed surface temperatures to those derived from the AIRS product, which is known to have some difficulty in extreme high latitude areas. The goals of this study are 1) understand computed surface downward irradiance and temperature error separated by surface type (e.g. sea ice or open water) and by season, and 2) error covariance in spatial and temporal space. We seek to use this information to extrapolate the error from the MOSAiC domain to a larger arctic region.

David A Rutan↗

Regime-based aerosol–cloud interactions from CALIPSO-MODIS and the Energy Exascale Earth System Model version 2 (E3SMv2) over the Eastern North Atlantic

This study investigates aerosol-cloud interactions in marine boundary layer (MBL) clouds using an advanced deep-learning-driven synoptic-regime-based framework, combining satellite data (CALIPSO vertically resolved aerosol extinction and MODIS cloud properties) with 1° nudged Energy Exascale Earth System Model version 2 (E3SMv2) simulation over the Eastern North Atlantic (ENA; ∼10°×10°, 2006–2014). The E3SMv2 captures observed seasonal variations in cloud droplet number concentrations (N d ) and liquid water path (LWP), though it systematically underestimates N d . We then partition ENA meteorology into four synoptic regimes (Pre-Trough, Post-Trough, Ridge, Trough) via a deep-learning clustering of ERA5 reanalysis fields, enabling regime-dependent aerosol-cloud interactions analyses. Both satellite and E3SMv2 exhibit an inverted-V LWP-N d relationship. In Post-Trough and Ridge regimes, the satellite shows stronger negative LWP-N d sensitivities than in Pre-Trough regime. The Trough regime displays a muted satellite LWP response. In comparison, the model predicts more exaggerated LWP responses across regimes, with LWP increasing too quickly at low N d and decreasing more sharply at high N d , especially in Pre-Trough and Trough regimes. These exaggerated model LWP sensitivities may stem from uncertainties in representing drizzle processes, entrainment, and turbulent mixing. As for N d susceptibility to aerosols, N d increases with MBL aerosol extinction in both datasets, but the simulated aerosol-cloud interactions appear oversensitive to meteorological conditions. Overall, E3SMv2 better captures aerosol effects under regimes that favor stratiform clouds (Post-Trough, Ridge), but performance deteriorates for regimes with deeper, dynamically complex clouds (Trough), highlighting the need for improved representations of those cloud processes in climate models.

Environmental sciences↗

CHUWD-H v1.0: a comprehensive historical hourly weather database for U.S. urban energy system modeling

Reliable and continuous meteorological data are crucial for modeling the responses of energy systems and their components to weather and climate conditions, particularly in densely populated urban areas. However, existing long-term datasets often suffer from spatial and temporal gaps and inconsistencies, posing great challenges for detailed urban energy system modeling and cross-city comparison under realistic weather conditions. Here we introduce the Historical Comprehensive Hourly Urban Weather Database (CHUWD-H) v1.0, a 23-year (1998-2020) gap-free and quality-controlled hourly weather dataset covering 550 weather station locations across all urban areas in the contiguous United States. CHUWD-H v1.0 synthesizes hourly weather observations from stations with outputs from a physics-based solar radiation model and a reanalysis dataset through a multi-step gap filling approach. A 10-fold Monte Carlo cross-validation suggests that the accuracy of this gap filling approach surpasses that of conventional gap filling methods. Designed primarily for urban energy system modeling, CHUWD-H v1.0 should also support historical urban meteorological and climate studies, including the validation and evaluation of urban climate modeling.

54 ENVIRONMENTAL SCIENCES↗

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↗

Evaluation of Regional Water and Energy Balance in Contemporary Reanalyses

Regional to continental water and energy balance in reanalyses provide crucial information on the climate and climate variations in that region and can contrast other regions. Reanalyses are a valuable tool in piecing together the complete balance owing to the incorporation of the budgets in the background model and the assimilation of observation. On the other hand, the background model has its own biases and observational assimilation can lead to non-negligible tendencies within the budgets, and these can vary over space and time. Here, we will compare the balance within several reanalyses and how they have evolved over recent decades. Several observational data products will provide reference points in the comparison. In addition, a recent project has integrated observations of the water and energy budgets in an optimization that aims to reduce the imbalance of regional budgets in using observational terms from disparate sources. In evaluating the regional balances, we will incorporate, when available, the water and energy convergence terms as well as their analysis increments (which may be diagnosed as residuals). The reanalyses considered here include: Japanese Reanalyses for Three Quarters Century (JRA3Q), ECMWF Reanalysis version 5 (ERA5), Modern Era Retrospective-analysis for Research and Applications version 2 (MERRA-2) and the NCEP Reanalysis II (NCEPR2). NCEPR2 is representative of the late 90s technology, while the others represent the contemporary climate reanalyses presently available. We will discuss the regional balances over key river basins (e.g. Mississippi and Amazon), and a range of continents including those well observed (North America and Europe) and those with fewer consistent observations (e.g. Africa). These contrasting regions show the degree of advancement in reanalyses and expose areas in need of improvement. The results also demonstrate the usefulness of producing complete atmospheric water and energy budgets as part of the diagnostic output.

Michael G Bosilovich↗

Evaluation of Regional Water and Energy Balance in Contemporary Reanalyses

Regional to continental water and energy balance in reanalyses provide crucial information on the climate and climate variations in that region and can contrast other regions. Reanalyses are a valuable tool in piecing together the complete balance owing to the incorporation of the budgets in the background model and the assimilation of observation. On the other hand, the background model has its own biases and observational assimilation can lead to non-negligible tendencies within the budgets, and these can vary over space and time. Here, we will compare the balance within several reanalyses and how they have evolved over recent decades. Several observational data products will provide reference points in the comparison. In addition, a recent project has integrated observations of the water and energy budgets in an optimization that aims to reduce the imbalance of regional budgets in using observational terms from disparate sources. In evaluating the regional balances, we will incorporate, when available, the water and energy convergence terms as well as their analysis increments (which may be diagnosed as residuals). The reanalyses considered here include: Japanese Reanalyses for Three Quarters Century (JRA3Q), ECMWF Reanalysis version 5 (ERA5), Modern Era Retrospective-analysis for Research and Applications version 2 (MERRA-2) and the NCEP Reanalysis II (NCEPR2). NCEPR2 is representative of the late 90s technology, while the others represent the contemporary climate reanalyses presently available. We will discuss the regional balances over key river basins (e.g. Mississippi and Amazon), and a range of continents including those well observed (North America and Europe) and those with fewer consistent observations (e.g. Africa). These contrasting regions show the degree of advancement in reanalyses and expose areas in need of improvement. The results also demonstrate the usefulness of producing complete atmospheric water and energy budgets as part of the diagnostic output.

Michael Bosilovich↗

Version 2.0 AERONET Dust Aerosol properties, Constraints and Application to Asian Dust Observations

In November 2006, AERONET released Version 2 of the Dubovik and King sky radiance and optical depth inversion. Reanalysis of the entire AERONET database revealed marked differences in aerosol properties in arid and semi arid regions with dust dominated aerosols. The change will be illustrated through sensitivity analysis and examples from the UAE2 (United Arab Emirates Unified Aerosol Experiment) field campaign. Properties of dust dominated aerosols will be presented from regional AERONET sites in China showing variations in dust aerosol properties. The constraints and limitations of the AERONET inversion will be presented that will facilitate analysis by the user community of these data.

Holben, Brent↗

TDCOSMO. I. An exploration of systematic uncertainties in the inference of $H_0$ from time-delay cosmography

Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, H0. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on H0 from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in H0 values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the “bulge-halo” conspiracy, H0 is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields H0 = 74.0−1.8+1.7 km s−1 Mpc−1, while the power-law model yields 74.2−1.6+1.6 km s−1 Mpc−1. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.

79 ASTRONOMY AND ASTROPHYSICS↗

Influence of Antarctic and Greenland Continental Shelf Circulation on High‐Latitude Oceans in E3SM

The science objectives of this project are to simulate and understand the impacts of both deep-basin warm-water intrusions and land-ice melt on the continental shelf circulations and sea-ice distributions around the margins of Greenland and Antarctica. As well, the role of subsurface ocean heat from the Atlantic on declining sea-ice cover in the Arctic is explored. Mesoscale processes and fine bathymetry are implicated in cross-shelf property transports around both Greenland and Antarctica. Therefore, we configured and ran an atmospheric reanalysis-forced global ocean/sea-ice simulation on a grid that reduces from 8 km at the Equator to 2 km at the poles (UH8to2) with 60 vertical levels. It was produced using the Energy Exascale Earth System Model ‘‘HiLAT’’ code (E3SMv0-HiLAT) that uses the Parallel Ocean Program (POP) and CICE5 as its ocean and sea-ice components, respectively. Two main UH8to2 simulations were carried out: one for 1975-2009 and the other for July 2016-2020 after it was initialized from a 1/25° data-assimilative ocean/sea-ice prediction system ocean/sea-ice state. The UH8to2 is not coupled to an active land-ice model. Rather, land-ice melt is represented by observationally informed freshwater fluxes (FWFs). Short (multi-year) UH8to2 simulations were conducted to understand sensitivities when Greenland ice sheet (GrIS) melt is released only at the ocean surface or when it is distributed over the upper water column in accordance with fjord melt plume behavior; these cases were compared with a no GrIS melt case. West Greenland continental shelf currents were fastest in the vertical distribution case and an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy was found relative to the surface release case. Further, salinity is lower and meltwater volume greater in the eastern Labrador Sea in the vertical distribution case. For the Arctic, the veracity of the UH8to2 was evaluated for 2017-2020 using available observations. Simulated seasonal sea-ice thickness and concentration are realistic, but the ice is unrealistically thin in the central and eastern Arctic in the fall. Comparisons of vertical sections of ocean temperature, salinity, and buoyancy collected from Ice-Tethered Profilers (ITPs) in the eastern Arctic in the fall and winter of 2019/2020 and co-located/concurrent UH8to2 fields show the stratification over the top 100 m of the water column is too low in the model, the simulated mixed layer too deep, and the simulated subsurface Atlantic Water (AW) too warm; these biases may contribute to the sea-ice biases. A model intercomparison study using the UH8to2 and a forced 1/25° regional Arctic ocean/sea-ice (uses the HYbrid Coordinate Ocean Model and CICE5) simulation further investigates the relationship between AW and sea-ice in the eastern Arctic. The models show a mesoscale-rich pulse of Atlantic Water extending into the eastern basin that reaches maximum intensity in late winter of 2018, after which it decreases in strength. Concurrent and co-located sea-ice melt or the inhibition of sea-ice growth is seen and is attributed to halocline mesoscale eddies doming into the mixed layer with convection bringing this heat into the vicinity of the sea-ice.

58 GEOSCIENCES↗

The Influence of Sea ice on Arctic Cloud Properties: What Can We Learn by Applying a In Situ Observational Strategy to Satellite Data?

Sea ice is declining because of anthropogenic climate change. This change alters many aspects of the Arctic climate system, including the way that the surface and atmosphere interact. Atmosphere-surface coupling processes represent an important cloud feedback mechanism that can alter the Arctic surface energy budget. For the Arctic, it has been hypothesized that a reduction in sea ice cover could lead to an increase in clouds. If this process were to occur as originally hypothesized, sea ice loss in all seasons would lead to an increase in clouds. Recent observational studies find a cloud response to sea ice loss in non-summer months and no cloud response in summer months. However, previous studies rely on inter-annual variability and reanalysis to control for the influence of meteorology, reducing the confidence in the resulting conclusions. We adopt a phenomenological, event-based approach that does not need to use meteorological reanalysis. The approach analyzes cloud properties derived from CALIPSO-CloudSat over sea ice and adjacent ice-free footprints by compositing individual satellite ground tracks that cross the Arctic sea ice edge. The underlying assumption, which we verify, is that footprints that are close to each other in space and time experience similar large-scale meteorological conditions. Our results show larger cloud fraction and more cloud liquid water over ice-free than over sea ice footprints and provide additional evidence for a seasonal dependence of cloud-sea ice coupling that is in line with previous work. We find a different result where the maximum cloud property differences between sea ice and ice-free ocean occurs in spring, not fall as earlier studies suggest. We argue that these cloud differences between sea ice and ice-free ocean are primarily caused by the influence of the surface type on the thermodynamic stability of the lower troposphere and not principally from an increase in surface evaporation. In addition, we explore the sensitivity of these results to marginal ice zone width, season, Atlantic vs. Pacific sector, and the cloud property dependence on the distance from the sea ice edge. Overall, our results provide further evidence that cloud-sea ice coupling processes are not offsetting the observed surface energy budget perturbation due to sea ice loss in summer; thus, the cloud response to declining sea ice appears to contribute to amplified Arctic warming.

Patrick C Taylor↗

Image-Subtraction Photometry of Variable Stars in the Globular Clusters NGC 6388 and NGC 6441

We have applied Alard's image subtraction method (ISIS v2.1) to the observations of the globular clusters NGC 6388 and NGC 6441 previously analyzed using standard photometric techniques (DAOPHOT, ALLFRAME). In this reanalysis of observations obtained at CTIO, besides recovering the variables previously detected on the basis of our ground-based images, we have also been able to recover most of the RR Lyrae variables previously detected only in the analysis of Hubble Space Telescope WFPC2 observations of the inner region of NGC 6441. In addition, we report five possible new variables not found in the analysis of the EST observations of NGC 6441. This dramatically illustrates the capabilities of image subtraction techniques applied to ground-based data to recover variables in extremely crowded fields. We have also detected twelve new variables and six possible variables in NGC 6388 not found in our previous groundbased studies. Revised mean periods for RRab stars in NGC 6388 and NGC 6441 are 0.676 day and 0.756 day, respectively. These values are among the largest known for any galactic globular cluster. Additional probable type II Cepheids were identified in NGC 6388, confirming its status as a metal-rich globular cluster rich in Cepheids.

Corwin, Michael T.↗

Wind Forcing of the North Sea Pole Tide

The Chandler wobble of the earth's rotation has a period near 14 months and sets up the 0.5 cm amplitude pole tide in the deep oceans, However, the pole tide is anomalously large in the North Sea, where the amplitude increases sharply up to 3 cm along the continental coast. It is shown here that the sea levels are well correlated with winds at the pole tide frequency. The Princeton Ocean Model is used to investigate the response of the North Sea to wind forcing. The barotropic numerical ocean model depicts realistic coastlines and bathymetry at 5 ft x 5 ft resolution, with 97 x 73 grid points. The monthly mean wind fields for the 40-year period (1958-1997) from the National Centers for Atmospheric Prediction (NCEP) reanalysis were used to force the model. The winds were converted to stress with a neutral drag coefficient that varied linearly with windspeed (instead of using the NCEP windstress). A 5-day simulation was made for each month until the resulting flow regime came into equilibrium, and model water levels at various station locations were saved for comparison with tidal-gauge observed sea levels from the Permanent Service for Mean Sea Level (PSMSL). The comparison is made for 10 North Sea stations with high quality tide gauge data. We find: (1) good agreement in annual and semi-annual phases and in the trend of amplitude w.r.t. latitude; (2) more importantly for this study, the model-predicted and observed non-seasonal sea level variations show a very significant temporal correlation as well as spectral coherence. However, a large amplitude difference exists between the two sets -- the overall amplitude variability of the observed is generally a factor of 2-3 larger than the model prediction (this same phenomenon has been reported in ocean circulation studies, although the cause is not yet clear.) Our results indicate that the wind forcing is the main cause of the observed large pole tide in the North Sea.

OConnor, W.↗

Long-Term Ozone Variability and Trends from Reanalyses: Can It Be Done?

Stratospheric ozone concentrations have begun to show early signs of recovery following the implementation of the Montreal Protocol and its amendments as well as in response to decreasing upper-stratospheric temperatures. Secular trends in stratospheric ozone are modulated by considerable interannual variability and systematic changes in transport patterns that are expected under increasing concentrations of greenhouse gases, especially in the lower stratosphere. These factors necessitate the continued close monitoring of stratospheric ozone in upcoming decades, with a special focus on the lower stratosphere.As highly resolved data sets combining a plethora of observations with model simulations atmospheric reanalyses are, in principle, well suited for the task. All major reanalyses generate ozone output. However, significant spurious discontinuities that arise from step changes in the observing systems prevent a straightforward analysis of ozone trends and long-term variability. Building on our recent work, in this presentation we will demonstrate that trend detection is nonetheless possible using the ozone record from NASA's MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, Version 2) reanalysis bias-corrected using a chemistry model simulation as a transfer function. Next, we will outline several strategies to reduce artificial discontinuities in the ozone record in future NASA reanalyses. This discussion will be illustrated by an example of joint assimilation of bias-corrected ozone profiles from the Microwave Limb Sounder (MLS) on the Aura satellite (2004 to present) and the Ozone Mapping Profiler Suite Limb Profiler (OMPS-LP) sensors that are expected to operate on future NOAA platforms.

Wargan, Krzysztof↗

The Reanalysis for Stratospheric Trace-gas Studies

In order to re-examine trace gas transport in the middle atmosphere for the period May 1991 until April 1995, a "reanalysis" is being performed using an up-to-date version of the DAO's "GEOS" assimilation system. The Reanalysis for Stratospheric Trace-gas Studies (ReSTS) is intended to provide state-of-the-art estimates of the atmosphere during a period when the Upper Atmospheric Research Satellite provided a high density of trace-gas observations, and when the aerosol loading from the eruption of Mount Pinatubo contaminated the lower stratosphere, at the same time performing a natural tracer transport experiment. This study will present the first results from ReSTS, focussing on the improvements over the meteorological analyses produced by the then-operational GEOS-1 data assimilation system; emphasis will be placed on the improved representations of physical processes between GEOS-1 and the current GEOS-4 systems, highlighting the transport properties of the datasets. Alongside the production of a comprehensive atmospheric dataset, important components of ReSTS include performing sensitivity studies to the formulation of the assimilation system (including the representation of physical processes in the GCM, such as feedbacks between ozone/aerosols and meteorology) and to the inclusion of additional data types (including limb-sounding temperature data alongside the TOVS observations). Impacts of some of these factors on the analyzed meteorology and transport will be discussed. Of particular interest are attempts to determine the relative importance of various steps in the assimilation process to the quality of the final analyses.

Pawson, Steven↗

TESS Revisits WASP-12: Updated Orbital Decay Rate and Constraints on Atmospheric Variability

After observing WASP-12 in the second year of the primary mission, the Transiting Exoplanet Survey Satellite (TESS) revisited the system in late 2021 during its extended mission. In this paper, we incorporate the new TESS photometry into a reanalysis of the transits, secondary eclipses, and phase curve. We also present a new Ks-band occultation observation of WASP-12b obtained with the Palomar/Wide-field Infrared Camera instrument. The latest TESS photometry spans three consecutive months, quadrupling the total length of the TESS WASP-12 light curve and extending the overall time baseline by almost two years. Based on the full set of available transit and occultation timings, we find that the orbital period is shrinking at a rate of −29.81 ± 0.94 ms yr−1. The additional data also increase the measurement precision of the transit depth, orbital parameters, and phase-curve amplitudes. We obtain a secondary eclipse depth of 466 ± 35 ppm, a 2σ upper limit on the nightside brightness of 70 ppm, and a marginal 6°.2 ± 2°.8 eastward offset between the dayside hotspot and the substellar point. The voluminous TESS data set allows us to assess the level of atmospheric variability on timescales of days, months, and years. We do not detect any statistically significant modulations in the secondary eclipse depth or day–night brightness contrast. Likewise, our measured Ks-band occultation depth of 2810 ± 390 ppm is consistent with most ∼2.2 μm observations in the literature.

Ian Wong↗

Assessing the Altitude and Dispersion of Volcanic Plumes Using MISR Multi-angle Imaging from Space: Sixteen Years of Volcanic Activity in the Kamchatka Peninsula, Russia

Volcanic eruptions represent a significant source of atmospheric aerosols and can display local, regional and global effects, impacting earth systems and human populations. In order to assess the relative impacts of these events, accurate plume injection altitude measurements are needed. In this work, volcanic plumes generated from seven Kamchatka Peninsula volcanoes (Shiveluch, Kliuchevskoi, Bezymianny, Tolbachik, Kizimen, Karymsky and Zhupanovsky), were identified using over 16 years of Multi-angle Imaging SpectroRadimeter (MISR) measurements. Eighty-eight volcanic plumes were observed by MISR, capturing 3-25% of reported events at individual volcanoes. Retrievals were most successful where high intensity events persisted over a period of weeks to months. Compared with existing ground and airborne observations, and alternative satellite-based reports compiled by the Global Volcanism Program (GVP), MISR plume height retrievals showed general consistency; the comparison reports appear to be skewed towards the region of highest concentration observed in MISR-constrained vertical plume extent. The report observations display less discrepancy with MISR toward the end of the analysis period, with improvements in the suborbital data likely the result of the deployment of new instrumentation. Conversely, the general consistency of MISR plume heights with conventionally reported observations supports the use of MISR in the ongoing assessment of volcanic activity globally, especially where other types of volcanic plume observations are unavailable. Differences between the northern (Shiveluch, Kliuchevskoi, Bezymianny and Tolbachik) and southern (Kizimen, Karymsky and Zhupanovsky) volcanoes broadly correspond to the Central Kamchatka Depression (CKD) and Eastern Volcanic Front (EVF), respectively, geological sub-regions of Kamchatka distinguished by varying magma composition. For example, by comparison with reanalysis-model simulations of local meteorological conditions, CKD plumes generally were less constrained by mid-tropospheric (< 6 km) layers of vertical stability above the boundary layer, suggesting that these eruptions were more energetic than those in the EVF

Kamchatka↗

Assessment of WRF (v 4.2.1) dynamically downscaled precipitation on subdaily and daily timescales over CONUS

Abstract. This study analyzes the quality of simulated historical precipitation across the contiguous United States (CONUS) in a 12 km Weather Research and Forecasting model version 4.2.1 (WRF v 4.2.1)-based dynamical downscaling of the fifth-generation ECMWF atmospheric reanalysis (ERA5). This work addresses the following questions. First, how well are the 3 and 24 h precipitation characteristics (diurnal and annual cycles, precipitation frequency, annual and seasonal mean and maximum precipitation, and distribution of seasonal maximum precipitation) represented in the downscaled simulation, compared to ERA5? And second, how does the performance of the simulated WRF precipitation vary across seasons, regions, and timescales? Performance is measured against the National Centers for Environmental Prediction/Environmental Modeling Center (NCEP/EMC) 4 km Stage IV and Oregon State University Parameter-Elevation Regressions on Independent Slopes Model (PRISM) data on 3 and 24 h timescales, respectively. Our analysis suggests that the 12 km WRF exhibits biases typically found in other WRF simulations, including those at convection-permitting scales. In particular, WRF simulates both the timing and magnitude of the summer diurnal precipitation peak as well as ERA5 over most of the CONUS, except for a delayed diurnal peak over the Great Plains. As compared to ERA5, both the month and the magnitude of the precipitation peak annual cycle are remarkably improved in the downscaled WRF simulation. WRF slightly overestimates 3 and 24 h precipitation maximum over the CONUS, in contrast to ERA5, which generally underestimates these quantities mainly over the eastern half of the CONUS. Notably, WRF better captures the probability density distribution (PDF) of 3 and 24 h annual and seasonal maximum precipitation. WRF exhibits seasonally dependent precipitation biases across the CONUS, while ERA5's biases are relatively consistent year round over most of the CONUS. These results suggest that dynamical downscaling to a higher resolution improves upon some precipitation metrics but is susceptible to common regional climate model biases. Consequently, if used as input data for domain-specific models, we suggest moderate bias correction be applied to the dynamically downscaled product.

54 ENVIRONMENTAL SCIENCES↗

Modeling Spatial Asymmetries in Teleconnected Extreme Temperatures

Abstract Combining strengths from deep learning and extreme value theory can help describe complex relationships between variables where extreme events have significant impacts (e.g., environmental or financial applications). Neural networks learn complicated nonlinear relationships from large datasets under limited parametric assumptions. By definition, the number of occurrences of extreme events is small, which limits the ability of the data-hungry, nonparametric neural network to describe rare events. Inspired by recent extreme cold winter weather events in North America caused by atmospheric blocking, we examine several probabilistic generative models for the entire multivariate probability distribution of daily boreal winter surface air temperature. We propose metrics to measure spatial asymmetries, such as long-range anticorrelated patterns that commonly appear in temperature fields during blocking events. Compared to vine copulas, the statistical standard for multivariate copula modeling, deep learning methods show improved ability to reproduce complicated asymmetries in the spatial distribution of ERA5 temperature reanalysis, including the spatial extent of in-sample extreme events.

Krock, Mitchell L.↗