Engineering Papers⌕ Search

Engineering topics

Schubert, Gerald

Publications and source records attributed to Schubert, Gerald.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

The Atmospheres of the Terrestrial Planets:Clues to the Origins and Early Evolution of Venus, Earth, and Mars

We review the current state of knowledge of the origin and early evolution of the three largest terrestrial planets - Venus, Earth, and Mars - setting the stage for the chapters on comparative climatological processes to follow. We summarize current models of planetary formation, as revealed by studies of solid materials from Earth and meteorites from Mars. For Venus, we emphasize the known differences and similarities in planetary bulk properties and composition with Earth and Mars, focusing on key properties indicative of planetary formation and early evolution, particularly of the atmospheres of all three planets. We review the need for future in situ measurements for improving our understanding of the origin and evolution of the atmospheres of our planetary neighbors and Earth, and suggest the accuracies required of such new in situ data. Finally, we discuss the role new measurements of Mars and Venus have in understanding the state and evolution of planets found in the habitable zones of other stars.

origin of mars↗
NASA NTRS2005

Amalthea's density is less than that of water

Radio Doppler data from the Galileo spacecraft's encounter with Amalthea, one of Jupiter's small inner moons, on 5 November 2002 yield a mass of (2.08 +/- 0.15) x 10(18) kilograms. Images of Amalthea from two Voyager spacecraft in 1979 and Galileo imaging between November 1996 and June 1997 yield a volume of (2.43 +/- 0.22) x 10(6) cubic kilometers. The satellite thus has a density of 857 +/- 99 kilograms per cubic meter. We suggest that Amalthea is porous and composed of water ice, as well as rocky material, and thus formed in a cold region of the solar system, possibly not at its present location near Jupiter.

Jupiter↗
NASA NTRS2001

Three-Dimensional Simulations of Mantle Convection in Io

Io has very high surface heat flow and an abundance of volcanic activity, which are thought to be driven by nonuniform tidal heating in its interior. This nonuniform heat is transported to the base of the lithosphere by very vigorous convection in Io's silicate mantle, the form of which is presumably responsible for the distribution of surface features such as volcanoes and mountains. We here present three-dimensional spherical calculations of mantle convection in Io, in order to ascertain the likely form of this convection and the resulting distribution of heat flow at the surface and core-mantle boundary. Different models of tidal dissipation are considered: the endmember scenarios (identified by M. N. Ross and G. Schubert) of dissipation in the entire mantle, or dissipation in a thin (approximately 100-km-thick) asthenosphere, as well as the 'preferred' distribution of M. N. Ross et al. comprising 1/3 mantle and 2/3 asthenosphere heating. The thermal structure of Io's mantle and asthenosphere is found to be strongly dependent on tidal heating mode, as well as whether the mantle-asthenosphere boundary is permeable or impermeable. Results indicate a large-scale flow pattern dominated by the distribution of tidal heating, with superimposed small-scale asthenospheric instabilities that become more pronounced with increasing Rayleigh number. These small-scale instabilities spread out the surface heat flux, resulting in smaller heat flux variations with increasing Rayleigh number. Scaled to Io's Rayleigh number of O(10(exp 12)) variations of order a few percent are expected. This small but significant variation in surface heat flux may be compatible with the observed distributions of volcanic centers and mountains, which appear fairly uniform at first sight but display a discernible distribution when suitably processed. The observed distribution of volcanic centers is similar to the asthenosphere heating distribution, implying that most of the tidal heating in Io occurs in an asthenosphere.

Tackley, Paul J.↗
NASA NTRS1997

Aeolian Processes and Features on Venus

Aeolian features on Venus include dune fields, eroded hills (yardangs), wind streaks, (miniature dunes of 10 to 30 cm wavelength). Although and possibly microdunes (in repetitive imaging by Magellan did show changes in the appearance of the surface, these changes are attributed to radar artifacts as a consequence of look direction rather than to physical changes of the surface. Nonetheless, measurements of wind speeds near the surface of Venus and wind tunnel simulations suggest that aeolian processes could be currently active on Venus. Study of radar images of terrestrial analogs shows that radar wavelength, polarization, and viewing geometry, including look direction and incidence angle, all influence the detection of dunes, yardangs, and wind streaks. For best detection, dune crests and yardangs should be oriented perpendicular to look direction. Longer wavelength systems can penetrate sand sheets a meter or more thick, rendering them invisible, especially in arid regions. For wind streaks to be visible, there must be a contrast in surface properties between the streak and the background on which it occurs. Nonetheless, more than 6000 aeolian features have been found on Magellan images of Venus, the most common of which are various wind streaks. Mapping wind streak orientations enables near-surface wind patterns to be inferred for the time of their formation. Type P streaks are associated with parabolic ejecta crater deposits and are considered to have formed in association with the impact event. Most Type P streaks are oriented westward, indicative of the upper altitude superrotation winds of Venus. Non Type P streaks have occurrences and orientations consistent with Hadley circulation. Some streaks in the southern hemisphere are oriented to the northeast, suggesting a Coriolis effect.

Greeley, Ronald↗
NASA NTRS1995

Spatial symmetry breaking in rapidly rotating convective spherical shells

Many problems in geophysical and astrophysical convection systems are characterized by fast rotation and spherical shell geometry. The combined effects of Coriolis forces and spherical shell geometry produce a unique spatial symmetry for the convection pattern in a rapidly rotating spherical shell. In this paper, we first discuss the general spatial symmetries for rotating spherical shell convection. A special model, a spherical shell heated from below, is then used to illustrate how and when the spatial symmetries are broken. Symmetry breaking occurs via a sequence of spatial transitions from the primary conducting state to the complex multiple-layered columnar structure. It is argued that, because of the dominant effects of rotation, the sequence of spatial transitions identified from this particular model is likely to be generally valid. Applications of the spatial symmetry breaking to planetary convection problems are also discussed.

Zhang, Keke↗
NASA NTRS1995

Mantle plume interaction with an endothermic phase change

High spatial resolution numerical simulations of mantle plumes impinging from below on the endothermic phase change at 660-km depth are used to investigate the effects of latent heat release on the plume-phase change interaction. Both axisymmetric and planar upflows are considered, and the strong temperature dependence of mantle viscosity is taken into account. For plume strengths considered, a Clapeyron slope of -4 MPa/K prevents plume penetration of the phase change. Plumes readily penetrate the phase change for a Clapeyron slope of -2 MPa/K and arrive in the upper mantle considerably hotter than if they had not traversed the phase change. For the same amount of thermal drive, i.e., the same excess basal temperature, axisymmetric plumes are hotter upon reaching the upper mantle than are planar upwellings. Heating of plumes by their passage through the spinel-perovskite endothermic phase change can have important consequences for the ability of the plume to thermally thin the lithosphere and cause melting and volcanism.

Schubert, Gerald↗
NASA NTRS1995

Lithospheric thickness and mantle/lithosphere density contrast beneath Beta Rigio, Venus

The spatial variation of the geoid/topography ratio over the large Venusian volcanic highland Beta Regio is suggestive of thermal compensation, i.e., support of the highland's topography by lithospheric thinning. Both the thickness of the lithosphere and the density contrast at its base can be inferred from a quadratic regression of suitably filtered (600 km less than wavelength less than 4000 km) geoid vs. topography data. The regression yields a mean lithospheric thickness of 270 km and a density contrast of magnitude 2.5% to 3.0%. Simple isostatic balance of the long-wavelength topography at Beta Regio requires thinning of the lithosphere by 50-60% beneath the rise.

Moore, William B.↗
NASA NTRS1995

Wind-Related Features and Processes on Venus Summary of Magellan Results

A search of Magellan synthetic aperture radar images covering about 98% of the venusian surface shows that aeolian features occur at all longitudes and latitudes. A global data base for wind streaks, the most common type of aeolian feature, was developed. For each of the 5970 streaks in the data base, information was compiled on location, streak type, radar backscatter, dimensions, azimuth, orientation with respect to local slope, and type of landform with which it is associated. In addition, streaks occurring in association with parabolic ejecta deposits were designated type P streaks, which constitute about 31% of the data base. Wind streak azimuths were analyzed to assess wind patterns at the time of their formation. Both hemispheres show strong westward and equatorward trends in azimuths, consistent with Hadley circulation and inferred upper atmospheric westward zonal winds. When type P streaks (those considered to result from transient impact events) were removed, the westward component was greatly reduced, suggesting that the upper zonal winds do not extend to the surface. The presence of equator-oriented streaks at high latitudes suggests that Hadley circulation extends to the poles. A field of possible yardangs found southwest of Mead Crater strikes NE-SW and occupies plains situated in a shallow topographic depression. Analysis of non-type P streaks in the area suggests that equatorward winds are funneled through the depression and are responsible for the erosion of the terrain to form the yardangs. Dune deposits are limited on Venus. Two dune fields were identified (Aglonice and Fortuna-Meshkenet) which total in area about 18,300 square km. Microdunes are proposed for some southern hemisphere areas which show distinctive radar reflectivities. Bragg scattering and/or subpixel reflections from the leeward faces of microdune bedforms could account for the unusual radar backscatter cross sections.

Greeley, Ronald↗
NASA NTRS1995

Wind-Related Features and Processes on Venus: Summary of Magellan Results

A search of Magellan synthetic aperture radar images covering approximately 98% of the venusian surface shows that aeolian features occur at all longitudes and latitudes. A global data base for wind streaks, the most common type of aeolian feature, was developed. For each of the 5970 streaks in the data base, information was compiled on location, streak type, radar backscatter, dimensions, azimuth, orientation with respect to local slope, and type of landform with which it is associated. In addition, streaks occurring in association with parabolic ejecta deposits were designated type P streaks, which constitute about 31% of the data base. Wind streak azimuths were analyzed to assess wind patterns at the time of their formation. Both hemispheres show strong westward and equatorward trends in azimuths, consistent with Hadley circulation and inferred upper atmospheric westward zonal winds. When type P streaks (those considered to result from transient impact events) were removed, the westward component was greatly reduced, suggesting that the upper zonal winds do not extend to the surface. The presence of equator-oriented streaks at high latitudes suggests that Hadley circulation extends to the poles. A field of possible yardangs found southwest of Mead Crater strikes NE-SW and occupies plains situated in a shallow topographic depression. Analysis of non-type P streaks in the area suggests that equatorward winds are funneled through the depression and are responsible for the erosion of the terrain to form the yardangs. Dune deposits are limited on Venus. Two dune fields were identified (Aglonice and Fortuna-Meshkenet) which total in area about 18,300 sq km. Microdunes are proposed for some southern hemisphere areas which show distinctive radar reflectivities. Bragg scattering and/or subpixel reflections from the leeward faces of microdune bedforms could account for the unusual radar backscatter cross sections.

Greeley, Ronald↗
NASA NTRS1994

Modeling of periodic great earthquakes on the San Andreas fault: Effects of nonlinear crustal rheology

We analyze the cycle of great earthquakes along the San Andreas fault with a finite element numerical model of deformation in a crust with a nonlinear viscoelastic rheology. The viscous component of deformation has an effective viscosity that depends exponentially on the inverse absolute temperature and nonlinearity on the shear stress; the elastic deformation is linear. Crustal thickness and temperature are constrained by seismic and heat flow data for California. The models are for anti plane strain in a 25-km-thick crustal layer having a very long, vertical strike-slip fault; the crustal block extends 250 km to either side of the fault. During the earthquake cycle that lasts 160 years, a constant plate velocity v(sub p)/2 = 17.5 mm yr is applied to the base of the crust and to the vertical end of the crustal block 250 km away from the fault. The upper half of the fault is locked during the interseismic period, while its lower half slips at the constant plate velocity. The locked part of the fault is moved abruptly 2.8 m every 160 years to simulate great earthquakes. The results are sensitive to crustal rheology. Models with quartzite-like rheology display profound transient stages in the velocity, displacement, and stress fields. The predicted transient zone extends about 3-4 times the crustal thickness on each side of the fault, significantly wider than the zone of deformation in elastic models. Models with diabase-like rheology behave similarly to elastic models and exhibit no transient stages. The model predictions are compared with geodetic observations of fault-parallel velocities in northern and central California and local rates of shear strain along the San Andreas fault. The observations are best fit by models which are 10-100 times less viscous than a quartzite-like rheology. Since the lower crust in California is composed of intermediate to mafic rocks, the present result suggests that the in situ viscosity of the crustal rock is orders of magnitude less the rock viscosity determined in the laboratory.

Reches, Ze'ev↗
NASA NTRS1994

Gravity over coronae and chasmata on Venus

The global spherical harmonic model of Venus' gravity field MGNP60FSAAP, with horizontal resolution of about 600 km, shows that most coronae have little or no signature in the gravity field. Nevertheless, some coronae and some segments of chasmata are associated with distinct positive gravity anomalies. No corona has been found to have a negative gravity anomaly. The spatial coincidence of the gravity highs over four closely spaced 300- to 400-km-diameter coronae in Eastern Eistla Regio with the structures themselves is remarkable and argues for a near-surface or lithospheric origin of the gravity signals over such relatively small features. Apparent depths of compensation (ADCs) of the prominent gravity anomalies at Artemis, Latona, and Heng-o Coronae are about 150 to 200 km. The geoid/topography ratios (GTRs) at Artemis, Latona, and Heng-o Coronae lie in the range 32 to 35 m/km. The large ADCs and GTRs of Artemis, Latona, and Heng-o Coronae are consistent with topographically related gravity and a thick Venus lithosphere or shallowly compensated topography and deep positive mass anomalies due to subduction of underthrusting at these coronae. At arcuate segments of Hecate and Parga Chasmata ADCs are about 125 to 150 km, while those at Fauta Corona, four coronae in Eastern Eistla Regio, and an arcuate segment of Wester Parga Chasmata are about 75 km. The GTRs at Fauta Corona, the four coronae in eastern Eistla Regio, and the accurate segments of Hecate, Parga, and Western Parga Chasmata are about 12 to 21 m/km. By analogy with gravity anomalies of similar horizontal scale (600 km-several thousand kilometers) on the concave sides of terrestrial subduction zone arcs, which are due in large part to subducted lithosphere, it is inferred that the gravity anomalies on Venus are consistent with retrograde subduction at Artemis Chasma, along the northern and southern margins of Latona Coronam, and elsewhere along Parga and Hecate Chasmata.

Schubert, Gerald↗
NASA NTRS1994

Venus' center of figure-center of mass offset

Magellan altimetry data reveal that the center of figure (CF) of Venus is displaced approximately 280 m from its center of mass (CM) toward 4.4 deg S, 135.8 deg E, a location in Aphrodite Terra. This offset is smaller than those of other terrestrial planets but larger than the estimated error, which is no more than a few tens of meters. We examine the possibility that the CF-CM offset is related to specific geologic provinces on Venus by deriving three simple models for the offset: a thick-crust model, a hotspot model, and a thick-lithosphere model. The offset caused by a region of thick crust depends upon the region's extent, the crust-mantle density contrast, and the thickness of excess crust. A hotspot-related offset depends on the extent of the thermally anomalous region and the magnitude of the thermal anomaly. Offset due to a region of thick lithosphere depends upon the extent of the region, the average temperature contrast across the lithosphere, and the amount of excess lithosphere. We apply the three models to Venus plateau-shaped highlands, volcanic rises, and lowlands, respectively, in an attempt to match the observed CF-CM offset location and magnitude. The influence of most volcanic rises and of Ishtar Terra on the CF-CM offset must be quite small if we are to explain the direction of the observed offset. The lack of influence of volcanic rises can be explained if the related thermal anomalies are limited to a few hundred degrees or less and are plume-shaped (i.e., characterized by a flattened sublithospheric `head' with a narrow cylindrical feeder `tail'). The unimportance of Ishtar Terra is most easily explained if it lies atop a significant mantle downwelling.

Bindschadler, Duane L.↗
NASA NTRS1994

Characteristics of finite amplitude stationary gravity waves in the atmosphere of Venus

This paper extends the study of stationary gravity waves generated near the surface of Venus reported previously by Young et al. to include finite amplitude effects associated with large amplitude waves. Waves are forced near the surface of Venus by periodic forcing. The height-dependent profiles of static stability and mean wind in the Venus atmosphere play a very important role in the evolution of the nonlinear behavior of the waves, just as they do in the linear wave solutions. Certain wave properties are qualitatively consistent with linear wave theory, such as wave trapping, resonance, and wave evanescence for short horizontal wavelenghts. However, the finite amplitude solutions also exhibit many other interesting features. In particular, for forcing amplitudes representative of those that could be expected in mountainous regions such as Aphrodite Terra, waves generated near the surface can reach large amplitudes at and above cloud levels, with clear signatures in the circulation pattern. At still higher levels, the waves can reach large enough amplitude to break, unless damping rates above the clouds are sufficient to limit wave amplitude growth. Well below cloud levels the waves develop complex flow patterns as the result of finite amplitude wave-wave interactions, and waves are generated having considerably shorter horizontal wavelenghts than that associated with the forcing near the surface. Nonlinear interactions can excite waves that are resonant with the background wind and static stability fields even when the primary surface forcing does not, and these waves can dominate the wave spectrum near cloud levels. A global map of Venus topographic slopes derived from Magellan altimetry data shows that slopes of magnitude comparable to or exceeding that used to force the model are ubiquitous over the surface.

Young, Richard E.↗
NASA NTRS1994

Spatial and temporal relations between coronae and extensional belts, northern Lada Terra, Venus

Preliminary studies of the distribution of coronae and volcanic rises on Venus show that many of these features tend to cluster along zones of rifting and extension. The plains north of Lada Terra are crossed by two such extensional belts. Each belt is composed of grabens, ridges, faults, volcanic flows, coronae and coronalike features. The longer and more prominent belt is the NW trending Alpha-Lada extensional belt, which is over 6000 km long and 50-200 km wide, and includes the coronae Eve, Tamfana, Carpo, Selu, Derceto, Otygen, and an unnamed corona south of Otygen. The second belt is the NNE trending Derceto-Quetzalpetlatl extensional belt, which is about 2000 km long and in places over 300 km wide, and includes the coronae Sarpanitum, Eithinoha, and Quetzalpetlatl. The two belts intersect at the 1600 x 600 km wide Derceto volcanic plateau. It is apparent that deformation along the two belts overlapped in time, though deformation along the Alpha-Lada extensional belt probably continued after the deformation along the Derceto-Quetzalpetlatl extensional belt terminated. In certain areas, volcanism originated in grabens within the extensional belts, whereas in other areas, such as in Eve, Selu, Derceto, and Quetzalpetlatl, volcanism originated in the coronae and flowed into the lower parts of the extensional belts. Regional extension has affected the evolution of all the coronae at some stage of their development. Regional deformation occurred before the initiation of Derceto and Eithinoha of their development. Regional deformation occurred before the initiation of Derceto and Eithinoha and after the initiation of Carpo, Tamfana, Otygen, and Sarpanitum. It is thus unlikely that coronae formation along the belts is solely a consequence of the regional extension, and it is also unlikely that regional extension has been caused solely by the coronae. No corona along the belts was formed subsequent to the cessation of the regional extension. We therefore suggest that the regional extension and the coronae are interrelated. Some of the coronae may have determined the location of the surface expression of the regional extension, whereas the locations of other coronae may have been influenced by the concentration of regional extensional stresses.

Baer, Gidon↗
NASA NTRS1994

Wind Streaks on Venus: Clues to Atmospheric Circulation

Magellan images reveal surface features on Venus attributed to wind processes. Sand dunes, wind-sculpted hills, and more than 5830 wind streaks have been identified. The streaks serve as local "wind vanes," representing wind direction at the time of streak formation and allowing the first global mapping of near-surface wind patterns on Venus. Wind streaks are oriented both toward the equator and toward the west. When streaks associated with local transient events, such as impact cratering, are deleted, the westward component is mostly lost but the equatorward component remains. This pattern is consistent with a Hadley circulation of the lower atmosphere.

Greeley, Ronald↗
NASA NTRS1993

The spatial distribution of coronae and related features on Venus

Coronae are large quasi-circular geologica features that are common on Venus. They appear to be the surface tectonic and volcanic expressions of mantle diapirs that have impinged on the underside of the venusian lithosphere. We have investigated the spatial distribution of 335 coronae and related features identified in Magellan radar data. It is more clustered than a Poisson distribution, with a statistical certainty of more than 99%. It is dominated by a single large cluster centered near the equator at about 245 deg longitude. The features are preferentially found at elevation and geoid values close to the planetary mean, with a paucity at both the highest and lowest levels of topgraphy and geoid. Some coronae appear aligned in quasi-linear chains. We attribute the clustering of coronae and related features to preferential formation of these features above regions of broad-scale mantle upwelling, and suggest that a major mantle upwelling underlies the one large cluster. We suggest that coronae are rare at the lowest elevations because these may be regions of mantle downwelling. The shortage of coronae at the highest elevations may result both from obscuration by other intense tectonism there and from suppression there of their formation by an unusually thick crust. Corona chains may be produced by enhanced passive mantle uplift below failed or incipient rifts.

Squyres, Steven W.↗
NASA NTRS1993

On the penetration of the 660 km phase change by mantle downflows

We present a simple analytic model of the interaction of cold convective downwelling currents with an endothermic phase change. The model describes the ponding and lateral spreading of downflows along the phase transition interface. A simple comparison of the vertical forces on the ponding material provides a necessary condition for a downflow to penetrate the phase boundary. This condition is fundamentally dependent on the geometry of the downflow. For planar downwellings, the model predicts a minimum ponding time before the structure can penetrate the phase boundary. For columnar (axisymmetric) downflows, there is no minimum time of spreading required before penetration can proceed. The model thus provides an explanation for the observation that in numerical models of three-dimensional convection with an endothermic phase change, cylindrical downflows penetrate the phase interface while planar ones do not. Since descending slabs in the Earth's mantle display a wide spectrum of geometries between planar and cylindrical (given various trench curvatures, as well as intersections of two or more subduction zones), this phenomenon may explain, in part, why some slabs appear to extend into the lower mantle while others are deflected at the 660 km discontinuity.

Bercovici, David↗