Search NASASearch

Engineering topics

Sojka, J. J.

Publications and source records attributed to Sojka, J. J..

At least 19 records

CEDAR-GEM Challenge for Systematic Assessment of Ionosphere/Thermosphere Models in Predicting TEC During the 2006 December Storm Event

In order to assess current modeling capability of reproducing storm impacts on total electron content (TEC), we considered quantities such as TEC, TEC changes compared to quiet time values, and the maximum value of the TEC and TEC changes during a storm. We compared the quantities obtained from ionospheric models against ground-based GPS TEC measurements during the 2006 AGU storm event (14-15 December 2006) in the selected eight longitude sectors. We used 15 simulations obtained from eight ionospheric models, including empirical, physics-based, coupled ionosphere-thermosphere, and data assimilation models. To quantitatively evaluate performance of the models in TEC prediction during the storm, we calculated skill scores such as RMS error, Normalized RMS error (NRMSE), ratio of the modeled to observed maximum increase (Yield), and the difference between the modeled peak time and observed peak time. Furthermore, to investigate latitudinal dependence of the performance of the models, the skill scores were calculated for five latitude regions. Our study shows that RMSE of TEC and TEC changes of the model simulations range from about 3 TECU (total electron content unit, 1 TECU = 1016 el m2) (in high latitudes) to about 13 TECU (in low latitudes), which is larger than latitudinal average GPS TEC error of about 2 TECU. Most model simulations predict TEC better than TEC changes in terms of NRMSE and the difference in peak time, while the opposite holds true in terms of Yield. Model performance strongly depends on the quantities considered, the type of metrics used, and the latitude considered.

Shim, J. S.

CEDAR Electrodynamics Thermosphere Ionosphere (ETI) Challenge for Systematic Assessment of Ionosphere/Thermosphere Models: NmF2, hmF2, and Vertical Drift Using Ground-Based Observations

Objective quantification of model performance based on metrics helps us evaluate the current state of space physics modeling capability, address differences among various modeling approaches, and track model improvements over time. The Coupling, Energetics, and Dynamics of Atmospheric Regions (CEDAR) Electrodynamics Thermosphere Ionosphere (ETI) Challenge was initiated in 2009 to assess accuracy of various ionosphere/thermosphere models in reproducing ionosphere and thermosphere parameters. A total of nine events and five physical parameters were selected to compare between model outputs and observations. The nine events included two strong and one moderate geomagnetic storm events from GEM Challenge events and three moderate storms and three quiet periods from the first half of the International Polar Year (IPY) campaign, which lasted for 2 years, from March 2007 to March 2009. The five physical parameters selected were NmF2 and hmF2 from ISRs and LEO satellites such as CHAMP and COSMIC, vertical drifts at Jicamarca, and electron and neutral densities along the track of the CHAMP satellite. For this study, four different metrics and up to 10 models were used. In this paper, we focus on preliminary results of the study using ground-based measurements, which include NmF2 and hmF2 from Incoherent Scatter Radars (ISRs), and vertical drifts at Jicamarca. The results show that the model performance strongly depends on the type of metrics used, and thus no model is ranked top for all used metrics. The analysis further indicates that performance of the model also varies with latitude and geomagnetic activity level.

Shim, J. S.

Investigation and Development of Data-Driven D-Region Model for HF Systems Impacts

Space Environment Corporation (SEC) and RP Consultants (RPC) are to develop and validate a weather-capable D region model for making High Frequency (HF) absorption predictions in support of the HF communications and radar communities. The weather-capable model will assimilate solar and earth space observations from NASA satellites. The model will account for solar-induced impacts on HF absorption, including X-rays, Solar Proton Events (SPE's), and auroral precipitation. The work plan includes: I . Optimize D-region model to quickly obtain ion and electron densities for proper HF absorption calculations. 2. Develop indices-driven modules for D-region ionization sources for low, mid, & high latitudes including X-rays, cosmic rays, auroral precipitation, & solar protons. (Note: solar spectrum & auroral modules already exist). 3. Setup low-cost monitors of existing HF beacons and add one single-frequency beacon. 4. Use PENEX HF-link database with HF monitor data to validate D-region/HF absorption model using climatological ionization drivers. 5. Develop algorithms to assimilate NASA satellite data of solar, interplanetary, and auroral observations into ionization source modules. 6. Use PENEX HF-link & HF-beacon data for skill score comparison of assimilation versus climatological D-region/HF absorption model. Only some satellites are available for the PENEX time period, thus, HF-beacon data is necessary. 7. Use HF beacon monitors to develop HF-link data assimilation algorithms for regional improvement to the D-region/HF absorption model.

Eccles, J. V.

What is the source of observed annual variations in plasmaspheric density?

Plasmaspheric densities have been observed previously to be higher in December than in June, with the ratio varying between 1.5 and 3.0 and with larger variations at lower L shells. In order to search for the cause of the observed annual variations, we have modeled plasmaspheric density, using a time-dependent hydrodynamic model. On an L = 2 field line with geomagnetic longitude equal to 300 deg, the modeled plasmaspheric densities were a factor of 1.5 times higher in December than in June. The modeled December to June density ratio was found to increase slightly with L shell, in contrast to observations; this discrepancy may be due to the fact that outer plasmaspheric flux tubes are never completely full. In addition, for an L = 2 field line with geomagnetic longitude equal to 120 deg, the modeled plasmaspheric density was higher in June than in December by a factor of about 1.2. Various numerical tests were also performed in order to examine the sensitivity of plasmaspheric density to various parameters. In particular, a large vertical neutral wind was applied in order to raise the O(+) profile which had the effect of raising plasmaspheric density by a factor of 6. This in conjunction with a theoretical analysis suggests that plasmaspheric density levels are very sensitive to O(+) levels in the upper ionosphere. We conclude that annual variations in plasmaspheric density are due to similar variations in ionospheric O(+).

Guiter, S. M.

A correlative comparison of the ring current and auroral electrojects usig geomagnetic indices

From a study of the 21 largest geomagnetic storms during solar cycle 21, a strong correlation is established between the ring current index Dst and the time-weighted accumulation of the 1-hour auroral electrojets indices, AE and AL. The time-weighted accumulation corresponds to convolution of the auroral electrojet indices with an exponential weighting function with an e-folding time of 9.4 hours. The weighted indices AE(sub w) and AL(sub w) have correltation coefficients against Dst ranging between 0.8 and 0.95 for 20 of the 21 storms. Correlation over the entire solar cycle 21 database is also strong but not as strong as for an individual storm. A set of simple Dst prediction functions provide a first approximation of the inferred dependence, but the specific functional relationship of Dst (AL(sub w)) or Dst (AL(sub w)) varies from one storm to the next in a systematic way. This variation reveals a missing parametric dependence in the transfer function. However, our results indicate that auroral electroject indices are potentially useful for predicting storm time enhancements of ring current intensity with a few hours lead time.

Cade, W. B., III

Ionospheric response to the sustained high geomagnetic activity during the March '89 great storm

A simulation was conducted to model the high-latitude ionospheric to the sustaied level of high geomagnetic activity for the great magnetic storm period of March 13-14, 1989. The geomagnetic and solar activity indices and the Defense Meterological Satellite Program (DMSP) F8 and F9 satellite data for particle precipitation and high-latitude convection were used as inputs to a time-dependent ionospheric model (TDIM). The results of the TDIM were compared to both DMSP plasma density data and ground-based total electron content (TEC) measurements for the great storm period as well as with earlier storm observations. The comparisons show that the overall structure of the high-latitude ionosphere was dominated by an increased convection speed within the polar cap that led to increased ion temperatures. In turn, this enhanced the NO(+) density, raised the atomic-to-molecular ion transition height to over 300 km, decreased N(sub m)F(sub 2), increased h(sub m)F(sub 2), and in places either increased n(sub e) at 800 km or slightly decreased it. The morphology of the ionosphere under these extreme conditions was considerably different than that modeled for less distributed intervals. These differences included the character of the dayside tongue of ionization that no longer extended deep into the polar cap. Instead, as a result of the ion heating and consequent reduction in N(sub m)F(sub 2), a large polar hole occupied much of the polar region. This polar hole extended beyond the auroral oval and merged with the night sector midatitude trough. The limitaions associated with the applicability of the TDIM to the geomagnetic conditions present on March 13 and 14 are discussed. The primary limitations of the TDIM derive from the limited temporal resolution of the model input parameters and the lack of suitably dynamic thermospheric specification for the great storm conditions. These limitations leads to midlatitude ionospheric storm phases that do no follow those observed.

Sojka, J. J.

Patches in the polar ionosphere: UT and seasonal dependence

The seasonal and Universal Time (UT) dependencies of patches in the polar ionosphere are simulated using the Utah State University time dependent ionospheric model (TDIM). Patch formation is achieved by changing the plasma convection pattern in response to temporal changes in the interplanetary magnetic field (IMF) B(sub y) component during periods of southwrd IMF. This mechanism redirects the plasma flow from the dayside high-density region, which is the source of the tongue of ionization (TOI) density feature, through the throat and leads to patches, rather than a continuous TOI. The model predicts that the patches are absent at winter solstice (northern hemisphere) between 0800 and 1200 UT and that they have their largest seasonal intensity at winter solstice between 2000 and 2400 UT. Between winter solstice and equinox, patches are strong and present all day. Patches are present in summer as well, although their intensity is only tens of percent above the background density. These winter-to-equinox findings are also shown to be consistent with observations. The model was also used to predict times at which patch observations could be performed to determine the contributions from other patch mechanisms. This observational window is +/- 20 days about winter solstice between 0800 and 1200 UT in the northern hemisphere. In this observational window the TOI is either absent or reduced to a very low density. Hence the time dependent electric field mechanism considered in this study does not produce patches, and if they are observed, then they must be due to some other mechanism.

Sojka, J. J.

Effect of high-latitude ionospheric convection on Sun-aligned polar caps

A coupled magnetospheric-ionospheric (M-I) magnetohydrodynamic (MHD) model has been used to simulate the formation of Sun-aligned polar cap arcs for a variety of interplanetary magnetic field (IMF) dependent polar cap convection fields. The formation process involves launching an Alfven shear wave from the magnetosphere to the ionosphere where the ionospheric conductance can react self-consistently to changes in the upward currents. We assume that the initial Alfven shear wave is the result of solar wind-magnetosphere interactions. The simulations show how the E region density is affected by the changes in the electron precipitation that are associated with the upward currents. These changes in conductance lead to both a modified Alfven wave reflection at the ionosphere and the generation of secondary Alfven waves in the ionosphere. The ensuing bouncing of the Alfven waves between the ionosphere and magnetosphere is followed until an asymptotic solution is obtained. At the magnetosphere the Alfven waves reflect at a fixed boundary. The coupled M-I Sun-aligned polar cap arc model of Zhu et al.(1993a) is used to carry out the simulations. This study focuses on the dependence of the polar cap arc formation on the background (global) convection pattern. Since the polar cap arcs occur for northward and strong B(sub y) IMF conditions, a variety of background convection patterns can exist when the arcs are present. The study shows that polar cap arcs can be formed for all these convection patterns; however, the arc features are dramatically different for the different patterns. For weak sunward convection a relatively confined single pair of current sheets is associated with the imposed Alfven shear wave structure. However, when the electric field exceeds a threshold, the arc structure intensifies, and the conductance increases as does the local Joule heating rate. These increases are faster than a linear dependence on the background electric field strength. Furthermore, above the threshold, the single current sheet pair splits into multiple current sheet pairs. For the fixed initial ionospheric and magnetospheric conditions used in this study, the separation distance between the current pairs was found to be almost independent of the background electric field strength. For either three-cell or distorted two-cell background convection patterns the arc formation favored the positive B(sub y) case in the northern hemisphere.

Sojka, J. J.

Model study of multiple polar cap arcs: Occurrence and spacing

A new scenario for the formation of multiple polar cap arcs is proposed based on the results from a time-dependent electrodynamic model of polar cap arcs developed by Zhu et al. (1993). The results suggest that the appearance of multiple polar cap arcs may not be due to multiple structures in the magnetospheric source region, but instead, may primarily be determined by the coupled magnetosphere-ionosphere system in which the ionosphere plays an active role. It was found that with the same magnetospheric driver, a strong ionospheric background convection and an ionospheric background Hall conductance in the range of from 0.5 to 1.5 mho favor the occurrence of multiple polar cap arcs. It was also found that the spacing between individual arcs is highly correlated with the ratio of the enhanced Hall and Pedersen conductances (i.e., the hardness of the primary magnetospheric precipitation), which is allowed to vary in the model runs. A harder precipitation causes a wider spacing. In the modelings, the spacing ranged from about 20 kilometers to about 60 kilometers when the precipitation hardness increased.

Zhu, L.

Modeling polar cap F-region patches using time varying convection

Creation of polar cap F-region patches are simulated for the first time using two independent physical models of the high latitude ionosphere. The patch formation is achieved by temporally varying the magnetospheric electric field (ionospheric convection) input to the models. The imposed convection variations are comparable to changes in the convection that result from changes in the B(y) IMF component for southward IMF. Solar maximum-winter simulations show that simple changes in the convection pattern lead to significant changes in the polar cap plasma structuring. Specifically, in winter, as enhanced dayside plasma convects into the polar cap to form the classic tongue-of-ionization the convection changes produce density structures that are indistinguishable from the observed patches.

Sojka, J. J.

A time-dependent model of polar cap arcs

In the present 2D time-dependent model of polar cap arcs, the electrodynamics of the polar cap arcs are treated self-consistently in the frame of the coupled magnetosphere-ionosphere system. The focus of this work is a description of the features of the spatial structure and temporal evolution of the polar cap arcs. Results indicate that the time constant for the formation of the polar cap arcs is around 10 min; an initial single-arc precipitation pattern associated with a polar cap arc tends to split into multiple precipitation regions, yielding a multiple polar cap arc structure. Strong downward field-aligned currents can develop near the intensive upward field-aligned currents and form a pair structure of the field-aligned current in the polar cap arcs. The model predicts plasma flow crossing the polar cap arcs, but its amplitude is small and its characteristic time scale is much larger than the time constant for the formation of the polar cap arcs. While polar cap arcs are developing, the associated current system tends to become further localized.

Zhu, L.

Parameterized study of the ionospheric modification associated with sun-aligned polar cap arcs

The local ionospheric modification that is due to a generalized steady state solar aligned (SA) arc structure is addressed. For a representative set of SA arc parameters which includes both convection and precipitation, emphasis is placed on the modification by SA polar cap arcs upon the F region electron density and the height integrated conductivity. At low fluxes and low characteristic energies, SA polar cap arcs have the most pronounced relative effect at F region altitudes in darkness for winter solar minimum conditions. The absolute enhancement in summer solar minimum and winter solar maximum is equivalent to that of winter solar minimum, but the higher ambient densities make the relative enhancement less. The TEC enhancement associated with an SA arc may be used to indicate the degree of plasma cross flow across the arc.

Crain, D. J.

Influence of horizontal inhomogeneity in the ionosphere on the reflection of Alfven waves

A study of the reflection of Alfven waves at a horizontally inhomogeneous ionosphere has been carried out. In this study, the Alfven speed above the ionosphere is assumed to be uniform and the ionosphere is treated as a height-integrated conducting slab. Analytical and numerical results indicate that the horizontal nonuniformity of the ionospheric conductivity can lead to a rotation of the reflected wave fields and cause field-aligned currents that originate in the ionosphere. A strong conductivity nonuniformity on the direction perpendicular to the incident wave field, large Hall to Pedersen conductivity ratios, and low conductivity values lead to a large rotation of the reflected wave field in the range from a few to 40 deg. The implications of the theoretical results for some ionospheric phenomena are presented, including comments on the Harang discontinuity and sun-aligned arcs.

Zhu, L.

Ionospheric simulation compared with Dynamics Explorer observations for November 22, 1981

The present study uses an extensive DE-2 data base to both constrain inputs to a time-dependent ionospheric model (TDIM) for a simulation of the ionosphere and then check the simulated densities. The investigation was carried out for both a quiet period and a substorm period. The quiet-day study produced very good agreement between modeled and observed electron densities in the topside ionosphere with two significant exceptions: First, across the polar region the DE-2 LANG densities showed fine structure in addition to the overall regional density morphology. Second, a surprising discrepancy arose in the presunrise and midlatitude trough. The TDIM densities were an order of magnitude lower than those observed by DE-2. The substorm study showed remarkably good agreement with the observed densities.

Sojka, J. J.

Latest developments in the display of large-scale ionospheric and thermospheric data sets

Over the past decade, data base sizes have continually increased and will continue to do so in the future. This problem of size is further compounded because the trend in present-day studies is to use data from many different locations and different instruments and then compare it with data from global scale physical models. The latter produce data bases of comparable if not even larger size. Much of the data can be viewed as 'image' time sequences and is most readily viewed on color display terminals. These data sets reside in national or owner-generated data bases linked together by computer networks. As the size increases, just moving this data around, taking 'quick-looks' at the data, or even storing it locally become severe problems compromising the scientific return from the data. Is the present-day technology with these analysis techniques being used in the best way? What are the prospects for reducing the storage and transmission size of the data sets? Examples of such problems and potential solutions are described in this paper.

Sojka, J. J.

Comparison of the USU ionospheric model with the UCL-Sheffield coupled thermospheric-ionospheric model

Several physical models of the high-latitude ionosphere have been developed that describe the time-dependent evolution of the E- and F-region plasma density. The models require a variety of inputs, including solar EUV fluxes, magnetospheric convection, auroral precipitation, and neutral atmosphere. Of specific relevance to this study is how the neutral atmosphere is incorporated into the ionospheric models. For the USU ionospheric model, the neutral atmosphere is the MSIS 1986 empirical model, while for the UCL-Sheffield coupled thermospheric-ionospheric model the neutral atmosphere is computed simultaneously with the ionosphere. Both models were run for similar solar and magnetospheric conditions (solar maximum, moderate geomagnetic activity, and winter solstice). Solar maximum conditions ensured a strong coupling between the ionosphere and thermosphere, which provided the possibility of a large ionospheric difference between the two physical models. This was further enhanced by choosing winter conditions so that the densities were not dominated by sunlight. The comparison of the two models indicated that both models predict the same morphological features with similar ionospheric densities, generally within about 30 percent.

Sojka, J. J.

Approaches to ionospheric modelling, simulation and prediction

The ionosphere is a complex, multispecies, anisotropic medium that exhibits a significant variation with time, space, season, solar cycle, and geomagnetic activity. In recent years, a wide range of models have been developed in an effort to describe ionospheric behavior. The modeling efforts include: (1) empirical models based on extensive worldwide data sets; (2) simple analytical models for a restricted number of ionospheric parameters; (3) comprehensive, 3D, time-dependent models that require supercomputers; (4) spherical harmonic models based on fits to output obtained from comprehensive numerical models; and (5) ionospheric models driven by real-time magnetospheric inputs. In an effort to achieve simplicity, some of the models have been restricted to certain altitude or latitude domains, while others have been restricted to certain ionospheric parameters, such as the F-region peak density, the auroral conductivity, and the plasma temperatures. The current status of the modeling efforts is reviewed.

Schunk, R. W.

Field-aligned current associated with a distorted two-cell convection pattern during northward interplanetary magnetic field

The influence of the ionospheric conductance on the field-aligned current associated with a distorted two-cell convection pattern during northward IMF was investigated using the Heppner-Maynard (1987) convection model and the Utah State University conductivity model described by Rasmussen and Schunk (1987). Results show that the variation of the ionospheric conductivity distribution can significantly affect the features of the field-aligned current for northward IMF, where matching or mismatching between the conductance gradient and the convection electric field plays a key role. It was found that the increase of the field-aligned current in the polar cap observed during summer is mainly due to the increasing contribution from the Pedersen current, and that the increase of the field-aligned current in both the oval region and the evening-midnight sector during the active aurora period is mainly due to the increasing contribution from the Hall current.

Zhu, L.