Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface data models”

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 19 records

El Nino and the Global Ocean Observing System

Until a decade ago, an often-quoted expression in oceanography is that very few observations are recorded throughout the ocean. Now, the sentiment is no longer valid in the uppermost 10% of the tropical Pacific Ocean nor at the surface of the global ocean. One of the remarkable legacies of the 1985-1994 Tropical Oceans Global Atmosphere (TOGA) Program is an in situ marine meteorological and upper oceanographic measurement array throughout the equatorial Pacific to monitor the development and maintenance of El Nino episodes. The TOGA Observing System, which initially consisted of moored- and drifting-buoy arrays, a network of commercial ships, and coastal and island stations, now includes a constellation of satellites and data-assimilating models to simulate subsurface oceanographic conditions. The El Nino and La Nina tropical Pacific Ocean observing system represents the initial phase of an integrated global ocean observing system. Remarkable improvements have been made in ocean model simulation of subsurface currents, but some problems persist. For example, the simulation of the South Equatorial Current (SEC) remains an important challenge in the 2S-2N Pacific equatorial wave guide. During El Nino the SEC at the equator is reduced and sometimes the direction is reversed, becoming eastward. Both conditions allow warm water stored in the western Pacific to invade the eastern region, creating an El Nino episode. Assimilation of data is a tenet of faith to correct simulation errors caused by deficiencies in surface fluxes (especially wind stress) and parameterizations of subgrid-scale physical processes. In the first of two numerical experiments, the Pacific SEC was simulated with and without assimilation of subsurface temperature data. Along the equator, a very weak SEC occurred throughout the eastern Pacific, independent of assimilation of data. However, as displayed in the diagram, in the western Pacific there was no satisfactory agreement between the two simulations. To help determine reliability of the simulated SEC in the western Pacific, current measurements recorded during the 9-19 October 1994 voyage of the French research vessel L'Atalante are also shown in the diagram. With data assimilation, the simulated SEC was in much better agreement with L'Atalante observations. The simulated SEC with data assimilation was far from perfect, in part because of the sparsity of subsurface temperature observations. In the next experiment, TOPEX/POSEIDON sea surface height data in combination with subsurface temperatures will be assimilated to assess further improvement of the simulation of the SEC.

Halpern, David↗

Simulation gravity modeling to spacecraft-tracking data - Analysis and application

It is proposed that line-of-sight gravity measurements derived from spacecraft-tracking data can be used for quantitative subsurface density modeling by suitable orbit simulation procedures. Such an approach avoids complex dynamic reductions and is analogous to the modeling of conventional surface gravity data. This procedure utilizes the vector calculations of a given gravity model in a simplified trajectory integration program that simulates the line-of-sight gravity. Solutions from an orbit simulation inversion and a dynamic inversion on Doppler observables compare well (within 1% in mass and size), and the error sources in the simulation approximation are shown to be quite small. An application of this technique is made to lunar crater gravity anomalies by simulating the complete Bouguer correction to several large young lunar craters. It is shown that the craters all have negative Bouguer anomalies.

Phillips, R. J.↗

OGCM Simulations of Equatorial Pacific Current and Temperature to ERS-1, FSU and NMC Surface Winds and to Assimilation of Subsurface Temperature Data

The relative accuracies of three surface wind data products for the tropical Pacific Ocean during April 1992 to March 1994 were examined by analyzing temperature and current fields along the equator, which were simulated with an ocean general circulation model. Simulations were made with and without assimilation of surface and subsurface temperature data. Simulated currents were compared with observations at three sites (170oW, 140oW, 110oW) at the equator. Model-generated currents and temperatures indicated that the ERS-1 westward wind speeds were low compared to the FSU and NMC winds. With data assimilation, the agreement between simulated and observed currents was highest at 170oW and lowest at 110oW.

ERS-1 FSU↗

Critique of 'Elastic thickness of the Venus lithosphere estimated from topography and gravity' by A. Cazenave and K. Dominh

A critique of a model of the near-surface structure of Venus by Cazenave and Dominh (1981) based on Pioneer Venus data on gravity and topography is presented. Two objections are raised, the first an assumption of Airy compensation at a depth of 6 km, and the second that a gravity-topography relation is sufficient to study the subsurface characteristics. It is shown that determination of the gravity field along the line of sight, and not along the vertical, is necessary for a valid geophysical model with the available data. It is concluded that the spacecraft data are sufficient to define constraints on models of the subsurface characteristics, but are not sufficient for actually describing the substrate.

Reasenberg, R. D.↗

Geosat Data Assimilation with Application to the Eastern North Atlantic

An attempt is made to determine the three-dimensional ocean circulation from satellite altimeter measurements by assimilating Geosat sea surface height data into an eddy-resolving QuasiGeostrophic (QG) model of the eastern North Atlantic Ocean. Results are tested against independent information from hydrographic field observations and moored current meter data collected during the Geosat ERM. The comparison supports the concept of inferring aspects of the three-dimensional flow field from sea surface height observations by combining altimetric measurements with the dynamics of ocean circulation models. A Holland-type QG model with open boundaries is set up on a 2000 km X 2000 km domain of the eastern North Atlantic between 25 deg. and 45 deg. N, 32 deg. and 8 deg. W. By using a simple nudging technique, about two years of Geosat altimeter data are assimilated into the model every five days as space-time objective analyses on the model grid. The error information resulting from the analysis is used during the assimilation procedure to account for data uncertainties. Results show an intense eddy field, which in the surface layer interacts with a meandering Azores Front. Compared to Geosat, the model leads to smoothed fields that follow the observations. Model simulations are significantly correlated with hydrographic data from March 1988 and June 1989, both close to the surface and in the subsurface. Good agreement is also found between the model velocity fields and moored current meter data in the top two model layers. The agreement is visually weak in the bottom layer, although a coherence analysis reveals an agreement between the model simulation and current meter data over the full water column at periods exceeding 80 days.

Stammer, Detlef↗

Repurposing Drilling Control Diagnostics for Subsurface Edge Detection and Boundary Advisement During Planetary Drilling

Informed decision-making during lunar drilling and sampling missions will require data monitoring tools and specialized ground data systems. Accurate and updated situational awareness, with ongoing data monitoring, is critical for timely responses by to incoming science data. Traverse plans and scheduled activities may need to be flexibly changed in order to react to unexpected data or situations. Unlike (for example) Mars missions, the relative lightspeed closeness of the Moon allows for near-real-time ground processing of incoming mission and instrument data. An Apollo-class lunar regolith drill will in a sense “travel” a meter or two vertically at a given subsurface characterization site. As the drill penetrates into lunar regolith, it is likely to encounter a range of material densities, orientations, fracture toughness, and (perhaps) ice percentages. Lunar drill telemetry can provide science teams with a valuable first look into the subsurface structure, the regolith bulk properties, and constituents at each drilled site. Real-time AI-based recognition and reaction to downhole situations has been developed for automated deeper drilling on Mars and beyond. We can leverage the same knowledge bases and pattern-matching as areal-time interpreter of the subsurface, a situational awareness tool during drilling operations. We recently (Sept. 2019) demonstrated this AI drilling monitoring and analysis capability, in control of in-situ drilling and sampling operations, mounted on a KREX-2 rover in Chile’s Atacama Desert. Terrestrial automated drilling log analyses in oil exploration have used similar machine learning techniques in classifying and identifying features in drilling logs –but these typically are designed assuming a drilling fluid influencing downhole measurements and data (permeability, resistivity). Drilling models and existing AI software designed to detect and respond to drilling faults and hard materials can be repurposed, for near-real-time (ground-based) interpretation of drilling telemetry –a potentially valuable advisory tool for strata boundaries and changes in drilling parameters. On the Moon, this approach could be used to study the structure and to some extent the composition of lunar regolith vs. borehole depth, based on recognizable variations in fracture hardness, drilling energy and penetration rates while actively drilling. Since the early 2000s, a series of increasingly-capable real-time drilling telemetry interpretation and characterization software tools have been developed. These subsurface models and software tools have monitored the real-time drilling data received, and automatically identified changes in drill behavior (e.g., encountering a harder target layer, bit inclusions, drill choking due to infall downhole, and others) correlating these with subsurface structures and features. We discuss the mappings between drill borehole parameters, faults or events detected, and modeled changes in rock layer boundaries, in examples drawn from field testing at analog sites in an Arctic impact crater, Rio Tinto, and Chile’s Atacama Desert. These demonstrate how subsurface structural boundaries led to fault detections and responses by the software.

drilling advisor↗

Applications of Data Assimilation to Analysis of the Ocean on Large Scales

It is commonplace to begin talks on this topic by noting that oceanographic data are too scarce and sparse to provide complete initial and boundary conditions for large-scale ocean models. Even considering the availability of remotely-sensed data such as radar altimetry from the TOPEX and ERS-1 satellites, a glance at a map of available subsurface data should convince most observers that this is still the case. Data are still too sparse for comprehensive treatment of interannual to interdecadal climate change through the use of models, since the new data sets have not been around for very long. In view of the dearth of data, we must note that the overall picture is changing rapidly. Recently, there have been a number of large scale ocean analysis and prediction efforts, some of which now run on an operational or at least quasi-operational basis, most notably the model based analyses of the tropical oceans. These programs are modeled on numerical weather prediction. Aside from the success of the global tide models, assimilation of data in the tropics, in support of prediction and analysis of seasonal to interannual climate change, is probably the area of large scale ocean modeling and data assimilation in which the most progress has been made. Climate change is a problem which is particularly suited to advanced data assimilation methods. Linear models are useful, and the linear theory can be exploited. For the most part, the data are sufficiently sparse that implementation of advanced methods is worthwhile. As an example of a large scale data assimilation experiment with a recent extensive data set, we present results of a tropical ocean experiment in which the Kalman filter was used to assimilate three years of altimetric data from Geosat into a coarsely resolved linearized long wave shallow water model. Since nonlinear processes dominate the local dynamic signal outside the tropics, subsurface dynamical quantities cannot be reliably inferred from surface height anomalies. Because of its potential for large scale synoptic coverage of the deep ocean, acoustic travel time data should be a natural complement to satellite altimetry. Satellite data give us vertical integrals associated with thermodynamic and dynamic processes.

Miller, Robert N.↗

Detection of Subsurface Material Separation in Shuttle Orbiter Slip-Side Joggle Region of the Wing Leading Edge using Infrared Imaging Data from Arc Jet Tests

The objective of the present study was to determine whether infrared imaging (IR) surface temperature data obtained during arc-jet tests of Space Shuttle Orbiter s reinforced carbon-carbon (RCC) wing leading edge panel slip-side joggle region could be used to detect presence of subsurface material separation, and if so, to determine when separation occurs during the simulated entry profile. Recent thermostructural studies have indicated thermally induced interlaminar normal stress concentrations at the substrate/coating interface in the curved joggle region can result in local subsurface material separation, with the separation predicted to occur during approach to peak heating during reentry. The present study was an attempt to determine experimentally when subsurface material separations occur. A simplified thermal model of a flat RCC panel with subsurface material separation was developed and used to infer general surface temperature trends due to the presence of subsurface material separation. IR data from previously conducted arc-jet tests on three test specimens were analyzed: one without subsurface material separation either pre or post test, one with pre test separation, and one with separation developing during test. The simplified thermal model trend predictions along with comparison of experimental IR data of the three test specimens were used to successfully infer material separation from the arc-jet test data. Furthermore, for the test specimen that had developed subsurface material separation during the arc-jet tests, the initiation of separation appeared to occur during the ramp up to the peak heating condition, where test specimen temperature went from 2500 to 2800 F.

Daryabeigi, Kamran↗

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

Determining the 3D Subsurface Density Structure of Taurus Littrow Valley Using Apollo 17 Gravity Data

Surface gravity surveys can detect subsurface density variations that can reveal subsurface geologic features. In 1972, the Apollo 17 (A17) mission conducted the Traverse Gravimeter Experiment (TGE) using a gravimeter that measured the local gravity field near Taurus Littrow Valley (TLV), located on the south-eastern rim of the Serenitatis basin. TLV is hypothesized to be a basaltfilled radial graben resulting from the impact that formed Mare Serenitatis. It is bounded by both the North and South Massifs (NM and SM) as well as other smaller mountains to the East that are thought to be mainly composed of brecciated highland material. The TGE is the first and only successful gravity survey on the surface of the Moon. Other more recent satellite surveys, such as NASA's Gravity Recovery and Interior Laboratory (GRAIL) mission (2011- 2012), have produced the best global gravity field to date (approx. 13km resolution). However, these satellite surveys are not sensitive enough to detect fine-scale (<1km) lunar subsurface structures. This underscores the value of the data collected at the surface by A17. In the original analysis of the data a 2D forward-modelling approach was used to derive a thickness of the subsurface basalt layer of 1.0 km by assuming a simple flat-faced rectangular geometry and using densities derived from Apollo lunar samples. We are investigating whether modern 3D modelling techniques in combination with high-resolution topographical and image datasets can reveal additional fine-scale subsurface structure in TLV.

TGE↗

Development of Great Lakes algorithms for the Nimbus-G coastal zone color scanner

A series of experiments in the Great Lakes designed to evaluate the application of the Nimbus G satellite Coastal Zone Color Scanner (CZCS) were conducted. Absorption and scattering measurement data were reduced to obtain a preliminary optical model for the Great Lakes. Available optical models were used in turn to calculate subsurface reflectances for expected concentrations of chlorophyll-a pigment and suspended minerals. Multiple nonlinear regression techniques were used to derive CZCS water quality prediction equations from Great Lakes simulation data. An existing atmospheric model was combined with a water model to provide the necessary simulation data for evaluation of the preliminary CZCS algorithms. A CZCS scanner model was developed which accounts for image distorting scanner and satellite motions. This model was used in turn to generate mapping polynomials that define the transformation from the original image to one configured in a polyconic projection. Four computer programs (FORTRAN IV) for image transformation are presented.

Tanis, F. J.↗

Aseismic crustal deformation in the Transverse Ranges of southern California

Using geodetic data at a plate boundary, the subsurface displacements that accompany plate motion were modeled as the sum of rigid block motion and the elastic effects of frictional drag between blocks in a model consisting of 12 blocks and 27 rectangular fault patches. The block motion inferred from the geodetic data has the same order of magnitude as the geologic estimates, and for many faults the agreement was excellent. The model implies about 20 mm/y of crustal shortening normal to the trend of the San Andreas fault.

Cheng, Abe↗

Mars Sulfate Formation Sourced in Sulfide-Enriched Subsurface Fluids: The Rio Tinto Model

The extensive evidence for sulfate deposits on Mars provided by analyses of MER and Mars Express data shows that the sulfur played an essential role in the geochemical cycles of the planet, including reservoirs in the atmosphere, hydro-sphere and geosphere. Overall the data are consistent with a fluvial/lacustrine-evaporative origin of at least some of the sulfate deposits, with mineral precipitation through oversaturation of salty acidic fluids enriched in sulfates. This scenario requires reservoirs of sulfur and associated cations, as well as an acidic and oxidizing hydrochemistry which could be provided by surface and subsurface catching of meteoric waters resulting in the presence of sulfur-bearing gases and steam photochemistry. In this work we suggest a new scenario for the extensive generation of sulfates in Mars based on the observation of seasonal changes in the redox and pH of subsurface waters enriched in sulfur that supply the acidic Mars process analog of Rio Tinto. This model considers the long-term subsurface storage of sulfur during most of Noachian and its release from the late Noachian to Hesperian time through weathering by meteoric fluids that would acidify and oxidize the sulfur bearing compounds stored in the subsurface.

Fernandez-Remolar, D. C.↗

Foreland Basin Structures and Remote Sensing

Rocky Mountain foreland basins are somewhat unique in that the basins may exhibit a variety of structural styles. It is generally agreed that shortening has occurred in the foreland basement but the cause is controversial: vertical vs compressional horizontal tectonics. Even when shortening is attributed to compression, the attitude (dip) of the fault plane and whether the horizontal or vertical component of movement is dominant is unconstrained. The controversy is difficult to resolve from surface data alone due to the variety of possible interpretations. Detailed surface mapping and geologic modeling are needed to constrain subsurface interpretations. In many areas of the Wind River and Bighorn basins detailed geologic maps do not exist. State-of-the-art remote sensing data could potentially provide an efficient means of mapping surface geology. State-of-the-art remote sensing systems now provide geometrically correct data at 30 meter pixel size and increased spectral coverage, capable of more detailed geologic analyses. These data can be photographically enlarged to 1:24,000 scale and combined with 7 1/2' uses topographic quads to provide an excellent base map for geologic interpretations.

Paylor, E. D.↗

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

The Evolution of Oblique Impact Flow Fields Using Maxwell's Z Model

Oblique impacts are the norm rather than the exception for impact craters on planetary surfaces. This work focuses on the excavation of experimental oblique impact craters using the NASA Ames Vertical Gun Range (AVGR). Three-dimensional particle image velocimetry (3D PIV) is used to obtain quantitative data on ejection positions, three dimensional velocities and angles. These data are then used to constrain Maxwell's Z Model and follow the subsurface evolution of the excavation-stage flow-field center during oblique impacts.

Anderson, J. L. B.↗

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↗

Use of Data to Improve El Nino Southern Oscillation Simulations and Forecasts

Various data over the period 1980-1998 have been used to improve El Nino Southern Oscillation (ENSO) simulations and forecasts performed by coupled ocean/atmosphere models. For such a study, oceanic and atmospheric satellite data are necessary. Models assume a baroclinic ocean with a mixed layer coupled to a baroclinic atmosphere to simulate the anomalous fields. The various data have been used to revisit the parameterizations of the subsurface temperature in the mixed layer, atmospheric convection, and friction in the baroclinic ocean. This new model, named Tsub.Conv, simulates 4-year oscillations like the model with the original parameterization, but with anomalous thermocline and wind anomalies that agree a lot better with reality. The atmospheric model still presents some deficiencies in the eastern Pacific and when it is replaced by a statistical atmosphere, the coupled model (Tsub.Astat) performs even better. It is found that the off-equatorial ocean and wind curl anomalies play a crucial role for all these models and that Tsub.Astat simulates an ocean heat content recharge coming from the north prior to a warm event like the ones observed during the 1985-1992 period. The ENSO forecasts over 1980-1998 are considerably improved with Tsub.Astat. One needs to add the simulation of westerly wind bursts to be able to forecast the 1997 event one year in advance. Additional information is contained in the original.

Perigaud, Claire↗