Results for “forward modeling” — Page 11

NASA NTRS2012

Limits on the Abundance and Burial Depth of Lunar Polar Ice

The Diviner imaging radiometer experiment aboard the Lunar Reconnaissance Orbiter has revealed that surface temperatures in parts of the lunar polar regions are among the lowest in the solar system. Moreover, modeling of these Diviner data using realistic thermal conductivity profiles for lunar regolith and topography-based illumination has been done, with surprising results. Large expanses of circum-polar terrain appear to have near-subsurface temperatures well below 110K, despite receiving episodic low-angle solar illumination [Paige et al., 2010]. These subsurface cold traps could provide areally extensive reservoirs of volatiles. Here we examine the limits to abundance and burial depth of putative volatiles, based on the signature they would create for orbital thermal and epithermal neutrons. Epithermals alone are not sufficient to break the abundance-depth ambiguity, while thermal neutrons provide an independent constraint on the problem. The subsurface cold traps are so large that even modest abundances, well below that inferred from LCROSS observations, would produce readily detectable signatures in the Lunar Prospector neutron spectrometer data [Colaprete et al., 2010]. Specifically, we forward-model the thermal and epithermal neutron leakage flux that would be observed for various ice concentrations, given the depth at which ice stability begins. The LCROSS results point to a water-equivalent hydrogen abundance (WEH) in excess of 10 wt%, when all hydrogenous species are added together (except for H2, detected by LAMP on LRO [Gladstone et al., 2010]). When such an ice abundance is placed in a layer below the stability depth of Paige et al., the epithermal and thermal neutron leakage fluxes are vastly reduced and very much at odds with orbital observations. So clearly an environment that is conducive to cold trapping is necessary but not sufficient for the presence of volatiles such as water. We present the limits on the abundances that are indeed consistent with orbital data. At the LCROSS impact site itself, the data are consistent with very high ice abundances at 50-100 cm depth. However, radar results rule out these high abundances.

Elphic, Richard C.↗

NASA NTRS2008

Meridional Circulation in the Sun

Measuring the depth variation of the meridional flows is important for understanding the solar cycle, at least according to a number of dynamo models. While attempting to extend the early observations of Giles (1999; Ph. D. thesis, Stanford Univ.) of time-distance measurements of flow, we have stumbled upon some systematic errors that can affect these measurements: 1) the additional distance traveled by radiation coming from points away from disk center causes an apparent 'shrinking' Sun, that is an apparent flow towards the disk center, 2) in measurements away from the central longitude, the rotation signal can leak into meridional flow signals, and 3) in measurements of the north-south mean travel times along the equator, a spurious error of 6 sec travel time is seen. That the signal is spurious is confirmed by observing half the time with the image rotated 180 degrees. Although this is an effect with mean travel times and not differences, it still seems useful to understand it. Attempts to understand and overcome these systematic problems will be presented. Forward modeling has been done using ray theory to test the sensitivity of travel times to various models.

Duvall, T. L., Jr.↗

NASA NTRSDate not supplied

Recent Developments in the Assimilation of Microwave and Radar Observations Into NWP Models

Microwave observations play a very important role in improving the weather forecasts. Although these observations are routinely assimilated into NWP models in clear-sky conditions, assimilation of all-sky microwave observations is very limited. Two main factors contributing to this limitation are inaccuracy in the input cloud and hydrometeor profiles used as input to the radiative transfer model and also error in scattering calculations performed by the radiative transfer model itself. The Mie theory is used by many fast RT models to estimate the optical properties of single particles. The Mie theory assumes spherical shapes for ice or snow particles with mixture of air and ice. However, hydrometeors scattering radiation at microwave frequencies have different shapes, sizes, and orientations. Therefore, using Mie theory to determine their optical properties leads to large uncertainties in all-sky radiative transfer calculations. The discrete dipole approximation (DDA) which approximates the optical properties of large objects in terms of discrete dipoles has shown promise in calculating the scattering properties of particles with different shapes in the microwave frequencies. This presentation focuses on recent advancements in the CRTM scattering calculations for frozen hydrometeors in the microwave frequencies using the DDA technique. In addition to using stand-alone CRTM calculations using collocated ATMS and reanalysis profiles, the data assimilation experiments conducted using the NOAA FV3GFS forecast system are used to evaluate the scattering improvements. Additionally, the backscattering information from the DDA database was used to implement a radar simulator into CRTM. The radar operator takes advantage of CRTM different modules to calculate clouds absorption and scattering properties. In addition to the forward model both adjoint and tangent linear of the radar simulator are implemented and evaluated as well. The radar simulator is currently being tested within the JEDI/GEOS data assimilation framework to facilitate the assimilation of radar measurements such as CloudSat CPR and GPM DPR into the NASA GEOS model.

Isaac Moradi↗

NASA NTRSDate not supplied

Enhancing CRTM All-Sky Simulations and Implementation of A New Active Sensor Module

Microwave observations play a very important role in improving the weather forecasts. Although these observations are routinely assimilated into NWP models in clear-sky conditions, assimilation of all-sky microwave observations is very limited. Two main factors contributing to this limitation are inaccuracy in the input cloud and hydrometeor profiles used as input to the radiative transfer model and also error in scattering calculations performed by the radiative transfer model itself. The Mie theory is used by many fast RT models to estimate the optical properties of single particles. The Mie theory assumes spherical shapes for ice or snow particles with mixture of air and ice. However, hydrometeors scattering radiation at microwave frequencies have different shapes, sizes, and orientations. Therefore, using Mie theory to determine their optical properties leads to large uncertainties in all-sky radiative transfer calculations. The discrete dipole approximation (DDA) which approximates the optical properties of large objects in terms of discrete dipoles has shown promise in calculating the scattering properties of particles with different shapes in the microwave frequencies. This presentation focuses on recent advancements in the CRTM scattering calculations for frozen hydrometeors in the microwave frequencies using the DDA technique. In addition to using stand-alone CRTM calculations using collocated ATMS and reanalysis profiles, the data assimilation experiments conducted using the NOAA FV3GFS forecast system are used to evaluate the scattering improvements. Additionally, the backscattering information from the DDA database was used to implement a radar simulator into CRTM. The radar operator takes advantage of CRTM different modules to calculate clouds absorption and scattering properties. In addition to the forward model both adjoint and tangent linear of the radar simulator are implemented and evaluated as well. The radar simulator is currently being tested within the JEDI/GEOS data assimilation framework to facilitate the assimilation of radar measurements such as CloudSat CPR and GPM DPR into the NASA GEOS model.

Isaac Moradi↗

NASA NTRSDate not supplied

CRTM Microwave Cloud Scattering Lookup Tables and Radar Simulator

Microwave observations play a very important role in improving the weather forecasts. Although these observations are routinely assimilated into NWP models in clear-sky conditions, assimilation of all-sky microwave observations is very limited. Two main factors contributing to this limitation are inaccuracy in the input cloud and hydrometeor profiles used as input to the radiative transfer model and also error in scattering calculations performed by the radiative transfer model itself. The Mie theory is used by many fast RT models to estimate the optical properties of single particles. The Mie theory assumes spherical shapes for ice or snow particles with mixture of air and ice. However, hydrometeors scattering radiation at microwave frequencies have different shapes, sizes, and orientations. Therefore, using Mie theory to determine their optical properties leads to large uncertainties in all-sky radiative transfer calculations. The discrete dipole approximation (DDA) which approximates the optical properties of large objects in terms of discrete dipoles has shown promise in calculating the scattering properties of particles with different shapes in the microwave frequencies. This presentation focuses on recent advancements in the CRTM scattering calculations for frozen hydrometeors in the microwave frequencies using the DDA technique. In addition to using stand-alone CRTM calculations using collocated ATMS and reanalysis profiles, the data assimilation experiments conducted using the NOAA FV3GFS forecast system are used to evaluate the scattering improvements. Additionally, the backscattering information from the DDA database was used to implement a radar simulator into CRTM. The radar operator takes advantage of CRTM different modules to calculate clouds absorption and scattering properties. In addition to the forward model both adjoint and tangent linear of the radar simulator are implemented and evaluated as well. The radar simulator is currently being tested within the JEDI/GEOS data assimilation framework to facilitate the assimilation of radar measurements such as CloudSat CPR and GPM DPR into the NASA GEOS model.

Isaac Moradi↗

NASA NTRS2012

Rainbow Fourier Transform

We present a novel technique for remote sensing of cloud droplet size distributions. Polarized reflectances in the scattering angle range between 135deg and 165deg exhibit a sharply defined rainbow structure, the shape of which is determined mostly by single scattering properties of cloud particles, and therefore, can be modeled using the Mie theory. Fitting the observed rainbow with such a model (computed for a parameterized family of particle size distributions) has been used for cloud droplet size retrievals. We discovered that the relationship between the rainbow structures and the corresponding particle size distributions is deeper than it had been commonly understood. In fact, the Mie theory-derived polarized reflectance as a function of reduced scattering angle (in the rainbow angular range) and the (monodisperse) particle radius appears to be a proxy to a kernel of an integral transform (similar to the sine Fourier transform on the positive semi-axis). This approach, called the rainbow Fourier transform (RFT), allows us to accurately retrieve the shape of the droplet size distribution by the application of the corresponding inverse transform to the observed polarized rainbow. While the basis functions of the proxy-transform are not exactly orthogonal in the finite angular range, this procedure needs to be complemented by a simple regression technique, which removes the retrieval artifacts. This non-parametric approach does not require any a priori knowledge of the droplet size distribution functional shape and is computationally fast (no look-up tables, no fitting, computations are the same as for the forward modeling).

Alexandrov, Mikhail D.↗

NASA NTRS2024

Spectral Correlation in MODIS Water-Leaving Reflectance Retrieval Uncertainty

Spectral remote sensing reflectance, Rrs(λ) (sr−1), is the fundamental quantity used to derive a host of bio-optical and biogeochemical properties of the water column from satellite ocean color measurements. Estimation of uncertainty in those derived geophysical products is therefore dependent on knowledge of the uncertainty in satellite-retrieved R rs . Furthermore, since the associated algorithms require R rs at multiple spectral bands, the spectral (i.e., band-to-band)error covariance in R rs is needed to accurately estimate the uncertainty in those derived properties. This study establishes a derivative-based approach for propagating instrument random noise, instrument systematic uncertainty, and forward model uncertainty into R rs as retrieved using NASA’s multiple-scattering epsilon (MSEPS) atmospheric correction algorithm, to generate pixel-level error covariance in R rs . The approach is applied to measurements from Moderate Resolution Imaging Spectroradiometer (MODIS) on the Aqua satellite and verified using Monte Carlo (MC) analysis. We also make use of this full spectral error covariance in R rs to calculate uncertainty in phytoplankton pigment chlorophyll-a concentration (chl a , mg/m 3 ) and diffuse attenuation coefficient of downwelling irradiance at 490 nm (K d (490), m -1 ). Accounting for the error covariance in R rs generally reduces the estimated relative uncertainty in chl a by ∼1-2% (absolute value) in waters with chl a < 0.25 mg/m 3 where the color index (CI) algorithm is used. The reduction is ∼5-10% in waters with chl a > 0.35 mg/m 3 where the blue-green ratio (OCX) algorithm is used. Such reduction can be higher than 30% in some regions. For K d (490), the reduction by error covariance is generally ∼2%, but can be higher than 20% in some regions. The error covariance in R rs is further verified through forward-calculating chl a from MODIS-retrieved and in situ R rs and comparing estimated uncertainty with observed differences. An 8-day global composite of propagated uncertainty shows that the goal of 35% uncertainty in chl a can be achieved over deep ocean waters (chl a ≤ 0.1 mg/m3). While the derivative-based approach generates reasonable error covariance in R rs some assumptions should be updated as our knowledge improves. These include the inter-band error correlation in top-of-atmosphere reflectance, and uncertainties in the calibration of MODIS 869 nm band, in ancillary data, and in the in situ data used for system vicarious calibration.

Ocean color↗

NASA NTRS2014

The Architecture of the LkCa 15 Transitional Disk Revealed By High-Contrast Imaging

We present four new epochs of Ks-band images of the young pre-transitional disk around LkCa 15 and perform extensive forward modeling to derive the physical parameters of the disk. We find indications of strongly anisotropic scattering (Henyey-Greenstein g = 0.67 (+0.18/−0.11)) and a significantly tapered gap edge ("round wall") but see no evidence that the inner disk, whose existence is predicted by the spectral energy distribution, shadows the outer regions of the disk visible in our images.We marginally confirm the existence of an offset between the disk center and the star along the line of nodes; however, the magnitude of this offset (x = 27 (+19/−20) mas) is notably lower than that found in our earlier H-band images. Intriguingly, we also find an offset of y = 69 (+49/−25) mas perpendicular to the line of nodes at high significance. If confirmed by future observations, this would imply a highly elliptical - or otherwise asymmetric - disk gap with an effective eccentricity of e ≈ 0.3. Such asymmetry would most likely be the result of dynamical sculpting by one or more unseen planets in the system. Finally, we find that the bright arc of scattered light we see in direct imaging observations originates from the near side of the disk and appears brighter than the far side because of strong forward scattering.

dynamical sculpting↗

NASA NTRS2014

The Architecture of the LkCa 15 Transitional Disk Revealed by High-contrast Imaging

We present four new epochs of Ks-band images of the young pre-transitional disk around LkCa 15, and perform extensive forward modeling to derive the physical parameters of the disk. We find indications of strongly anisotropic scattering (Henyey-Greenstein g = 0.67+0.18 −0.11) and a significantly tapered gap edge ('round wall'), but see no evidence that the inner disk, whose existence is predicted by the spectral energy distribution, shadows the outer regions of the disk visible in our images. We marginally confirm the existence of an offset between the disk center and the star along the line of nodes; however, the magnitude of this offset (x = 27+19 −20 mas) is notably lower than that found in our earlier H-band images (Thalmann et al. 2010). Intriguingly, we also find, at high significance, an offset of y = 69+49 −25 mas perpendicular to the line of nodes. If confirmed by future observations, this would imply a highly elliptical- or otherwise asymmetric-disk gap with an effective eccentricity of e ≈ 0.3. Such asymmetry would most likely be the result of dynamical sculpting by one or more unseen planets in the system. Finally, we find that the bright arc of scattered light we see in direct imaging observations originates from the near side of the disk, and appears brighter than the far side because of strong forward scattering.

anisotropic scattering↗

NASA NTRS1982

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.↗

NASA NTRS1983

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.↗

NASA NTRS1990

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.↗

NASA NTRS1991

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.↗

NASA NTRS1992

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.↗

NASA NTRS1991

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.↗

NASA NTRS1993

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.↗

NASA NTRS1995

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.↗

NASA NTRS1994

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.↗