Search NASA⌕ Search

SEARCH · Search NASA

Results for “zonal”

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 703 records · Page 39

Temperature and Chemical Species Distributions in the Middle Atmosphere Observed during Titan's Late Northern Spring to Early Summer

We present a study of the seasonal evolution of Titan’s thermal field and distributions of haze, C2H2, C2H4, C2H6, CH3C2H, C3H8, C4H2, C6H6, HCN, and HC3N from March 2015 (Ls = 66°) to September 2017 (Ls = 93°) (i.e., from the last third of northern spring to early summer). We analyzed thermal emission of Titan’s atmosphere acquired by the Cassini Composite Infrared Spectrometer with limb and nadir geometry to retrieve the stratospheric and mesospheric temperature and mixing ratios pole-to-pole meridional cross sections from 5 mbar to 50 μbar (120–650 km). The southern stratopause varied in a complex way and showed a global temperature increase from 2015 to 2017 at high-southern latitudes. Stratospheric southern polar temperatures, which were observed to be as low as 120 K in early 2015 due to the polar night, showed a 30 K increase (at 0.5 mbar) from March 2015 to May 2017 due to adiabatic heating in the subsiding branch of the global overturning circulation. All photochemical compounds were enriched at the south pole by this subsidence. Polar cross sections of these enhanced species, which are good tracers of the global dynamics, highlighted changes in the structure of the southern polar vortex. These high enhancements combined with the unusually low temperatures (<120 K) of the deep stratosphere resulted in condensation at the south pole between 0.1 and 0.03 mbar (240–280 km) of HCN, HC3N, C6H6 and possibly C4H2 in March 2015 (Ls = 66°). These molecules were observed to condense deeper with increasing distance from the south pole. At high-northern latitudes, stratospheric enrichments remaining from the winter were observed below 300 km between 2015 and May 2017 (Ls = 90°) for all chemical compounds and up to September 2017 (Ls = 93°) for C2H2, C2H4, CH3C2H, C3H8, and C4H2. In September 2017, these local enhancements were less pronounced than earlier for C2H2, C4H2, CH3C2H, HC3N, and HCN, and were no longer observed for C2H6 and C6 H6, which suggests a change in the northern polar dynamics near the summer solstice. These enhancements observed during the entire spring may be due to confinement of this enriched air by a small remaining winter circulation cell that persisted in the low stratosphere up to the northern summer solstice, according to predictions of the Institut Pierre Simon Laplace Titan Global Climate Model (IPSL Titan GCM). In the mesosphere we derived a depleted layer in C2H2, HCN, and C2H6 from the north pole to mid-southern latitudes, while C4H2, C3H4, C2H4, and HC3N seem to have been enriched in the same region. In the deep stratosphere, all molecules except C2H4 were depleted due to their condensation sink located deeper than 5 mbar outside the southern polar vortex. HCN, C4H2, and CH3C2H volume mixing ratio cross section contours showed steep slopes near the mid-latitudes or close to the equator, which can be explained by upwelling air in this region. Upwelling is also supported by the cross section of the C2H4 (the only molecule not condensing among those studied here) volume mixing ratio observed in the northern hemisphere. We derived the zonal wind velocity up to mesospheric levels from the retrieved thermal field. We show that zonal winds were faster and more confined around the south pole in 2015 (Ls = 67−72°) than later. In 2016, the polar zonal wind speed decreased while the fastest winds had migrated toward low-southern latitudes.

S Vinatier↗

Interaction of Saturn’s Hexagon with Convective Storms

In March 2020 a convective storm erupted at planetographic latitude 76°N in the southern flank of Saturn’s long-lived hexagonal wave. The storm reached a zonal size of 4,500 km and developed a tail extending zonally 33,000 km. Two new short-lived storms erupted in May in the hexagon edge. These storms formed after the convective storms that took place in 2018 in nearby latitudes. There were no noticeable changes in the zonal profile of Saturn's polar winds in 2018-2020. Measurements of the longitude position of the vertices of the hexagon throughout this period yield a value for its period of rotation equal to that of System III of radio-rotation measured at the time of Voyagers. We report changes in the hexagon clouds related to the activity of the storms. Our study reinforces the idea that Saturn’s hexagon is a well rooted structure with a possible direct relationship with the bulk rotation of the planet.

Convective storms↗

The Salience of Nonlinearities in the Boreal Winter Response to ENSO: Arctic Stratosphere and Europe

The Arctic stratospheric response to El Niño (EN) and La Niña (LN) is evaluated in a 41 member ensemble of the period 1980 to 2009 in the Goddard Earth Observing System Chemistry-Climate Model. We consider whether the responses to EN and LN are equal in magnitude and opposite in sign, whether the responses to moderate and extreme events are proportion-ate, and if the response depends on whether sea surface temperature anomalies (SSTs) peak in the Eastern Pacific (EP) or Central Pacific (CP). There is no indication of any nonlinearities between EN and LN, though in ~ 15% of the ensemble members the stratospheric sudden warming (SSW) frequencies for EN and LN are similar, suggesting that a similar SSW frequency for EN and LN, as has occurred over the past ~ 60 years, can occur by chance. The response to extreme EN events is not proportionate to the amplitude of the underlying SST anomalies in spring. EP EN events preferentially increase zonal wavenumber 1 and decrease zonal wavenumber 2 as compared to CP EN events, however the zonal-mean Arctic stratospheric and subpolar surface response is generally little different between EP EN and CP EN once one accounts for the relative weakness of CP events. These differences between EP and CP events and between moderate and extreme EN events only emerge if at least 25 events are composited, however, due to the small signal-to-noise ratio, and hence these differences may be of little practical benefit.

Sea surface temperature anomalies (SSTs) peak in t↗

Intercomparison of Middle Atmospheric Meteorological Analyses for the Northern Hemisphere Winter 2009-2010

Detailed meteorological analyses based on observations extending through the middle atmosphere (~15 to 100 km altitude) can provide key information to whole atmosphere modelling systems regarding the physical mechanisms linking day-to-day changes in ionospheric electron density to meteorological variability near the Earth’s surface. However, the extent to which independent middle atmosphere analyses differ in their representation of wave-induced coupling to the ionosphere is unclear. To begin to address this issue, we present the first intercomparison among four such analyses, JAGUAR-DAS, MERRA-2, NAVGEM-HA, and WACCMX+DART, focusing on the Northern Hemisphere (NH) 2009-2010 winter, which includes a major sudden stratospheric warming (SSW). This intercomparison examines the altitude, latitude, and time dependences of zonal mean zonal winds and temperatures among these four analyses over the 1 December 2009 – 31 March 2010 period, as well as latitude and altitude dependences of monthly mean amplitudes of the diurnal and semidiurnal migrating solar tides, the eastward propagating diurnal zonal wave number 3 nonmigrating tide, and traveling planetary waves associated with the quasi-5-day and quasi-2-day Rossby modes. Our results show generally good agreement among the four analyses up to the stratopause (~50 km altitude). Large discrepancies begin to emerge in the mesosphere and lower thermosphere owing to (1) differences in the types of satellite data assimilated by each system and (2) differences in the details of the global atmospheric models used by each analysis system. The results of this intercomparison provide initial estimates of uncertainty in analyses commonly used to constrain middle atmospheric meteorological variability in whole atmosphere model simulations.

John P Mccormack↗

Key Observations and Discoveries from the Vector Electric Field Investigation on the C/NOFS Satellite

Launched in 2008 and operating for 7.5 years, the Air Force Communication /Navigation Outage Forecasting System (C/NOFS) satellite included the Vector Electric Field Investigation (VEFI) designed and built at NASA’s Goddard Space Flight Center. VEFI successfully met its objectives: 1) investigate the role of ambient electric fields in initiating nighttime density depletions and turbulence; 2) determine the quasi-DC electric fields associated with abrupt density depletions, and 3) quantify the spectrum of the irregularities associated with density depletions, providing many key observations and discoveries including: -- Global (low latitude) vector DC electric fields at 16 s/sec revealing variations with local time, longitude, and season between extremely low and moderate solar activity-- Reversed zonal E x B drifts below the F-region ledge at sunset and simultaneous observations of large scale Kelvin-Helmholtz instabilities as seeds of spread-F -- Large expanses of quasi-coherent kilometer-scale vector wave observations (electric field and density) below 450 km and their discovery as a source of scintillations using the C/NOFS GPS -- Reversed zonal DC electric fields and simultaneous observations of afternoon counter electrojet-- Enhanced zonal DC electric fields at sunrise and full vector plasma flow continuity at the terminator-- Vector electric field and density irregularities extending to meter-scales within equatorial plasma depletions -- Intense electric field structure within the equatorial ionosphere at night in the absence of density depletions -- First observations of the westward equatorial electrojet in post-midnight ionosphere and possible association with downward meridional drifts-- Spaceborne Dst observations and implications for ring current asymmetries-- Vector electric and magnetic field Poynting flux within depletions and TIDs-- Measurements of ionospheric reflectance and Alfvenic waves associated with TIDs-- Observations of Alfven resonators -- ULF magnetic field structure within the nightime equatorial ionosphere-- Discovery of Schumman resonances in space and implications of ELF radiowave propagation-- Vector observations of 50-60 Hz powerline radiation without harmonics-- Ion cyclotron resonance absorption lines associated with ELF hiss and the identification of ambient ions-- Parallel electric fields associated with lower hybrid waves driven by thunderstorm lightning-related sferics-- Discovery of Z-mode radiation in the equatorial ionosphere including its associated with density depletions-- Optical lightning detector waveforms and electric field sferics observed up to altitudes of 800 km

Robert F Pfaff↗

Ecosystem age-class dynamics and distribution in the LPJ-wsl v2.0 global ecosystem model

Forest ecosystem processes follow classic responses with age, peaking production around canopy closure and declining thereafter. Although age dynamics might be more dominant in certain regions over others, demographic effects on net primary production (NPP) and heterotrophic respiration (Rh) are bound to exist. Yet, explicit representation of ecosystem demography is notably absent in many global ecosystem models. This is concerning because the global community relies on these models to regularly update our collective understanding of the global carbon cycle. This paper aims to present the technical developments of a computationally efficient approach for representing age-class dynamics within a global ecosystem model, the Lund–Potsdam–Jena – Wald, Schnee, Landschaft version 2.0 (LPJ-wsl v2.0) dynamic global vegetation model and to determine if explicit representation of demography influenced ecosystem stocks and fluxes at global scales or at the level of a grid cell. The modeled age classes are initially created by simulated fire and prescribed wood harvesting or abandonment of managed land, otherwise aging naturally until an additional disturbance is simulated or prescribed. In this paper, we show that the age module can capture classic demographic patterns in stem density and tree height compared to inventory data, and that simulated patterns of ecosystem function follow classic responses with age. We also present two scientific applications of the model to assess the modeled age-class distribution over time and to determine the demographic effect on ecosystem fluxes relative to climate. Simulations show that, between 1860 and 2016, zonal age distribution on Earth was driven predominately by fire, causing a 45- to 60-year difference in ages between older boreal (50–90° N) and younger tropical (23° S–23° N) ecosystems. Between simulation years 1860 and 2016, land-use change and land management were responsible for a decrease in zonal age by −6 years in boreal and by −21 years in both temperate (23–50° N) and tropical latitudes, with the anthropogenic effect on zonal age distribution increasing over time. A statistical model helped to reduce LPJ-wsl v2.0 complexity by predicting per-grid-cell annual NPP and Rh fluxes by three terms: precipitation, temperature, and age class; at global scales, R2 was between 0.95 and 0.98. As determined by the statistical model, the demographic effect on ecosystem function was often less than 0.10 kg C/sq. myr but as high as 0.60 kg C/sq. myr where the effect was greatest. In the eastern forests of North America, the simulated demographic effect was of similar magnitude, or greater than, the effects of climate; simulated demographic effects were similarly important in large regions of every vegetated continent. Simulated spatial datasets are provided for global ecosystem ages and the estimated coefficients for effects of precipitation, temperature and demography on ecosystem function. The discussion focuses on our finding of an increasing role of demography in the global carbon cycle, the effect of demography on relaxation times (resilience) following a disturbance event and its implications at global scales, and a finding of a 40 Pg C increase in biomass turnover when including age dynamics at global scales. Whereas time is the only mechanism that increases ecosystem age, any additional disturbance not explicitly modeled will decrease age. The LPJ-wsl v2.0 age module represents another step forward towards understanding the role of demography in global ecosystems.

LPJ-wsl↗

Low-Level Jets and the Convergence of Mars Data Assimilation Algorithms

Data assimilation is an increasingly popular technique in Mars atmospheric science, but its effect on the mean states of the underlying atmosphere models has not been thoroughly examined. The robustness of results to the choice of model and assimilation algorithm also warrants further study. We investigate these issues using two Mars general circulation models (MGCMs), with particular emphasis on zonal wind and temperature fields. When temperature retrievals from the Mars Global Surveyor Thermal Emission Spectrometer (TES) are assimilated into the U.K.-Laboratoire de M ́et ́eorologie Dynamique (UK-LMD) MGCM to create the Mars Analysis Correction Data Assimilation (MACDA) reanalysis, low-level zonal jets in the winter northern hemisphere shift equatorward and weaken relative to a free-running control simulation from the same MGCM. The Ensemble Mars Atmosphere Reanalysis System (EMARS) reanalysis, which is also based on TES temperature retrievals, also shows jet weakening (but less if any shifting) relative to a control simulation performed with the underlying Geophysical Fluid Dynamics Laboratory (GFDL) MGCM. Examining higher levels of the atmosphere, monthly mean three-dimensional temperature and zonal wind fields are in generally better agreement between the two reanalyses than between the two control simulations. In conjunction with information about the MGCMs’ physical parametrizations, intercomparisons between the various reanalyses and control simulations suggest that overall the EMARS control run is plausibly less biased (relative to the true state of the Martian atmosphere) than the MACDA control run. Implications for future observational studies are discussed. Plain Language Summary An increasingly popular way to study Martian weather and climate is to combine atmospheric temperature observations with a computer model (specifically, a global climate model). The process of combining model and observations is called “data assimilation”, and the resulting merged data set is called a “reanalysis”. One advantage of reanalyses is that they include variables (such as wind) that are not directly observed. For scientific and practical applications we want these variables to be reasonably accurate— however, it is not clear how well data assimilation algorithms compute them. Our study investigates this issue using two Mars reanalyses and two model simulations that do not assimilate temperature data. We focus on slowly-varying atmospheric phenomena (timescales from 10 Mars days to a season). Assimilating temperature data into two different global climate models changes the strength and/or spatial pattern of east-west winds at low altitudes. Furthermore, monthly mean three-dimensional temperature and east-west wind fields agree better between reanalyses than between non-assimilating model simulations. This suggests that the data assimilation process is basically successful. One non-assimilating model simulation has less realistic representations of atmospheric physical processes than the other—we argue that this plausibly gives it larger biases relative to the true state of the atmosphere.

Todd A Mooring↗

Global Total Ozone Recovery Trends Attributed to Ozone-Depleting Substance (ODS) Changes Derived From Five Merged Ozone Datasets

We report on updated trends using different merged zonal mean total ozone datasets from satellite and ground-based observations for the period from 1979 to 2020. This work is an update of the trends reported in Weber et al. (2018) using the same datasets up to 2016. Merged datasets used in this study include NASA MOD v8.7 and NOAA Cohesive Data (COH) v8.6, both based on data from the series of Solar Backscatter Ultraviolet (SBUV), SBUV-2, and Ozone Mapping and Profiler Suite (OMPS) satellite instruments (1978–present), as well as the Global Ozone Monitoring Experiment (GOME)-type Total Ozone – Essential Climate Variable (GTO-ECV) and GOME-SCIAMACHY-GOME-2 (GSG) merged datasets (both 1995–present), mainly comprising satellite data from GOME, SCIAMACHY, OMI, GOME-2A, GOME-2B, and TROPOMI. The fifth dataset consists of the annual mean zonal mean data from ground-based measurements collected at the World Ozone and Ultraviolet Radiation Data Centre (WOUDC). Trends were determined by applying a multiple linear regression (MLR) to annual mean zonal mean data. The addition of 4 more years consolidated the fact that total ozone is indeed slowly recovering in both hemispheres as a result of phasing out ozone-depleting substances (ODSs) as mandated by the Montreal Protocol. The near-global (60° S–60° N) ODS-related ozone trend of the median of all datasets after 1995 was 0.4 ± 0.2 (2σ) %/decade, which is roughly a third of the decreasing rate of 1.5 ± 0.6 %/decade from 1978 until 1995. The ratio of decline and increase is nearly identical to that of the EESC (equivalent effective stratospheric chlorine or stratospheric halogen) change rates before and after 1995, confirming the success of the Montreal Protocol. The observed total ozone time series are also in very good agreement with the median of 17 chemistry climate models from CCMI-1 (Chemistry-Climate Model Initiative Phase 1) with current ODS and GHG (greenhouse gas) scenarios (REF-C2 scenario). The positive ODS-related trends in the Northern Hemisphere (NH) after 1995 are only obtained with a sufficient number of terms in the MLR accounting properly for dynamical ozone changes (Brewer–Dobson circulation, Arctic Oscillation (AO), and Antarctic Oscillation (AAO)). A standard MLR (limited to solar, Quasi-Biennial Oscillation (QBO), volcanic, and El Niño–Southern Oscillation (ENSO)) leads to zero trends, showing that the small positive ODS-related trends have been balanced by negative trend contributions from atmospheric dynamics, resulting in nearly constant total ozone levels since 2000.

Total Column Ozone Trends↗

Diurnal Waves Forced By Horizontal Convergence of Near-Surface Winds on Mars

Boundary-layer winds near steeply sloped topography excite a distinct type of atmospheric wave. Basic properties of the waves were determined through analysis of 7 Mars years of temperature profiles from the Mars Climate Sounder on Mars Reconnaissance Orbiter. A more complete understanding of the waves was derived from a simulation with the NASA Ames Mars General Circulation Model. This type of wave appears in the tropics in the vicinity of conspicuous surface features, such as large volcanoes and the rim of Isidis Planitia. The amplitude is typically 5–10 K, the vertical wavelength is about 40 km, the zonal wavelength is about 1400 km, and the zonal phase speed relative to the surface is about 25 m s−1. The wave forcing comes from horizontal convergence of the near-surface winds, which generates a daily surge in the vertical winds. These convergence zones were not previously recognized as an important source of atmospheric waves on Mars. Strong forcing occurs only within a limited range of local time, producing wave packets (diurnal pulses) rather than steady oscillations. Waves are present throughout the year and their properties vary with season. Near the solstices, the wave amplitude is largest in the tropics of the summer hemisphere, where the convergence zones are stronger than in the winter hemisphere. The diurnal-mean zonal winds, which vary with latitude and season, control the direction of propagation of the predominant waves through filtering by critical levels.

Mars↗

Grid Refinement Techniques for the 𝜸-Re𝜽𝒕 Transition Model in FUN3D

There has been an increased focus on the overall accuracy and grid convergence of Reynolds-averaged Navier-Stokes (RANS)-based transition models from the recent AIAA and NATO-AVT workshops. Even though satisfactory grid convergence could be achieved for simple two-dimensional flow configurations, it required mesh counts that are substantially larger than those used in typical applications. In this paper, we focus our efforts on understanding how the grid resolution and topology influences the accuracy and convergence of results by studying the Schubauer and Skramstad flat-plate configuration and the NLF-0416 airfoil at an angle of attack equal to five degrees using FUN3D, a second-order finite-volume code. By focusing on these cases, we can analyze both natural and separation-induced transition scenarios. Multiple grid refinement strategies are investigated. First, we determine the relative effectiveness of zonal streamwise refinement in the transition region within structured grids as an alternative to the costly option of globally uniform refinement of a baseline grid. We also complement this zonal technique by globally varying the wall-normal resolution keeping the streamwise resolution fixed. The zonal streamwise refinement can accurately model natural transition in a flat-plate boundary layer and separation-induced transition, but struggles to accurately model natural transition in airfoil flows. A series of unstructured prismatic grids that have similar node counts and viscous wall spacings as the structured hexahedral grids are also tested, and they do not achieve grid convergence until an extremely fine resolution. Last, we employ adjoint-based unstructured grid adaptation in FUN3D to natural and separation-induced transition on the NLF-0416 airfoil. The adjoint-based refinement process converges to the same solution as the baseline family of structured grids, but leads to smaller errors on coarser grids.

Grid Adaptation↗

Grid Refinement Techniques for the 𝜸-Re 𝜽 𝒕 Transition Model in FUN3D

There has been an increased focus on the overall accuracy and grid convergence of Reynolds-averaged Navier-Stokes (RANS)-based transition models from the recent AIAA and NATO-AVT workshops. Even though satisfactory grid convergence could be achieved for simple two-dimensional flow configurations, it required mesh counts that are substantially larger than those used in typical applications. In this paper, we focus our efforts on understanding how the grid resolution and topology influences the accuracy and convergence of results by studying the Schubauer and Skramstad flat-plate configuration and the NLF-0416 airfoil at an angle of attack equal to five degrees using FUN3D, a second-order finite-volume code. By focusing on these cases, we can analyze both natural and separation-induced transition scenarios. Multiple grid refinement strategies are investigated. First, we determine the relative effectiveness of zonal streamwise refinement in the transition region within structured grids as an alternative to the costly option of globally uniform refinement of a baseline grid. We also complement this zonal technique by globally varying the wall-normal resolution keeping the streamwise resolution fixed. The zonal streamwise refinement can accurately model natural transition in a flat-plate boundary layer and separation-induced transition, but struggles to accurately model natural transition in airfoil flows. A series of unstructured prismatic grids that have similar node counts and viscous wall spacings as the structured hexahedral grids are also tested, and they do not achieve grid convergence until an extremely fine resolution. Last, we employ adjoint-based unstructured grid adaptation in FUN3D to natural and separation-induced transition on the NLF-0416 airfoil. The adjoint-based refinement process converges to the same solution as the baseline family of structured grids, but leads to smaller errors on coarser grids.

Grid Adaptation↗

Three-Dimensional Atmospheric Dynamics of Jupiter from Ground-Based Doppler Imaging Spectroscopy in the Visible

We present three-dimensional (3D) maps of Jupiter’s atmospheric circulation at cloud-top level from Doppler-imaging data obtained in the visible domain with JIVE, the second node of the JOVIAL network, which is mounted on the Dunn Solar Telescope at Sunspot, New Mexico. We report on 12 nights of observations between May 4 and May 30, 2018, representing a total of about 80 hours. Firstly, the average zonal wind profile derived from our data is compatible with that derived from cloud-tracking measurements performed on Hubble Space Telescope images obtained in April 2018 from the Outer Planet Atmospheres Legacy (OPAL) program. Secondly, we present the first ever two-dimensional maps of Jupiter’s atmospheric circulation from Doppler measurements. The zonal velocity map highlights well-known atmospheric features, such as the equatorial hot spots and the Great Red Spot (GRS). In addition to zonal winds, we derive meridional and vertical velocity fields from the Doppler data. The motions attributed to vertical flows are mainly located at the boundary between the equatorial belts and tropical zones, which could indicate active motion in these regions. Qualitatively, these results compare well to recent Juno data that have unveiled the three-dimensional structure of Jupiter’s wind field. To the contrary, the motions attributed to meridional circulation are very different from what is obtained by cloud tracking, except at the GRS. Because of limitations with data resolution and processing techniques, we acknowledge that our measurement of vertical or meridional flows of Jupiter are still to be confirmed.

Francois-Xavier Schmider↗

Evidence for Coupling Between the Subseasonal Oscillations in the Southern Hemisphere Midlatitude Ocean and Atmosphere

Abstract A regional perspective is taken to investigate the intriguing ∼25‐day periodicity in the Southern Hemisphere (SH) baroclinic annular mode (BAM) and the possible underlying mechanisms. Results show that the 25‐day periodicity exists not only in the zonal average but also in the regional average of eddy activity as long as the domain of the average is wide enough to accommodate a wave packet. The favorable domain for the ∼25 days periodicity is located downstream of the largest baroclinicity in midlatitudes. The wave activity developed there can travel circumglobally several times with a typical group velocity of 25 m s −1 , and reenter the domain every 12.5 days, setting the base periodicity. A slower oceanic oscillation at a time scale of ∼100 days in the same region appears to select and reinforce the multiples of 12.5 days, giving rise to a pronounced ∼25‐day periodicity in the wave activity.

Xue, Daokai↗

Learning how structures form in drift-wave turbulence

Drift-wave turbulence produces anomalous transport via cross-correlations between fluctuations. This transport has profound implications for confinement, structure formation, and virtually all aspects of the non-linear turbulent dynamics. Herein, we use a data-driven method based on deep learning in order to study turbulent transport in the 2D Hasegawa–Wakatani system and infer a reduced mean-field model from numerical solution. In addition to the usual turbulent diffusion, we find an effect which couples the particle flux to the local gradient of vorticity, which tends to modulate the density profile. The direct coupling to the shear is relatively weak. In addition, the deep learning method finds a model for spontaneous zonal flow generation by negative viscosity, stabilized by non-linear and hyperviscous terms. We compare these results to analytic calculations using quasilinear theory and wave kinetics, finding qualitative agreement, though the calculations miss certain higher-order effects. A simplified, 1-D model for the evolution of the profile, flow, and intensity based on the deep learning results is solved numerically and compared to previous models for staircasing based on bistability. We see that the physics uncovered by the deep learning method provided simple explanations for the formation of zonal structures in the density, flow, and turbulence fields. We highlight the important role of symmetry in the deep learning method and speculate on the portability of the method to other applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of resonant magnetic perturbations on radial electric fields in DIII-D tokamak

Here, gyrokinetic simulations of DIII-D tokamak equilibrium find that resonant magnetic perturbation (RMP) drives a neoclassical non-ambipolar electron particle flux, which causes a rapid change of equilibrium radial electric fields consistent with experimental observations during the suppression of the edge localized mode (ELM). The simulation results provide a support for the conjecture that RMP-induced changes of radial electric fields lead to the enhanced turbulent transport at the pedestal top during the ELM suppression (Taimourzadeh et al 2019 Nucl. Fusion 59 046005). Furthermore, gyrokinetic simulations of collisionless damping of zonal flows show that resonant responses to the RMP decrease the residual level of the zonal flows and damp the geodesic acoustic mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theory of the tertiary instability and the Dimits shift from reduced drift-wave models

Tertiary modes in electrostatic drift-wave turbulence are localized near extrema of the zonal velocity $U(x)$ with respect to the radial coordinate $x$. We argue that these modes can be described as quantum harmonic oscillators with complex frequencies, so their spectrum can be readily calculated. The corresponding growth rate $\gamma_{\rm TI}$ is derived within the modified Hasegawa--Wakatani model. We show that $\gamma_{\rm TI}$ equals the primary-instability growth rate plus a term that depends on the local $U''$; hence, the instability threshold is shifted compared to that in homogeneous turbulence. This provides a generic explanation of the well-known yet elusive Dimits shift, which we find explicitly in the Terry--Horton limit. Linearly unstable tertiary modes either saturate due to the evolution of the zonal density or generate radially propagating structures when the shear $|U'|$ is sufficiently weakened by viscosity. The Dimits regime ends when such structures are generated continuously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

(abstract) Role of Long Equatorial Wave Reflection in the Low-Frequency Variability Observed During the TOPEX/POSEIDON Period: Data Analysis and Coupled Model Study

The analysis of TOPEX/POSEIDON sea level and ERS-1 zonal wind stress data gives evidence that first mode meridonal long Rossby waves fully reflected into Kelvin waves at the equatorial western Pacific boundary at different periods including the January-June 1994 period. The evolution of the conditions (zonal wind stress, sea surface temperature and sea level anomalies) in the central Pacific suggests that this reflection may have played a role in the reversal of cold to warm anomalies observed in the central Pacific in June-July 1994. To investigate the actual role of reflected Kelvin waves observed by TOPEX/POSEIDON, a simple ocean-atmosphere coupled model is designed. Starting from these initial conditions, several simulations are run for testing the role western bundary reflection evidenced in TOPEX/POSEIDON data. Coupled mechanisms are suggested to explain the low-frequency variability observed during this period.

Rossby waves Kelvin waves zonal wind stress TOPEX/↗

Advantages of a variable‐resolution global climate model in reproducing the seasonal evolution of East Asian summer monsoon

Abstract The East Asian summer monsoon (EASM) is unique among monsoon systems that it features meridional evolution of the summer monsoon. In this study, we evaluate the performances of a Variable‐Resolution Community Earth System Model (VR‐CESM) regionally refined over eastern China (14 km) in reproducing the seasonal evolution of EASM precipitation over China. Compared with reference datasets, VR‐CESM shows better performance than the corresponding globally uniform coarse‐resolution model CESM (quasi‐uniform 1°), especially over western China where complex local topography exists. The northward monsoon migration is closely related to low‐level southerly flows and vertical moisture advection, which are more reasonably simulated in VR‐CESM. The four critical timings of the EASM (monsoon onset, withdrawal, peak, and duration) are also better captured in VR‐CESM than in CESM. The corresponding spatial Pearson correlation coefficients of the four critical timings with respect to reference datasets are about 0.1 higher in VR‐CESM than those in CESM. Both models are most accurate in simulating monsoon onset and least accurate at simulating the monsoon peak. The overestimated zonal thermal contrast in CESM is responsible for the earlier monsoon onset and excessive precipitation in September over the Yangtze River valley. Finer resolution in VR‐CESM, especially over the Tibetan Plateau (TP), appears to be a main factor in simulating better zonal thermal contrast and seasonal evolution of the EASM.

54 ENVIRONMENTAL SCIENCES↗