Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface data”

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 145 records · Page 8

Ocean color spectrum calculations

The development is considered of procedures for measuring a number of subsurface oceanographic parameters using remotely sensed ocean color data. It is proposed that the first step in this effort should be the development of adequate theoretical models relating the desired oceanographic parameters to the upwelling radiances to be observed. A portion of a contributory theoretical model is shown to be described by a modified single scattering approach based upon a simple treatment of multiple scattering. The resulting quasi-single scattering model can be used to predict the upwelling distribution of spectral radiance emerging from the sea. The shape of the radiance spectrum predicted by this model for clear ocean water shows encouraging agreement with measurments made at the edge of the Sargasso Sea off Cape Hatteras.

Mccluney, W. R.↗

Ocean color spectrum calculations

There is obvious value in developing the means for measuring a number of subsurface oceanographic parameters using remotely sensed ocean color data. The first step in this effort should be the development of adequate theoretical models relating the desired oceanographic parameters to the upwelling radiances to be observed. A portion of a contributory theoretical model can be described by a modified single scattering approach based on a simple treatment of multiple scattering. The resulting quasisingle scattering model can be used to predict the upwelling distribution of spectral radiance emerging from the sea. The shape of the radiance spectrum predicted by this model for clear ocean water shows encouraging agreement with measurements made at the edge of the Sargasso Sea off Cape Hatteras.

Mccluney, W. R.↗

The gravity method and interpretive techniques for detecting vertical crustal movements

Observations of temporal variations of gravity can be used as an inexpensive and rapid means of detecting, monitoring, and studying crustal deformation associated with many active geological processes. Such observations can yield qualitative or semiquantitative information on elevation changes and, when combined with independent elevation data, can yield information about changes of the subsurface density field arising from both subsurface displacements and temporal variations of the density of materials in the subsurface. The effectiveness of gravity techniques in specific cases of crustal deformation depends on the configuration of the local gravity field, the physical processes involved in the deformation, and the accuracy with which temporal gravity variations can be measured.

Jachens, R. C.↗

Fiber Optic Thermal Health Monitoring of Aerospace Structures and Materials

A new technique is presented for thermographic detection of flaws in materials and structures by performing temperature measurements with fiber Bragg gratings. Individual optical fibers with multiple Bragg gratings employed as surface temperature sensors were bonded to the surfaces of structures with subsurface defects or thickness variations. Both during and following the application of a thermal heat flux to the surface, the individual Bragg grating sensors measured the temporal and spatial temperature variations. The investigated structures included a 10-ply composite specimen with subsurface delaminations of various sizes and depths. The data obtained from grating sensors were further analyzed with thermal modeling to reveal particular characteristics of the interested areas. These results were found to be consistent with those from conventional thermography techniques. Limitations of the technique were investigated using both experimental and numerical simulation techniques. Methods for performing in-situ structural health monitoring are discussed.

Wu, Meng-Chou↗

Preliminary Regional Analysis of the Kaguya Lunar Radar Sounder (LRS) Data through Eastern Mare Imbrium

The Lunar Radar Sounder (LRS) experiment on board the Kaguya spacecraft is observing the subsurface structure of the Moon, using ground-penetrating radar operating in the frequency range of 5 MHz [1]. Because LRS data provides in-formation about lunar features below the surface, it allows us to improve our understanding of the processes that formed the Moon, and the post-formation changes that have occurred (such as basin formation and volcanism). We look at a swath of preliminary LRS data, that spans from 7 to 72 N, and from 2 to 10 W, passing through the eastern portion of Mare Imbrium (Figure 1). Using software, designed for the mineral exploration industry, we produce a preliminary, coarse 3D model, showing the regional structure beneath the study area. Future research will involve smaller subsets of the data in regions of interest, where finer structures, such as those identified in [2], can be studied.

Cooper, B.L.↗

Lessons Learned from Assimilating Altimeter Data into a Coupled General Circulation Model with the GMAO Augmented Ensemble Kalman Filter

Satellite altimetry measurements have provided global, evenly distributed observations of the ocean surface since 1993. However, the difficulties introduced by the presence of model biases and the requirement that data assimilation systems extrapolate the sea surface height (SSH) information to the subsurface in order to estimate the temperature, salinity and currents make it difficult to optimally exploit these measurements. This talk investigates the potential of the altimetry data assimilation once the biases are accounted for with an ad hoc bias estimation scheme. Either steady-state or state-dependent multivariate background-error covariances from an ensemble of model integrations are used to address the problem of extrapolating the information to the sub-surface. The GMAO ocean data assimilation system applied to an ensemble of coupled model instances using the GEOS-5 AGCM coupled to MOM4 is used in the investigation. To model the background error covariances, the system relies on a hybrid ensemble approach in which a small number of dynamically evolved model trajectories is augmented on the one hand with past instances of the state vector along each trajectory and, on the other, with a steady state ensemble of error estimates from a time series of short-term model forecasts. A state-dependent adaptive error-covariance localization and inflation algorithm controls how the SSH information is extrapolated to the sub-surface. A two-step predictor corrector approach is used to assimilate future information. Independent (not-assimilated) temperature and salinity observations from Argo floats are used to validate the assimilation. A two-step projection method in which the system first calculates a SSH increment and then projects this increment vertically onto the temperature, salt and current fields is found to be most effective in reconstructing the sub-surface information. The performance of the system in reconstructing the sub-surface fields is particularly impressive for temperature, but not as satisfactory for salt.

Keppenne, Christian↗

GRAIL Refinements to Lunar Seismic Structure

Joint interpretation of disparate geophysical datasets helps to reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co‐located at one corner. Due to this limited geographical extent, near‐surface ray coverage from moonquakes is low, but increases with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefore well suited to joint inversion, since the complementary information reduces inherent model ambiguity. Previous joint inversions of the Apollo seismic data (seismic phase arrival times) and Clementine‐ or Lunar Prospector‐derived gravity data (mass and moment of inertia) attempted to recover the subsurface structure of the Moon by focusing on hypothetical lunar compositions that explore the density/velocity relationship. These efforts typically search for the best fitting thermodynamically calculated velocity/density model, allowing variables like core size, velocity, and/or composition to vary freely. Seismic velocity profiles previously derived from the Apollo seismic data through inversion of travel times vary both in the depth of the crust and mantle layers, and the seismic velocities and densities assigned to those layers. The lunar mass and moment of inertia likewise only constrain gross variations in the density profile beyond that of a uniform density sphere. As a result, composition and structure models previously obtained by jointly inverting these data retain the original uncertainties inherent in the input data sets. We will perform a joint inversion of Apollo seismic delay times and gravity data collected by the GRAIL lunar gravity mission, in order to recover seismic velocities and density as a function of latitude, longitude, and depth within the Moon. We will relate density to seismic velocity using a linear relationship that is allowed to be depth‐dependent. The corresponding coefficient (B) can reflect a variety of material properties that vary with depth, including temperature and composition. The inversion seeks to recover the set of density, velocity, and B‐coefficient perturbations that minimize (in a least‐squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗

GRAIL Refinements to Lunar Seismic Structure

Joint interpretation of disparate geophysical datasets helps reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co-located at one corner. Due to this limited geographical extent, near-surface ray coverage from moonquakes is low, but increases with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefore well suited to joint inversion, since the complementary information reduces inherent model ambiguity. Previous joint inversions of the Apollo seismic data (seismic phase arrival times) and Clementine- or Lunar Prospector-derived gravity data (mass and moment of inertia) attempted to recover the subsurface structure of the Moon by focusing on hypothetical lunar compositions that explored the density/velocity relationship. These efforts typically searched for the best fitting thermodynamically calculated velocity/density model, and allowed variables like core size, velocity, and/or composition to vary freely. Seismic velocity profiles derived from the Apollo seismic data through travel time inversion vary both in the depth of the crust and mantle layers, and the seismic velocities and densities assigned to those layers. The lunar mass and moment of inertia likewise only constrain gross variations in the density profile beyond that of a uniform density sphere. As a result, composition and structure models previously obtained by jointly inverting these data retain the original uncertainties inherent in the input data sets. We perform a joint inversion of Apollo seismic delay times and gravity data collected by the GRAIL lunar gravity mission, in order to recover seismic velocity and density as a function of latitude, longitude, and depth within the Moon. We relate density (p) to seismic velocity (v) using a depth-dependent linear relationship. The corresponding coefficient (B) can reflect a variety of material properties, including temperature and composition. The inversion seeks to recover the set of p, v, and B perturbations that minimize (in a least-squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗

Multi-Wavelength Dielectrometer (MWD) Sensor For Planetary Subsurface Exploration

NASA's mission to Moon, Mars, and Beyond envisions landing of a light weight measurement platform on the planetary surface. The Multi-Wavelength Dielectrometer (MWD) on-board consists of essential electronics and metallic plates acting as electrodes attached to the body of such platform. An electric signal applied to one of the electrodes acting as a cathode sets up electric field pattern (in the soil medium) between the cathode and other electrodes acting as anodes. The electrodes are swept through multiple wavelengths (1Hz-1MHz) and the electric current drawn by the electrodes is measured at each frequency. The measured current whose amplitude and phase depend upon electrode spacing, dielectric constant of the subsurface soil, and the frequency is then used to estimate electrical properties of the soil. In this paper the MWD sensor that will measure the dielectric properties of Moon/Mars s soil is presented. A procedure to process the MWD measured data for extracting the soil properties is also described. Assuming the subsurface soil structure as multilayer strata having varying electric properties, an electric equivalent circuit of the multiple electrodes configuration placed on a multi-layer soil sample is obtained. The current drawn by the equivalent circuit from the low frequency signal generator is then calculated. By minimizing the difference between the model s estimated current and measured MWD data the electric properties of soil samples are extracted. Experimental and simulated results will be presented to validate the proposed procedure for extracting soil properties.

Deshpande, Manohar↗

Dislocation model for aseismic fault slip in the transverse ranges of Southern California

Geodetic data at a plate boundary can reveal the pattern of subsurface displacements that accompany plate motion. These displacements are modelled as the sum of rigid block motion and the elastic effects of frictional interaction between blocks. The frictional interactions are represented by uniform dislocation on each of several rectangular fault patches. The block velocities and fault parameters are then estimated from geodetic data. Bayesian inversion procedure employs prior estimates based on geological and seismological data. The method is applied to the Transverse Ranges, using prior geological and seismological data and geodetic data from the USGS trilateration networks. Geodetic data imply a displacement rate of about 20 mm/yr across the San Andreas Fault, while the geologic estimates exceed 30 mm/yr. The prior model and the final estimates both imply about 10 mm/yr crustal shortening normal to the trend of the San Andreas Fault. Aseismic fault motion is a major contributor to plate motion. The geodetic data can help to identify faults that are suffering rapid stress accumulation; in the Transverse Ranges those faults are the San Andreas and the Santa Susana.

Cheng, A.↗

GRAIL-Identified Gravity Anomalies in Oceanus Procellarum: Insight into Subsurface Impact and Magmatic Structures on the Moon

Four, quasi-circular, positive Bouguer gravity anomalies (PBGAs) that are similar in diameter (~90-190 km) and gravitational amplitude (> 140 mGal contrast) are identified within the central Oceanus Procellarum region of the Moon. These spatially associated PBGAs are located south of Aristarchus Plateau, north of Flamsteed crater, and two are within the Marius Hills volcanic complex (north and south). Each is characterized by distinct surface geologic features suggestive of ancient impact craters and/or volcanic/plutonic activity. Here, we combine geologic analyses with forward modeling of high-resolution gravity data from the Gravity Recovery and Interior Laboratory (GRAIL) mission in order to constrain the subsurface structures that contribute to these four PBGAs. The GRAIL data presented here, at spherical harmonic degrees 6–660, permit higher resolution analyses of these anomalies than previously reported, and reveal new information about subsurface structures. Specifically, we find that the amplitudes of the four PBGAs cannot be explained solely by mare-flooded craters, as suggested in previous work; an additional density contrast is required to explain the high-amplitude of the PBGAs. For Northern Flamsteed (190 km diameter), the additional density contrast may be provided by impact-related mantle uplift. If the local crust has a density ~2800 kg/cu.m, then ~7 km of uplift is required for this anomaly, although less uplift is required if the local crust has a lower mean density of ~2500 kg/cu.m. For the Northern and Southern Marius Hills anomalies, the additional density contrast is consistent with the presence of a crustal complex of vertical dikes that occupies up to ~50% of the regionally thin crust. The structure of Southern Aristarchus Plateau (90 km diameter), an anomaly with crater-related topographic structures, remains ambiguous. Based on the relatively small size of the anomaly, we do not favor mantle uplift; however, understanding mantle response in a region of especially thin crust needs to be better resolved. It is more likely that this anomaly is due to subsurface magmatic material given the abundance of volcanic material in the surrounding region. Overall, the four PBGAs analyzed here are important in understanding the impact and volcanic/plutonic history of the Moon, specifically in a region of thin crust and elevated temperatures characteristic of the Procellarum KREEP Terrane.

Deutsch, Ariel N.↗

Radar Sounding of Mars: A Focus on MARSIS

Radar has the unique capability of looking under the dry and cold surfaces of Mars. The depth of penetration of radio waves depends on a number of surface and subsurface parameters such as surface topography, subsurface geological structure and surface and subsurface electromagnetic properties. Among these parameters, the surface topography is known best largely due to valuable data provided by Mars Global Surveyor's MOLA instrument. However, little information is available on the electromagnetic properties and subsurface characteristics of Mars.

radar sounding↗

A Polarimetric Scattering Model for the 2-Layer Problem

In this paper, polarimetric signatures from two layers are examined. In the two layer problem a low-loss dielectric layer sits on top of a shallow subsurface layer of scatterers. Volume scattering within the top surface layer is not considered. We show that for polarimetric backscatter data acquired from just one incidence angle, no parameter is proportional to the depth of the subsurface layer. We also show that separation of the two components and inversion of the polarimetric backscatter measurements is not possible, in the latter case due to the mismatch between the number of variables that determine the backscatter behavior, versus the degrees of freedom in the measured values. The interesting results obtained at a sand dune site in SW France courtesy of the group led by Dr. Phillippe Paillou are examined through the prism of this model. It is shown that the model agrees well with results obtained for a subsurface layer at 3 m depth, where the subsurface returns are significant, and that surface scattering alone can account for observations made over a subsurface layer at a depth of 5.8 m, provided the surface slope is adjusted by approx. 10 degrees. The implications for exploration of the near-surface of Mars, a significant proportion of which is obscured by a dust mantle, are discussed.

Freeman, A.↗

Improving Ocean Reanalyses and ENSO Forecasts By Assimilation of Rain Corrected Satellite Sea Surface Salinity Using the GMAO S2S Forecast System

During the past years, we have seen that the La Nina to El Nino transition has had a significant global socio-economic impact and so has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of satellite altimetry and subsurface temperature and salinity from (mostly) Argo help to improve the initialization of the thermocline, while satellite Sea Surface Temperature (SST) aids in constraining surface heat-fluxes, leading to improved subseasonal to seasonal forecasts of the coupled system. However, few studies have focused on improving the fresh-water flux and near-surface density and mixing through assimilation of satellite sea surface salinity (SSS). For expediency, the few ocean models that do assimilate SSS bias-correct the satellite SSS data to normalize towards the near-surface Argo data. However, in rainy regions, where buoyant water sits as a fresh lens at the surface, this assumption is likely inadequate. In previous work, we have shown that adjusting SSS data to bulk salinity (Sb) using the Rain Impact Model (RIM) of Santos-Garcia et al., 2014 has improved the near-surface density and mixed layer depth, leading to deeper thermocline and improved the NINO3.4 SST forecasts. Now we utilize the Soil Moisture/Ocean Salinity, Rain Corrected (SMOS_RC) SSS product provided by the Centre Aval de Traitement des données SMOS (CATDS CPDC) to represent the Sb more accurately at first model layer (in our case 5 m). Rather than using a diffusivity model as with RIM, SMOS_RC relies on an observed relationship between the spatial heterogeneity of SMOS SSS and instantaneous rain rate (RR) (Supply et al., 2020). In order to test the impact of SMOS_RC versus SMOS, we compare two reanalyses over the period 2014 to 2021. For both reanalysis experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF scheme (Penny et al., 2013). One reanalysis additionally assimilates SMOS SSS data as is (i.e., with the fresh bias), and a separate reanalysis is performed assimilating the SMOS_RC data. We assess the impact for near-surface and subsurface dynamics within ocean reanalyses by validating against observations and explore how SSS assimilation (SMOS versus SMOS_RC) impacts dynamical ENSO forecasts using the NASA GMAO Sub-seasonal to Seasonal coupled forecast system (GEOS S2S-3, Molod et al., 2020, Hackert et al., 2023). We will show that improved SSS estimates and near-surface density and mixing led to more accurate coupled air/sea interaction and better ENSO forecasts.

Eric Hackert↗

Comparison of Retracking Algorithms Using Airborne Radar and Laser Altimeter Measurements of the Greenland Ice Sheet

This paper compares four continental ice sheet radar altimeter retracking algorithms using airborne radar and laser altimeter data taken over the Greenland ice sheet in 1991. The refurbished Advanced Application Flight Experiment (AAFE) airborne radar altimeter has a large range window and stores the entire return waveform during flight. Once the return waveforms are retracked, or post-processed to obtain the most accurate altitude measurement possible, they are compared with the high-precision Airborne Oceanographic Lidar (AOL) altimeter measurements. The AAFE waveforms show evidence of varying degrees of both surface and volume scattering from different regions of the Greenland ice sheet. The AOL laser altimeter, however, obtains a return only from the surface of the ice sheet. Retracking altimeter waveforms with a surface scattering model results in a good correlation with the laser measurements in the wet and dry-snow zones, but in the percolation region of the ice sheet, the deviation between the two data sets is large due to the effects of subsurface and volume scattering. The Martin et al model results in a lower bias than the surface scattering model, but still shows an increase in the noise level in the percolation zone. Using an Offset Center of Gravity algorithm to retrack altimeter waveforms results in measurements that are only slightly affected by subsurface and volume scattering and, despite a higher bias, this algorithm works well in all regions of the ice sheet. A cubic spline provides retracked altitudes that agree with AOL measurements over all regions of Greenland. This method is not sensitive to changes in the scattering mechanisms of the ice sheet and it has the lowest noise level and bias of all the retracking methods presented.

Ferraro, Ellen J.↗

Radar Reflectivity of the Martian Polar Regions

Radar experiments provide a unique method of probing the surfaces and subsurfaces of planetary bodies. Information on surface and subsurface structures and properties can be extracted from radar data. There is a well developed history of radar investigations of the planet Mars, beginning with the first reports of variations of scattering properties as a function of martian longitude in the mid-1960's. Because of the rapid rotation of Mars, the standard technique of delay-doppler mapping cannot be used as effectively as it is on Mercury and Venus, making it a more difficult radar target. Techniques have been developed to overcome this difficulty, including an inversion technique to combine many doppler-only (CW) experiments from different viewing geometries into a map of surface reflectivity, and new random long-code techniques. These techniques still suffer from ambiguities, however. A third technique which does not suffer from these ambiguities is the combination of the powerful Goldstone transmitter with the VLA as the receiving instrument to create a combined radar imaging instrument. We have used this combined radar instrument to image the surface of Mars in 3.5-cm radar reflectivity during the 1988, 1992/93, and 1999 oppositions. During the 1988 experiments, the residual south polar ice cap (RSPIC) was the brightest radar reflector on the planet - intrinsically brighter than even the Tharsis lava flows. This was quite an unexpected result. In contrast, during the 1992/93 experiments, the residual north polar ice cap (RNPIC) was not nearly so bright, and in fact showed no enhancement at all, This was puzzling, given the 1988 results for the RSPIC. We attributed the lack of a radar reflection enhancement to a combination of three effects: the geometry was different; the season was different; and the intrinsic scattering from the residual ice caps was different. The 1999 experiments provided a chance to test the relative importance of these three effects, since the RNPIC would be visible in a geometry very similar to the RSPIC in 1988, and the season would also be very similar. Preliminary reduction of the data taken in 1999 shows that the RNPIC did in fact show a radar reflectivity enhancement, but that it was still not as bright as the RSPIC in 1988. Details regarding the radar reflectivity of the residual ice caps, as well as the surrounding polar layered terrains (in light of our radar reflectivity data) will be discussed. Additional information is contained in the original extended abstract.

Butler, B. J.↗

Numerical Simulations of Silverpit Crater Collapse

The Silverpit crater is a recently discovered, 60-65 Myr old complex crater, which lies buried beneath the North Sea, about 150 km east of Britain. High-resolution images of Silverpit's subsurface structure, provided by three-dimensional seismic reflection data, reveal an inner-crater morphology similar to that expected for a 5-8 km diameter terrestrial crater. The crater walls show evidence of terrace-style slumping and there is a distinct central uplift, which may have produced a central peak in the pristine crater morphology. However, Silverpit is not a typical 5-km diameter terrestrial crater, because it exhibits multiple, concentric rings outside the main cavity. External concentric rings are normally associated with much larger impact structures, for example Chicxulub on Earth, or Orientale on the Moon. Furthermore, external rings associated with large impacts on the terrestrial planets and moons are widely-spaced, predominantly inwardly-facing, asymmetric scarps. However, the seismic data show that the external rings at Silverpit represent closely-spaced, concentric faultbound graben, with both inwardly and outwardly facing fault-scarps. This type of multi-ring structure directly analogous to the Valhalla-type multi-ring basins found on the icy satellites. Thus, the presence and style of the multiple rings at Silverpit is surprising given both the size of the crater and its planetary setting. A further curiosity of the Silverpit structure is that the external concentric rings appear to be extensional features on the West side of the crater and compressional features on the East side. The crater also lies in a local depression, thought to be created by postimpact movement of a salt layer buried beneath the crater.

Collins, G. S.↗

Numerical Simulations of Silverpit Crater Collapse

The Silverpit crater is a recently discovered, 60-65 Myr old complex crater, which lies buried beneath the North Sea, about 150 km east of Britain. High-resolution images of Silverpit's subsurface structure, provided by three-dimensional seismic reflection data, reveal an inner-crater morphology similar to that expected for a 5-8 km diameter terrestrial crater. The crater walls show evidence of terracestyle slumping and there is a distinct central uplift, which may have produced a central peak in the pristine crater morphology. However, Silverpit is not a typical 5-km diameter terrestrial crater, because it exhibits multiple, concentric rings outside the main cavity. External concentric rings are normally associated with much larger impact structures, for example Chicxulub on Earth, or Orientale on the Moon. Furthermore, external rings associated with large impacts on the terrestrial planets and moons are widely-spaced, predominantly inwardly-facing, asymmetric scarps. However, the seismic data show that the external rings at Silverpit represent closely-spaced, concentric fault-bound graben, with both inwardly and outwardly facing faults-carps. This type of multi-ring structure is directly analogous to the Valhalla-type multi-ring basins found on the icy satellites. Thus, the presence and style of the multiple rings at Silverpit is surprising given both the size of the crater and its planetary setting.

Collins, G. S.↗