Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geophysical constraints”

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 181 records · Page 10

Serendipitous Geodesy from Bennu's Short-Lived Moonlets

The Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer (OSIRIS-REx; or OREx) spacecraft arrived at its target, near-Earth asteroid (101955) Bennu, on December 3, 2018. The OSIRIS-REx spacecraft has since collected a wealth of scientific information in order to select a suitable site for sampling. Shortly after insertion into orbit on December 31, 2018, particles were identified in starfield images taken by the navigation camera (NavCam 1). Several groups within the OSlRlS-REx team analyzed the particle data in an effort to better understand this newfound activity of Bennu and to investigate the potential sensitivity of the particles to Bennu's geophysical parameters. A number of particles were identified through automatic and manual methods in multiple images, which could be turned into short sequences of optical tracking observations. Here, we discuss the precision orbit determination (OD) effort focused on these particles at NASA GSFC, which involved members of the Independent Navigation Team (INT) in particular. The particle data are combined with other OSIRIS-REx tracking data (radiometric from OSN and optical landmark data) using the NASA GSFC GEODYN orbit determination and geodetic parameter estimation software. We present the results of our study, particularly those pertaining to the gravity field of Bennu. We describe the force modeling improvements made to GEODYN specifically for this work, e.g., with a raytracing-based modeling of solar radiation pressure. The short-lived, low-flying moonlets enable us to determine a gravity field model up to a relatively high degree and order: at least degree 6 without constraints, and up to degree 10 when applying Kaula-like regularization. We can backward- and forward-integrate the trajectory of these particles to the ejection and landing sites on Bennu. We assess the recovered field by its impact on the OSIRIS-REx trajectory reconstruction and prediction quality in the various mission phases (e.g., Orbital A, Detailed Survey, and Orbital B).

Mazarico, Erwan↗

Inversion of Time-Lapse Seismic Reservoir Monitoring Data Using CycleGAN: A Deep Learning-Based Approach for Estimating Dynamic Reservoir Property Changes

Carbon capture and storage is being pursued globally as a geoengineering measure for reducing the emission of anthropogenic CO 2 the atmosphere. Comprehensive monitoring, verification, and accounting programs must be established for demonstrating the safe storage of injected CO 2 . One of the most commonly deployed monitoring techniques is time-lapse seismic reservoir monitoring (also known as 4-D seismic), which involves comparing 3-D seismic survey data taken at the same study site but over different times. Analyses of 4-D seismic data volumes can help improve the quality of storage reservoir characterization, track the movement of injected CO 2 plume, and identify potential CO 2 spillover/leakage from the storage reservoirblue. However, the derivation of high-resolution CO 2 saturation maps from 4-D seismic data is a highly nonlinear and ill-posed inverse problem, often requiring significant computational effort. In this research, we apply a physics-based deep learning method to facilitate the solution of both the forward and inverse problems in seismic inversion while honoring physical constraints. A cycle generative adversarial neural network (CycleGAN) model is trained to learn the bidirectional functional mappings between the reservoir dynamic property changes and seismic attribute changes, such that both forward and inverse solutions can be obtained efficiently from the trained model. We show that our CycleGAN-based approach not only improves the reliability of 4-D seismic inversion but also expedites the quantitative interpretation. Our deep learning-based workflow is generic and can be readily used for reservoir characterization and reservoir model updates involving the use of 4-D seismic data.

58 GEOSCIENCES↗

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

An End-To-End Earthquake Detection Method for Joint Phase Picking and Association Using Deep Learning

Earthquake monitoring by seismic networks typically involves a workflow consisting of phase detection/picking, association, and location tasks. In recent years, the accuracy of these individual stages has been improved through the use of machine learning techniques. Here, in this study, we introduce a new end-to-end approach that improves overall earthquake detection accuracy by jointly optimizing each stage of the detection pipeline. We propose a neural network architecture for the task of multi-station processing of seismic waveforms recorded over a seismic network. This end-to-end architecture consists of three sub-networks: a backbone network that extracts features from raw waveforms, a phase picking sub-network that picks P- and S-wave arrivals based on these features, and an event detection sub-network that aggregates the features from multiple stations to associate and detect earthquakes across a seismic network. We use these sub-networks together with a shift-and-stack module based on back-projection that introduces kinematic constraints on arrival times, allowing the neural network model to generalize to different velocity models and to variable station geometry in seismic networks. We evaluate our proposed method on the STanford EArthquake Dataset (STEAD) and on the 2019 Ridgecrest, CA earthquake sequence. The results demonstrate that our end-to-end approach can effectively pick P- and S-wave arrivals and achieve earthquake detection accuracy rivaling that of other state-of-the-art approaches. Because our approach preserves information across tasks in the detection pipeline, it has the potential to outperform approaches that do not.

58 GEOSCIENCES↗

Flashline Mars Arctic Research Station (FMARS) 2009 Expedition Crew Perspectives

The Flashline Mars Arctic Research Station (FMARS), located on the rim of the Haughton Crater on Devon Island in the Canadian Arctic, is a simulated Mars habitat that provides operational constraints similar to those which will be faced by future human explorers on Mars. In July 2009, a six-member crew inhabited the isolated habitation module and conducted the twelfth FMARS mission. The crew members conducted frequent EVA operations wearing mock space suits to conduct field experiments under realistic Mars-like conditions. Their scientific campaign spanned a wide range of disciplines and included many firsts for Mars analog research. Among these are the first use of a Class IV medical laser during a Mars simulation, helping to relieve crew stress injuries during the mission. Also employed for the first time in a Mars simulation at FMARS, a UAV (Unmanned Aerial Vehicle) was used by the space-suited explorers, aiding them in their search for mineral resources. Sites identified by the UAV were then visited by geologists who conducted physical geologic sampling. For the first time, explorers in spacesuits deployed passive seismic equipment to monitor earthquake activity and characterize the planet's interior. They also conducted the first geophysical electromagnetic survey as analog Mars pioneers to search for water and characterize geological features under the surface. The crew collected hydrated minerals and attempted to produce drinkable water from the rocks. A variety of equipment was field tested as well, including new cameras that automatically geotag photos, data-recording GPS units, a tele-presence rover (operated from Florida), as well as MIT-developed mission planning software. As plans develop to return to the Moon and go on to Mars, analog facilities like FMARS can provide significant benefit to NASA and other organizations as they prepare for robust human space exploration. The authors will present preliminary results from these studies as well as their perspectives on topics including human factors, logistics, EVA operations, and the use of social media throughout the mission.

Cusack, Stacy↗

The Cretaceous-Tertiary Impact Crater and the Cosmic Projectile that Produced it

Evidence gathered to date from topographic data, geophysical data, well logs, and drill-core samples indicates that the buried Chicxulub basin, the source crater for the approximately 65 Ma Cretaceous-Tertiary (K/T) boundary deposits, is approximately 300 km in diameter. A prominent topographic ridge and a ring of gravity anomalies mark the position of the basin rim at approximately 150 km from the center. Wells in this region recovered thick sequences of impact-generated breccias at 200-300 m below present sea level. Inside the rim, which has been severely modified by erosion following impact, the subsurface basin continues to deepen until near the center it is approximately 1 km deep. The best planetary analog for this crater appears to be the 270 km-diameter Mead basin on Venus. Seismic reflection data indicate that the central zone of downward displacement and excavation (the transient crater is approximately 130 km in diameter, consistent with previous studies of gravity anomaly data). Our analysis of projectile characteristics utilizes this information, coupled with conventional scaling relationships, and geochemical constraints on the mass of extraterrestrial material deposited within the K/T boundary layer. Results indicate that the Chicxulub crater would most likely be formed by a long-period comet composed primarily of nonsilicate materials (ice, hydrocarbons, etc.) and subordinate amounts (less than or equal to 50 percent) primitive chondritic material. This collision would have released the energy equivalent to between 4 x 10(exp 8) and 4 x 10(exp 9) megatons of TNT. Studies of terrestrial impact rates suggest that such an event would have a mean production rate of approximately 1.25 x 10(exp -9) y(exp -1). This rate is considerably lower than that of the major mass extinctions over the last 250 million years (approximately 5 x 10(exp -7) y(exp -1). Consequently, while there is substantial circumstantial evidence establishing the cause-effect link between the Chicxulub basin forming event and the K/T biological extinctions, the results of our analysis do not support models of impact as a common or singular causative agent of mass extinctions on Earth.

Sharpton, Virgil L.↗

A Common Calibration Source Framework for Fully-Polarimetric and Interferometric Radiometers

Two types of microwave radiometry--synthetic thinned array radiometry (STAR) and fully-polarimetric (FP) radiometry--have received increasing attention during the last several years. STAR radiometers offer a technological solution to achieving high spatial resolution imaging from orbit without requiring a filled aperture or a moving antenna, and FP radiometers measure extra polarization state information upon which entirely new or more robust geophysical retrieval algorithms can be based. Radiometer configurations used for both STAR and FP instruments share one fundamental feature that distinguishes them from more 'standard' radiometers, namely, they measure correlations between pairs of microwave signals. The calibration requirements for correlation radiometers are broader than those for standard radiometers. Quantities of interest include total powers, complex correlation coefficients, various offsets, and possible nonlinearities. A candidate for an ideal calibration source would be one that injects test signals with precisely controllable correlation coefficients and absolute powers simultaneously into a pair of receivers, permitting all of these calibration quantities to be measured. The complex nature of correlation radiometer calibration, coupled with certain inherent similarities between STAR and FP instruments, suggests significant leverage in addressing both problems together. Recognizing this, a project was recently begun at NASA Goddard Space Flight Center to develop a compact low-power subsystem for spaceflight STAR or FP receiver calibration. We present a common theoretical framework for the design of signals for a controlled correlation calibration source. A statistical model is described, along with temporal and spectral constraints on such signals. Finally, a method for realizing these signals is demonstrated using a Matlab-based implementation.

Kim, Edward J.↗

Source Characterization of the Declared North Korean Nuclear Tests From Regional Distance Coda Wave Spectral Ratios

Abstract Seismic observations of underground nuclear explosions provide crucial data on source yield and depth that cannot easily be estimated from other geophysical methods. However, it is difficult to obtain reliable yield estimates for test sites for which we do not have direct seismic calibration experiments. To obtain source information from uncalibrated sites and paths, local and regional seismic records of six, proximal, declared underground nuclear explosions in North Korea are used to compute spectral ratios of narrow‐band waveform envelopes of body‐wave coda that remove path and site effects to reveal precise, relative source moment. The yields of these explosions are obtained from the observed source ratios by simultaneously fitting the classical source model of Mueller and Murphy (1971), https://doi.org/10.1785/bssa0610061675 to all event pairs. The source model provides an impressive fit to the observations considering that the P phase coda derived source spectral ratios did not require a prior knowledge of the source or regionally calibrated corrections to be applied to the data. However, the observed corner frequencies from S wave coda spectral ratios are lower than the source model predictions, but are well fit by models calculated using the corner frequency consistent with the Fisk conjecture. The results presented here provide novel constraints on the spectral distributions of the energy radiated by the sources of the DPRK test series, and allows for an independent evaluation of existing estimated source model parameters.

58 GEOSCIENCES↗

Geodetic Investigations of the Mission Concept MAGIC to Reveal Callisto’s Internal Structure

Geodetic and geophysical investigations of the Galilean moon Callisto address fundamental questions regarding the formation and evolution of the Jovian system. Callisto's evolution and internal structure appear to significantly differ from the other Jovian satellites. Similarly-sized Ganymede is a highly evolved ice-rock moon with a differentiated interior, intrinsic magnetic field, and abundant surface evidence of internal activity. In contrast, Callisto's surface is ancient, and Galileo spacecraft data suggest its interior is only incompletely differentiated, despite the presumed presence of a sub-surface ocean. These properties make Callisto uniquely able to constrain the timing and nature of the Jovian system formation. The Magnetics, Altimetry, Gravity, and Imaging of Callisto (MAGIC) mission concept is conceived to fully characterize the properties of this enigmatic moon from its deep interior to the icy shell. Three main instruments are included as a scientific payload. Highly accurate measurements of Callisto's topography, magnetic field, and morphology are obtained by the onboard laser altimeter, magnetometer, and camera, respectively. The telecommunication system supports an additional gravity and radio science investigation. Long- and short-wavelength gravity anomalies afford powerful constraints on internal differentiation and the properties of the hydrosphere (water and ice). Comprehensive numerical simulations and covariance analyses of MAGIC mission scenarios presented in this paper show that the gravitational degree-2 normalized coefficients and the pole obliquity enable the determination of the moment of inertia with an accuracy better than 0.015%. The combination of gravity and altimetry measurements acquired by MAGIC are essential to the characterization of Callisto's interior if – as is likely – the degree-2 gravity includes non-hydrostatic terms. MAGIC's radio science data yield the estimation of Callisto's gravity field with spatial resolutions of <100 km. The combination of gravitational and deformation tides that are retrieved by the radio science and altimetry investigations, respectively, leads to the recovery of the rigid ice shell thickness to within ∼3 km. Together these datasets would resolve ambiguities inherent in Galileo flyby data, revealing Callisto's interior structure as well as the existence and properties of its postulated internal ocean.

Antonio Genova↗

Planetary Geophysics and Tectonics

Research supported by grant NAGW-1928 has addressed a variety of problems related to planetary evolution. One important focus has been on questions related to the role of chemical buoyancy in planetary evolution with application to both Venus and the Moon. We have developed a model for the evolution of the Moon (Hess and Parmentier, 1995) in which dense, highly radioactive, late stage magma ocean cumulates sink forming a core. This core heats the overlying, chemically layered mantle giving rise to a heated, chemically well-mixed layer that thickens with time. This Mixed layer eventually becomes hot enough and thick enough that its top begins to melt at a pressure low enough that melt is buoyant, thus creating mare basalts from a high pressure source of the correct composition and at an appropriate time in lunar evolution. In work completed during the last year, numerical experiments on convection in a chemically stably stratified fluid layer heated from below have been completed. These results show us how to calculate the evolution of a mixed layer in the Moon, depending on the heat production in the ilmenite- cumulate core and the chemical stratification of the overlying mantle. Chemical stratification of the mantle after its initial differentiation is would trap heat in the deep interior and prevent the rapid rise of plumes with accompanying volcanism. This trapping of heat in the interior can explain the thickness of the lunar lithosphere as a function of time as well as the magmatic evolution. We show that heat transported to the base of the lithosphere at a rate determined by current estimates of radioactivity in the Moon would not satisfy constraints on elastic lithosphere thickness from tectonic feature associated with basin loading. Trapping heat at depth by a chemically stratified mantle may also explain the absence of global compressional features on the surface that previous models predict for an initially hot lunar interior. For Venus, we developed a model in which the chemical buoyancy of crust and a depleted mantle layer stabilizes the lithosphere for long periods of time and provides a mechanism of episodic planetary evolution (Parmentier and Hess, 1992). Continued thickening of a residual depleted mantle layer eventually suppresses pressure release melting and the creation of depleted mantle. Continued cooling then allows the lithosphere to become heavier than the underlying hotter, undepleted mantle. This repeated instability can occur on time scales appropriate for episodic global resurfacing on Venus. We have also examined the role of the gabbro-eclogite phase transformation on crust and lithosphere stability and as a mechanism of crustal recycling in the absence of plate tectonics. Our work thus far concentrates on the scale of instability that would occur due to cooling or crustal thickening associated with horizontal shortening. Whether repeated overturn can explain the evolution of Venus depends in part on whether sufficient heat transfer can occur between overturns and on constraints provided by understanding observed surface features and evolution.

Parmentier, E. M.↗

Connectivity-informed drainage network generation using deep convolution generative adversarial networks

Abstract Stochastic network modeling is often limited by high computational costs to generate a large number of networks enough for meaningful statistical evaluation. In this study, Deep Convolutional Generative Adversarial Networks (DCGANs) were applied to quickly reproduce drainage networks from the already generated network samples without repetitive long modeling of the stochastic network model, Gibb’s model. In particular, we developed a novel connectivity-informed method that converts the drainage network images to the directional information of flow on each node of the drainage network, and then transforms it into multiple binary layers where the connectivity constraints between nodes in the drainage network are stored. DCGANs trained with three different types of training samples were compared; (1) original drainage network images, (2) their corresponding directional information only, and (3) the connectivity-informed directional information. A comparison of generated images demonstrated that the novel connectivity-informed method outperformed the other two methods by training DCGANs more effectively and better reproducing accurate drainage networks due to its compact representation of the network complexity and connectivity. This work highlights that DCGANs can be applicable for high contrast images common in earth and material sciences where the network, fractures, and other high contrast features are important.

42 ENGINEERING↗

Measurements of Acidic Gases and Aerosol Species Aboard the NASA DC-8 Aircraft During the Pacific Exploratory Mission in the Tropics (PEM-Tropics A)

We received funding to provide measurements of nitric acid (HNO3), formic acid (HCOOH), acetic acid (CH3COOH), and the chemical composition of aerosols aboard the NASA Ames DC-8 research aircraft during the PEM-Tropics A mission. These measurements were successfully completed and the final data resides in the electronic archive (ftp-gte.larc.nasa.gov) at NASA Langley Research Center. For the PEM-Tropics A mission the University of New Hampshire group was first author of four different manuscripts. Three of these have now appeared in the Journal of Geophysical Research-Atmospheres, included in the two section sections on PEM-Tropics A. The fourth manuscript has just recently been submitted to this same journal as a stand alone paper. All four of these papers are included in this report. The first paper (Influence of biomass combustion emissions on the distribution of acidic trace gases over the Southern Pacific basin during austral springtime) describes the large-scale distributions of HNO3, HCOOH, and CH3COOH. Arguments were presented to show, particularly in the middle tropospheric region, that biomass burning emissions from South America and Africa were a major source of acidic gases over the South Pacific basin. The second paper (Aerosol chemical composition and distribution during the Pacific Exploratory Mission (PEM) Tropics) covers the aerosol aspects of our measurement package. Compared to acidic gases, O3, and selected hydrocarbons, the aerosol chemistry showed little influence from biomass burning emissions. The data collected in the marine boundary layer showed a possible marine source of NH3 to the troposphere in equatorial areas. This source had been speculated on previously, but our data was the first collected from an airborne platform to show its large-scale features. The third paper (Constraints on the age and dilution of Pacific Exploratory Mission-Tropics biomass burning plumes from the natural radionuclide tracer Pb-210) utilized the unexpectedly high concentrations of Pb-210 in the combustion plumes to estimate their ages and mixing along the transport route to the South Pacific basin. The final paper in the group (Tropospheric reactive-odd nitrogen over the South Pacific in austral springtime). This paper provides a summary of reactive nitrogen during PEM-Tropics A, with HNO3 and PAM showing the most impact from combustion emission.

Talbot, Robert W.↗

Improving multiwell petrophysical interpretation from well logs via machine learning and statistical models

Well-log interpretation estimates in situ rock properties along well trajectory, such as porosity, water saturation, and permeability, to support reserve-volume estimation, production forecasts, and decision making in reservoir development. However, due to measurement errors, variability of well logs caused by multiple measurement vendors, different borehole tools, and nonuniform drilling/borehole conditions, estimations of rock properties with original well logs without proper preprocessing may not be accurate, especially in the context of multiwell estimation. Well-log normalization techniques such as two-point scaling and mean-variance normalization are commonly used to improve the robustness of multiwell rock-property estimation. However, these techniques do not consider the correlation between well logs and require subjective knowledge for their effective implementation. To reduce uncertainties and processing time associated with multiwell rock-property estimation from well logs, we develop discriminative adversarial (DA) and linear constraint models for well-log normalization and rock-property estimation. The DA neural network model developed for well-log normalization and interpretation can perform linear and nonlinear well-log normalization while considering the joint distribution of each well log and rock properties. However, the linear constraint model uses an ensemble of predictions from linear models to constrain well-log normalization and rock-property estimation. We also develop a divergence-based type well identification method to select type (training) wells for a test well based on the statistical similarity of associated well-log distributions instead of the interwell distance. We apply the DA model to perform well-log normalization and prediction of permeability for the Seminole San Andres Unit carbonate reservoir. Compared with the permeability predicted with the classical machine learning model without well-log normalization and models with two-point scaling normalization, the DA model yields the most accurate permeability prediction by decreasing the mean-squared error of permeability prediction by 20%–50%.

Geochemistry & Geophysics↗

Hayabusa2#’s Exploration to Asteroids 2001 CC21 and 1998 KY26 Provides Key Insights Into Planetary Defense

The Hayabusa2 extended mission, nicknamed Hayabusa2# (SHARP: Small Hazardous Asteroid Reconnaissance Probe), started its mission after the Hayabusa2 spacecraft successfully returned in December 2020 to the Earth with Ryugu's samples and released the Sample Return Capsule containing it. Hayabusa2#’s key operations include flying by Near Earth Asteroid (98943) 2001 CC21 and rendezvousing with Near-Earth Asteroid 1998 KY26. The flyby is planned to be in 2026, and the rendezvous will be in 2031. Detailed investigations of these targets enable maturing key knowledge and technologies necessary for planetary defense. This paper reviews the characteristics of these targets and strategies to maximize outcomes from the mission. The spacecraft will fly by 2001 CC21 at a speed of ~5 km/s in 2026. Given constraints on the spacecraft condition, not designed to perform a flyby operation, operational plans require careful assessments to maximize the proximity observations of this asteroid. The ~700 m diameter asteroid's shape is currently unknown, though lightcurve observations suggest the asteroid is elongated. While the spin pole is not constrained well, the spin period is about 5.02 h with uncertainties of 0.01 h [1]. The taxonomic class may be L-type, which implies the presence of Calcium-aluminum-rich inclusions (CAI), one of the primitive materials in the solar system [2], although there are studies suggesting an S-type. Therefore, this must be carefully examined. Photometric and spectroscopic observations will strongly constrain this asteroid’s properties. Our rendezvous target, 1998 KY26, is a ~30 m diameter roundish object spinning at a spin period of 10.7 min, which may represent one of the common groups of Near Earth Asteroids in size and composition that give a higher likelihood of threatening the Earth. Radar and optical observations reconstructed this asteroid's shape with high uncertainties and inferred that its surface composition might be similar to that of carbonaceous asteroids [3]. No spacecraft has ever visited such small bodies. Hayabusa2# visiting 1998 KY26 will be the first mission to document its geophysical properties, offering strong insights into what most frequent invaders look like [4].

Hayabusa2#↗

A simulation to study the feasibility of improving the temporal resolution of LAGEOS geodynamic solutions by using a sequential process noise filter

A key drawback to estimating geodetic and geodynamic parameters over time based on satellite laser ranging (SLR) observations is the inability to accurately model all the forces acting on the satellite. Errors associated with the observations and the measurement model can detract from the estimates as well. These 'model errors' corrupt the solutions obtained from the satellite orbit determination process. Dynamical models for satellite motion utilize known geophysical parameters to mathematically detail the forces acting on the satellite. However, these parameters, while estimated as constants, vary over time. These temporal variations must be accounted for in some fashion to maintain meaningful solutions. The primary goal of this study is to analyze the feasibility of using a sequential process noise filter for estimating geodynamic parameters over time from the Laser Geodynamics Satellite (LAGEOS) SLR data. This evaluation is achieved by first simulating a sequence of realistic LAGEOS laser ranging observations. These observations are generated using models with known temporal variations in several geodynamic parameters (along track drag and the J(sub 2), J(sub 3), J(sub 4), and J(sub 5) geopotential coefficients). A standard (non-stochastic) filter and a stochastic process noise filter are then utilized to estimate the model parameters from the simulated observations. The standard non-stochastic filter estimates these parameters as constants over consecutive fixed time intervals. Thus, the resulting solutions contain constant estimates of parameters that vary in time which limits the temporal resolution and accuracy of the solution. The stochastic process noise filter estimates these parameters as correlated process noise variables. As a result, the stochastic process noise filter has the potential to estimate the temporal variations more accurately since the constraint of estimating the parameters as constants is eliminated. A comparison of the temporal resolution of solutions obtained from standard sequential filtering methods and process noise sequential filtering methods shows that the accuracy is significantly improved using process noise. The results show that the positional accuracy of the orbit is improved as well. The temporal resolution of the resulting solutions are detailed, and conclusions drawn about the results. Benefits and drawbacks of using process noise filtering in this type of scenario are also identified.

Hartman, Brian Davis↗

Soil Moisture Active Passive Mission L4_C Data Product Assessment (Version 2 Validated Release)

The SMAP satellite was successfully launched January 31st 2015, and began acquiring Earth observation data following in-orbit sensor calibration. Global data products derived from the SMAP L-band microwave measurements include Level 1 calibrated and geolocated radiometric brightness temperatures, Level 23 surface soil moisture and freezethaw geophysical retrievals mapped to a fixed Earth grid, and model enhanced Level 4 data products for surface to root zone soil moisture and terrestrial carbon (CO2) fluxes. The post-launch SMAP mission CalVal Phase had two primary objectives for each science product team: 1) calibrate, verify, and improve the performance of the science algorithms, and 2) validate accuracies of the science data products as specified in the L1 science requirements. This report provides analysis and assessment of the SMAP Level 4 Carbon (L4_C) product pertaining to the validated release. The L4_C validated product release effectively replaces an earlier L4_C beta-product release (Kimball et al. 2015). The validated release described in this report incorporates a longer data record and benefits from algorithm and CalVal refinements acquired during the SMAP post-launch CalVal intensive period. The SMAP L4_C algorithms utilize a terrestrial carbon flux model informed by SMAP soil moisture inputs along with optical remote sensing (e.g. MODIS) vegetation indices and other ancillary biophysical data to estimate global daily net ecosystem CO2 exchange (NEE) and component carbon fluxes for vegetation gross primary production (GPP) and ecosystem respiration (Reco). Other L4_C product elements include surface (10 cm depth) soil organic carbon (SOC) stocks and associated environmental constraints to these processes, including soil moisture and landscape freeze/thaw (FT) controls on GPP and respiration (Kimball et al. 2012). The L4_C product encapsulates SMAP carbon cycle science objectives by: 1) providing a direct link between terrestrial carbon fluxes and underlying FT and soil moisture constraints to these processes, 2) documenting primary connections between terrestrial water, energy and carbon cycles, and 3) improving understanding of terrestrial carbon sink activity in northern ecosystems. There are no L1 science requirements for the L4_C product; however self-imposed requirements have been established focusing on NEE as the primary product field for validation, and on demonstrating L4_C accuracy and success in meeting product science requirements (Jackson et al. 2012). The other L4_C product fields also have strong utility for carbon science applications; however, analysis of these other fields is considered secondary relative to primary validation activities focusing on NEE. The L4_C targeted accuracy requirements are to meet or exceed a mean unbiased accuracy (ubRMSE) for NEE of 1.6 g C/sq m/d or 30 g C/sq m/yr, emphasizing northern (45N) boreal and arctic ecosystems; this is similar to the estimated accuracy level of in situ tower eddy covariance measurement-based observations (Baldocchi 2008).

Cal/Val↗

Wasatch Fault Structure from Machine Learning Arrival Times and High-Precision Earthquake Locations

Abstract On 18 March 2020, a magnitude 5.7 earthquake hit the Salt Lake valley in the state of Utah, United States. Using a dense geophone deployment and machine learning (ML), an additional several thousand events were detected and located. Currently, both the mainshock and the majority of the aftershocks are suspected to have occurred on or near a deeper portion of the Salt Lake segment of the Wasatch fault—part of a large range-bounding fault system thought to be capable of generating an Mw 7.2 earthquake. However, a small subset of aftershocks may have occurred on a portion of the more steeply, eastward dipping, and poorly understood West Valley fault. Unfortunately, the catalog locations and lack of focal mechanisms for this subset of aftershocks provide only a crude constraint on the true fault structure. To better illuminate fault structure, we relocate the ML-generated catalog with a range of magnitudes from −2 to 4.6, using: (1) NonLinLoc, a nonlinear location algorithm, (2) source-specific station terms, and (3) waveform coherence. We further compute first-motion focal mechanisms for 68 events. Results of the relocation suggest a simpler, minimally listric Wasatch fault geometry, contrary to what has been previously proposed. We also find that analysis of the focal mechanisms and waveform similarity indicates minimal event similarity throughout the Magna sequence, suggesting a highly complex and heterogeneous rupture zone, as opposed to rupture on a single plane. These findings suggest an increased seismic hazard due to the overall shallowness of the earthquake sequence and highly varied rupture mechanisms.

Geochemistry & Geophysics↗

GPS Monitoring of Surface Change During and Following the Fortuitous Occurrence of the M(sub w) = 7.3 Landers Earthquake in our Network

Accomplishments: (1) Continues GPS monitoring of surface change during and following the fortuitous occurrence of the M(sub w) = 7.3 Landers earthquake in our network, in order to characterize earthquake dynamics and accelerated activity of related faults as far as 100's of kilometers along strike. (2) Integrates the geodetic constraints into consistent kinematic descriptions of the deformation field that can in turn be used to characterize the processes that drive geodynamics, including seismic cycle dynamics. In 1991, we installed and occupied a high precision GPS geodetic network to measure transform-related deformation that is partitioned from the Pacific - North America plate boundary northeastward through the Mojave Desert, via the Eastern California shear zone to the Walker Lane. The onset of the M(sub w) = 7.3 June 28, 1992, Landers, California, earthquake sequence within this network poses unique opportunities for continued monitoring of regional surface deformation related to the culmination of a major seismic cycle, characterization of the dynamic behavior of continental lithosphere during the seismic sequence, and post-seismic transient deformation. During the last year, we have reprocessed all three previous epochs for which JPL fiducial free point positioning products available and are queued for the remaining needed products, completed two field campaigns monitoring approx. 20 sites (October 1995 and September 1996), begun modeling by development of a finite element mesh based on network station locations, and developed manuscripts dealing with both the Landers-related transient deformation at the latitude of Lone Pine and the velocity field of the whole experiment. We are currently deploying a 1997 observation campaign (June 1997). We use GPS geodetic studies to characterize deformation in the Mojave Desert region and related structural domains to the north, and geophysical modeling of lithospheric behavior. The modeling is constrained by our existing and continued GPS measurements, which will provide much needed data on far-field strain accumulation across the region and on the deformational response of continental lithosphere during and following a large earthquake, forming the basis for kinematic and dynamic modeling of secular and seismic-cycle deformation. GPS geodesy affords both regional coverage and high precision that uniquely bear on these problems.

Miller, M. Meghan↗