Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forward modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

Deriving cloud droplet number concentration from surface-based remote sensors with an emphasis on lidar measurements

Abstract. Given the importance of constraining cloud droplet number concentrations (Nd) in low-level clouds, we explore two methods for retrieving Nd from surface-based remote sensing that emphasize the information content in lidar measurements. Because Nd is the zeroth moment of the droplet size distribution (DSD), and all remote sensing approaches respond to DSD moments that are at least 2 orders of magnitude greater than the zeroth moment, deriving Nd from remote sensing measurements has significant uncertainty. At minimum, such algorithms require the extrapolation of information from two other measurements that respond to different moments of the DSD. Lidar, for instance, is sensitive to the second moment (cross-sectional area) of the DSD, while other measures from microwave sensors respond to higher-order moments. We develop methods using a simple lidar forward model that demonstrates that the depth to the maximum in lidar-attenuated backscatter (Rmax⁡) is strongly sensitive to Nd when some measure of the liquid water content vertical profile is given or assumed. Knowledge of Rmax⁡ to within 5 m can constrain Nd to within several tens of percent. However, operational lidar networks provide vertical resolutions of > 15 m, making a direct calculation of Nd from Rmax⁡ very uncertain. Therefore, we develop a Bayesian optimal estimation algorithm that brings additional information to the inversion such as lidar-derived extinction and radar reflectivity near the cloud top. This statistical approach provides reasonable characterizations of Nd and effective radius (re) to within approximately a factor of 2 and 30 %, respectively. By comparing surface-derived cloud properties with MODIS satellite and aircraft data collected during the MARCUS and CAPRICORN II campaigns, we demonstrate the utility of the methodology.

54 ENVIRONMENTAL SCIENCES↗

Radar-based Snow Property Profile Retrievals for NSA

The product contains estimated snow microphysical properties retrieved primarily from vertically-profiling radar observations using a priori constraints that are dependent on atmospheric state. The product provides data at 1-minute time resolution, including snowfall rate at the surface and vertical profiles at 300-meter resolution of snow particle size distribution parameters and snow water contents. Estimated uncertainties for the snowfall rates and snow water contents are also provided. The product is obtained using an optimal estimation retrieval in which radar forward model and a priori state characteristics are based on expected snow particle microphysical and radar scattering properties evaluated from field experiment data.

54 ENVIRONMENTAL SCIENCES↗

MLP-NN vs Gauss-Newton files

-Input files for randomly selected subset data (30 instances) from primary dipole-dipole forward modeling data.-Inversion results in surfer grid format as well as in .DAT format.-Scatterplots for MLP-NN vs Gauss-Newton present in the excel sheet.

58 GEOSCIENCES↗

The deep lunar electrical conductivity profile - Structural and thermal inferences

Simultaneous lunar surface and orbital magnetometer records are reexamined, to ascertain intervals which may be suitable for measuring lunar inductive response in the solar wind and terrestrial magnetosheath. Power spectral estimates of the response tangent to the lunar surface, defined in terms of transfer and gain functions, are obtained for the 0.0001-0.01 Hz frequency range. The maximum consistency of estimates from different time intervals is found when the initial analysis is limited to the tangential direction of maximum incident power, or that direction in which the ratio of signal to background noise is greatest. Spherically symmetric plasma confinement theory is used in the interpretation of transfer function data, by way of forward model calculations, under the assumption of continuous electrical conductivity increase with depth. Results are presented for internal electrical conductivity profile, metallic core radius, and selenotherm limits.

Hood, L. L.↗

An upper bound on the radius of a highly electrically conducting lunar core

Parker's (1980) nonlinear inverse theory for the electromagnetic sounding problem is converted to a form suitable for analysis of lunar day-side transfer function data by: (1) transforming the solution in plane geometry to that in spherical geometry; and (2) transforming the theoretical lunar transfer function in the dipole limit to an apparent resistivity function. The theory is applied to the revised lunar transfer function data set of Hood et al. (1982), which extends in frequency from 10 to the -5th to 10 to the -3rd Hz. On the assumption that an iron-rich lunar core, whether molten or solid, can be represented by a perfect conductor at the minimum sampled frequency, an upper bound of 435 km on the maximum radius of such a core is calculated. This bound is somewhat larger than values of 360-375 km previously estimated from the same data set via forward model calculations because the prior work did not consider all possible mantle conductivity functions.

Hobbs, B. A.↗

Study of the crust and mantle using magnetic surveys by Magsat and other satellites

This article summarizes some of the methodology and results for studies of the earth's mantle and crust. Mantle conductivity studies can be made either by studying signals impressed on the earth from outside, e.g., the ionosphere or magnetosphere, or by studing signals originating in the core and transmitted through the mantle. Crustal field studies begin with a careful selection of the data and subsequent removal of core and external fields by some sort of filtering. Average maps from different local times sometimes differ, presumably due to the remaining presence of fields of external origin. Several techniques for further filtering are discussed. Where large-area aeromagnetic maps are available, crustal maps derived from satellite data can be compared with upward continued data. In general, the comparisons show agreement, with some differences, particularly in and near the auroral belts. The satellite data are further reduced by various methods of inverse and forward modeling, sometimes including reduction to the pole. These techniques are generally unstable at the equator. Common methods of stabilizing the inversions include principle components analysis and ridge regression. Because of the presence of the core field, the entire crustal contribution from the field is not known. Also, there is a basic nonuniqueness to the inverse solutions. Nevertheless, magnetizations that are interpretable can be derived.

Langel, Robert A.↗

Magsat magnetic anomaly contrast across Labrador Sea passive margins

Many passive margins not complicated by nearby anomalous crustal structure have satellite elevation crustal magnetic anomaly contrasts across them that are recognizable in reduced-to-pole versions of the Magsat and POGO data. In the Labrador Sea region this contrast is particularly well developed with strong positive anomalies overlying the continental crust of Greenland and eastern Canada and prominent negative anomalies situated over the Labrador Sea and Baffin Bay. In this work, forward modeling of the large-scale crustal bodies in this region (continental, oceanic, passive margin, several anomalous structures) was used to show that the Magsat anomaly contrast is due simply to the change in crustal susceptibility and thickness at the continental/oceanic crustal transition. Because the thickness varies more than the average susceptibility from continental to oceanic crust, the strong anomaly contrast is essentially an edge effect due mostly to the change in crustal structure.

Bradley, Lauren M.↗

Ground-based microwave monitoring of stratospheric ozone

A microwave instrument developed for operational measurements of ozone for the Network for Detection of Stratospheric Change is discussed. The instrument observes two spectral lines near 3-mm wavelength with a bandwidth of 630 MHz, allowing profile retrieval from 20 to 70 km. The observing technique and calibration procedures are described. The measurement forward model and retrieval algorithm are formulated. Preliminary comparisons with a colocated ground-based lidar and the SAGE II instrument are presented. The measurements are shown to typically agree to within 5 to 10 percent.

Parrish, A.↗

Detection of stratospheric ozone trends by ground-based microwave observations

We present an assessment of the trend detection capabilities of ground-based microwave measurements of ozone, focusing on the 20-60 km region. Uncertainty due to variable forward model and calibration errors is 3-5 percent; due to temperature profile errors, 1-2 percent; and due to spectral measurement errors, 1-3 percent. The measurement resolution is 8-10 km below 40 km, and increases to 17 km at 60 km. The net error for trend detection is 4-6 percent. Comparison with other measurements suggests interannual variations can be measured to 3-5 percent. We conclude a trend of 0.5 percent per year would be detected in about 10 years.

Connor, Brian J.↗

Chicxulub impact basin: Gravity characteristics and implications for basin morphology and deep structure

The K-T-aged Chicxulub Impact Structure is buried beneath the Tertiary carbonate rocks of the Northern Yucatan Platform. Consequently its morphology and structure are poorly understood. Reprocessed Bouguer (onshore) and Free Air (offshore) gravity data over Northern Yucatan reveal that Chicxulub may be a 200-km-diameter multi-ring impact basin with at least three concentric basin rings. The positions of these rings follow the square root of 2 spacing rule derived empirically from analysis of multi-ring basins on other planets indicating that these rings probably correspond to now-buried topographic basin rings. A forward model of the gravity data along a radial transect from the southwest margin of the structure indicates that the Chicxulub gravity signature is compatible with this interpretation. We estimate the basin rim diameter to be 204 +/- 16 km and the central peak ring diameter (D) is 104 +/- 6 km.

Sharpton, Virgil L.↗

Ozone vertical distribution retrieval from ground-based high resolution infrared solar spectra

A practical procedure for the retrieval of ozone vertical profiles from ground-based high resolution Fourier transform infrared solar spectra has been developed. The analysis is based on a multilayer line-by-line forward model and a semi-empirical version of the optimal estimation inversion method of Rodgers. The 1002.6-1003.2 cm(exp -1) spectral interval has been selected for the analysis on the basis of synthetic spectrum calculations. This interval contains numerous ozone lines covering a range of intensities and providing retrieval sensitivity from ground level to about 35 km. Characterization of the method and an error analysis have been performed. For a spectral resolution of 0.05-0.01 cm(exp -1) and a signal-to-noise ratio greater than or equal to 100 the retrieval is stable with a vertical resolution of approximately 5 km attainable near the surface degrading to approximately 10 km in the stratosphere. Synthetic spectra studies show that the a priori profile and weak constraints selected for the retrievals introduce no significant biases for a wide range of ozone profiles.

Pougatchev, N. S.↗

Spectroscopy of Solid State Laser Materials

We retrieved the vertical distribution of ozone from a series 0.005-0.013/cm resolution infrared solar spectra recorded with the McMath Fourier Transform spectrometer at the Kitt Peak National Solar Observatory. The analysis is based on a multi-layer line-by-line forward model and a semi-empirical version of the optimal estimation inversion method by Rodgers. The 1002.6-1003.2/cm spectral interval has been selected for the analysis on the basis of synthetic spectrum calculations. The characterization and error analysis of the method have been performed. It was shown that for the Kitt Peak spectral resolution and typical signal-to-noise ratio (greater than or equal to 100) the retrieval is stable, with the vertical resolution of approximately 5 km attainable near the surface degrading to approximately 10 km in the stratosphere. Spectra recorded from 1980 through 1993 have been analyzed. The retrieved total ozone and vertical profiles have been compared with total ozone mapping spectrometer (TOMS) satellite total columns for the location and dates of the Kitt Peak Measurements and about 100 ozone ozonesoundings and Brewer total column measurements from Palestine, Texas, from 1979 to 1985. The total ozone measurements agree to +/- 2%. The retrieved profiles reproduce the seasonally averaged variations with altitude, including the ozone spring maximum and fall minimum measured by Palestine sondes, but up to 15% differences in the absolute values are obtained.

Buoncristiani, A. M.↗

Magnetospheric Image Unfolding

The Grant was a three year grant funded under the Space Physics Supporting Research and Technology and Suborbital Program. Our objective was to develop automated techniques needed to unfold or "invert" global images of the magnetospheric ion populations obtained by the new magnetospheric imaging techniques (ENA, EUV) in anticipation of future missions such as the Magnetospheric Imager and, now, IMAGE. Our focus on the present three year grant is to determine the degree to which such images can quantitatively constrain the global electromagnetic properties of the magnetosphere. In a previous three year grant period we successfully automated a forward modeling inversion algorithm, demonstrated that these inversions are robust in the face of realistic instrumental considerations such as counting statistics and backgrounds, applied error analysis techniques to the extracted parameters using variational procedures, implemented very realistic magnetospheric test images to test the inversion algorithms using the Rice University Magnetospheric Specification Model, and began the process of generating parametric models with the flexibility to handle the realistic magnetospheric images (e.g. Roelof et al, 1992; 1993). Our plan for the present 3 year grant period was to complete the development of the inversion tools needed to handle realistic magnetospheric images, assess the degree to which global electrodynamics is quantitatively constrained by ENA images of the magnetosphere, and bring the inversion of EUV images up to the maturity that we will have achieved for the ENA imaging. Below the accomplishments of our three year effort are present followed by a list of our presentations and publications. The accomplishments of all three years are presented here, and thus some of these items appeared on interim progress reports.

Source record↗

Global Geodetic Signatures of the Antarctic Ice Sheet

Four scenarios of present day Antarctic ice sheet mass change are developed from comprehensive reviews of the available glaciological and oceanographic evidence. The gridded scenarios predict widely varying contributions to secular sea level change xi ranging from -1.1 to 0.45 mm/yr, and predict polar motion m and time-varying low-degree gravitational coefficients J(sub l) that differ significantly from earlier estimates. A reasonably linear relationship between the rate of sea level change from Antarctica xi(sub A) and the predicted Antarctic J(sub l) is found for the four scenarios. This linearity permits a series of forward models to be constructed that incorporate the effects of ice mass changes in Antarctica, Greenland, and distributed smaller glaciers, as well as postglacial rebound (assuming the ICE-3G deglaciation history), with the goal of obtaining optimum reconciliation between observed constraints on J(sub l) and sea level rise xi. Numerous viable combinations of lower mantle viscosity and hydrologic sources are found that safely "observed" in the range of 1 to 2-2.5 mm/yr and observed J(sub l) for degrees 2, 3, and 4. In contrast, rates of global sea level rise above 2.5 mm/yr are inconsistent with available J(sub l) observations. The successful composite models feature a pair of lower mantle viscosity solutions arising from the sensitivity of J(sub l) to glacial rebound. The paired values are well separated at xi = 1 mm/yr, but move closer together as xi is increased, and, in fact, merge around xi = 2 - 2.5 mm/yr, revealing an intimate relation between xi and preferred lower mantle viscosity. This general pattern is quite robust and persists for different J(sub l) solutions, for variations in source assumptions, and for different styles of lower mantle viscosity stratification. Tighter J(sub l) constraints for l greater than 2 may allow some viscosity stratification schemes and source assumptions to be excluded in the future. For a given total observed xi the sea level rise from Antarctica xi(sub A) is tightly constrained and ranges from 0 to + 1 mr/yr (corresponding to an ablating ice sheet) as estimates of xi are raised from 1 to 2.5 mm/yr. However, when the degree 3 zonal harmonic constraint is removed, the solutions show little sensitivity to Antarctic mass balance, emphasizing the need for a well determined odd-degree secular zonal harmonic for determining polar ice mass balance.

James, Thomas S.↗

The Surface Displacement Field of the November 8, 1997, Mw7.6 Manyi (Tibet) Earthquake Observed with ERS InSAR Data

ERS2 radar data acquired before and after the Mw7.6, Manyi (Tibet) earthquake of November 8, 1997, provide geodetic information about the surface displacement produced by the earthquake in two ways. (1) The sub-pixel geometric adjustment of the before and after images provides a two dimensional offset field with a resolution of approx, 1m in both the range (radar line of sight) and azimuth (satellite track) directions. Comparison of offsets in azimuth and range indicates that the displacement along the fault is essentially strike-slip and in a left-lateral sense. The offset map reveals a relatively smooth and straight, N78E surface rupture that exceeds 150 km in length, consistent with the EW plane of the Harvard CMT solution. The rupture follows the trace of a quaternary fault visible on satellite imagery (Tapponnier and Molnar, 1978; Wan Der Woerd, pers. comm.). (2) Interferometric processing of the SAR data provides a range displacement map with a precision of a few millimeters. The slip distribution along the rupture reconstructed from the range change map is a bell-shaped curve in the 100-km long central section of the fault with smaller, local maxima near both ends. The curve shows that the fault slip exceeds 2.2 m in range, or 6.2 in strike-slip, along a 30-km long section of the fault and remains above 1 m in range, approx. 3 m strike-slip, along most of its length. Preliminary forward modeling of the central section of the rupture, assuming a uniform slip distribution with depth, indicates that the slip occur-red essentially between 0 and the depth of 10 km, consistent with a relatively shallow event (Velasco et al., 1998).

Peltzer, G.↗

A Systematic Error Correction Method for TOVS Radiances

Treatment of systematic errors is crucial for the successful use of satellite data in a data assimilation system. Systematic errors in TOVS radiance measurements and radiative transfer calculations can be as large or larger than random instrument errors. The usual assumption in data assimilation is that observational errors are unbiased. If biases are not effectively removed prior to assimilation, the impact of satellite data will be lessened and can even be detrimental. Treatment of systematic errors is important for short-term forecast skill as well as the creation of climate data sets. A systematic error correction algorithm has been developed as part of a 1D radiance assimilation. This scheme corrects for spectroscopic errors, errors in the instrument response function, and other biases in the forward radiance calculation for TOVS. Such algorithms are often referred to as tuning of the radiances. The scheme is able to account for the complex, air-mass dependent biases that are seen in the differences between TOVS radiance observations and forward model calculations. We will show results of systematic error correction applied to the NOAA 15 Advanced TOVS as well as its predecessors. We will also discuss the ramifications of inter-instrument bias with a focus on stratospheric measurements.

Joiner, Joanna↗

Charge Exchange Contribution to the Decay of the Ring Current, Measured by Energetic Neutral Atoms (ENAs)

In this paper we calculate the contribution of charge exchange to the decay of the ring current. Past works have suggested that charge exchange of ring current protons is primarily responsible for the decay of the ring current during the late recovery phase, but there is still much debate about the fast decay of the early recovery phase. We use energetic neutral atom (ENA) measurements from Polar to calculate the total ENA energy escape. To get the total ENA escape we apply a forward modeling technique, and to estimate the total ring current energy escape we use the Dessler-Parker-Sckopke relationship. We find that during the late recovery phase of the March 10, 1998 storm ENAs with energies greater than 17.5 keV can account for 75% of the estimated energy loss from the ring current. During the fast recovery the measured ENAs can only account for a small portion of the total energy loss. We also find that the lifetime of the trapped ions is significantly shorter during the fast recovery phase than during the late recovery phase, suggesting that different processes are operating during the two phases.

Jorgensen, A. M.↗