Ray tracing studies for the topside ionosphere.
Nonvertical propagation effects on high altitude topside ionosphere sounder data reduction, using ray tracing
SEARCH · Search NASA
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.
Nonvertical propagation effects on high altitude topside ionosphere sounder data reduction, using ray tracing
Two-dimensional ray tracing of Z-mode radiation from sources in the auroral zone has been performed. The Z-mode waves are assumed to be excited by the cyclotron maser mechanism at frequencies near the electron gyrofrequency and at wave normal angles near 90 deg. The emission is found to propagate primarily perpendicular to the magnetic field line (horizontal propagation) for all sources considered, and thus not to substantially lower altitudes. The frequency bandwidth estimated from the ray paths is much less than that observed by DE-1. It is suggested that the wave source needs to be in the region of field-aligned current in order for the cyclotron maser mechanism to explain the broad frequency bandwidth of Z-mode emissions.
Using a model of Jupiter's magnetized plasma environment, radiation raypaths were calculated with a three-dimension ray tracing program. It is assumed that energetic particles produce the emission in the planet's auroral zone at frequencies just above the electron gyrofrequencies. This radiation is generated in narrow sheets defined by the angle of a ray with respect to the magnetic field line. By specifying the source position: latitude, longitude, and radial distance from the planet, signatures in the spectrum of frequency versus time seen by Voyager 1 and 2 were duplicated. The frequency range and the curvature of the decametric arcs in these dynamic spectra are the result of the geometry of the radiation sheets (imposed by the plasma and by the B-field) and illumination of Voyager 1 and 2 as the rotating magnetosphere mimics a pulsar.
We use an analytical fit to an emission lobe profile together with three-dimensional ray tracing to model the broad-banded smooth Uranian kilometric radiation (UKR). We assume the radiation is gyroemission from sources along magnetic field lines. Using an iterative technique that modifies the lobe function and source region, the results are compared to observations at a frequency of 481 kHz. The best-fit calculations are compared to previously published models and to recent ultraviolet (UV) observations.
Radio wave propagation predictions are described in which modern comprehensive theoretical ionospheric models are coupled with ray-tracing programs. In the computer code described, a network of electron density and collision frequency parameters along a band about the great circle path is calculated by specifying the transmitter and receiver geographic coordinates, time, the day number, and the 2800-MHz solar flux. The ray paths are calculated on specifying the frequency, mode, range of elevation angles, and range of azimuth angles from the great circle direction. The current program uses a combination of the Penn State MKI E and F region models and the Mitra-Rowe D and E region model. Application of the technique to the prediction of satellite to ground propagation and calculation of oblique incidence propagation paths and absorption are described. The implications of the study to the development of the next generation of ionospheric models are discussed.
Raypaths for decametric wavelength radiation in Jupiter's magnetosphere were calculated. The model-dependent raypaths with the Voyager observations were compared. Characteristics of the source regions and the influence of propagation effects were deduced. A three dimensional ray tracing program was employed to calculate the raypaths. Families of rays were launched at particular angles with respect to the magnetic field lines to generate conical sheets of radiation for various frequencies and various source locations. As the planet's magnetic field rotates, these warped sheets of radiation sweep past the observer, producing signatures in frequency versus time plots. These signatures match some of those found in the Voyager data. The greatest propagation effects occur in and around the source regions in the Io auroral oval.
Errors may be introduced in satellite laser ranging data by atmospheric refractivity. Ray tracing data have indicated that horizontal refractivity gradients may introduce nearly 3-cm rms error when satellites are near 10-degree elevation. A correction formula to compensate for the horizontal gradients has been developed. Its accuracy is evaluated by comparing it to refractivity profiles. It is found that if both spherical and gradient correction formulas are employed in conjunction with meteorological measurements, a range resolution of one cm or less is feasible for satellite elevation angles above 10 degrees.
Generating graphics to faithfully represent information can be a computationally intensive task. A way of using the Massively Parallel Processor to generate images by ray tracing is presented. This technique uses sort computation, a method of performing generalized routing interspersed with computation on a single-instruction-multiple-data (SIMD) computer.
The physical arc length difference between the red and blue beam of a dual-color laser generates an error term in determining Delta-n as well as n (the refractive index) and D (the line length). Numerical ray trace examples and theoretical approximations show that the resulting relative error Delta-D/D increases as D-squared. Error reduction by a factor of nearly 10 is possible by using the pressure, temperature, and humidity of both endpoints to calculate n. The present study results from a desire to improve the precision of monitoring the local and regional position of the Lunar Laser Ranging Observatory atop Haleakala, Maui, Hawaii; the monitoring of the earth's crustal deformation and strain is a primary concern of the study.
A fundamental problem in remote sensing and radiative transfer simulations involving ice clouds is the ability to compute accurate optical properties for individual ice particles. While relatively simple and intuitively appealing, the conventional geometric-optics method (CGOM) is used frequently for the solution of light scattering by ice crystals. Due to the approximations in the ray-tracing technique, the CGOM accuracy is not well quantified. The result is that the uncertainties are introduced that can impact many applications. Improvements in the Invariant Imbedding T-matrix method (II-TM) and the Improved Geometric-Optics Method (IGOM) provide a mechanism to assess the aforementioned uncertainties. The results computed by the II-TMþIGOM are considered as a benchmark because the IITM solves Maxwell's equations from first principles and is applicable to particle size parameters ranging into the domain at which the IGOM has reasonable accuracy. To assess the uncertainties with the CGOM in remote sensing and radiative transfer simulations, two independent optical property datasets of hexagonal columns are developed for sensitivity studies by using the CGOM and the II-TMþIGOM, respectively. Ice cloud bulk optical properties obtained from the two datasets are compared and subsequently applied to retrieve the optical thickness and effective diameter from Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. Additionally, the bulk optical properties are tested in broadband radiative transfer (RT) simulations using the general circulation model (GCM) version of the Rapid Radiative Transfer Model (RRTMG) that is adopted in the National Center for Atmospheric Research (NCAR) Community Atmosphere Model (CAM, version 5.1). For MODIS retrievals, the mean bias of uncertainties of applying the CGOM in shortwave bands (0.86 and 2.13 micrometers) can be up to 5% in the optical thickness and as high as 20% in the effective diameter, depending on cloud optical thickness and effective diameter. In the MODIS infrared window bands centered at 8.5, 11, and 12 micrometers biases in the optical thickness and effective diameter are up to 12% and 10%, respectively. The CGOM-based simulation errors in ice cloud radiative forcing calculations are on the order of 10Wm(exp 2).
Certain Z-mode wave emissions in the earth's magnetosphere have been identified using the wave spectra and polarization measurements of the DE 1 satellite. Although such emissions accompany the aurora, and thus presumably originate from the evening-sector auroral zone, they are found to occur over much wider ranges of latitude and longitude. Since the predicted cyclotron maser emission at the cyclotron frequency could not have produced waves which travel such great distances, as shown by three-dimensional ray tracing, it is proposed instead that these emissions must originate from lower altitudes within the auroral zone and probably from near the plasma frequency inside the auroral plasma cavity.
Atmospheric turbulence has significant impairments on the operation of Free-Space Optical (FSO) communication systems, in particular temporal and spatial intensity fluctuations at the receiving aperture resulting in power surges and fades, changes in angle of arrival, spatial coherence degradation, etc. The refractive index structure parameter Cn2 is a statistical measure of the strength of turbulence in the atmosphere and is highly dependent upon vertical height. Therefore to understand atmospheric turbulence effects on vertical FSO communication links such as space-to-ground links, it is necessary to specify Cn2 profiles along the atmospheric propagation path. To avoid the limitations on the applicability of classical approaches, propagation simulation through geometrical ray tracing is applied. This is achieved by considering the atmosphere along the optical propagation path as a spatial distribution of spherical bubbles with varying relative refractive index deviations representing turbulent eddies. The relative deviations of the refractive index are statistically determined from altitude-dependent and time-varying temperature fluctuations, as measured by a microwave profiling radiometer. For each representative atmosphere ray paths are analyzed using geometrical optics, which is particularly advantageous in situations of strong turbulence where there is severe wavefront distortion and discontinuity. The refractive index structure parameter is then determined as a function of height and time.
Gravity waves with wave periods of 13 to 15 min and horizontal phase velocities of 90 to 220 m/sec were present in ground-based observations of the upper atmosphere during time periods when tornadoes were occurring and gravity waves with wave periods of 20 to 25 min and horizontal phase velocities of 100 to 200 m/sec were detected when a hurricane was present. Combinations of available neutral atmosphere data and model parameter values were used with a group ray tracing technique in an attempt to locate the sources of these waves. Computed sources of the waves with periods of 13 to 15 min were located within 50 km of the locations where tornadoes touched down from 2 to 4 h later. In the case of the waves with periods of 20 to 25 min it was found that the computed location of the source was roughly where the hurricane would be located 3 h after the waves were excited. The applicability of the present study to a tornado and hurricane warning system is noted.
Sunlight reflected from the earth is, to a certain extent, polarized. Radiometers, such as the MODIS instrument on board the TERRA and AQUA spacecraft, are to a certain extent polarizers. Accurate radiometric measurements must take into account both the polarization state of the scene and the polarization sensitivity of the measuring instrument. The measured polarization characteristics of the MODIS instruments are contained in various radiometric models. Continued use of these radiometric math models, over a number of years, have shown where these models can be improved. Currently a MODIS polarization ray trace model has been created which models the thin film structure on the optical elements. This approach is described and modeled and measured instrument polarization sensitivity results presented.
Five computational methods for solution of the radiative transfer equation in an absorbing-emitting and non-scattering gray medium were compared on a 2 m JP-8 pool fire. The temperature and absorption coefficient fields were taken from a synthetic fire due to the lack of a complete set of experimental data for fires of this size. These quantities were generated by a code that has been shown to agree well with the limited quantity of relevant data in the literature. Reference solutions to the governing equation were determined using the Monte Carlo method and a ray tracing scheme with high angular resolution. Solutions using the discrete transfer method, the discrete ordinate method (DOM) with both S(sub 4) and LC(sub 11) quadratures, and moment model using the M(sub 1) closure were compared to the reference solutions in both isotropic and anisotropic regions of the computational domain. DOM LC(sub 11) is shown to be the more accurate than the commonly used S(sub 4) quadrature technique, especially in anisotropic regions of the fire domain. This represents the first study where the M(sub 1) method was applied to a combustion problem occurring in a complex three-dimensional geometry. The M(sub 1) results agree well with other solution techniques, which is encouraging for future applications to similar problems since it is computationally the least expensive solution technique. Moreover, M(sub 1) results are comparable to DOM S(sub 4).
Lindzen (1981, 1984) has considered the effects of monochromatic, steady gravity waves propagating into the upper atmosphere. These waves reach such large amplitudes in the mesosphere that they become convectively unstable. The wave is effectively dissipated by the convection, and the pseudomomentum carried by the wave is transferred to the basic flow. The net effect is to accelerate the background flow to the phase speed of the breaking gravity wave. Schoeberl et al. (1983) have discussed modifications to Lindzen's parameterization. The present paper has the objective to develop a type of parameterization scheme for wave breaking that, in principle, can handle lateral wave propagation. Ray tracing is used to follow the gravity wave packets through varying wind conditions. The theory considered uses both the propagation properties of the wave packet and the concept of the conservation of wave action density to determine the local amplitude and position of the wave packet.
During the 1986 Project FIRE (First International Satellite Cloud Climatology Project Regional Experiment) field campaign, four 22 deg halo-producing cirrus clouds were studied jointly from a ground-based polarization lidar and an instrumented aircraft. The lidar data show the vertical cloud structure and the relative position of the aircraft, which collected a total of 84 slides by impaction, preserving the ice crystals for later microscopic examination. Although many particles were too fragile to survive impaction intact, a large fraction of the identifiable crystals were columns and radial bullet rosettes, with both displaying internal cavitations and radial plate-column combinations. Particles that were solid or displayed only a slight amount of internal structure were relatively rare, which shows that the usual model postulated by halo theorists, i.e., the randomly oriented, solid hexagonal crystal, is inappropriate for typical cirrus clouds. With the aid of new ray-tracing simulations for hexagonal hollow-ended column and bullet-rosette models, we evaluate the effects of more realistic ice-crystal structures on halo formation and lidar depolarization and consider why the common halo is not more common in cirrus clouds.
Results of an ongoing effort to quantify the role turbulence in scattering sound in jets are reported. Using a direct numerical simulation database to provide the flow data, ray paths traced through the mean flow are compared with those traced through the actual time evolving turbulent flow. Significant scattering by the turbulence is observed. The most notable effect is that upstream traveling waves that are trapped in the potential core by the mean flow, which acts as a wave guide, easily escape in the turbulent flow. A crude statistical estimate based on ray number density suggests that directivity is modified by the turbulence, but no rigorous treatment of non-uniformities in the high-frequency approximation is attempted.