Search NASA⌕ Search

Engineering topics

Richards, P. G.

Publications and source records attributed to Richards, P. G..

At least 19 records

Achievements and Prospects of Global Broadband Seismographic Networks After 30 Years of Continuous Geophysical Observations

Global seismographic networks (GSNs) emerged during the late nineteenth and early twentieth centuries, facilitated by seminal international developments in theory, technology, instrumentation, and data exchange. The mid- to late-twentieth century saw the creation of the World-Wide Standardized Seismographic Network (1961) and International Deployment of Accelerometers (1976), which advanced global geographic coverage as seismometer bandwidth increased greatly allowing for the recording of the Earth's principal seismic spectrum. The modern era of global observations and rapid data access began during the 1980s, and notably included the inception of the GEOSCOPE initiative (1982) and GSN (1988). Through continual improvements, GEOSCOPE and the GSN have realized near-real time recording of ground motion with state-of-art data quality, dynamic range, and timing precision to encompass 180 seismic stations, many in very remote locations. Data from GSNs are increasingly integrated with other geophysical data (e.g., space geodesy, infrasound and Interferometric Synthetic Aperture Radar). Globally distributed seismic data are critical to resolving crust, mantle, and core structure; illuminating features of the plate tectonic and mantle convection system; rapid characterization of earthquakes; identification of potential tsunamis; global nuclear test verification; and provide sensitive proxies for environmental changes. As the global geosciences community continues to advance our understanding of Earth structure and processes controlling elastic wave propagation, GSN infrastructure offers a springboard to realize increasingly multi-instrument geophysical observatories. Here, we review the historical, scientific, and monitoring heritage of GSNs, summarize key discoveries, and discuss future associated opportunities for Earth Science.

58 GEOSCIENCES↗

On the Relative Importance of Convection and Temperature on the Behavior of the Ionosphere in North American during January 6-12, 1997

Measurements from a network of digisondes and an incoherent scatter radar In Eastern North American For January 6-12, 1997 have been compared with the Field Line Interhemispheric Plasma (FLIP) model which now includes the effects of electric field convective. With the exception of Bermuda, the model reproduces the daytime electron density very well most of the time. As is typical behavior for winter solar minimum on magnetically undisturbed nights, the measurements at Millstone Hill show high electron temperatures before midnight followed by a rapid decay, which is accompanied by a pronounced density enhancement in the early morning hours. The FLIP model reproduces the nighttime density enhancement well, provided the model is constrained to follow the topside electron temperature and the flux tube is full. Similar density enhancements are seen at Goose Bay, Wallops Island and Bermuda. However, the peak height variation and auroral images indicate the density enhancements at Goose Bay are most likely due to particle precipitation. Contrary to previously published work we find that the nighttime density variation at Millstone Hill is driven by the temperature behavior and not the other way around. Thus, in both the data and model, the overall nighttime density is lowered and the enhancement does not occur if the temperature remains high all night. Our calculations show that convections of plasma from higher magnetic latitudes does not cause the observed density maximum but it may enhance the density maximum if over-full flux tubes are convected over the station. On the other had, convection of flux tubes with high temperatures and depleted densities may prevent the density maximum from occurring. Despite the success in modeling the nighttime density enhancements, there remain two unresolved problems. First, the measured density decays much faster than the modeled density near sunset at Millstone Hill and Goose Bay though not at lower latitude stations. Second, we cannot fully explain the large temperatures before midnight nor the sudden decay near midnight.

Richards, P. G.↗

Issues in Quantitative Analysis of Ultraviolet Imager (UV) Data: Airglow

The GGS Ultraviolet Imager (UVI) has proven to be especially valuable in correlative substorm, auroral morphology, and extended statistical studies of the auroral regions. Such studies are based on knowledge of the location, spatial, and temporal behavior of auroral emissions. More quantitative studies, based on absolute radiometric intensities from UVI images, require a more intimate knowledge of the instrument behavior and data processing requirements and are inherently more difficult than studies based on relative knowledge of the oval location. In this study, UVI airglow observations are analyzed and compared with model predictions to illustrate issues that arise in quantitative analysis of UVI images. These issues include instrument calibration, long term changes in sensitivity, and imager flat field response as well as proper background correction. Airglow emissions are chosen for this study because of their relatively straightforward modeling requirements and because of their implications for thermospheric compositional studies. The analysis issues discussed here, however, are identical to those faced in quantitative auroral studies.

Germany, G. A.↗

Analysis of Auroral Morphology: Substorm Precursor and Onset on January 10, 1997

The solar wind interaction with the geomagnetic field is studied using global auroral images obtained by the Ultraviolet Imager (UVI) on Polar. We study the dynam,cs of the poleward and equatorward boundaries of the auroral oval in response to the solar wind IMF on January 10, 1997 using a neural network algorithm to perform an automated morphological analysis. Poleward and equatorward boundaries identified by the algorithm demonstrate a clear growth motion with the southward turning of the IMF and growth and poleward expansion at substorm onset. The area poleward of the oval (polar cap) is found to increase in size coincident with the'southward turning of the IMF Bz component at 0220 UT and peaks at substorm onset at 0334 UT. The area of the oval, however, decreases continuously throughout the period of the polar cap area increase with a slight recovery observed during the substorm onset. These observations are consistent with the concept that magnetospheric dynamics are directly driven by the solar wind-geomagnetic field interactions.

Germany, G. A.↗

Global Auroral Imaging as a Remote Diagnostic of Geospace

Images of the Earth's aurora, taken from space, can be used to examine plasma behavior throughout the magnetospheric regions surrounding the earth. The coupling of the magnetospheric plasmas through the ionosphere are discussed. A summary of past and current imaging technology is given and then specific examples of remote sensing are given using images from the Ultraviolet Imager aboard the POLAR satellite.

Germany, G. A.↗

Compensation for Spherical Geometric and Absorption Effects on Lower Thermospheric Emission Intensities Derived from High Earth Orbit Images

Remote sensing of the atmosphere from high earth orbit is very attractive due to the large field of view obtained and a true global perspective. This viewpoint is complicated by earth curvature effects so that slant path enhancement and absorption effects, small from low earth orbit, become dominant even at small nadir view angles. The effect is further complicated by the large range of local times and solar zenith angles in a single image leading to a modulation of the image intensity by a significant portion of the diurnal height variation of the absorbing layer. The latter effect is significant in particular for mesospheric, stratospheric and auroral emissions due to their depth in the atmosphere. As a particular case, the emissions from atomic oxygen (130.4 and 135.6 nm) and molecular nitrogen (two LBH bands, LBHS from 140 to 160 nm and LBHL from 160 to 180 nm) as viewed from the Ultraviolet Imager (UVI) are examined. The LBH emissions are of particular interest since LBHS has significant 02 absorption while LBHL does not, In the case of auroral emissions this differential absorption, well examined in the nadir, gives information about the height of the emission and therefore the energy of the precipitating particles. Using simulations of the viewing geometry and images from the UVI we examine these effects and obtain correction factors to adjust to the nadir case with a significant improvement of the derived characteristic energy. There is a surprisingly large effect on the images from the 02 diurnal layer height changes. An empirical compensation to the nadir case is explored based on the local nadir and local zenith angles for each portion of the image. These compensations are demonstrated as applied to the above emissions in both auroral and dayglow images and compared to models. The extension of these findings to other instruments, emissions and spectral regions is examined.

Swift, W.↗

A modified thermal conductivity for low density plasma magnetic flux tubes

In response to inconsistencies which have arisen in results from a hydrodynamic model in simulation of high ion temperature (1-2 eV) observed in low density, outer plasmasphere flux tubes, we postulate a reduced thermal conductivity coefficient in which only particles in the loss cone of the quasi-collisionless plasma contribute to the thermal conduction. Other particles are assumed to magnetically mirror before they reach the topside ionosphere and therefore not to remove thermal energy from the plasmasphere. This concept is used to formulate a mathematically simple, but physically limiting model for a modified thermal conductivity coefficient. When this modified coefficient is employed in the hydrodynamic model in a case study, the inconsistencies between simulation results and observations are largely resolved. The high simulated ion temperatures are achieved with significantly lower ion temperatures in the topside ionosphere. We suggest that this mechanism may be operative under the limited low density, refilling conditions in which high ion temperatures are observed.

Comfort, R. H.↗

Sensitivity of the 6300 A twilight airglow to neutral composition

The field line interhemispheric plasma (FLIP) model is used to study the 6300 A line intensity measured during three morning twilights from the McDonald Observatory in Texas. The Imaging Spectrometric Observatory (ISO) measured the 6300 A intensity during the winter of 1987 and the spring and summer of 1988. The FLIP model reproduces the measured intensity and its variation through the twilight well on each day using neutral densities from the MSIS-86 empirical model. This is in spite of the fact that different component sources dominate the integrated volume emission rate on each of the days analyzed. The sensitivity of the intensity to neutral composition is computed by varying the N2, O2, and O densities in the FLIP model and comparing to the intensity computed with the unmodified MSIS-86 densities. The ion densities change self-consistently. Thus the change in neutral composition also changes the electron density. The F2 peak height is unchanged in the model runs for a given day. The intensity changes near 100 deg SZA are comparable to within 10% when either (O2), (N2), or (O) is changed, regardless of which component source is dominant. There is strong sensitivity to changes in (N2) when dissociative recombination is dominant, virtually no change in the nighttime (SZA greater than or equal to 108 deg) intensity with (O2) doubled, and sensitivity of over 50% to doubling or halving (O) at night. When excitation by conjugate photoelectrons is the dominant nighttime component source, the relative intensity change with (O) doubled or halved is very small. This study shows the strong need for simultaneous measurements of electron density and of emissions proportional to photoelectron fluxes if the 6300 A twilight airglow is to be used to retrieve neutral densities.

Melendez-Alvira, D. J.↗

Comparisions of modeled N(+), O(+), H(+) and He(+) in the midlatitude ionospherewith mean densities and temperatures from Atmosphere Explorer

In this study, Atmosphere Explorer data and model results for the ion and electron temperature and the density of N(+), O(+), H(+), and He(+) between 120 and 1400 km altitude are compared for two midlatitude ranges (L=2 and L=4), noon and midnight local time, winter and summer, at solar minimum. The data for the heavy atomic ions (O(+) and N(+)) show that their densities are greater at noon than at midnight for a given season and greater in summer than winter for a given local time. There is only a weak latitudinal variation in the density of these ions. The data show that the light ion (H(+) and He(+)) densities are greater at midnight than at noon and are generally greater in winter than summer. There is a strong latitudinal variation of the light ion densities, with the densities decreasing with increasing latitude. The model densities are in good agreement with the AE densities for N(+), O(+), and H(+). Model He(+) densities are lower, by a factor of 2 or more, than the measured densities. Model ion and electron temperatures agree well with the measured temperatures with only a modest increase in plasmapheric heating.

Craven, P. D.↗

Determination of ionospheric conductivities from FUV auroral emissions

The purpose of this paper is to examine the viability of using Lyman-Birge-Hopfeld (LBH) emission ratios to infer auroral conductances and to quantify the strengths and weaknesses of this technique. We show that column-integrated Hall and Pedersen conductances may be determined from a single remote measurement of a pair of auroral LBH emissions, one in the region of strong O2 absorption (1464 A) and one lying outside of this region (1838 A). The dependence of the determined conductivities on incident average energy, total energy flux, and changes in solar and magnetic activity levels is examined. We show that, for eneriges above 5 keV, auroral conductances may be scaled by the square root of the incident energy flux with errors less than about 20%. For lower energies, however, the scaling may deviate significantly from a power of 0.5. We provide appropriate scaling factors as a function of average energy. We also note that the choice of either a Gaussian or Maxwellian distribution can significantly affect low-energy (less than keV) conductance calculations. Finally, we quantify the conditions under which the mean energy of a Gaussian energy distribution may differ significantly from the characteristic energy.

Germany, G. A.↗

Ionoshperic effects of the March 1990 magnetic storm: Comparison of theory and measurement

This paper presents a comparison of the measured and modeled ionospheric response to magnetic storms at Millstone Hill and Arecibo during March 16-23, 1990. Magnetic activity was low until midday Universal Time (UT) on day 18 when Kp reached 6, days 19 and 20 were quiet, but a large storm occurred around midnight UT on day 20 (Kp=7) ad it was a moderately disturbed (Kp=4) for the remainder of the study period. At Millstone Hill, the daytime peak density (Nm F2) showed only a modest 30% decrease in response to the first storm and recovered to prestorm values before the onset of the second storm. The model reproduces the daytime peak electron density well for this period. However, the severe storm on March 20 caused a factor of 4 deplection in electron density, while the model densities were not greatly affected. the inclusion of vibrationally excited nitrogen N2(sup *) in the model was unable to account for the observed large electron density depletions afterward March 20. The storm did not appear to affect the overall magnitude of the electron density at Arecibo very much, but did cause unusual wavelike structure in the peak density and peak height following the storm. The model reproduces the daytime Nm F2 very well for Arecibo, but after sunset the model densities decay too rapidly. This study indicates that successful modeling of severe ionospheric stroms will require better definition of the storm time inputs, especially of the neutral atmosphere.

Richards, P. G.↗

Global observations and modeling of the ionosphere, thermosphere and mesosphere

The Imaging Spectrometric Observatory (ISO) flown on the Atmospheric Laboratory for Applications and Science (ATLAS) 1 mission between 24 Mar 1992 and 2 Apr 1992, acquired a database designed to study several outstanding problems in the ionosphere, thermosphere and mesosphere. In this paper we discuss the goals and preliminary results from three of these studies. To support these studies, the ISO acquired a database of: (1) emissions for the retrieval of neutral and ion densities to test global models of the ionosphere and thermosphere; (2) emissions for the retrieval of mesospheric composition of major and minor constituents needed to test models of the oxygen-hydrogen photochemistry, (3) emissions of the bands of the metastable states of O2, and O(S-1) produced by three-body recombination of O in the mesosphere.

Torr, D. G.↗

F2 peak electron density at Millstone Hill and Hobart: Comparsion of theory and measurement at solar maximum

This paper compares the observed behavior of the (F2) layer of the ionosphere at Millstone Hill and Hobart with calculations from the field line interhemispheric plasma (FLIP) model for solar maximum, solstice conditions in 1990. During the study period the daily F(sub 10.7) index varied by more than a factor of 2 (123 to 280), but the 81-day mean F(sub 10.7) (F(sub 10.7 A)) was almost constant near 190. Calculations were performed with and without the effects of vibrationally excited N2 (N(sup *)(sub 2) which affects the loss rate of atomic oxygen ions. In the case without N(sup *)(sub 2) there is generally good agreement between the model and measurement for the daytime, peak density of the F region (NmF2). Both the model and the measurement show a strong seasonal anomaly with the winter noon densities a factor of 3 to 4 greater than the summer noon densities at Millstone Hill and a factor of 2 greater at Hobart. The seasonal anomaly in the model is caused by changes in the neutral composition as given by the mass spectrometer and incoherent scatter (MSIS) 86 neutral density model. There is generally little or no increase in the observed noon NmF2 as a function of daily F(sub 10.7) except at Millstone Hill in winter. In contrast to the generally good agreement between model and data at noon, the model badly underestimates the density at night at Millstone Hill at all seasons. At Hobart the model reproduces the nighttime density variations well in both winter and summer. The international reference ionosphere (IRI) model generally provides a good representation of the average behavior of noon NmF2 and hmF2 but because the data show a lot of day-to-day variability, there are often large differences. The FLIP model is able to reproduce this variability when hmF2 is specified. The IRI model peak densities are better than the FLIP densities at night, but the IRI model does not represent the Millstone Hill summer data very well at night in 1990.

Richards, P. G.↗

EUVAC: A solar EUV flux model for aeronomic calculations

This paper presents a new solar Extreme Ultraviolet (EUV) flux model for aeronomic calulations (EUVAC), which is based on the measured F74113 solar EUV reference spectrum. The model provides fluxes in the 37 wavelength bins that are in widespread use. This paper also presents cross sections to be used with the EUVAC flux model to calculate photoionization rates. The flux scaling for solar activity is accomplished using a proxy-based on the F10.7 index and its 81-day average together with the measured solar flux variation from the EUVS instrument on the Atmosphere Explorer E satellite. This new model produces 50-575 A integrated EUV fluxes in good agreement with rocket observations. The solar cycle variation of the chromospheric fluxes agrees well with the measured variation of the Lyman alpha flux between 1982 and 1988. In addition, the theoretical photoelectron fluxes, calculated using the new EUV flux model, are in good agreement with the solar minimum photoelectron fluxes from the Atmosphere Explorer E satellite and also with the solar maximum photoelectron fluxes from the Dynamics Explorer satellite. Its relative simplicity coupled with its ability to reproduce the 50-575 A solar EUV flux as well as the measured photoelectron spectrum makes the model well suited for aeronomic applications. However, EUVAC is not designed to accurately predict the solar flux variability for numerous individual lines.

Richards, P. G.↗

New sources for the hot oxygen geocorona

This paper investigates new sources of thermospheric non thermal (hot) oxygen due to exothermic reactions involving numerous minor (ion and neutral) and metastable species. Numerical calculations are performed for low altitude, daytime, winter conditions, with moderately high solar activity and low magnetic activity. Under these conditions we find that the quenching of metastable species are a significant source of hot oxygen, with kinetic energy production rates a factor of ten higher than those due to previously considered O2(+) and NO(+) dissociative recombination reactions. Some of the most significant new sources of hot oxygen are reactions involving quenching of O(+)((sup 2)D), O((sup 1)D), N((sup 2)D), O(+)((sup 2)P) and vibrationally excited N2 by atomic oxygen.

Richards, P. G.↗

Thermospheric O I 844.6-nm emission in twilight

The thermospheric O I 844.6-nm column emission rate was measured over Toronto, a midlatitude station, in the autumn of 1991 using an imaging Fabry-Perot spectrometer. Twilight decay curves were measured on four clear evenings when the solar zenith angle was between 95 degs and 104 degs, giving corresponding column emission rates between 874 R and 130 R at 20 degs elevation angle in the azimuth of the Sun. The expected decay curves were calculated from the field line interhemisperic plasma model assuming only photoelectron impact excitation as the production mechanism with a cross section appropriate to an optically thin atmosphere. The agreement was good when the solar and geomagnetic activity levels were low to moderate, but the emission rate was overestimated during high activity periods. The comparison indicates that the photoelectron impact mechanism with a thin-atmosphere cross section is sufficient to explain the twilight decay of the thermospheric O I 844.6-nm emission.

Bahsoun-Hamade, F.↗

Simultaneous retrieval of the solar EUV flux and neutral thermospheric O, O2, N2, and temperature from twilight airglow

We present a method to retrieve neutral thermospheric composition and the solar EUV flux from ground-based twilight optical measurements of the O(+) ((exp 2)P) 7320 A and O((exp 1)D) 6300 A airglow emissions. The parameters retrieved are the neutral temperature, the O, O2, N2 density profiles, and a scaling factor for the solar EUV flux spectrum. The temperature, solar EUV flux scaling factor, and atomic oxygen density are first retrieved from the 7320-A emission, which are then used with the 6300-A emission to retrieve the O2 and N2 densities. The retrieval techniques have been verified by computer simulations. We have shown that the retrieval technique is able to statistically retrieve values, between 200 and 400 km, within an average error of 3.1 + or - 0.6% for thermospheric temperature, 3.3 + or - 2.0% for atomic oxygen, 2.3 + or - 1.3% for molecular oxygen, and 2.4 + or - 1.3% for molecular nitrogen. The solar EUV flux scaling factor was found to have a retrieval error of 5.1 + or - 2.3%. All the above errors have a confidence level of 95%. The purpose of this paper is to prove the viability and usefulness of the retrieval technique by demonstrating the ability to retrieve known quantities under a realistic simulation of the measurement process, excluding systematic effects.

Fennelly, J. A.↗

Use of FUV auroral emissions as diagnostic indicators

In an earlier study we modeled selected Far Ultraviolet (FUV) auroral emissions (O I(1356 A), N2 Lyman-Birge-Hopfield (LBH) (1464 A), and LBH (1838 A)) to examine the sensitivity of these emissions and their ratios to likely changes in the neutral atmosphere. In this paper we extend that study to examine the dependence of these same emissions and their ratios on the shape of the energy distribution of the auroral electrons. In particular, we wish to determine whether changes in energy spectra might interfere with our determination of the characteristic energy. Modeled column-integrated emissions show relatively small (less than 30%) dependences on the shape and width of the incident energy spectrum, provided the average energy and total energy flux of the energy distribution are held constant. Long-wavelength FUV emissions, which are relatively unaffected by O2 absorption losses, exhibit virtually no dependence on the shape of the incident energy distribution. Changes in ratios of FUV short- to long-wavelength emissions as a function of characteristic energy are much larger than those due to changes in energy distribution. As a result, the determination of characteristic energy using these emission ratios is relatively unambiguous. We also examine the relative intensities of the aurora and the dayglow for various conditions. The intensities of modeled FUV auroral emissions relative to the dayglow emissions are presented as a function of solar zenith angle and incident energy flux. Under certain conditions (energy flux less than or equal 1 erg/sq cm/s and solar zenith angle less than or equal 50 deg) the dayglow will be limiting factor in the detection of weak auroras.

Germany, G. A.↗