Search NASA⌕ Search

SEARCH · Search NASA

Results for “model reanalysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 91 records · Page 5

A new metrics framework for quantifying and intercomparing atmospheric rivers in observations, reanalyses, and climate models

We present a new atmospheric river (AR) analysis and benchmarking tool, namely Atmospheric River Metrics Package (ARMP). It includes a suite of new AR metrics that are designed for quick analysis of AR characteristics via statistics in gridded climate datasets such as model output and reanalysis. This package can be used for climate model evaluation in comparison with reanalysis and observational products. Integrated metrics such as mean bias and spatial pattern correlation are efficient for diagnosing systematic AR biases in climate models. For example, the package identifies the fact that, in CMIP5 and CMIP6 (Coupled Model Intercomparison Project Phases 5 and 6) models, AR tracks in the South Atlantic are positioned farther poleward compared to ERA5 reanalysis, while in the South Pacific, tracks are generally biased towards the Equator. For the landfalling AR peak season, we find that most climate models simulate a completely opposite seasonal cycle over western Africa. This tool can also be used for identifying and characterizing structural differences among different AR detectors (ARDTs). For example, ARs detected with the Mundhenk algorithm exhibit systematically larger size, width, and length compared to the TempestExtremes (TE) method. The AR metrics developed from this work can be routinely applied for model benchmarking and during the development cycle to trace performance evolution across model versions or generations and set objective targets for the improvement of models. They can also be used by operational centers to perform near-real-time climate and extreme event impact assessments as part of their forecast cycle.

58 GEOSCIENCES↗

Land Surface Verification Toolkit (LVT) - A Generalized Framework for Land Surface Model Evaluation

Model evaluation and verification are key in improving the usage and applicability of simulation models for real-world applications. In this article, the development and capabilities of a formal system for land surface model evaluation called the Land surface Verification Toolkit (LVT) is described. LVT is designed to provide an integrated environment for systematic land model evaluation and facilitates a range of verification approaches and analysis capabilities. LVT operates across multiple temporal and spatial scales and employs a large suite of in-situ, remotely sensed and other model and reanalysis datasets in their native formats. In addition to the traditional accuracy-based measures, LVT also includes uncertainty and ensemble diagnostics, information theory measures, spatial similarity metrics and scale decomposition techniques that provide novel ways for performing diagnostic model evaluations. Though LVT was originally designed to support the land surface modeling and data assimilation framework known as the Land Information System (LIS), it also supports hydrological data products from other, non-LIS environments. In addition, the analysis of diagnostics from various computational subsystems of LIS including data assimilation, optimization and uncertainty estimation are supported within LVT. Together, LIS and LVT provide a robust end-to-end environment for enabling the concepts of model data fusion for hydrological applications. The evolving capabilities of LVT framework are expected to facilitate rapid model evaluation efforts and aid the definition and refinement of formal evaluation procedures for the land surface modeling community.

Kumar, Sujay V.↗

Computationally Efficient Modeling and Data Assimilation of Near-Surface Variability

Near-surface (< 20m) ocean exhibits high variability due to coupled interactions, for e.g., with the atmosphere, sea ice, land, etc. Here we focus on atmospheric heat and momentum (wind) forcing, which are known to cause diurnal variability within the mixed layer. Only recently with a combination of sufficiently high vertical/horizontal resolution (75L, 1/4deg) and sub-daily atmospheric forcing fields, ocean models are starting to resolve this diurnal variability. However, the computation expense of such a high vertical resolution is burdensome in the context of coupled modeling and data assimilation. An alternative approach is to parameterize this diurnal variability with a prognostic model, that is embedded into the ocean model.In the first part of this presentation, we will demonstrate results with the above two approaches, by comparing them to profiles of near-surface temperature and salinity. In the context of data assimilation and reanalysis, this modeling capability opens the door to re-examine and perhaps improve specification of background (or, ensemble) error characteristics. The second half of this talk will focus on illustrating diurnally varying errors within an ensemble DA, and possible approaches to improve localization (horizontal/vertical) to extract maximum possible observational information content from in-situ and satellite observations of sea surface temperature.

Akella, Santha↗

Identifying snowfall elevation patterns by assimilating satellite-based snow depth retrievals

Precipitation in mountain regions is highly variable and poorly measured, posing important challenges to water resource management. Traditional methods to estimate precipitation include in-situ gauges, Doppler weather radars, satellite radars and radiometers, numerical modeling and reanalysis products. Each of these methods is unable to adequately capture complex orographic precipitation. Here, we propose a novel approach to characterize orographic snowfall over mountain regions. We use a particle batch smoother to leverage satellite information from Sentinel-1 derived snow depth retrievals and to correct various gridded precipitation products. This novel approach is tested using a simple snow model for an alpine basin located in Trentino Alto Adige, Italy. Here, we quantify the precipitation biases across the basin and found that the assimilation method (i) corrects for snowfall biases and uncertainties, (ii) leads to cumulative snowfall elevation patterns that are consistent across precipitation products, and (iii) results in overall improved basin-wide snow variables (snow depth and snow cover area) and basin streamflow estimates.

54 ENVIRONMENTAL SCIENCES↗

Clifford Neural Operators on Atmospheric Data Influenced Partial Differential Equations

Mathematical representations of the atmosphere are key to forecasting and research tasks across Earth science. Numerically solving the underlying partial differential equations(PDEs) of the atmosphere, however, can be difficult and computationally expensive with numerous trade-offs between computing efficiency and accuracy. Utilizing neural net-works to learn approximations of the PDE solutions from the data can help us model complex phenomena more efficiently than traditional numerical schemes. Here, we have applied Clifford algebra-based neural operators for predicting atmospheric variables. Clifford Fourier neural operators are used with two different backbone architectures, ResNet and UNet, on custom data of U10, V10 and surface pressure as well as U500, V500 and Z500. Clifford Fourier neural operators, coupled with ResNet and UNet architectures, are applied to a key reanalysis dataset. Model performance is initially strong, but we observe increasing errors, resulting in the model becoming highly unstable.

Sujit Roy↗

Characterizing Ice Nucleating Particles Over the Southern Ocean Using Simultaneous Aircraft and Ship Observations

Abstract Supercooled liquid clouds are ubiquitous over the Southern Ocean (SO), even to temperatures below −20°C, and comprise a large fraction of the marine boundary layer (MBL) clouds. Earth system models and reanalysis products have struggled to reproduce the observed cloud phase distribution and occurrence of cloud ice in the region. Recent simulations found the microphysical representation of ice nucleation and growth has a large impact on these properties, however, measurements of SO ice nucleating particles (INPs) to validate simulations are sparse. This study presents measurements of INPs from simultaneous aircraft and ship campaigns conducted over the SO in austral summer 2018, which include the first in situ observations in and above cloud in the region. Our results confirm recent observations that INP concentrations are uniformly lower than measurements made in the late 1960s. While INP concentrations below and above cloud are similar, higher ice nucleation efficiency above cloud supports model simulations that the dominant INP composition varies with height. Model parameterizations based solely on aerosol properties capture the mean relationship between INP concentration and temperature but not the observed variability, which is likely related to the only modest correlations observed between INPs and environmental or aerosol metrics. Including wind speed in addition to activation temperature in a marine INP parameterization reduces bias but does not explain the large range of observed INP concentrations. Direct and indirect inference of marine INP size suggests MBL INPs, at least during Austral summer, are dominated by particles with diameters smaller than 500 nm.

Meteorology & Atmospheric Sciences↗

Land surface Verification Toolkit (LVT)

LVT is a framework developed to provide an automated, consolidated environment for systematic land surface model evaluation Includes support for a range of in-situ, remote-sensing and other model and reanalysis products. Supports the analysis of outputs from various LIS subsystems, including LIS-DA, LIS-OPT, LIS-UE. Note: The Land Information System Verification Toolkit (LVT) is a NASA software tool designed to enable the evaluation, analysis and comparison of outputs generated by the Land Information System (LIS). The LVT software is released under the terms and conditions of the NASA Open Source Agreement (NOSA) Version 1.1 or later. Land Information System Verification Toolkit (LVT) NOSA.

surface model↗

An Evaluation of Biomass Burning Aerosol Mass, Extinction, and Size Distribution in Geos Using Observations From Camp2EX

Biomass burning aerosol impacts aspects of the atmosphere and Earth system through radiative forcing, serving as cloud condensation nuclei, and air quality. Despite its importance, the representation of biomass burning aerosol is not always accurate in numerical weather prediction and climate models or reanalysis products. Using observations collected as part of the Cloud, Aerosol and Monsoon Processes Philippines Experiment (CAMP2Ex) in August through October of 2019, aerosol concentration and optical properties are evaluated within the Goddard Earth Observing System (GEOS) and its underlying aerosol module, GOCART. In the operational configuration, GEOS assimilates aerosol optical depth observations at 550 nm to constrain aerosol fields. Particularly for biomass burning aerosol, without the assimilation of aerosol optical depth, aerosol extinction is underestimated compared to observations collected in the Philippines region during the CAMP2Ex campaign. The assimilation process adds excessive amounts of carbon to account for the underestimated extinction, resulting in positive biases in the mass of black and organic carbon, especially within the boundary layer, relative to in situ observations from the Langley Aerosol Research Group Experiment. Counteracting this, GEOS is deficient in sulphate and nitrate aerosol just above the boundary layer. Aside from aerosol mass, extinction within GEOS is a function of ambient relative humidity and an assumed particle size distribution. The relationship between dry and ambient extinction in GEOS reveals that hygroscopic growth is too aggressive within the model for biomass burning aerosol. An additional concern lies in the assumed particle size distribution for GEOS, which has a mode radius that is too small for organic carbon. Variability in the observed particle size distribution for biomass burning aerosol within a single flight also illuminates the fact that a single assumed particle size distribution is not sufficient and that for a proper representation, a more advanced aerosol module with GEOS may be necessary.

Aerosol Mass↗

Evaluation of GEOS Aerosols Using CAMP2Ex Observations

Biomass burning aerosol impacts aspects of the atmosphere and Earth system through direct and semi-direct effects, as well as influencing air quality. Despite its importance, the representation of biomass burning aerosol is not always accurate in numerical weather prediction and climate models or reanalysis products. Using observations collected as part of the Cloud, Aerosol and Monsoon Processes Philippines Experiment (CAMP2Ex) in August through October of 2019, aerosol concentration and optical properties are evaluated within the Goddard Earth Observing System (GEOS) and its underlying aerosol module, GOCART. In the operational configuration, GEOS assimilates aerosol optical depth observations at 550 nm from AERONET and MODIS to constrain aerosol fields. Particularly for biomass burning aerosol, without the assimilation of aerosol optical depth, aerosol extinction is underestimated compared to observations collected in the Philippines region during the CAMP2Ex campaign. The assimilation process adds excessive amounts of carbon to account for the underestimated extinction, resulting in positive biases in the mass of black and organic carbon, especially within the boundary layer, relative to in situ observations from the Langley Aerosol Research Group Experiment. Counteracting this, GEOS is deficient in sulfate and nitrate aerosol just above the boundary layer. Aerosol extinction within GEOS is a function of the mass of different aerosol species, the ambient relative humidity, the assumed spectral optical properties, and particle size distribution per species. The relationship between dry and ambient extinction in GEOS reveals that hygroscopic growth is too high within the model for biomass burning aerosol. An additional concern lies in the assumed particle size distribution for GEOS, which has a single mode radius that is too small for organic carbon. Variability in the observed particle size distribution for biomass burning aerosol within a single flight also illuminates the fact that a single assumed particle size distribution is not sufficient and that for a proper representation, a more advanced aerosol module within GEOS may be necessary.

Allison Collow↗

Assessing Radiative Feedbacks and Their Contribution to the Arctic Amplification Measured by Various Metrics

Arctic amplification (AA), characterized by a more rapid surface air temperature (SAT) warming in the Arctic than the global average, is a major feature of global climate warming. Various metrics have been used to quantify AA based on SAT anomalies, trends, or variability, and they can yield quite different conclusions regarding the magnitude and temporal patterns of AA. This study examines and compares various AA metrics for their temporal consistency in the region north of 70°N from the early twentieth to the early 21st century using observational data and reanalysis products. We also quantify contributions of different radiative feedback mechanisms to AA based on short-term climate variability in reanalysis and model data using the Kernel-Gregory approach. Albedo and lapse rate feedbacks are positive and comparable, with albedo feedback being the leading contributor for all AA metrics. The net cloud feedback, which has large uncertainties, depends strongly on the data sets and AA metrics used. By quantifying the influence of internal variability on AA and related feedbacks based on global climate model ensemble simulations, we find that water vapor and cloud feedbacks are most heavily affected by internal variability.

54 ENVIRONMENTAL SCIENCES↗

Earth System Reanalysis in Support of Climate Model Improvements

Recent climate model developments, established through increased model resolution, have led to substantial improvements in model simulations of the time-evolving, coupled Earth system and its subcomponents. However, regardless of resolution, climate models will always produce climate features and variability that differ from the real world and will be prone to biases. This is due to many remaining uncertainties, such as in parametric and structural model uncertainty, in the initial conditions prescribed, and in the prescribed (scenario) forcing which varies on decadal to centennial timescales. Further model improvements are expected to arise specifically from improved representation of physical processes realized through model-data fusion. This will create an unprecedented opportunity to better exploit a large array of Earth observations, from in situ measurements to weather radars and satellite observations, as the resolved scales of the models approach those of the observations. For this, climate DA will be the central tool to bring models and observations into consistency, by improving initial conditions, inferring uncertain model parameters and structure, and quantifying uncertainty. Generally, there will be advantages and complementarities of adjoint-based smoother approaches, ensemble-based filter approaches, or new ML-inspired approaches. Yet, the ever-increasing model resolution will present growing challenges arising from computational cost, calling for new ways of performing data assimilation and model optimization. Using the complementarity in a hybrid approach, blending tools and concepts from variational, ensemble and ML methods might be what is required in the future. In this context ML could be important to handle non-linear responses, and to better approximate non-Gaussian distributions.

54 ENVIRONMENTAL SCIENCES↗

Metrics tool for for evaluating atmospheric rivers in climate data

The metric tool code is designed for evaluating atmospheric rivers (AR) in climate models and reanalysis. It is python based, with built in metrics of AR frequency, AR precipitation, AR peak day, AR characteristics (width, length, area, latitude and longitude), and output of diagnostic plots. Is there

Dong, Bo↗

Quantifying Spatial Separation Error in Tropospheric Wind Measurements

Wind measurements from the surface through the troposphere are an important asset for both range and Day of Launch (DOL) operations. Weather balloons have long been used to fulfill this operational need. Decker (2017) found balloons at Kennedy Space Center (KSC) can drift as far as 200 km from the launch site. Curtis et al. (2019) found Root Mean Square (RMS) wind deltas between 0.5 m/s and 6.3 m/s for separations between 30 km and 200 km from KSC using the North American Regional Reanalysis (NARR) model. The purpose of this study is to validate the model-based results found in Curtis et al. (2019) with observed wind and separation deltas calculated between the NASA Tropospheric Doppler Radar Wind Profiler (TDRWP) and weather balloons.

Curtis, Nathan↗

Unprecedented Atmospheric Conditions (1948–2019) Drive the 2019 Exceptional Melting Season Over the Greenland Ice Sheet

Understanding the role of atmospheric circulation anomalies on the surface mass balance of the Greenland ice sheet (GrIS) is fundamental for improving estimates of its current and future contributions to sea level rise. Here, we show, using a combination of remote sensing observations, regional climate model outputs, reanalysis data, and artificial neural networks, that unprecedented atmospheric conditions (1948–2019) occurring in the summer of 2019 over Greenland promoted new record or close-to-record values of surface mass balance (SMB), runoff, and snowfall. Specifically, runoff in 2019 ranked second within the 1948–2019 period (after 2012) and first in terms of surface mass balance negative anomaly for the hydrological year 1 September 2018–31 August 2019. The summer of 2019 was characterized by an exceptional persistence of anticyclonic conditions that, in conjunction with low albedo associated with reduced snowfall in summer, enhanced the melt–albedo feedback by promoting the absorption of solar radiation and favored advection of warm, moist air along the western portion of the ice sheet towards the north, where the surface melt has been the highest since 1948. The analysis of the frequency of daily 500 hPa geopotential heights obtained from artificial neural networks shows that the total number of days with the five most frequent atmospheric patterns that characterized the summer of 2019 was 5 standard deviations above the 1981–2010 mean, confirming the exceptional nature of the 2019 season over Greenland.

Greenland ice sheet↗

Diagnosis of Antarctic Blowing Snow Properties Using MERRA-2 Reanalysis with a Machine Learning Model

This paper presents the work on using a machine learning model to diagnose Antarctic blowing snow (BLSN) properties with the Modern Era Retrospective analysis for Research and Applications v2 (MERRA-2) data. We adopt the random forest classifier for BLSN identification and the random forest regressor for BLSN optical depth and height diagnosis. BLSN properties observed from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) are used as the truth for training the model. Using MERRA-2 fields such as snow age, surface elevation and pressure, temperature, specific humidity, and temperature gradient at the 2m level, and wind speed at the 10m level as input, reasonable results are achieved. Hourly blowing snow property diagnostics are generated with the trained model. Using the year 2010 as an example, it is shown that the Antarctic BLSN frequency is much higher over East than West Antarctica. High frequency months are from April to September, during which BLSN frequency exceeds 20% over East Antarctica. For May 2010, the BLSN snow frequency in the region is as high as 37%. Due to the suppression by strong surface-based inversions, larger values of BLSN height and optical depth are usually limited to the coastal regions, wherein the strength of surface-based inversions is weaker.

Antarctic↗

Subseasonal Forecasting and MJO Teleconnections in Machine Learning Weather Prediction Models

Abstract In recent years, machine‐learning (ML) models trained on reanalysis data have rivaled physics‐based forecast models in terms of performance skill for global weather forecasting. With increased rollout stability, the question of how these models perform for subseasonal to seasonal (S2S, week 3–8) forecasting has emerged. In this study we run a large set of subseasonal hindcasts over 2004–2023 to evaluate two ML weather forecast models at the S2S time scale, SFNO‐HENS (Nvidia, fully ML) and NeuralGCM (Google Research, hybrid). Corresponding hindcasts from the European Centre for Medium‐Range Weather Forecasts (ECMWF) are used as a baseline for comparison to a physics‐based model. Because our focus is on predicting moisture transport over the Western United States between October and March, we evaluate the models' prediction skill for the Madden‐Julian Oscillation (MJO) and its associated teleconnections in the North Pacific. We find that both ML models are competitive with the ECWMF model, with comparable skill in predicting the North Pacific large‐scale circulation and the MJO at week 3 and beyond. Even though overall the mid‐latitude subseasonal prediction skill remains low, the ML models exhibit interesting behavior such as a realistic propagation of the MJO across the Maritime Continent and realistic teleconnections. A SFNO‐HENS sensitivity experiment with altered initial conditions in the tropics demonstrates the stability of the model, and it illustrates the capability of ML models to represent important physical processes of the atmosphere at the S2S time scale. Plain Language Summary Predicting weather patterns and precipitation a few weeks in advance (subseasonal time scale) is of great interest for stakeholders such as water managers in the Southwest United States (US), where arid conditions prevail. Subseasonal forecasts from traditional weather forecast models exhibit low skill in the region, limiting their applicability. Here we examine whether the recent breakthrough in weather forecasting made with machine learning/artificial intelligence models can translate to improved subseasonal forecasts. Recently‐developed machine learning models exhibit comparable skill to a state‐of‐the‐art physics‐based model for predicting weather patterns in the North Pacific/North America region, and associated moisture transport. The same applies to their skill in predicting the tropical pattern, the Madden‐Julian Oscillation, and its important remote perturbations over the midlatitude East Pacific and Southwest US. Additionally, a perturbation experiment carried out with one of the machine learning models illustrates their ability to not only predict the evolution of atmospheric fields, but also to learn and represent physical processes such as tropics‐extratropics Rossby wave propagation. Key Points Two machine learning weather forecast models exhibit state‐of‐the‐art prediction skill at the subseasonal time scale in the Pacific sector The models equal ECWMF in terms of Madden‐Julian oscillation (MJO) prediction skill, and they accurately predict the MJO propagation and associated teleconnections The two machine‐learning models represent key physical processes for subseasonal prediction, despite being trained for weather forecasting

Peings, Yannick↗

Simulated convective systems using a cloud resolving model: Impact of large-scale temperature and moisture forcing using observations and GEOS-3 reanalysis

The GCE (Goddard Cumulus Ensemble) model, which has been developed and improved at NASA Goddard Space Flight Center over the past two decades, is considered as one of the finer and state-of-the-art CRMs (Cloud Resolving Models) in the research community. As the chosen CRM for a NASA Interdisciplinary Science (IDS) Project, GCE has recently been successfully upgraded into an MPI (Message Passing Interface) version with which great improvement has been achieved in computational efficiency, scalability, and portability. By basically using the large-scale temperature and moisture advective forcing, as well as the temperature, water vapor and wind fields obtained from TRMM (Tropical Rainfall Measuring Mission) field experiments such as SCSMEX (South China Sea Monsoon Experiment) and KWAJEX (Kwajalein Experiment), our recent 2-D and 3-D GCE simulations were able to capture detailed convective systems typical of the targeted (simulated) regions. The GEOS-3 [Goddard EOS (Earth Observing System) Version-3] reanalysis data have also been proposed and successfully implemented for usage in the proposed/performed GCE long-term simulations (i.e., aiming at producing massive simulated cloud data -- Cloud Library) in compensating the scarcity of real field experimental data in both time and space (location). Preliminary 2-D or 3-D pilot results using GEOS-3 data have generally showed good qualitative agreement (yet some quantitative difference) with the respective numerical results using the SCSMEX observations. The first objective of this paper is to ensure the GEOS-3 data quality by comparing the model results obtained from several pairs of simulations using the real observations and GEOS-3 reanalysis data. The different large-scale advective forcing obtained from these two kinds of resources (i.e., sounding observations and GEOS-3 reanalysis) has been considered as a major critical factor in producing various model results. The second objective of this paper is therefore to investigate and present such an impact of large-scale forcing on various modeled quantities (such as hydrometeors, rainfall, and etc.). A third objective is to validate the overall GCE 3-D model performance by comparing the numerical results with sounding observations, as well as available satellite retrievals.

Shie, C.-L.↗

Estimates of Southern Hemispheric Gravity Wave Momentum Fluxes across Observations, Reanalyses, and Kilometer-Scale Numerical Weather Prediction Model

Abstract Gravity waves (GWs) are among the key drivers of the meridional overturning circulation in the mesosphere and upper stratosphere. Their representation in climate models suffers from insufficient resolution and limited observational constraints on their parameterizations. This obscures assessments of middle atmospheric circulation changes in a changing climate. This study presents a comprehensive analysis of stratospheric GW activity above and downstream of the Andes from 1 to 15 August 2019, with special focus on GW representation ranging from an unprecedented kilometer-scale global forecast model (1.4 km ECMWF IFS), ground-based Rayleigh lidar (CORAL) observations, modern reanalysis (ERA5), to a coarse-resolution climate model (EMAC). Resolved vertical flux of zonal GW momentum (GWMF) is found to be stronger by a factor of at least 2–2.5 in IFS compared to ERA5. Compared to resolved GWMF in IFS, parameterizations in ERA5 and EMAC continue to inaccurately generate excessive GWMF poleward of 60°S, yielding prominent differences between resolved and parameterized GWMFs. A like-to-like validation of GW profiles in IFS and ERA5 reveals similar wave structures. Still, even at ∼1 km resolution, the resolved waves in IFS are weaker than those observed by lidar. Further, GWMF estimates across datasets reveal that temperature-based proxies, based on midfrequency approximations for linear GWs, overestimate GWMF due to simplifications and uncertainties in GW wavelength estimation from data. Overall, the analysis provides GWMF benchmarks for parameterization validation and calls for three-dimensional GW parameterizations, better upper-boundary treatment, and vertical resolution increases commensurate with increases in horizontal resolution in models, for a more realistic GW analysis. Significance Statement Gravity wave–induced momentum forcing forms a key component of the middle atmospheric circulation. However, complete knowledge of gravity waves, their atmospheric effects, and their long-term trends are obscured due to limited global observations, and the inability of current climate models to fully resolve them. This study combines a kilometer-scale forecast model, modern reanalysis, and a coarse-resolution climate model to first compare the resolved and parameterized momentum fluxes by gravity waves generated over the Andes, and then evaluate the fluxes using a state-of-the-art ground-based Rayleigh lidar. Our analysis reveals shortcomings in current model parameterizations of gravity waves in the middle atmosphere and highlights the sensitivity of the estimated flux to the formulation used.

Meteorology & Atmospheric Sciences↗