Search NASA⌕ Search

Engineering topics

Schubert, Gerald

Publications and source records attributed to Schubert, Gerald.

At least 19 records

A fully self-consistent multi-layered model of jupiter

We construct a three-dimensional, fully self-consistent, multi-layered, non-spheroidal model of Jupiter consisting of an inner core, a metallic electrically conducting dynamo region, and an outer molecular electrically insulating envelope. We assume that the Jovian zonal winds are on cylinders parallel to the rotation axis but, due to the effect of magnetic braking, are confined within the outer molecular envelope. We also assume that the location of the molecular-metallic interface is characterized by its equatorial radius HR{sub e}, where R {sub e} is the equatorial radius of Jupiter at the 1 bar pressure level and H is treated as a parameter of the model. We solve the relevant mathematical problem via a perturbation approach. The leading-order problem determines the density, size, and shape of the inner core, the irregular shape of the 1 bar pressure level, and the internal structure of Jupiter that accounts for the full effect of rotational distortion, but without the influence of the zonal winds; the next-order problem determines the variation of the gravitational field solely caused by the effect of the zonal winds on the rotationally distorted non-spheroidal Jupiter. The leading-order solution produces the known mass, the known equatorial and polar radii, and the known zonal gravitational coefficient J {sub 2} of Jupiter within their error bars; it also yields the coefficients J {sub 4} and J {sub 6} within about 5% accuracy, the core equatorial radius 0.09R{sub e} and the core density ρ{sub c}=2.0×10{sup 4} kg m{sup −3} corresponding to 3.73 Earth masses; the next-order solution yields the wind-induced variation of the zonal gravitational coefficients of Jupiter.

79 ASTRONOMY AND ASTROPHYSICS↗

THERMAL-GRAVITATIONAL WIND EQUATION FOR THE WIND-INDUCED GRAVITATIONAL SIGNATURE OF GIANT GASEOUS PLANETS: MATHEMATICAL DERIVATION, NUMERICAL METHOD, AND ILLUSTRATIVE SOLUTIONS

The standard thermal wind equation (TWE) relating the vertical shear of a flow to the horizontal density gradient in an atmosphere has been used to calculate the external gravitational signature produced by zonal winds in the interiors of giant gaseous planets. We show, however, that in this application the TWE needs to be generalized to account for an associated gravitational perturbation. We refer to the generalized equation as the thermal-gravitational wind equation (TGWE). The generalized equation represents a two-dimensional kernel integral equation with the Green’s function in its integrand and is hence much more difficult to solve than the standard TWE. We develop an extended spectral method for solving the TGWE in spherical geometry. We then apply the method to a generic gaseous Jupiter-like object with idealized zonal winds. We demonstrate that solutions of the TGWE are substantially different from those of the standard TWE. We conclude that the TGWE must be used to estimate the gravitational signature of zonal winds in giant gaseous planets.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

A THREE-DIMENSIONAL NUMERICAL SOLUTION FOR THE SHAPE OF A ROTATIONALLY DISTORTED POLYTROPE OF INDEX UNITY

We present a new three-dimensional numerical method for calculating the non-spherical shape and internal structure of a model of a rapidly rotating gaseous body with a polytropic index of unity. The calculation is based on a finite-element method and accounts for the full effects of rotation. After validating the numerical approach against the asymptotic solution of Chandrasekhar that is valid only for a slowly rotating gaseous body, we apply it to models of Jupiter and a rapidly rotating, highly flattened star ({alpha} Eridani). In the case of Jupiter, the two-dimensional distributions of density and pressure are determined via a hybrid inverse approach by adjusting an a priori unknown coefficient in the equation of state until the model shape matches the observed shape of Jupiter. After obtaining the two-dimensional distribution of density, we then compute the zonal gravity coefficients and the total mass from the non-spherical model that takes full account of rotation-induced shape change. Our non-spherical model with a polytropic index of unity is able to produce the known mass of Jupiter with about 4% accuracy and the zonal gravitational coefficient J {sub 2} of Jupiter with better than 2% accuracy, a reasonable result considering that there is only one parameter in the model. For {alpha} Eridani, we calculate its rotationally distorted shape and internal structure based on the observationally deduced rotation rate and size of the star by using a similar hybrid inverse approach. Our model of the star closely approximates the observed flattening.

79 ASTRONOMY AND ASTROPHYSICS↗

ON THE GRAVITATIONAL FIELDS OF MACLAURIN SPHEROID MODELS OF ROTATING FLUID PLANETS

Hubbard recently derived an important iterative equation for calculating the gravitational coefficients of a Maclaurin spheroid that does not require an expansion in a small distortion parameter. We show that this iterative equation, which is based on an incomplete solution of the Poisson equation, diverges when the distortion parameter is not sufficiently small. We derive a new iterative equation that is based on a complete solution of the Poisson equation and, hence, always converges when calculating the gravitational coefficients of a Maclaurin spheroid.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

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

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↗

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↗

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↗

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

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↗

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↗

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↗

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↗

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

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