Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Source and propagation modeling”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

CODEX Optical Design and Alignment

The COronal Diagnostic EXperiment (CODEX) is a Heliophysics mission to measure the density, temperature, and velocity of the electrons in the solar corona with the primary goal of improving our understanding of the physical conditions of the solar wind in the acceleration region. The temperature and velocity measurement requires much higher signal-to-noise ratio than the density measurements. In solar coronagraphs, the diffraction of the solar disk light due to the occulting element is the dominant source of noise. Therefore, to further suppress the diffracted sun light with respect to the existing coronagraphs is a critical element of the CODEX design. To minimize the stray light due to diffraction, the selected optical design is a two-stage standard coronagraph with an external occulter, an internal occulter, and a Lyot stop. What is unique for this design is that a focal mask was inserted at the telescope focal plane. It works together with the field lens suppressing the stray light down by ~ another order of magnitude as compared to a traditional three-stage approach. During the optical design, a Fourier Transform based beam propagation software, i.e., GLAD, was used to model the beam path through the full coronagraph, from the external occulter to the detector array. All diffraction sensitive elements: external occulter, internal occulter, focal mask, and Lyot stop were carefully modeled and optimized. As a result, the requirement of achieving a stray light level which is one order of magnitude lower than F-corona was satisfied. On the other hand, to achieve the final suppression, a precision optical alignment is another must. This paper also presents our creative alignment procedure: using the combination of metrology, precision alignment equipment, and real time diffraction ring monitoring to minimize the diffraction. The final test results show that the suppression ratio (B/B 0 ) reaches 10 -11 level, which is equivalent to one order of magnitude lower than F-corona.

CODEX optical design and alignment↗
DOE OSTI2021

A temperature profile diagnostic for radiation waves on OMEGA-60

Predicting and matching radiation wave propagation with computational models has proven difficult. Information provided by experiments studying radiation flow has been limited when only radiation breakout is measured. We have developed the COAX (co-axial) diagnostic platform to provide spatial temperature profiles of a radiation wave through low density foams as a more detailed constraint for simulations. COAX uses a standard, laser-driven OMEGA-60 halfraum to drive radiation down a titanium-laden silicon oxide foam. Point-projection X-ray absorption spectroscopy perpendicular to the radiation flow measures the spatial profile of titanium ionization. The spectroscopic measurement utilizes a broadband capsule backlighter. Imaging and streak spectroscopy are used to characterize the size and spectrum of this source. Radiography provides an additional constraint by capturing the developing shock as the radiation flow becomes subsonic. The DANTE diagnostic is used to measure the halfraum temperature. Here, we provide a spectroscopic analysis of COAX data to determine temperature, and we describe experimental sources of uncertainty. The temperature is obtained by comparison to multi-temperature synthetic spectra post-processed from radiation-hydrodynamics simulations. Quantitative comparison between data and synthetic spectra generated from temperature profiles at relevant simulation times enable determination of a peak temperature of 114 ± 8 eV at 265 ± 22.4 μm from the halfraum. This represents an improvement over the temperature uncertainties of previous radiation flow experiments. Further refinements to the spectroscopic analysis could achieve ±4 eV. The combination between space-resolved spectroscopy and radiography enables us to determine the distance from the halfraum of both the radiation front and the shock front at the time of measurement. For the example shown in this paper the radiation front position is 600–630 μm at 3.43 ± 0.16 ns and the shock front position is 633 μm at 3.3 ± 0.24 ns.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1977

Experiments on whistler wave filamentation and VLF hiss in a laboratory plasma

With the development of a large magnetized plasma source it has become possible to investigate space plasma physics problems in the laboratory. First, the nonlinear effects associated with the excitation of a large amplitude whistler wave have been explored. It is found that the radiation pressure of the wave and thermal effects give rise to a field-aligned density depression in which the wave becomes completely trapped. Hyperfine filaments with diameters small compared with the parallel wavelength are observed. Second, the stability of oblique whistler waves in the presence of an electron beam has been studied. A broadband whistler instability is observed and identified as a Cherenkov interaction between beam electrons and whistlers propagating near the resonance cone. These observations confirm the present model for the generation of VLF hiss in the aurora.

Stenzel, R. L.↗
NASA NTRS1987

Two dimensional imaging observations of meter-decameter bursts associated with the February 1986 flare activity

An analysis is presented of the two dimensional imaging observations of a flare observed on 3 Feb. l986 using the Clark Lake Multifrequency Radioheliograph. The flare produced almost all types of Meter-decimeter radio emission: enhanced storm radiation, type III/V bursts, II and IV and flare continuum. The flare continuum had early (FCE) and late (FC II) components and the type II occurred during the period between these two components. Comparing the source positions of type III/V and FCE it was found that these bursts must have occurred along adjacent open and closed field lines, respectively. The positional analysis of type II and FC II implies that the nonthermal electrons responsible for FC II need not be accelerated by type II shock and this conclusion is further supported by the close association of FC II with a microwave peak. Using the positional and temporal analysis of all these bursts and the associated hard X-ray and microwave emissions, a schematic model is developed for the magnetic field configuration in the flaring region in which the nonthermal particles responsible for these bursts are confined or along which they propagate.

Gopalswamy, N.↗
NASA NTRS2011

Observation of Schumann Resonances in the Earth's Ionosphere

The surface of the Earth and the lower edge of the ionosphere define a cavity in which electromagnetic waves propagate. When the cavity is excited by broadband electromagnetic sources, e.g., lightning, a resonant state can develop provided the average equatorial circumference is approximately equal to an integral number of wavelengths of the electromagnetic waves. This phenomenon, known as Schumann resonance, corresponds to electromagnetic oscillations of the surface-ionosphere cavity, and has been used extensively to investigate atmospheric electricity. Using measurements from the Communications/Navigation Outage Forecasting System (C/NOFS) satellite, we report, for the first time, Schumann resonance signatures detected well beyond the upper boundary of the cavity. These results offer new means for investigating atmospheric electricity, tropospheric-ionospheric coupling mechanisms related to lightning activity, and wave propagation in the ionosphere. The detection of Schumann resonances in the ionosphere calls for revisions to the existing models of extremely low frequency wave propagation in the surface-ionosphere cavity. Additionally, these measurements suggest new remote sensing capabilities for investigating atmospheric electricity at other planets.

Simoes, Fernando↗
NASA NTRSDate not supplied

Flux Comparison of Master-8 and Ordem 3.1 Modelled Space Debris Population

With ESA’s Meteoroid And Space debris Terrestrial Environment Reference (MASTER-8) model and NASA’s Orbital Debris Engineering Model(ORDEM) 3.1, the two premier orbital debris engineering models have been officially released. The two models come with significant enhancements and now represent state-of-the-art orbital debris modelling for their respective agencies. Both models provide the community with estimates of the space debris environment from low Earth orbit (LEO) up to at least geostationary altitude. The MASTER population is an event-based simulation of all known events that generate debris and objects that are part of the U.S. Space Surveillance Network (SSN) catalog, which provides coverage of objects with diameters down to approximately 10cm in LEO and 1 m in geosynchronous Earth orbit (GEO). Different models are used to simulate the artificial objects and their orbit evolution over time. These models are called “sources” since they assign an origin to each individual object and consist of fragments, solid rocket motor (SRM) remainders, sodium-potassium(NaK) droplets, paintflakes, ejecta, and multi-layer insulation (MLI) fragments. The objects from each source are characterized by having individual release mechanisms, as well as orbital distributions, material composition, size, and mass distributions. Dedicated radar and telescope observation data is used to calibrate the model for objects larger than 1 cm in LEO and larger than 10 cm in GEO. For calibrating the small-sized objects, below 1 cm, impact data from returned surfaces are analyzed. Because the >1cm object population is dominated by fragments, the fragmentation event database was updated to include new events, as well as re-evaluate past events. Special attention was drawn to re-evaluating theFengyun-1C anti-satellite test from 2007 and Cosmos-Iridium collision event from 2009 since these events shape the fragment population because of their severity. After 2009, the two largest fragmentations in terms of number of tracked debris are the Briz-M explosion in 2012 and the NOAA-16explosion in 2015. In total, there are 261 confirmed fragmentations in the database up to November 2016.The baseline population for ORDEM 3.1 is based on the U.S. SSN catalog, and observational datasets from radar, in situ, and optical sources provide a foundation from which the model populations are statistically extrapolated to smaller size regions. These regions are not well-covered by the SSN catalog yet may pose the greatest threat to operational spacecraft. The NASA Standard Satellite Breakup Model is used to generate fragments greater than 1mm from collisions and explosions, and these fragment populations are scaled using ground-based radar data. Specific major debris-producing events, including the Fengyun-1C, Iridium 33, and Cosmos2251 debris clouds, and unique populations, such as NaK droplets, were re-examined, modelled, and added to the ORDEM environment separately. Optical measurement data is used to model the GEO population down to 10 cm. The debris environment is propagated using NASA’s LEO-to-GEO Environment Debris model, and future explosions of intact objects and collisions involving objects greater than 10 cm are assessed statistically. The environment from a few millimetres down to 10 𝜇m is modelled using a special degradation model where small particles are generated from intact spacecraft and rocket bodies, then the populations are scaled to fit in situ cratering data from Space Shuttle returned surfaces. Fragments smaller than 10 cm are differentiated based on material density categories, i.e., high-, medium-, and low-density, to better characterize the potential debris risk posed to upper stages and spacecraft. This paper will discuss the MASTER and ORDEM approaches for modelling populations and compare fluxes for specific orbits, including sun-synchronous, ISS-altitude, geosynchronous transfer, and GEO. In the end, a conclusion is drawn towards the importance of having multiple fundamentally different, yet validated, models to estimate the space debris population.

Andre Horstmann↗
DOE OSTI2025

A Numerical Model for the Dynamics of Pickup Ions Outside the Heliopause and IBEX “Ribbon” Observation

The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1975

Spinning mode sound propagation in ducts with acoustic treatment

Recent acoustic data show larger noise attenuations than predicted for acoustically treated aircraft engine inlets without splitter rings. A theoretical study of the acoustic propagation of spinning modes in acoustically treated open circular ducts is presented, and a suppressor with splitter rings was modeled by using the rectangular approximation to the annular duct. Theoretical models were used to determine optimum impedance and maximum attenuation for several spinning lobe numbers from 0 to 50. Results of the analysis indicate that for circular ducts the maximum possible attenuation and the optimum wall impedance are strong functions of the lobe number, and for annular ducts the attenuation and optimum wall impedance are insensitive to the spinning lobe number for well cut-on modes. These results explain why suppressors with splitter rings were quite effective in spite of the lack of detailed information on the noise source modal structure. Conversely, effective use of outer wall treatment alone will require expanded knowledge of the noise source structure. Approximate solutions are presented to help interpret the more exact theoretical results.

Rice, E. J.↗
NASA NTRS1974

Spinning mode sound propagation in ducts with acoustic treatment

Recent acoustic data have shown larger noise attenuations than predicted for acoustically treated aircraft engine inlets without splitter rings. These data have stimulated a more detailed theoretical study of the acoustic propagation of spinning modes in acoustically treated open circular ducts. In addition, the suppressor with splitter rings was modeled by using the rectangular approximation to the annular duct. The theoretical models were used to determine optimum impedance and maximum attenuation for several spinning lobe numbers from 0 to 50. It is found that for circular ducts the maximum possible attenuation and the optimum wall impedance are strong functions of the lobe number. For annular ducts the attenuation and optimum wall impedance are insensitive to the spinning lobe number for well cut-on modes. The results help explain why suppressors with splitter rings have been quite effective in spite of the lack of detailed information on the noise-source modal structure. Conversely, effective use of outer-wall treatment alone will require expanded knowledge of the noise-source structure. Approximate solutions are presented to help interpret the more exact theoretical results.

Rice, E. J.↗
NASA NTRS1988

An analysis of the Kalman filter in the Gamma Ray Observatory (GRO) onboard attitude determination subsystem

The Gamma Ray Observatory (GRO) spacecraft needs a highly accurate attitude knowledge to achieve its mission objectives. Utilizing the fixed-head star trackers (FHSTs) for observations and gyroscopes for attitude propagation, the discrete Kalman Filter processes the attitude data to obtain an onboard accuracy of 86 arc seconds (3 sigma). A combination of linear analysis and simulations using the GRO Software Simulator (GROSS) are employed to investigate the Kalman filter for stability and the effects of corrupted observations (misalignment, noise), incomplete dynamic modeling, and nonlinear errors on Kalman filter. In the simulations, on-board attitude is compared with true attitude, the sensitivity of attitude error to model errors is graphed, and a statistical analysis is performed on the residuals of the Kalman Filter. In this paper, the modeling and sensor errors that degrade the Kalman filter solution beyond mission requirements are studied, and methods are offered to identify the source of these errors.

Snow, Frank↗
NASA NTRS2008

Modeling the Diffuse Cloud-Top Optical Emissions from Ground and Cloud Flashes

A number of studies have indicated that the diffuse cloud-top optical emissions from intra-cloud (IC) lightning are brighter than that from normal negative cloud-to-ground (CG) lightning, and hence would be easier to detect from a space-based sensor. The primary reason provided to substantiate this claim has been that the IC is at a higher altitude within the cloud and therefore is less obscured by the cloud multiple scattering medium. CGs at lower altitudes embedded deep within the cloud are more obscured, so CG detection is thought to be more difficult. However, other authors claim that because the CG source current (and hence luminosity) is typically substantially larger than IC currents, the greater CG source luminosity is large enough to overcome the effects of multiple scattering. These investigators suggest that the diffuse cloud top emissions from CGs are brighter than from ICs, and hence are easier to detect from space. Still other investigators claim that the detection efficiency of CGs and ICs is about the same because modern detector sensitivity is good enough to "see" either flash type no matter which produces a brighter cloud top emission. To better assess which of these opinions should be accepted, we introduce an extension of a Boltzmann lightning radiative transfer model previously developed. It considers characteristics of the cloud (geometry, dimensions, scattering properties) and specific lightning channel properties (length, geometry, location, current, optical wave front propagation speed/direction). As such, it represents the most detailed modeling effort to date. At least in the few cases studied thus far, it was found that IC flashes appear brighter at cloud top than the lower altitude negative ground flashes, but additional model runs are to be examined before finalizing our general conclusions.

Solakiewicz, Richard↗
NASA NTRS1992

Outer heliospheric radio emissions. II - Foreshock source models

Observations of LF radio emissions in the range 2-3 kHz by the Voyager spacecraft during the intervals 1983-1987 and 1989 to the present while at heliocentric distances greater than 11 AU are reported. New analyses of the wave data are presented, and the characteristics of the radiation are reviewed and discussed. Two classes of events are distinguished: transient events with varying starting frequencies that drift upward in frequency and a relatively continuous component that remains near 2 kHz. Evidence for multiple transient sources and for extension of the 2-kHz component above the 2.4-kHz interference signal is presented. The transient emissions are interpreted in terms of radiation generated at multiples of the plasma frequency when solar wind density enhancements enter one or more regions of a foreshock sunward of the inner heliospheric shock. Solar wind density enhancements by factors of 4-10 are observed. Propagation effects, the number of radiation sources, and the time variability, frequency drift, and varying starting frequencies of the transient events are discussed in terms of foreshock sources.

Cairns, Iver H.↗
NASA NTRS1988

Wave/current interaction model

The wave-current interaction for the application to remote sensing data via numerical simulations and data comparison is modelled. Using the field data of surface current shear, wind condition and ambient wave spectrum, the numerical simulations of directional wave spectrum evolution were used to interpret and to compare with the aircraft data from Radar Ocean Wave Spectrometer (ROWS) and Surface Contour Radar (SCR) across the front during Frontal Air Sea Interaction Experiment (FASINEX). The wave-ice interaction was inspired by the observation of large amplitude waves hundreds of kms inside the ice pack in the Weddell Sea, resulting in breakup of the ice pack. The developed analysis of processes includes the refraction of waves at the pack edge, the effects of pack compression on wave propagation, wave train stability and buckling stability in the ice pack. Sources of pack compression and interaction between wave momentum and pack compression are investigated. Viscous camping of propagating waves in the marginal ice zone are also studied. The analysis suggests an explanation for the change in wave dispersion observed from the ship and the sequence of processes that cause ice pack breakup, pressure ridge formation and the formation of open bands of water.

Liu, A. K.↗
NASA NTRS2000

Evolution of Bromoform in a Global Chemistry and Transport Model

It is well known that many chlorine and bromine compounds that are inert in the troposphere are destroyed in the stratosphere and contribute to the stratospheric burden of reactive chlorine and bromine species. But the contribution from those chlorine and bromine compounds which are reactive in the troposphere is less certain because it is not known whether convection can transport these gases to the upper troposphere rapidly enough to overcome their short tropospheric lifetimes. We examine this issue using a three-dimensional chemistry and transport model to simulate the evolution of three gases which have surface sources, bromoform (CHBr3), methyl chloroform (CH3CCl3), and carbon dioxide (CO2). Our objective is to determine if CHBr3 might enhance the lower stratospheric burden of reactive bromine. The other two gases provide tests of the quality of the simulation. Both CHBr3 and CH3CCl3 are destroyed in the troposphere by reaction with hydroxyl (OH), whose latitudinal and monthly variation is provided by a two-dimensional model and upon which a diurnal variation is imposed. Comparison of the lifetime of CH3CCl3 computed from observations (5 years) with the lifetime computed from the simulation provides an integrated test of the model's transport and photochemistry. Observations also show that CO2 exhibits a strong seasonal cycle in the northern hemisphere troposphere that is not propagated directly across the tropopause into the lower stratosphere. Thus, maintenance of the observed troposphere-stratosphere distinctness of CO2 in the presence of convection is a critical benchmark for meeting our objective.

Douglass, Anne R.↗
DOE OSTI2024

Inverse design for waveguide dispersion with a differentiable mode solver

Inverse design of optical components based on adjoint sensitivity analysis has the potential to address the most challenging photonic engineering problems. However, existing inverse design tools based on finite-difference-time-domain (FDTD) models are poorly suited for optimizing waveguide modes for adiabatic transformation or perturbative coupling, which lies at the heart of many important photonic devices. Among these, dispersion engineering of optical waveguides is especially challenging in ultrafast and nonlinear optical applications involving broad optical bandwidths and frequency-dependent anisotropic dielectric material response. In this work, we develop gradient back-propagation through a general-purpose electromagnetic eigenmode solver and use it to demonstrate waveguide dispersion optimization for second harmonic generation with maximized phase-matching bandwidth. This optimization of three design parameters converges in eight steps, reducing the computational cost of optimization by ∼100x compared to exhaustive search and identifying new designs for broadband optical frequency doubling of laser sources in the 1.3–1.4 µm wavelength range. Furthermore, we demonstrate that the computational cost of gradient back-propagation is independent of the number of parameters, as required for optimization of complex geometries. This technique enables practical inverse design for a broad range of previously intractable photonic devices.

Gray, Dodd (ORCID:000000030469599X)↗
NASA NTRS1984

Interaction of aerodynamic noise with laminar boundary layers in supersonic wind tunnels

The interaction between incoming aerodynamic noise and the supersonic laminar boundary layer is studied. The noise field is modeled as a Mach wave radiation field consisting of discrete waves emanating from coherent turbulent entities moving downstream within the supersonic turbulent boundary layer. The individual disturbances are likened to miniature sonic booms and the laminar boundary layer is staffed by the waves as the sources move downstream. The mean, autocorrelation, and power spectral density of the field are expressed in terms of the wave shapes and their average arrival rates. Some consideration is given to the possible appreciable thickness of the weak shock fronts. The emphasis in the interaction analysis is on the behavior of the shocklets in the noise field. The shocklets are shown to be focused by the laminar boundary layer in its outer region. Borrowing wave propagation terminology, this region is termed the caustic region. Using scaling laws from sonic boom work, focus factors at the caustic are estimated to vary from 2 to 6 for incoming shocklet strengths of 1 to .01 percent of the free stream pressure level. The situation regarding experimental evidence of the caustic region is reviewed.

Schopper, M. R.↗
NASA NTRS1998

Inner Magnetospheric Superthermal Electron Transport: Photoelectron and Plasma Sheet Electron Sources

Two time-dependent kinetic models of superthermal electron transport are combined to conduct global calculations of the nonthermal electron distribution function throughout the inner magnetosphere. It is shown that the energy range of validity for this combined model extends down to the superthermal-thermal intersection at a few eV, allowing for the calculation of the entire distribution function and thus an accurate heating rate to the thermal plasma. Because of the linearity of the formulas, the source terms are separated to calculate the distributions from the various populations, namely photoelectrons (PEs) and plasma sheet electrons (PSEs). These distributions are discussed in detail, examining the processes responsible for their formation in the various regions of the inner magnetosphere. It is shown that convection, corotation, and Coulomb collisions are the dominant processes in the formation of the PE distribution function, and that PSEs are dominated by the interplay between the drift terms. Of note is that the PEs propagate around the nightside in a narrow channel at the edge of the plasmasphere as Coulomb collisions reduce the fluxes inside of this and convection compresses the flux tubes inward. These distributions are then recombined to show the development of the total superthermal electron distribution function in the inner magnetosphere and their influence on the thermal plasma. PEs usually dominate the dayside heating, with integral energy fluxes to the ionosphere reaching 10(exp 10) eV/sq cm/s in the plasmasphere, while heating from the PSEs typically does not exceed 10(exp 8)eV/sq cm/s. On the nightside, the inner plasmasphere is usually unheated by superthermal electrons. A feature of these combined spectra is that the distribution often has upward slopes with energy, particularly at the crossover from PE to PSE dominance, indicating that instabilities are possible.

Khazanov, G. V.↗
NASA NTRS2003

Non-Migrating Tides, with Zonally Symmetric Component, Generated in the Mesosphere

For comparison with measurements from the TIMED satellite and coordinated ground based observations, we discuss results from our Numerical Spectral Model (NSM) that incorporates the Doppler Spread Parameterization (Hines, 1997) for small-scale gravity waves (GWs). The NSM extends from the ground into the thermosphere and describes the major dynamical features of the atmosphere including the wave driven equatorial oscillations (QBO and SAO), and the seasonal variations of tides and planetary waves. With emphasis on the non-migrating tides, having periods of 24 and 12 hours, we discuss our modeling results that account for the classical migrating solar excitation sources only. As reported earlier, the NSM reproduces the observed seasonal variations and in particular the large equinoctial maxima in the amplitude of the migrating diurnal tide at altitudes around 90 km. Filtering of the tide by the zonal circulation and GW momentum deposition was identified as the cause. The GWs were also shown to produce a strong non-linear interaction between the diurnal and semi-diurnal tides. Confined largely to the mesosphere, the NSM produces through dynamical interactions a relatively large contribution of non-migrating tides. A striking feature is seen in the diurnal and semi-diurnal oscillations of the zonal mean (m = 0). Eastward propagating tides are also generated for zonal wave numbers m = 1 to 4. When the NSM is run without GWs, the amplitudes for the non-migrating tides, including m = 0, are generally small. Planetary wave interaction and non-linear coupling that involves the filtering of GWs and related height integration of dynamical features are discussed as possible mechanisms for generating these non-migrating tides in the NSM. As is the case for the solar migrating tides, the non-migrating tides reveal persistent seasonal variations. Under the influence of the QBO and SAO, interannual variations are produced.

Mayr, H. G.↗