Search NASA⌕ Search

SEARCH · Search NASA

Results for “geostrophic”

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 127 records · Page 7

Moist wind relationships

The equations describing the temporal and spatial behavior of the kinematic moisture and heat flux are described. In these nonlinear equations, the contribution by diabatic processes to the large-scale flux is composed of two parts. One part is associated with a Rayleigh damping term while the other arises from temporal and spatial changes in the pressure gradient term. It was found that the influence of diabatic processes on large-scale moisture fluxes depends greatly on the degree of balance between forcing and damping terms in the governing kinematic flux equations. The existence of a near balance requires a reduction in the large-scale horizontal geostrophic wind speed. Based on an examination of the moisture flux equations, it is argued that reductions in the large-scale horizontal wind speed observed within major cumulus cloud systems help conserve large-scale moisture fluxes. The deviation of the wind from geostrophic conditions is easily estimated for a near balanced state. This wind modification induces secondary vertical circulations that contribute to convergence, creating or supporting long-lived mesoscale flows. We believe this process to be a major supporter of the mesoscale circulations observed in severe storms and squall lines. In the tropics the wind modification has an antitriptic relationship. These diagnostic findings suggest possible modifications to the wind field in the application of a cumulus parameterization, and may be important in diabatic initialization of numerical weather prediction models.

Raymond, William H.↗

Jupiter's winds and Arnol'd's second stability theorem: Slowly moving waves and neutral stability

Since the Voyager encounters in 1979, it has been known that Jupiter's cloud-top zonal winds violate the barotropic stability criterion. A vortex-tube stretching analysis of the Voyager wind data indicates that the more general Charney-Stern stability criterion is also violated. On the other hand, the zonal winds determined by tracking cloud features in Hubble Space Telescope images taken in 1991 precisely match the zonal winds determined by tracking cloud features in Voyager images, and it is hard to understand how a complicated zonal wind profile like Jupiter's could be unstable and yet not change at all in 12 years. In fact, there are at least two unknown ways to violate the Charney-Stern stability criterion and still have a stable flow. The better known of these is called Fjortoft's theorem, or Arnol'd's 1st theorem for the case of large-amplitude perturbations. Although the Fjortoft-Arnol'd theorem has been extended from the quasi-geostrophic equations to the primitive equations, the basic requirement that the potential vorticity be an increasing function of streamfunction is opposite to the case found in Jupiter, where the Voyager data indicate that the potential vorticity is a decreasing function of streamfunction. But this second case is precisely that which is covered by Arnol'd's 2nd stability theorem. In fact, the Voyager data suggest that Jupiter's zonal winds are neutrally stable with respect to Arnol'd's 2nd stability theorem. Here, we analyze the linear stability problem of a one-parameter family of sinusoidal zonal wind profiles that are close to neutral stability with respect to Arnol'd's 2nd stability theorem. We find numerically that the most unstable mode is always stationary, which may help to explain the slowly moving mode 10 waves observed on Jupiter. We find that violation of Arnol'd's 2nd stability theorem is both necessary and sufficient for instability of sinusoidal profiles. However, there appears to be no simple extension of Arnol'd's 2nd stability theorem to the primitive equations. Nevertheless, the primitive growth rates are small, and the primitive system is still governed by the quasi-geostrophic neutral-stability configuration.

Stamp, Andrew P.↗

On the Linearly-Balanced Kinetic Energy Spectrum

It is well known that the earth's atmospheric motion can generally be characterized by the two dimensional quasi-geostrophic approximation, in which the constraints on global integrals of kinetic energy, entrophy and potential vorticity play very important roles in redistributing the wave energy among different scales of motion. Assuming the hypothesis of Kolmogrov's local isotropy, derived a -3 power law of the equilibrium two-dimensional kinetic energy spectrum that entails constant vorticity and zero energy flows from the energy-containing wave number up to the viscous cutoff. In his three dimensional quasi-geostrophic theory, showed that the spectrum function of the vertical scale turbulence - expressible in terms of the available potential energy - possesses the same power law as the two dimensional kinetic energy spectrum. As the slope of kinetic energy spectrum in the inertial range is theoretically related to the predictability of the synoptic scales (Lorenz, 1969), many general circulation models includes a horizontal diffusion to provide reasonable kinetic energy spectra, although the actual power law exhibited in the atmospheric general circulation is controversial. Note that in either the atmospheric modeling or the observational analyses, the proper choice of wave number Index to represent the turbulence scale Is the degree of the Legendre polynomial.

Lu, Huei,-Iin↗

Can We Infer Ocean Dynamics from Altimeter Wavenumber Spectra?

The wavenumber spectra of sea surface height (SSH) and kinetic energy (KE) have been used to infer the dynamics of the ocean. When quasi-geostrophic dynamics (QG) or surface quasi-geostrophic (SQG) turbulence dominate and an inertial subrange exists, a steep SSH wavenumber spectrum is expected with k-5 for QG turbulence and a flatter k-11/3 for SQG turbulence. However, inspection of the spectral slopes in the mesoscale band of 70 to 250 km shows that the altimeter wavenumber slopes typically are much flatter than the QG or SQG predictions over most of the ocean. Comparison of the altimeter wavenumber spectra with the spectra estimated from the output of an eddy resolving global ocean circulation model (the Hybrid Coordinate Ocean Model, HYCOM, at 1/25 resolution), which is forced by high frequency winds and includes the astronomical forcing of the sun and the moon, suggests that the flatter slopes of the altimeter may arise from three possible sources, the presence of internal waves, the lack of an inertial subrange in the 70 to 250 km band and noise or submesoscales at small scales. When the wavenumber spectra of SSH and KE are estimated near the internal tide generating regions, the resulting spectra are much flatter than the expectations of QG or SQG theory. If the height and velocity variability are separated into low frequency (periods greater than 2 days) and high frequency (periods less than a day), then a different pattern emerges with a relatively flat wavenumber spectrum at high frequency and a steeper wavenumber spectrum at low frequency. The stationary internal tides can be removed from the altimeter spectrum, which steepens the spectral slopes in the energetic internal wave regions. Away from generating regions where the internal waves

Richman, James↗

Robustness of the Ensemble Score Filter to the Type of Assimilated Observation Networks

Recent advances in data assimilation (DA) have focused on developing more flexible approaches that can better accommodate nonlinearities in models and observations. However, it remains unclear how the performance of these advanced methods depends on the observation network characteristics. In this study, we present initial experiments with the surface quasi‐geostrophic model, in which we compare a recently developed ensemble filter using score‐based diffusion models with the standard Local Ensemble Transform Kalman Filter (LETKF). Our results show that the analysis solutions respond differently to the number, spatial distribution, and nonlinear fraction of assimilated observations. We also find notable changes in the multiscale characteristics of the analysis errors. Given that standard DA techniques will eventually be replaced by more advanced methods, we hope this study sets the ground for future efforts to reassess the value of Earth observing systems in the context of newly emerging algorithms.

97 MATHEMATICS AND COMPUTING↗

Conservative Numerical Schemes with Optimal Dispersive Wave Relations: Part II. Numerical Evaluations

A new energy and enstrophy conserving scheme (EEC) for the shallow water equations is proposed and evaluated using a suite of test cases over the global spherical or bounded domain. The evaluation is organized around a set of pre-defined properties: accuracy of individual operators, accuracy of the whole scheme, conservation of key quantities, control of the divergence variable, representation of the energy and enstrophy spectra, and simulation of nonlinear dynamics. The results confirm that the scheme is between the first and second order accurate, and conserves the total energy and potential enstrophy up to the time truncation errors. Here, the scheme is capable of producing more physically realistic energy and enstrophy spectra, indicating that it can help prevent the unphysical energy cascade towards the finest resolvable scales. With an optimal representation of the dispersive wave relations, the scheme is able to keep the flow close to being non-divergent, and maintain the geostrophically balanced structures with large-scale geophysical flows over long-term simulations.

54 ENVIRONMENTAL SCIENCES↗

Operator inference with roll outs for learning reduced models from scarce and low-quality data

Data-driven modeling has become a key building block in computational science and engineering. However, data that are available in science and engineering are typically scarce, often polluted with noise and affected by measurement errors and other perturbations, which makes learning the dynamics of systems challenging. Here, in this work, we propose to combine data-driven modeling via operator inference with the dynamic training via roll outs of neural ordinary differential equations. Operator inference with roll outs inherits interpretability, scalability, and structure preservation of traditional operator inference while leveraging the dynamic training via roll outs over multiple time steps to increase stability and robustness for learning from low-quality and noisy data. Numerical experiments with data describing shallow water waves and surface quasi-geostrophic dynamics demonstrate that operator inference with roll outs provides predictive models from training trajectories even if data are sampled sparsely in time and polluted with noise of up to 10%.

97 MATHEMATICS AND COMPUTING↗

Data-driven variational multiscale reduced order models

We propose a new data-driven reduced order model (ROM) framework that centers around the hierarchical structure of the variational multiscale (VMS) methodology and utilizes data to increase the ROM accuracy at a modest computational cost. The VMS methodology is a natural fit for the hierarchical structure of the ROM basis: In the first step, we use the ROM projection to separate the scales into three categories: (i) resolved large scales, (ii) resolved small scales, and (iii) unresolved scales. In the second step, we explicitly identify the VMS–ROM closure terms, i.e., the terms representing the interactions among the three types of scales. In the third step, we use available data to model the VMS–ROM closure terms. Thus, instead of phenomenological models used in VMS for standard numerical discretizations (e.g., eddy viscosity models), we utilize available data to construct new structural VMS–ROM closure models. Specifically, we build ROM operators (vectors, matrices, and tensors) that are closest to the true ROM closure terms evaluated with the available data. We test the new data-driven VMS–ROM in the numerical simulation of four test cases: (i) the 1D Burgers equation with viscosity coefficient v = 10 -3 ; (ii) a 2D flow past a circular cylinder at Reynolds numbers, Re = 100, Re = 500, and Re= 1000; (iii) the quasi-geostrophic equations at Reynolds number Re = 450 and Rossby number Ro = 0.0036; and (iv) a 2D flow over a backward facing step at Reynolds number Re = 1000. The numerical results show that the data-driven VMS–ROM is significantly more accurate than standard ROMs.

42 ENGINEERING↗

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗

Observations and Modeling of Ocean Circulation in the Seychelles Plateau Region

The ocean circulation around and over the Seychelles Plateau is characterized using 35 months of temperature and velocity measurements along with a numerical model. The results here provide the first documented description of the ocean circulation atop the Seychelles Plateau. The Seychelles Plateau is an unusually broad (~200 km), shallow (~50 m) plateau, dropping off steeply to the abyss. It is situated in a dynamic location (3.5-5.5S, 54-57E) in the south-western tropical Indian Ocean where northwesterly winds are present during austral summer and become southeasterly in austral winter, following the reversal of the western Indian Ocean monsoon winds. Measurements around the Inner Islands, on the Seychelles Plateau, have been carried out since 2015. Velocity measurements show that most of the depth-averaged current variance on the Seychelles Plateau arises from near-inertial oscillations and lower-frequency variability. Lower-frequency variability encompasses seasonal and intraseasonal variability, the latter of which includes the effects of mixed Rossby-gravity waves and mesoscale eddies. Here, a global 0.1-degree numerical ocean simulation is used in conjunction with these observations to describe the regional circulation around the Seychelles Plateau. Atop the Seychelles Plateau, circulation is dominated by ageostrophic processes consistent with Ekman dynamics, while around the Seychelles Plateau, both geostrophic and ageostrophic processes are important and vary seasonally. Stratification responds to the sea surface height semi-annual signal which is due to Ekman pumping-driven upwelling (related to the Seychelles-Chagos Thermocline Ridge) and the arrival of an annual downwelling Rossby wave.

58 GEOSCIENCES↗

Assessing Historical Variability of South Asian Monsoon Lows and Depressions With an Optimized Tracking Algorithm

Abstract Cyclonic low‐pressure systems (LPS) produce abundant rainfall in South Asia, where they are traditionally categorized as monsoon lows, monsoon depressions, and more intense cyclonic storms. The India Meteorological Department (IMD) has tracked monsoon depressions for over a century, finding a large decline in their number in recent decades, but their methods have changed over time and do not include monsoon lows. This study presents a fast, objective algorithm for identifying monsoon LPS and uses it to assess interannual variability and trends in reanalyses. Variables and thresholds used in the algorithm are selected to best match a subjectively analyzed LPS data set while minimizing disagreement between four reanalyses in a training period. The stream function of 850 hPa horizontal wind is found to be optimal in this sense; it is less noisy than vorticity and represents the complete nondivergent wind, even when flow is not geostrophic. Using this algorithm, LPS statistics are computed for five reanalyses, and none show a detectable trend in monsoon depression counts since 1979. Both the Japanese 55‐year Reanalysis (JRA‐55) and the IMD data set show a step‐like reduction in depression counts when they began using geostationary satellite data, in 1979 and 1982, respectively; the 1958–2018 linear trend in JRA‐55, however, is smaller than in the IMD data set, and its error bar includes 0. There are more LPS in seasons with above‐average monsoon rainfall and in La Niña years, but few other large‐scale modes of interannual variability are found to modulate LPS counts, lifetimes, or track length consistently across reanalyses.

54 ENVIRONMENTAL SCIENCES↗

Upper-Tropospheric Troughs and North American Monsoon Rainfall in a Long-Term Track Dataset

The North American monsoon is frequently affected by transient, propagating upper tropospheric vorticity anomalies. Sometimes called Tropical Upper-Tropospheric Troughs (TUTTs), these features have been claimed to episodically enhance monsoon rainfall. Here, we track long-lived TUTTs in 40 years of reanalysis data, producing composites and case studies from 340 TUTTs which last, on average, 7 days as they move westward across the North American monsoon region. TUTTs are thought to form from midlatitude Rossby wave breaking; case studies from our dataset support this theory. TUTTs move westward within the easterly upper-level flow in which they are embedded. In vortex-centered composites along the full tracks of long-lived TUTTs, here, we find no detectable increase in rainfall within the main TUTT circulation. Instead, negative precipitation anomalies lie within about 500 km of the TUTT center. Quasi-geostrophic ascent occurs in the southeast quadrant of TUTTs but is confined to the upper troposphere and does not appear to interact with precipitation. Positive anomalies of ascent and rainfall occur south and southeast of TUTTs but lie outside the main TUTT vortex, perhaps indicating concurrent variations in nearby climatological precipitation maxima. In contrast with previous case studies and subjective analyses that showed TUTTs enhance precipitation in parts of northwestern Mexico, our composites along the tracks of long-lived TUTTs portray these systems, to first order, as strong vorticity anomalies trapped in the upper troposphere that interact only weakly and indirectly with precipitation.

54 ENVIRONMENTAL SCIENCES↗

Synoptic Drivers of Atmospheric River Induced Precipitation Near Dronning Maud Land, Antarctica

Atmospheric rivers (ARs) that reach the Antarctic Ice Sheet (AIS) transport anomalous moisture from lower latitudes and can impact the AIS via extreme precipitation and increased downward longwave radiation. ARs contribute significantly to the interannual variability of precipitation over the AIS and thus are likely to play a key role in understanding future changes in the surface mass balance of the AIS. Dronning Maud Land (DML) is one of four maxima in AR frequency over coastal East Antarctica, with AR precipitation explaining 77% of the interannual variability in precipitation for this region. We employ a 16-node self-organizing map (SOM) trained with MERRA-2 sea-level pressure anomalies to identify synoptic-scale environments associated with landfalling ARs in and around DML. Node composites of atmospheric variables reveal common drivers of precipitation associated with ARs reaching DML including anomalous high-low surface pressure couplets, anomalously high integrated water vapor, and coastal barrier jets. Using a quasi-geostrophic framework, we find that upward vertical motion associated with the occlusion process of attendant cyclones dominates atmospheric lift in AR environments. We further identify mechanisms that explain the variability in AR precipitation intensity across nodes, such as the lift associated with the occlusion process of attendant cyclones and the spatial coincidence of ascent induced by the occlusion process and frontogenesis. Here, the latter suggests that ARs making landfall during the mature phase of cyclogenesis result in higher precipitation intensity compared to landfalling ARs that occur during the occluded phase.

54 ENVIRONMENTAL SCIENCES↗

Evaluating the Representations of Atmospheric Rivers and Their Associated Precipitation in Reanalyses With Satellite Observations

Atmospheric rivers (ARs) are filaments of enhanced horizontal moisture transport in the atmosphere. Due to their prominent role in the meridional moisture transport and regional weather extremes, ARs have been studied extensively in recent years. Yet, the representations of ARs and their associated precipitation on a global scale remains largely unknown. In this study, we developed an AR detection algorithm specifically for satellite observations using moisture and the geostrophic winds derived from 3D geopotential height field from the combined retrievals of the Atmospheric Infrared Sounder and the Advanced Microwave Sounding Unit on NASA Aqua satellite. This algorithm enables us to develop the first global AR catalog based solely on satellite observations. The satellite–based AR catalog is then combined with the satellite–based precipitation (Integrated Muti–SatellitE Retrievals for GPM) to evaluate the representations of ARs and AR–induced precipitation in reanalysis products. Here our results show that the spreads in AR frequency and AR length distribution are generally small across data sets, while the spread in AR width is relatively larger. Reanalysis products are found to consistently underestimate both mean and extreme AR–related precipitation. However, all reanalyses tend to precipitate too often under AR conditions, especially over low latitude regions. This finding is consistent with the “drizzling” bias which has plagued generations of climate models. Overall, the findings of this study can help to improve the representations of ARs and associated precipitation in reanalyses and climate models.

54 ENVIRONMENTAL SCIENCES↗

Temporal Subsampling Diminishes Small Spatial Scales in Recurrent Neural Network Emulators of Geophysical Turbulence

The immense computational cost of traditional numerical weather and climate models has sparked the development of machine learning (ML) based emulators. Because ML methods benefit from long records of training data, it is common to use data sets that are temporally subsampled relative to the time steps required for the numerical integration of differential equations. Here, we investigate how this often overlooked processing step affects the quality of an emulator's predictions. We implement two ML architectures from a class of methods called reservoir computing: (a) a form of Nonlinear Vector Autoregression (NVAR), and (b) an Echo State Network (ESN). Despite their simplicity, it is well documented that these architectures excel at predicting low dimensional chaotic dynamics. We are therefore motivated to test these architectures in an idealized setting of predicting high dimensional geophysical turbulence as represented by Surface Quasi-Geostrophic dynamics. In all cases, subsampling the training data consistently leads to an increased bias at small spatial scales that resembles numerical diffusion. Interestingly, the NVAR architecture becomes unstable when the temporal resolution is increased, indicating that the polynomial based interactions are insufficient at capturing the detailed nonlinearities of the turbulent flow. The ESN architecture is found to be more robust, suggesting a benefit to the more expensive but more general structure. Spectral errors are reduced by including a penalty on the kinetic energy density spectrum during training, although the subsampling related errors persist. Future work is warranted to understand how the temporal resolution of training data affects other ML architectures.

58 GEOSCIENCES↗

A Non‐Intrusive Machine Learning Framework for Debiasing Long‐Time Coarse Resolution Climate Simulations and Quantifying Rare Events Statistics

Abstract Due to the rapidly changing climate, the frequency and severity of extreme weather is expected to increase over the coming decades. As fully‐resolved climate simulations remain computationally intractable, policy makers must rely on coarse‐models to quantify risk for extremes. However, coarse models suffer from inherent bias due to the ignored “sub‐grid” scales. We propose a framework to non‐intrusively debias coarse‐resolution climate predictions using neural‐network (NN) correction operators. Previous efforts have attempted to train such operators using loss functions that match statistics. However, this approach falls short with events that have longer return period than that of the training data, since the reference statistics have not converged. Here, the scope is to formulate a learning method that allows for correction of dynamics and quantification of extreme events with longer return period than the training data. The key obstacle is the chaotic nature of the underlying dynamics. To overcome this challenge, we introduce a dynamical systems approach where the correction operator is trained using reference data and a coarse model simulation nudged toward that reference. The method is demonstrated on debiasing an under‐resolved quasi‐geostrophic model and the Energy Exascale Earth System Model (E3SM). For the former, our method enables the quantification of events that have return period two orders longer than the training data. For the latter, when trained on 8 years of ERA5 data, our approach is able to correct the coarse E3SM output to closely reflect the 36‐year ERA5 statistics for all prognostic variables and significantly reduce their spatial biases.

Barthel Sorensen, B.↗

Observations of Subduction, Downward Heat Flux and Dense Filament Collapse in the Northern Gulf of Mexico

Submesoscale processes are important contributors to the global heat budget and generally support upward heat transport through restratification. However, in salinity‐stratified regions, such as the northern Gulf of Mexico with its influx of freshwater from the Mississippi‐Atchafalaya river system, temperature can act like a passive tracer and submesoscale processes can contribute to downward heat transport. Oceanic heat content is a factor in many environmental risks the region faces, for example, hurricane intensification, and marine heatwaves. During the 2022 field campaign of the Submesoscales Under Near‐Resonant Inertial Shear Experiment, a sampling plan was developed to study such submesoscale processes in high resolution. Over 31 hr, four assets (two research ships and two remotely controlled boats) drove in parallel across a dense filament, capturing its evolution in time and space. The observations show that surface waters, warmed by daytime solar radiation, were subducted and that the associated overturning circulation transported heat below the surface layer where it was later irreversibly mixed away. The estimated downward heat flux was as strong as the concurrent net air‐sea heat flux into the ocean. The filament was then observed to rapidly collapse which we attribute to boundary layer turbulence and the breakdown of geostrophic balance. The collapsing fronts display behaviors indicative of gravity currents. These observations highlight how in salinity‐stratified regions, frontal dynamics can be associated with downward heat flux and how the submesoscale can play an important role in the oceanic heat budget.

58 GEOSCIENCES↗

Pollen antigens and atmospheric circulation driven seasonal respiratory viral outbreak and its implication to the Covid-19 pandemic

Abstract The patterns of respiratory virus illness are expressed differently between temperate and tropical climates. Tropical outbreaks often peak in wet seasons. Temperate outbreaks typically peak during the winter. The prevailing causal hypotheses focus on sunlight, temperature and humidity variations. Yet no consistent factors have been identified to sufficiently explain seasonal virus emergence and decline at any latitude. Here we demonstrate close connections among global-scale atmospheric circulations, IgE antibody enhancement through seasonal pollen inhalation, and respiratory virus patterns at any populated latitude, with a focus on the US. Pollens emerge each Spring, and the renewed IgE titers in the population are argued to terminate each winter peak of respiratory illness. Globally circulated airborne viruses are postulated to subsequently deposit across the Southern US during lower zonal geostrophic winds each late Summer. This seasonally refreshed viral load is postulated to trigger a new influenza outbreak, once the existing IgE antibodies diminish to a critical value each Fall. Our study offers a new and consistent explanation for the seasonal diminishment of respiratory viral illnesses in temperate climates, the subdued seasonal signature in the tropics, the annually circulated virus phenotypes, and the northerly migration of influenza across the US every year. Our integrated geospatial and IgE hypothesis provides a new perspective for prediction, mitigation and prevention of the outbreak and spread of seasonal respiratory viruses including Covid-19 pandemic.

60 APPLIED LIFE SCIENCES↗