SEARCH · Engineering Papers
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.
The International GPS Service for Geodynamics (IGS) is an organization which operates under the auspices of the International Association of Geodesy (IAG) and has been operational since January 1994. The primary objective of the IGS is to provide precise GPS data and data products to support geodetic and geophysical research activities.
The International GPS Service for Geodynamics (IGS) is an approved activity of the International.
Beginning June 1992 and continuing indefinitely as part of our contribution to FLINN (Fiducial Laboratories for an International Natural Science Network), DOSE (NASA's Dynamics of the Solid Earth Program), and the IGS (International GPS Geodynamics Service), analysts at the Jet Propulsion Laboratory (JPL) have routinely been reducing data from a globally-distributed network of Rogue Global Positioning System (GPS) receivers.
Global geodesy has improved dramatically over the past decade starting with the GPS for IERS and Geodynamics demonstration campaign in 1991 (GIG 91). At the time it took over a week of CPU time to process a network solution based on 21 global receivers and orbit overlaps were in the 40 cm range. Today it is possible to process a network solution based on 80 global receivers in less then one day of CPU time and orbit overlaps are in the 4 cm range. Special methods are under development for efficient processing of increasingly large regional networks which may contain hundreds or thousands of GPS receivers. Along the way there have been many lessons learned about GPS satellites, receivers, monuments, antennas, radomes, analysis, reference frames, error sources, and interpretation. A wide range of scientific disciplines have been impacted including studies of plate motion, post-glacial rebound, seasonal loading, deformation in plate boundary zones, coseismic displacements due to major earthquakes, postseiemic relaxation, and interseismic strain accumulation related to assessment of seismic hazards. Lessons learned will be presented in the context of new dense networks such as the Plate Boundary Observation (PBO).
We model the interior of Rhea on the basis of observational constraints and the results from geodynamical models available in the literature. Ten main types of models are defined, depending on the presence or absence of a high-pressure ice layer (ice II), and the extent of separation of the rock component from the volatiles. The degree-two gravity coefficients are computed for each of these models in order to assess which properties of the interior are likely to be inferred from Cassini radio science measurements scheduled on 26 November 2005. C22 greater than 2.5 x 10(exp -4) indicates that the satellite is undifferentiated, except for a slight increase in density with depth resulting from material self-compression. C22 between 1.67 x 10(exp -4)(lower bound) and 1.90 x 10(exp -4) indicates the presence of a rocky core, whose radius can be determined from the satellite's mass and ices densities, for a given temperature profile. For other values, most of the ten models cannot be distinguished from each other. However, assumptions on the density of the rock phase, presence or absence of ice II, and the degree of differentiation could allow a unique model to be determined in many cases. While the calculation presented in this work assumes that Rhea is in hydrostatic equilibrium, it is likely that Rhea' gravity field is partly affected by nonhydrostatic anomalies.
Precision mean element (PME) satellite theories play a key role in orbit dynamics analyses. These theories employ: nonsingular orbital elements comprehensive force models Generalized Method of Averaging Numerical interpolation concepts The Draper Semianalytical Satellite Theory (DSST) (Refs. 1 - 6), whose development was led by the author, and the independently-developed Universal Semianalytical Method (USM) (Ref. 7) are examples of such theories. These theories provide the capability to tailor the force modeling to meet the desired computational speed vs. accuracy trade-off. The flexibility of such theories is demonstrated by their ability to include complicated atmosphere density models and spacecraft models in the perturbation theory context. The value of high speed satellite theories, in this era of computational plenty, is that they allow new ways of looking at astrodynamical problems such as orbit design (Refs. 8, 9) and atmosphere density updating (Refs. 10, 11). In the mid to late-1980 s, the geodynamics community led the development of very precise geopotential models such as GEM T2 and GEM T3 (Ref. 12), and with the subsequent analysis of the TOPEX flight data, JGM-2 and JGM-3 (Ref. 13). These were high degree and order geopotentials, at least 50 x 50. In 1993, the DSST implementation in the GTDS program was extended to include the 50 x 50 geopotential models (Ref. 14). The 50 x 50 geopotential, J2000 integration coordinate system, and solid Earth tide capabilities were integrated in GTDS by Scott Carter (Ref. 15). This capability demonstrated 1 m accuracy versus the TOPEX Precise Orbit Ephemerides. Subsequently the DSST Standalone program was also extended to include high degree and order geopotential models (Ref. 5). More recently GTDS has been hosted in the Linux PC environment. However, all of these efforts have been limited to modeling the motion of an artificial Earth satellite. They did not consider the additional complexities associated with lunar, planetary, or other natural satellite orbiters. Such complexities include: additional coordinate systems (associated with the direction of the north pole of rotation and the prime meridian of the new central bodies) (Ref. 16) normalized gravity model coefficients (desirable for high degree and order fields) (Ref. 17) indirect oblateness
Lakshmi Planum is a high-standing plateau (3.5-4.5 km above MPR) surrounded by the highest mountain ranges on Venus. Lakshmi represents a unique type of elevated region different from dome-shaped and rifted rises and tessera-bearing crustal plateaus. The unique characteristics of Lakshmi suggest that it formed by an unusual combination of processes and played an important role in Venus geologic history. Lakshmi was studied with Venera-15/16 and Magellan data, resulting in two classes of models, divergent and convergent, to explain its unusual topographic and morphologic characteristics. Divergent models explain Lakshmi as a site of mantle upwelling due to rising and subsequent collapse of a mantle diapir; such models explain emplacement of a lava plateau inside Lakshmi and, in some circumstances, formation of the mountain ranges. The convergent models consider Lakshmi as a locus of mantle downwelling, convergence, underthrusting, and possible subduction. Key features in these models are the mountain ranges, high topography of Lakshmi interior, and the large volcanic centers in the plateau center. These divergent and convergent models entail principally different mechanisms of formation and suggest different geodynamic regimes on Venus. Almost all models make either explicit or implicit predictions about the type and sequence of major events during formation and evolution of Lakshmi and thus detailed geological mapping can be used to test them. Here we present the results of such geological mapping (the V-7 quadrangle, 50-75degN, 300-360degE; scale 1:5M) that allows testing the proposed models for Lakshmi.
Carbon is present in various forms in the Earth s upper mantle (carbonate- or diamond-bearing mantle xenoliths, carbonatite magmas, CO2 emissions from volcanoes...). Moreover, there is enough carbon in chondritic material to stabilize carbonates into the mantles of Mars or Venus as well as in the Earth. However, the interactions with iron have to be constrained, because Fe is commonly thought to buffer oxygen fugacity into planetary mantles. [1] and [2] show evidences of the stability of clinopyroxene Ca(Mg,Fe)Si2O6 + magnesite (Mg,Fe)CO3 in the Earth s mantle around 6GPa (about 180km). The stability of oxidized forms of carbon (like magnesite) depends on the oxygen fugacity of the system. In the Earth s mantle, the maximum carbon content is 10000 ppm [3]. The fO2 parameter varies vertically as a function of pressure, but also laterally because of geodynamic processes like subduction. Thus, carbonates, graphite, diamond, C-rich gases and melts are all stable forms of carbon in the Earth s mantle. [4] show that the fO2 variations observed in SNC meteorites can be explained by polybaric graphite-CO-CO2 equilibria in the Martian mantle. [5] inferred from thermodynamic calculations that the stable form of carbon in the source regions of the Martian basalts should be graphite (and/or diamond). After [6], a metasomatizing agent like a CO2-rich melt may infiltrate the mantle source of nakhlites. However, according to [7] and [8], the FeO wt% value in the Martian bulk mantle is more than twice that of the Earth s mantle (KLB-1 composition by [9]). As iron and carbon are two elements with various oxidation states, Fe/C interaction mechanisms must be considered.
Surface geomagnetic observation can determine up to degree L = 14 time-varying spherical harmonic coefficients of the poloidal magnetic field. Assimilation of these coefficients to numerical dynamo simulation could help us understand better the dynamical processes in the Earth's outer core, and to provide more accurate forecast of geomagnetic secular variations (SV). In our previous assimilation studies, only the poloidal magnetic field in the core is corrected by the observations in the analysis. Unobservable core state variables (the toroidal magnetic field and the core velocity field) are corrected via the dynamical equations of the geodynamo. Our assimilation experiments show that the assimilated core state converges near the CMB, implying that the dynamo state is strongly constrained by surface geomagnetic observations, and is pulled closer to the truth by the data. We are now carrying out an ensemble of assimilation runs with 1000 years of geomagnetic and archeo/paleo magnetic record. In these runs the cross correlation between the toroidal and the poloidal magnetic fields is incorporated into the analysis. This correlation is derived from the physical boundary conditions of the toroidal field at the core-mantle boundary (CMB). The assimilation results are then compared with those of the ensemble runs without the cross-correlation, aiming at understanding two fundamental issues: the effect of the crosscorrelation on (1) the convergence of the core state, and (2) the SV prediction accuracies. The constrained dynamo solutions will provide valuable insights on interpreting the observed SV, e.g. the near-equator magnetic flux patches, the core-mantle interactions, and possibly other geodynamic observables.
This report updates the description of the space geodesy facilities of the Spanish National Geographic Institute (IGN). The current 40-meter radio telescope at Yebes, a network station for IVS, has performed geodetic VLBI observations regularly since September 2008. In addition to this, the project to establish an Atlantic Network of Geodynamical and Space Stations (RAEGE) is progressing with the construction of the first antenna, which is being erected at Yebes.
The main activities of the VLBI group at the Department of Geodesy and Geoinformation of the Vienna University of Technology were related to the development of the Vienna VLBI Software VieVS (http://vievs.hg.tuwien.ac.at/) and its application for various studies. For example, we dealt with scheduling, satellite tracking, and the estimation of geodynamical and astronomical parameters from VLBI observations. One highlight was the release of VieVS 2.0 just before the third VieVS User Workshop in September 2012.
Determining the volatile budget of the interior of Mars is crucial for our understanding of that planet's formation, geodynamics, cooling history and the origin of its volcanism and atmosphere as well as its potential for life. Surficial water is evident from spacecraft and rover data in polar caps and the atmosphere, in the presence of river channels, and in the detection of water-bearing minerals. Meteorites, however, are our best candidates for estimating the amount of water present at depth, even if all are crustal samples. The last 10 years have seen a blooming of studies measuring water and halogens in Martian meteorites. The bulk of these studies target phosphate, a typically late-stage phase in the igneous Martian meteorites that potentially would concentrate incompatible element hydrogen (H quantified traditionally as "water", i.e., H2O concentrations in weight) near the end of the crystallization sequence. However, determining the amount of water, F, and Cl in the magma from which a phosphate crystallized from is not straightforward and in most instances not possible. On the other hand, phosphates have turned out to be very useful in identifying hydrothermal processes that could have added water while or after the magma flowed and crystallized. Another caveat of analyzing Martian meteorite phases for water is that shocked phases such as maskelynite and impact melts appear to have incorporated water from the Martian atmosphere, as evidenced by high H isotope ((delta)D) signatures, and therefore their water concentrations cannot be interpreted in terms of deep planetary processes. The best candidates for estimating the water content of the Martian interior have been melt inclusions (glass or amphibole-bearing) which the enclosing mineral (usually olivine) would have prevented from exchanging volatiles with the surroundings after crystallization. Even some of these, however, have high (delta)D, meaning they were affected by H exchange via impact events or with crustal reservoirs or hydrothermal fluids. Here, nominally anhydrous minerals (pyroxene, olivine, plagioclase, or maskelynite) in orthopyroxenite ALH 84001 and selected nakhlites are analyzed for water and major elements, in order to determine 1) whether they contain any water; 2) if they do, what controls its distribution (crystallization, degassing, hydrothermal or impact processes); and 3) if any of these measurements can be used to infer the water contents of the parent magma and their mantle sources. A shock-reverberation experiment was also performed on terrestrial orthopyroxenes (opx) to simulate the heavily shocked conditions of ALH 84001 (> 31 GPa [17]).
Most volatile studies of Mars have targeted its surface via spacecraft and rover data, and have evidenced surficial water in polar caps and the atmosphere, in the presence of river channels, and in the detection of water bearing minerals. The other focus of Martian volatile studies has been on Martian meteorites which are all from its crust. Most of these studies are on hydrous phases like apatite, a late-stage phase, i.e. crystallizing near the end of the differentiation sequence of Martian basalts and cumulates. Moreover, calculating the water content of the magma a phosphate crystallized from is not always possible, and yet is an essential step to estimate how much water was present in a parent magma and its source. Water, however, is primarily dissolved in the interiors of differentiated planets as hydrogen in lattice defects of nominally anhydrous minerals (olivine, pyroxene, feldspar) of the crust and mantle. This hydrogen has tremendous influence, even in trace quantities, on a planet's formation, geodynamics, cooling history and the origin of its volcanism and atmosphere as well as its potential for life. Studies of hydrogen in nominally anhydrous phases of Martian meteorites are rare. Measuring water contents and hydrogen isotopes in well-characterized nominally anhydrous minerals of Martian meteorites is the goal of our study. Our work aims at deciphering what influences the distribution and origin of hydrogen in Martian minerals, such as source, differentiation, degassing and shock.
The history of the Greenland Ice Sheet has been influenced by the geodynamic response to ice sheet fluctuations, and this interaction may help explain past deglaciations under modest climate forcing. We hypothesize that when the Iceland hot spot passed beneath north‐central Greenland, it thinned the lithosphere and left anomalous heat likely with partially melted rock; however, it did not break through the crust to supply voluminous flood basalts. Subsequent Plio‐Pleistocene glacial‐interglacial cycles caused large and rapidly migrating stresses, driving dike formation and other processes that shifted melted rock toward the surface. The resulting increase in surface geothermal flux favored a thinner, faster‐responding ice sheet that was more prone to deglaciation. If this hypothesis of control through changes in geothermal flux is correct, then the long‐term (10 (sup 5) to 10 (sup 6) years) trend now is toward lower geothermal flux, but with higher‐frequency (less than or equal to 10( sup 4) to 10 (sup 5) years) oscillations linked to glacial‐interglacial cycles. Whether the geothermal flux is increasing or decreasing now is not known but is of societal relevance due to its possible impact on ice flow. We infer that projections of the future of the ice sheet and its effect on sea level must integrate geologic and geophysical data as well as glaciological, atmospheric, oceanic, and paleoclimatic information.
A fully automated near real-time GPS tracking system has been developed around JPL's GIPSY/OASIS II software. The system produces <25 cm (3D rms) GPS orbits and one-half nanosecond (15 cm) clock estimates. The process starts automatically when a favorable global distribution of ground data from the IGS network (International GPS Service for Geodynamics) becomes available. Ionospherically corrected phase and pseudorange data are optimally combined to remove satellite and ground receiver clock errors, including selective availability. After the GPS orbits are determined within the data arc, they are then propagated with empirically determined dynamic force models. Real-time <2 meter (3D rms) GPS orbits are always available. As a by-product of this process, other calibration estimates such as station clocks, troposphere estimates, and Earth orientation parameters are also produced. For daily arc fits, the process requires 7-8 hours of CPU time on an HP9000/735 workstation.
VATMOS-SR (Venus ATMOSpheric - Sample Return) is a small spacecraft mission concept which would return a gas sample from the upper atmosphere of Venus to Earth for scientific analysis. This could be the first sample return mission for an extra-terrestrial atmosphere, and potentially the first sample return from an Earth-sized planet. The VATMOS-SR mission concept consists of a SmallSat atmospheric sampling probe (45 deg. sphere/cone geometry, <1 m diameter) that is designed to skim through the Venus upper atmosphere and acquire gas samples below the homopause altitude (around ~110 km altitude), where the different atmospheric gases are mixed. The velocity of the spacecraft where sampling would occur is expected to be between ~10.5 km/s and ~13.1 km/s, depending on the trajectory chosen. This presentation will discuss hypervelocity sampling in the upper atmosphere of Venus, with respect to the VATMOS-SR mission concept. VATMOS-SR would enable critical atmospheric measurements to form a full picture of how, why, and when Venus evolved to be so different from Earth and Mars. The abundances and isotopic compositions of volatile elements (such as N, C, S, O and the noble gases) in planetary atmospheres record volatile delivery during accretion, outgassing from planetary interiors, and atmospheric loss to space. Precise and accurate determinations of atmospheric volatile signatures are the key to understanding the origins and geodynamical evolution of Venus compared to the other terrestrial planets. Hypersonic sampling poses unique technical and scientific challenges. To ensure it is possible to relate the composition of the sampled gases to the free stream atmospheric composition, large-scale numerical simulations are employed to model the flow through the VATMOS-SR sampling system. In particular, an emphasis is placed on quantifying noble gas isotopic fractionation that occurs during the sample acquisition and transfer process in order to determine how measured isotopic ratios of noble gases in the sample compare to actual isotopic ratios in the Venusian atmosphere. The Direct Simulation Monte Carlo (DSMC) code SPARTA, an open source software package developed by Sandia National Laboratories, is used in this work. SPARTA, based on Bird’s DSMC method, is a molecular-level gas-kinetic technique. As SPARTA is able to model hypervelocity reacting flows in strong chemical and thermal non-equilibrium, this software package is well suited to determine relevant flow properties for the VATMOS-SR mission concept, and to numerically quantify the expected level of elemental and/or isotopic fractionation in the sample acquired by VATMOS-SR. This presentation will show results from 3D simulations correlating the noble gas isotopic fractionation in the gas acquired at hypervelocity speeds to its ambient atmosphere value. In particular, emphasis will be placed at Xenon isotopes of masses 128 and 130, as precise measurements of that ratio would yield a comparison to Earth’s atmosphere. Additionally, sensitivity studies that quantify the uncertainties due to the freestream parameters as well as the modeling parameters will be performed.
Probes to penetrate the thick ice shells of our solar system’s Ocean Worlds have been studied for nearly 20 years, since scientific evidence strongly suggested a subsurface ocean on the Jupiter moon called Europa. There is keen scientific interest in exploring four significant themes on such proposed missions: 1) Geodynamics, 2) Geochemistry, 3) Habitability, and 4) Life Detection. The ice shells of Ocean Worlds are predicted to be up to 40 km thick; they exhibit extreme thermal environments, with ice temperatures from 100 K to 270 K, and extreme pressure environments from vacuum to 53 MPa. Jet Propulsion Laboratory has conducted a broad-look investigation of proposed mission concepts to Europa to identify the significant technology and operational challenges of Europa ice-penetration. The thermal-mechanical system (TMS) of an ice penetration probe (IPP) designed to access the ocean of an icy moon using radioisotope thermoelectric generators for heat and power faces technological hurdles exacerbated by severe thermal and volume constraints. This study identified thermal management and control (TMC) challenges that are strongly linked to: ice penetration start-up, mobility and navigation in the ice, communications while in the ice sheet, and detecting and avoiding in-ice hazards. The major objectives of the TMC system are: 1) Absorb internal thermal energy from the IPP radioisotope power source, 2) Maintain liquid water conditions around the IPP at all times, 3) Manage and control thermal flows from probe nose to tail, and 4) Provide pressure containment for all internal probe components. This work discusses the baseline TMC system architecture and design developed to accomplish these objectives, and survive and transit the extreme ice thicknesses in pursuit of Icy/Ocean Worlds science goals. The proposed TMC system consisting of an internal pumped two-phase fluid loop “thermal bus” for thermal energy capture, variable conductance heat pipe system for passively adaptive thermal energy transport around the probe, and water jetting system for ice cutting is described and discussed. Critical testing performed to date is described.
The White Mountain Magma Series (WMMS) is the largest Mesozoic felsic igneous province on the Eastern North American Margin (ENAM). Existing geochronology has suggested that magmatism occurred over 50 Myr with published ages for the oldest units apparently coeval with the ca. 201 Ma Central Atlantic Magmatic Province (CAMP), the flood basalt province associated with the end-Triassic mass extinction and the opening of the Atlantic Ocean. We use zircon U-Pb geochronology to show that emplacement of WMMS plutons was already underway at 207.5 Ma. The largest volcanic-plutonic complex, the White Mountain Batholith, was emplaced episodically from ca. 198.5 Ma to ca. 180 Ma, ca. 25 Myr older than published ages suggest, and all samples we dated from the Moat Volcanics are between ca. 185 to 180 Ma. This result shows that the Moat Volcanics and the White Mountain Batholith are broadly comagmatic, and it also constrains the age of a key Jurassic paleomagnetic pole. Our data suggest that a regional mantle thermal anomaly in eastern North America developed at least ca. 5 Myr prior to the main stage of CAMP flood basalt volcanism and suggests a geodynamic link between the WMMS and the CAMP.