Search NASA⌕ Search

SEARCH · Search NASA

Results for “ensemble dynamics”

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 163 records · Page 9

Synergy of Observations and Dynamo Models to Understand and Predict Solar Activity Cycles

The long-standing problem of understanding the evolution of the global magnetic fields that drive solar activity through different temporal scales is becoming more tractable because, in addition to 400 years of sunspot records, we now have almost 4 solar cycles of magnetic field observations. These observations allow us to discern physical connections between dynamo model variables and observations using data assimilation analysis. In particular, the Ensemble Kalman Filter approach takes into account uncertainties in both observations and modeling and allows us to make reliable forecasts of solar cycle activity cycles by using a relatively simple non-linear dynamical model of the solar dynamo. To expand this approach for more complex 2D and 3D dynamo modeling, it is necessary to decompose the observed synoptic magnetograms into poloidal and toroidal field components. In this presentation I will present initial results on magnetogram decomposition and assimilation of magnetogram data into dynamo modeling.

Synergy↗

Gas Giants

The gas giants (Jupiter and Saturn) and icy giants (Uranus and Neptune) are fluid planets with atmospheres primarily made of hydrogen and helium. The part of their atmospheres accessible to remote sensing occupies only a small fraction of their radii (0.05%). Clouds and hazes form around the 1 bar altitude pressure level and extend vertically, according to the thermochemical models, in a layer with a thickness of 200_500 km where temperature increases with depth (usually known as the "weather layer"). Clouds made of NH3, NH4SH, H2O (in Jupiter and Saturn), with the addition of CH4 (in Uranus and Neptune), cover the planet in stratified layers that are mixed with unknown hromophore agents. Dynamical phenomena in the weather layer shape different cloud patterns that define the visible appearance of these planets. In the thermal part of the spectrum clouds act as opacity sources providing brightness contrasts. The ensemble of cloud morphologies in terms of shapes, sizes and albedos allows their use as tracers of the atmospheric motions in the weather layer (Fig. 4.1). This is the main tool employed so far to study the winds on these fourplanets.

Sanchez-Lavega, A.↗

Version 4 of the SMAP Level-4 Soil Moisture Algorithm and Data Product

The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with a mean latency of ~2.5 days. The underlying L4_SM algorithm assimilates SMAP radiometer brightness temperature (Tb) observations into the NASA Catchment land surface model using a spatially-distributed ensemble Kalman filter. Version 4 of the L4_SM modeling system includes a reduction in the upward recharge of surface soil moisture from below under non-equilibrium conditions, resulting in reduced bias and improved dynamic range of L4_SM surface soil moisture compared to earlier versions. This change and additional technical modifications to the system reduce the mean and standard deviation of the observation-minus-forecast Tb residuals and overall soil moisture analysis increments while maintaining the skill of the L4_SM soil moisture estimates versus independent in situ measurements; the average, bias-adjusted RMSE in Version 4 is 0.039 m(exp 3) m(exp -3) for surface and 0.026 m(exp 3) m(exp -3) for root-zone soil moisture. Moreover, the coverage of assimilated SMAP observations in Version 4 is near-global owing to the use of additional satellite Tb records for algorithm calibration. L4_SM soil moisture uncertainty estimates are biased low (by 0.01-0.02 m(exp 3) m(exp -3)) against actual errors (computed versus in situ measurements). L4_SM runoff estimates, an additional product of the L4_SM algorithm, are biased low (by 35 mm year (exp -1)) against streamflow measurements. Compared to Version 3, bias in Version 4 is reduced by 46% for surface soil moisture uncertainty estimates and by 33% for runoff estimates.

RMSE↗

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↗

Solar wind three; Proceedings of the Third Conference, Pacific Grove, Calif., March 25-29, 1974

Papers are presented relating to the study of solar abundances, the evolution of the solar wind, solar corona dynamics, solar wind interactions, and stellar winds. Some of the topics covered include the nuclear composition of cosmic rays, the solar wind as deduced from lunar samples, spectral characteristics of flares, numerical simulation of interplanetary shock ensembles, radial gradient of solar wind velocity from 1 to 5 AU, the enhancement of solar wind fluctuations with scale size near the proton gyroradius, solar wind-Venus interaction observed from magnetic field experiment on Mariner 10, and binary stellar winds. Individual items are announced in this issue.

Russell, C. T.↗

Multivariate Error Covariance Estimates by Monte-Carlo Simulation for Assimilation Studies in the Pacific Ocean

One of the most difficult aspects of ocean state estimation is the prescription of the model forecast error covariances. The paucity of ocean observations limits our ability to estimate the covariance structures from model-observation differences. In most practical applications, simple covariances are usually prescribed. Rarely are cross-covariances between different model variables used. Here a comparison is made between a univariate Optimal Interpolation (UOI) scheme and a multivariate OI algorithm (MvOI) in the assimilation of ocean temperature. In the UOI case only temperature is updated using a Gaussian covariance function and in the MvOI salinity, zonal and meridional velocities as well as temperature, are updated using an empirically estimated multivariate covariance matrix. Earlier studies have shown that a univariate OI has a detrimental effect on the salinity and velocity fields of the model. Apparently, in a sequential framework it is important to analyze temperature and salinity together. For the MvOI an estimation of the model error statistics is made by Monte-Carlo techniques from an ensemble of model integrations. An important advantage of using an ensemble of ocean states is that it provides a natural way to estimate cross-covariances between the fields of different physical variables constituting the model state vector, at the same time incorporating the model's dynamical and thermodynamical constraints as well as the effects of physical boundaries. Only temperature observations from the Tropical Atmosphere-Ocean array have been assimilated in this study. In order to investigate the efficacy of the multivariate scheme two data assimilation experiments are validated with a large independent set of recently published subsurface observations of salinity, zonal velocity and temperature. For reference, a third control run with no data assimilation is used to check how the data assimilation affects systematic model errors. While the performance of the UOI and MvOI is similar with respect to the temperature field, the salinity and velocity fields are greatly improved when multivariate correction is used, as evident from the analyses of the rms differences of these fields and independent observations. The MvOI assimilation is found to improve upon the control run in generating the water masses with properties close to the observed, while the UOI failed to maintain the temperature and salinity structure.

Borovikov, Anna↗

Solar Activity Modeling: From Subgranular Dynamical Scales to the Solar Cycles

The dynamical effects of solar magnetoconvection span a wide range spatial and temporal scales that extend from the interior to the corona and from fast turbulent motions to global magnetic activity. To study the solar activity on short temporal scales (from minutes to hours), we use 3D radiative MHD simulations that allow us to investigate complex turbulent interactions that drive various phenomena, such as plasma eruptions, spontaneous formation of magnetic structures, funnel-like structures and magnetic loops in the corona, and others. In particular, we focus on multi-scale processes of energy exchange across layers of the solar interior and atmosphere, which contribute to coronal heating and eruptive dynamics. For modeling global-scale activity, we use a data assimilation approach that has demonstrated great potential for building reliable long-term forecasts of solar activity. In particular, it has been shown that the Ensemble Kalman Filter (EnKF) method applied to the Parker-Kleeorin-Ruzmakin dynamo model is capable of predicting solar activity up to one sunspot cycle ahead in time, as well as estimating the properties of the next cycle a few years before it begins. In this presentation, using the available magnetogram data, we discuss development of the methodology and forecast quality criteria (including forecast uncertainties and sources of errors). We demonstrate the influence of observational limitations on prediction accuracy, and we present the EnKF predictions of the upcoming Solar Cycle (25) based on both the sunspot number series and observed magnetic fields and discuss the uncertainties and potential of the data assimilation approach for modeling and forecasting solar activity.

Kitiashvili, I. N.↗

Uncertainty in Soil Moisture Retrievals: an Ensemble Approach Using SMOS L-Band Microwave Data

The uncertainty of soil moisture (SM) retrievals from satellite brightness temperature (TB) observations depends primarily on the choice of radiative transfer model (RTM) parameters, prior SM information and TB inputs. This paper studies the sensitivity of several (quasi-)operational and experimental SM retrieval products from the Soil Moisture Ocean Salinity (SMOS) mission to these choices at 11 reference sites, located in 7 watersheds across the United States (US). Different literature-based RTM parameter sets cause large biases between retrievals. Whereas typical RTM parameter sets are calibrated for SM retrievals, it is shown that a parameter set carefully optimized for TB forward modeling can also be used for retrieving SM. It is also shown that the inclusion of dynamic prior SM estimates in a Bayesian retrieval scheme can strongly improve SM retrievals, regardless of the choice of RTM parameters, and that the use of multi-angular and multi-polarization TB does not necessarily lead to superior retrievals compared to retrievals based on TB data at a single incidence angle and polarization. The second part of this paper evaluates ensemble uncertainty metrics for SM retrievals obtained by propagating a wide range of RTM parameters through the RTM. As expected for bounded variables, the spread in the ensemble SM retrievals is smallest for wet and dry SM values and highest for intermediate SM values. After removal of the strong long-term mean bias associated with the RTM parameter values for individual ensemble members, the remaining anomaly ensemble SM spread of 0.037 cu m/cu m approximates the actual time series unbiased root-mean-square-difference of 0.042 cu m/cu m between ensemble mean retrievals and in situ reference data across the reference sites. However, the temporal variability in the anomaly ensemble spread reveals higher-order biases in the retrieval error, which should be accounted for when characterizing retrieval error.

Jan Quets↗

Computation of launch vehicle system requirements using hybrid computer.

Formulating an adequate statistical statement concerning space vehicle dynamic states requires the combination of the statistics of the environment and the vehicle's basic parameters. The basic ingredient of the environment for the Space Shuttle launch phase is the winds, which are represented by an ensemble of measured winds (150/month), which today constitute the best statistical representation. The problem treated in this paper then becomes twofold: (1) how can the vehicle response be analyzed using wind ensembles, and (2) how can the vehicle parameter variations be treated in conjunction with wind ensembles.

Ryan, R. S.↗

Statistical mechanics of negative temperature states

The dynamics of two-dimensional interacting line vortices is identical to that of the two-dimensional electrostatic guiding center plasma. Both are Hamiltonian systems and are therefore susceptible to statistical mechanical treatments. The predictions of the microcanonical ensemble are explored for this system. Interest focuses primarily on the regime of total positive interaction energy, which should be above the Onsager negative temperature threshold. Calculations of the probability distribution for a component by means of the central limit theorem are carried out in the manner of Khinchin. The probability distribution of a component reduced to the usual Gibbs distribution in the regime of positive temperatures, and is still explicitly calculable for negative temperatures. The negative temperature states are neither quiescent nor spatially uniform. Expressions for the temperature are explicitly provided in terms of the total particle energy and particle number. A BBGKY hierarchy can be derived for both temperature regimes. Numerical simulations involving solutions of the equations of motion of 4008 particles are presented.

Montgomery, D. C.↗

An Intercomparison of Changes Associated with Earth's Lower Tropospheric Temperature Using Traditional and AMIP-Style Reanalyses

Reanalyses have become an integral tool for evaluating regional and global climate variations, and an important component of this is modifications to the energy budget. Reductions in Arctic Sea ice extent has induced an albedo feedback, causing the Arctic to warm more rapidly than anywhere else in the world, referred to as "Arctic Amplification." This has been demonstrated by observations and numerous reanalyses, including the Modern Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2). However, the Arctic Amplification signal is non-existent in a ten member ensemble of the MERRA-2 Atmospheric Model Intercomparison Project (M2AMIP) simulations, using the same prescribed climate forcing, including Sea Surface Temperature (SST) and ice. An evaluation of the temperature tendency within the lower troposphere due to radiation, moisture, and dynamics as well as the surface energy budget in MERRA-2 and M2AMIP will demonstrate that despite identical prescribed SSTs and sea ice in both versions, enhanced warming in the Arctic in MERRA-2 is in response to the analysis increment tendency due to temperature observations. Furthermore, the role of boundary conditions, model biases and changes in observing systems on the Arctic Amplification signal will be assessed. Literature on the topic of Arctic Amplification demonstrates that the enhanced warming begins in the mid-1990s. Anomalously warm Arctic SST in the early time period of MERRA-2 can mute the trend in Arctic lower troposphere temperature without the constraint of observations in M2AMIP. Additionally, MERRA-2 uses three distinct datasets of SST and sea ice concentration, which also plays a role.

Arcti↗

Dynamical Forecasts of Tropical Terrestrial Carbon Fluxes with the NASA S2S Retrospective Forecast System

Recent advances in the ability to predict climate anomalies at sub-seasonal to seasonal (S2S) timescales allow us to explore the possibility of forecasting carbon flux anomalies. Although carbon flux forecasting is a relatively new concept, it is potentially beneficial as it can help us better understand global and regional land-atmosphere carbon feedbacks associated with climate variations and can provide guidance for future field mission design. Here we evaluate the skill of forecasted terrestrial carbon anomalies generated from meteorological anomalies produced with the NASA Global Modeling and Assimilation Office (GMAO) S2S forecast system. We focus here on three representative time periods (the most recent 2015-2016 El Nino, 2011 La Nina, and 2014 as a neutral year), with each corresponding 9-month forecast comprising four ensemble members initialized in the preceding December. The meteorological variables produced by the GMAO forecast system were bias-corrected using a climatology derived from the Modern Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) before being used to drive a suite of offline simulations with the NASA Catchment-CN terrestrial biosphere model, a model that computes water-energy-carbon dynamics. Forecasts are evaluated by comparing against satellite-driven estimates of gross primary production (GPP) and inverse model estimates of net carbon flux that incorporate satellite carbon dioxide measurements. We find that the restrospectively predicted carbon fluxes in the tropics reasonably reproduce the signs and magnitudes of the observed anomalies between the 2015-2016 El Nino and the 2011 La Nina for both net flux and GPP. For instance, for the El Nino period, the magnitude of the forecasted negative GPP anomaly in the South American tropics (which undergoes anomalously warm and dry conditions) agrees with the observed GPP anomaly at leads of up to three or four months. Overall, this study demonstrates potential skill in the forecast of biospheric carbon fluxes a few months in advance, a capability that could contribute to attribution studies focusing on carbon flux variations and support innovative observation strategies in the future.

Lee, Eunjee↗

The Dynamics of Coupled Oscillator Phase Control

Arrays of coupled oscillators have been proposed as means of realizing high power rf sources via coherent spatial power combining. In such applications, a uniform phase distribution over the aperture is usually desired. However, it has been shown that by detuning some of the oscillators away from the oscillation frequency of the ensemble of oscillators, one may achieve other useful aperture phase distributions. Of particular interest among those achievable are linear phase distributions because these result in steering of the output rf beam away from the broadside direction. The theory describing the behavior of such arrays of coupled oscillators is quite complicated since the phenomena involved are inherently nonlinear. However, a simplified theory has been developed which facilitates intuitive understanding. This simplified theory is based on a "continuum model" in which the aperture phase is represented by a continuous function of the aperture coordinates. A challenging aspect of the development of this theory is the derivation of appropriate boundary conditions at the edges or ends of the array.

Pogorzelski, R. J.↗

Inertial measurement unit pre-processors and post-flight STS-1 comparisons

The flight results show that the relative tri-redundant Inertial Measurement Unit IMU performance throughout the entire entry flight was within the expected accuracy. Comparisons are presented which show differences in the accumulated sensed velocity changes as measured by the tri-redundant IMUs (in Mean Equator and Equinox of 1950.0), differences in the equivalent inertial Euler angles as measured with respect to the M50 system, and finally, preliminary instrument calibrations determined relative to the ensemble average measurement set. Also, differences in the derived body axes rates and accelerations are presented. Because of the excellent performance of the IMUs during the STS-1 entry, the selection as to which particular IMU would best serve as the dynamic data source for entry reconstruction is arbitrary.

Findlay, J. T.↗

Multidecadal Changes in Lower Stratospheric Ozone: Variability Vs. Trends

As upper stratospheric ozone appears to be recovering as a result of decreasing chlorine loading following the implementation of the Montreal Protocol and its amendments and in agreement with model projections, several recent studies report an apparent decline of ozone concentrations in the lower stratosphere in the last two decades, particularly in the extratropics. Our previous work as well as at least two other studies provide evidence that this decline results from transport changes rather than an intensification of chemical depletion. It remains unclear whether these changes represent long-term internal variability or are a consequence of a climate forcing. Here we perform free-running ensembles of the recent past (1980-2016) using the Goddard Earth Observing System Model (GEOS) at the cubed sphere C180 (approximately half degree) resolution. Two suites of 10-member ensembles are performed, one in which observed sea surface temperature (SSTs) are fully prescribed, and the other in which the linear SST trend over the recent past is removed so as to only retain internal variability. We evaluate the trends in both ozone as well as two idealized tracers with prescribed uniform loss that are used to isolate the role of transport from chemistry and emissions. Probability-distribution-functions of the trends in both ozone and idealized tracers are compared among ensemble members and with observed trends in order to evaluate the likelihood of recent observed declines in lower stratospheric ozone, relative to large internal variability. Moreover, comparisons among simulations with and without imposed SST trends indicate the extent to which dynamically-driven ozone trends reflect forced trends or internal variability in lower stratospheric dynamics.

Wargan, Krzysztof↗

Assimilation of DAWN Doppler Wind Lidar Data During the 2017 Convective Processes Experiment (CPEX): Impact on Precipitation and Flow Structure

An improved representation of 3-D air motion and precipitation structure through forecast models and assimilation of observations is vital for improvements in weather forecasting capabilities. However, there are few independent data to properly validate a model forecast of precipitation structure when the underlying dynamics are evolving on short convective timescales. Using data from the JPL Ku/Ka-band Airborne Precipitation Radar (APR-2) and the 2 μmDoppler Aerosol Wind (DAWN) lidar collected during the2017 Convective Processes Experiment (CPEX), the NASA Unified Weather Research and Forecasting (WRF) Ensemble Data Assimilation System (EDAS) modeling system was used to quantify the impact of high-resolution sparsely sampled DAWN measurements on the analyzed variables and on the forecast when the DAWN winds were assimilated. Over-all, the assimilation of the DAWN wind profiles had a discernible impact on the wind field as well as the evolution and timing of the 3-D precipitation structure. Analysis of individual variables revealed that the assimilation of the DAWN winds resulted in important and coherent modifications of the environment. It led to an increase in the near-surface convergence, temperature, and water vapor, creating more favorable conditions for the development of convection exactly where it was observed (but not present in the control run). Comparison to APR-2 and observations by the Global Precipitation Measurement (GPM) satellite shows a much-improved forecast after the assimilation of the DAWN winds – development of precipitation where there was none, more organized precipitation where there was some, and a much more intense and organized cold pool, similar to the analysis of the dropsonde data. The onset of the vertical evolution of the precipitation showed similar radar-derived cloud-top heights, but delayed in time. While this investigation was limited to a single CPEX flight date, the investigation design is appropriate for further investigation of the impact of airborne Doppler wind lidar observations upon short-term convective precipitation forecasts

DAWN↗

Impact of Satellite Sea Surface Salinity Observations on ENSO Predictions from the GMAO S2S Forecast System

Sea surface salinity (SSS) observations from space allow us to investigate if improved estimates of near-surface density stratification and associated mixing can positively impact seasonal to interannual variability of tropical Pacific Ocean dynamics as well as dynamical ENSO forecasts. For the first part of the presentation, we utilize our intermediate-complexity coupled model. Baseline experiments assimilate satellite sea level (multi-satellite gridded AVISO, 2013), SST (Reynolds et al., 2004), and in situ subsurface temperature and salinity observations (GTSPP NODC, 2006). These baseline experiments are then compared with experiments that additionally assimilate Aquarius (V5.0 Lilly and Lagerloef, 2008) and SMAP (V4.0 Fore et al., 2016) SSS. Twelve-month forecasts are initialized for each month from September 2011 to September 2017. For initialization of the coupled forecast, the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10S-5N. This pattern enhances air/sea interaction and amplifies the equatorial Kelvin wave signal. We find that including satellite SSS significantly improves NINO3.4 sea surface temperature anomaly validation over most forecast lead times. We next assess how different satellite SSS products impact the validation of ENSO forecasts. SMAP V4 reduces the salty bias in the western Pacific and so is an improvement upon SMAP V2 and SMOS V2 (Boutin et al., 2017) has similar validation characteristics as a combination of Aquarius and SMAP V4. Next we shift to present results from the NASA GMAO Sub-seasonal to seasonal (S2S_v2.1) production coupled model (i.e. the same model that contributes ENSO forecasts to the North American Multi-Model Ensemble Experiment). From March to June 2015, the availability of two overlapping satellite SSS instruments, Aquarius and SMAP, allows a unique opportunity to compare and contrast forecasts initialized with the benefit of these two satellite SSS observation types. We assess the impact of satellite sea surface salinity (SSS) observations on dynamical ENSO forecasts for the big 2015 El Nino event. We will present distinct experiments for the overlap period that include 1) freely evolving SSS (i.e. no satellite SSS as the production system), 2) Aquarius, and 3) SMAP initialization. Our results show that using Aquarius slightly improves validation of the reanalysis (including sea level and temperature statistics). Our production system without SSS assimilation generated too warm forecasts for the 2015 El Nino from March initial conditions. Incorporating Aquarius into initialization of the coupled system leads to a deeper, more realistic MLD that acts to damp the downwelling Kelvin signal and slightly cool NINO3.4 SST. With Aquarius the forecasts better match the observed amplitude of the 2015 event. On the other hand, SMAP V2 relaxation generally degrades validation statistics. At forecast initialization, SMAP is much too salty within 10o of the equator, leading to deeper MLD east of 165W. This deeper MLD leads to over-damping of the downwelling signal (i.e. relative upwelling), in turn leading to relatively too cool ENSO forecasts.

Hackert, E.↗

Self-Consistent Magnetosphere-Ionosphere Coupling and Associated Plasma Energization Processes

Magnetosphere-Ionosphere (MI) coupling and associated with this process electron and ion energization processes have interested scientists for decades and, in spite of experimental and theoretical research efforts, are still ones of the least well known dynamic processes in space plasma physics. The reason for this is that the numerous physical processes associated with MI coupling occur over multiple spatial lengths and temporal scales. One typical example of MI coupling is large scale ring current (RC) electrodynamic coupling that includes calculation of the magnetospheric electric field that is consistent with the ring current (RC) distribution. A general scheme for numerical simulation of such large-scale magnetosphere-ionosphere coupling processes has been presented earlier in many works. The mathematical formulation of these models are based on "modified frozen-in flux theorem" for an ensemble of adiabatically drifting particles in the magnetosphere. By tracking the flow of particles through the inner magnetosphere, the bounce-averaged phase space density of the hot ions and electrons can be reconstructed and the magnetospheric electric field can be calculated such that it is consistent with the particle distribution in the magnetosphere. The new a self-consistent ring current model has been developed that couples electron and ion magnetospheric dynamics with calculation of electric field. Two new features were taken into account in addition to the RC ions, we solve an electron kinetic equation in our model, self-consistently including these results in the solution. Second, using different analytical relationships, we calculate the height integrated ionospheric conductances as the function of precipitated high energy magnetospheric electrons and ions as produced by our model. This results in fundamental changes to the electric potential pattern in the inner magnetosphere, with a smaller Alfven boundary than previous potential formulations would predict but one consistent with recent satellite observations. This leads to deeper penetration of the plasma sheet ions and electrons into the inner magnetosphere and more effective ring current ions and electron energization.

Khazanov, G. V.↗