Engineering topics
Grimm, Robert E.
Publications and source records attributed to Grimm, Robert E..
Lithospheric Dynamics of Mars: Water, Flow, and Failure
Some of the largest Martian erosive features were influenced by groundwater, and include valley networks, outflow channels, and possibly landslides. We argue that hydrothermal systems attending crustal formation processes were able to drive sufficient groundwater to the surface to form the Noachian southern highlands valley networks, which show a spatial correlation to crustal magnetic anomalies, also results of crustal formation. Hydrothermal activity is quantified through numerical simulations of convection in a porous medium due to the presence of a hot intruded magma chamber. The parameter space includes magma chamber depth, volume, aspect ratio, and host rock permeability and porosity. For permeabilities as low as l0(exp -17) sq m and intrusion volumes as low as 50 km , the total discharge due to intrusions building that part of the southern highlands crust associated with magnetic anomalies spans a comparable range as the inferred discharge from the overlying valley networks. The Hesperian circum-Chryse outflow channels are further manifestations of groundwater discharge and Clifford and Parker (2001) suggest that the large volumes of water required for their formation flows beneath a confining cryosphere from the South Pole where meltwater beneath an ice cap recharges a global aquifer. We argue that recharge occurs instead over the nearby Tharsis aquifer at high obliquity, assisted by cryosphere melting due to volcanic activity. Numerical simulations quantify the strength and duration of outflow discharge given either South Polar or Tharsis recharge. The contribution of South Pole recharge given Clifford and Parker aquifer properties is negligible compared to that of the initial Tharsis inventory. Tharsis recharge, despite the restrictions of improved aquifer properties, makes a significant contribution and, unlike South Pole recharge under the same conditions, fulfills discharge requirements. Groundwater may have influenced long run-out landslide formation in the Valles Marineris. We present simulations of Martian, terrestrial, and lunar landslides that gauge the role of pore fluid pressure in reproducing accurate geometries and run-out with frictional, Bingham, and fluidization rheologies. The results indicate that pore fluid is a necessary component of Martian landslide formation and we suggest scenarios that might explain its presence.
Origin of Magnetic Lineations on Mars
The magnetic lineations discovered by MGS (Mars Global Surveyor) have been considered to be evidence of early plate tectonics on Mars. However, the lineations approximately follow lines of latitude, i.e., small circles. This presents significant geometrical problems for plate-like spreading, particularly at high latitudes. However, the sublatitudinal orientation of the lineations is consistent with meridianal extension and perhaps limited crustal spreading due to a stress event centered near the geographic pole. We hypothesize that this event was the early formation of the crustal dichotomy through mantle-convective processes. This could have taken the form of a southern megaplume that formed the thick highlands crust or as subduction or downwelling in the north. Both would result in tensional stresses in the south that would form extensional fractures perpendicular to the CM-CF (center-of-mass/center-of-gravity) offset. The observed magnitude and distribution of magnetization indicates that crustal intrusion associated with this major mantle-convective event resulted in approximately 1000 km of extension in the Southern highlands. Subsequent spin-axis reorientation due to loss of crust in the north or gain of crust in the south brought the CM-CF offset into its present N-S alignment. A portion of the ancient valley networks observed in the southern highlands are spatially associated with crustal magnetism and are quantitatively shown to be consistent with hydrothermal discharge over crustal intrusions.
Solid-State Multimission Magnetometer (SSM(3)): Application to Groundwater Exploration on Mars
This report describes work to develop solid-state magnetometers using magnetoresistive thin films, low-frequency electric-field measurements, and methods for electromagnetic detection of water and ice in the subsurface of Mars.
Low-Frequency Electromagnetic Exploration for Groundwater on Mars
Water with even a small amount of dissolved solids has an electrical conductivity orders of magnitude higher than dry rock and is therefore a near-ideal exploration target on Mars for low frequency, diffusive electromagnetic methods. Models of the temperature- and frequency-dependent electrical properties of rock-ice-water mixtures are used to predict the electromagnetic response of the Martian subsurface. Detection of ice is difficult unless it is massively segregated. In contrast, liquid water profoundly affects soundings, and even a small amount of adsorbed water in the cryosphere can be detected. Subcryospheric water is readily distinguishable at frequencies as low as 100 Hz for fresh water to 10 mHz for brines. These responses can be measured using either natural or artificial sources. Ultra low frequency signals from solar wind and diurnal-heating perturbations of the ionosphere are likely, and disturbances of regional crustal magnetic fields may also be observable. Spherics, or extremely to very low frequency signals from lightning discharge, would provide optimal soundings; however, lightning may be the least likely of the possible natural sources. Among the active techniques, only the time-domain electromagnetic (TDEM) method can accommodate a closely spaced transmitter and receiver and sound to depths of hundreds of meters or more. A ground- or aircraft-based TDEM system of several kilograms can detect water to a depth of several hundred meters, and a system of tens of kilograms featuring a large, fixed, rover- or ballistically deployed loop can detect water to several kilometers depth.
Controls on Martian Hydrothermal Systems: Application to Valley Network and Magnetic Anomaly Formation
Models of hydrothermal groundwater circulation can quantify limits to the role of hydrothermal activity in Martian crustal processes. We present here the results of numerical simulations of convection in a porous medium due to the presence of a hot intruded magma chamber. The parameter space includes magma chamber depth, volume, aspect ratio, and host rock permeability and porosity. A primary goal of the models is the computation of surface discharge. Discharge increases approximately linearly with chamber volume, decreases weakly with depth (at low geothermal gradients), and is maximized for equant-shaped chambers. Discharge increases linearly with permeability until limited by the energy available from the intrusion. Changes in the average porosity are balanced by changes in flow velocity and therefore have little effect. Water/rock ratios of approximately 0.1, obtained by other workers from models based on the mineralogy of the Shergotty meteorite, imply minimum permeabilities of 10(exp -16) sq m2 during hydrothermal alteration. If substantial vapor volumes are required for soil alteration, the permeability must exceed 10(exp -15) sq m. The principal application of our model is to test the viability of hydrothermal circulation as the primary process responsible for the broad spatial correlation of Martian valley networks with magnetic anomalies. For host rock permeabilities as low as 10(exp -17) sq m and intrusion volumes as low as 50 cu km, the total discharge due to intrusions building that part of the southern highlands crust associated with magnetic anomalies spans a comparable range as the inferred discharge from the overlying valley networks.
Recent deformation rates on Venus
Constraints on the recent geological evolution of Venus may be provided by quantitative estimates of the rates of the principal resurfacing processes, volcanism and tectonism. This paper focuses on the latter, using impact craters as strain indicators. The total postimpact tectonic strain lies in the range 0.5-6.5%, which defines a recent mean strain rate of 10(exp -18)-10(exp -17)/s when divided by the mean surface age. Interpretation of the cratering record as one of pure production requires a decline in resurfacing rates at about 500 Ma (catastrophic resurfacing model). If distributed tectonic resurfacing contributed strongly before that time, as suggested by the widespread occurrence of tessera as inliers, the mean global strain rate must have been at least approximately 10(exp -15)/s, which is also typical of terrestrial active margins. Numerical calculations of the response of the lithosphere to inferred mantle convective forces were performed to test the hypothesis that a decrease in surface strain rate by at least two orders of magnitude could be caused by a steady decline in heat flow over the last billion years. Parameterized convection models predict that the mean global thermal gradient decreases by only about 5 K/km over this time; even with the exponential dependence of viscosity upon temperature, the surface strain rate drops by little more than one order of magnitude. Strongly unsteady cooling and very low thermal gradients today are necessary to satisfy the catastrophic model. An alternative, uniformitarian resurfacing hypothesis holds that Venus is resurfaced in quasi-random 'patches' several hundred kilometers in size that occur in response to changing mantle convection patterns.
The deep structure of Venusian plateau highlands
Magellan gravity data confirm that several of the large, tectonically deformed, plateau-like highlands on venus are shallowly compensated, most likely by crustal thickness variations. Apparent depths of isostatic compensation, computed in the spatial domain, range from 30 to 50 km for Alpha, Tellus, Ovda, and Thetis Regiones. Using a two-layer model for isostatic compensation, Alpha, Tellus, and Ovda are best represented as nearly completely compensated in crust that is regionally 20-40 km thick around these highlands, with little contribution from deeper mantle sources. In contrast to these three areas, a stronger regional gravity high associated with Thetis requires a significant upper mantle component to compensation. This is evident in the spectral admittance as a pronounced deep, long-wavelength anomaly. In the two-layer isostatic model, a broad, deeply compensated upland underlies a shallowly compensated central block of Thetis. If this deep component is interpreted as a thermal anomaly, the loci of maximum upwelling agree well with sites of recent extension. The plateau highlands are thus physiographically and isostatically equivalent to terrestrial continents, though probably not compositionally. They also share the record of a long tectonic history. The large regional gravity anomaly of Thetis indicates that active mantle proceses continue even beneath some areas (tessera) thought to be a relic of a former geological regime. The excellent agreement of modeled crustal thicknesses around Alpha, Tellus, and Ovda Regiones suggests that 20-40 km is a representative global value for the plains. Such a crust is thicker than previously estimated and about twice as thick as the expected thickness of crust produced at venusian spreading centers
Shapes of Venusian 'pancake' domes imply episodic emplacement and silicic composition
The main evidence available for constraining the composition of the large circular 'pancake' domes on Venus is their gross morphology. Laboratory simulations using polyethylene glycol show that the height to diameter (aspect) ratios of domes of a given total volume depend critically on whether their extrusion was continuous or episodic, with more episodes leading to greater cooling and taller domes. Thus without observations of their emplacement, the compositions of Venusian domes cannot be uniquely constrained by their morphology. However, by considering a population of 51 Venusian domes to represent a sampling of many stages during the growth of domes with comparable histories, and by plotting aspect ratio versus total volume, we find that the shapes of the domes are most consistent with episodic emplacement. On Earth this mode of dome growth is found almost exclusively in lavas of dacite to rhyolite composition, strengthening earlier inferences about the presence of evolved magmas on Venus.
Flexure and the role of inplane force around coronae on Venus
Large coronae on Venus, such as Artemis and Latona, are rimmed by conspicuous trenches and associated outer rises. Sandwell and Schubert have observed that these systems resemble terrestrial subduction zones in planform and have succeeded in fitting an elastic plate bending equation to the inferred flexural topography. However, the first zero crossing bending moments required are -2.5 x 10(exp 17) N for Artemis and -5.0 x 10(exp 16) N for Latona. Since these moments are similar in magnitude to those of subducting slabs on Earth, a rollback subduction mechanism was proposed to explain the flexure around the largest coronae, although a differential thermal subsidence model is sufficient to account for the topography around some coronae. The purpose is to investigate the effect of inplane force as a possible alternative to large applied moments in producing flexure at Artemis and Latona. The close correlation of gravity to topography on Venus implies the absence of a low viscosity zone and the strong coupling of the lithosphere to mantle convection. If coronae are the surface manifestations of mantle plumes, they may be the sites of active convective stress coupling. As the upwelling reaches the lithosphere, it spreads radially outward, inducing shear tractions on the base of the plate. In addition, the hot, expanding corona may load the surrounding plate horizonally. Both the basal shear stresses and radial loading can be treated as an equivalent compressive inplane force in the mechanical lithosphere, which contributes to the bending of the outlying plate. Using a model that relates inplane force to the measured gravity anomalies, a rough value of the inplane force at Artemis was calculated. Recent Pioneer Venus spherical harmonic gravity models indicate a geoid anomaly of about 75 m over Artemis, which corresponds to an estimated inplane force on the order of -1x10(exp 13) N/m. The gravity model is unable to resolve Latona, but an inplane force of similar dimensions is assumed. The maximum possible inplane force based on the expected rheology can be constrained by using the approximate 5 K/km thermal gradient inferred from the best fit 30 km elastic plate at Artemis and Latona. For a dry olivine flow law in the upper mantle, the compressional load limit of the 60 km thick mechanical lithosphere is -4 x 10(exp 13) N/m. This value is equivalent to a load of -8 x 10(exp 13) N/m on a 30 km thick elastic plate.
Viscous relaxation of the Moho under large lunar basins
Viscously relaxed topography on the Moon is evidence of a period in lunar history of higher internal temperatures and greater surface activity. Previous work has demonstrated the viscous relaxation of the Tranquilitatis basin surface. Profiles of the lunar Moho under nine basins were constructed from an inversion of lunar gravity data. These profiles show a pattern of increasingly subdued relief with age, for which two explanations have been proposed. First, ancient basins may have initially had extreme Moho relief like that of younger basins like Orientale, but, due to higher internal temperatures in early lunar history, this relief viscously relaxed to that observed today. Second, ductile flow in the crust immediately after basin formation resulted in an initially shallow basin and subdued mantle uplift. The intent is to test the first hypothesis.
Tectonic resurfacing of Venus
Impact crater distributions and morphologies have traditionally played an important role in unraveling the geologic histories of terrestrial objects, and Venus has proved no exception. The key observations are: mean crater retention age about 500 Ma; apparently random spatial distribution; modest proportion (17 percent) of modified craters; and preferential association of modified craters with areas of low crater density. The simplest interpretation of these data alone is that Venus experienced global resurfacing (assumed to be largely volcanic) prior to 500 Ma, after which time resurfacing rates decreased dramatically. This scenario does not totally exclude present geological activity: some resurfacing and crater obliteration is occurring on part of the planet, but at rates much smaller than on Earth. An alternative endmember model holds that resurfacing is also spatially randomly distributed. Resurfacing of about 1 sq km/yr eliminates craters such that a typical portion of the surface has an age of 500 Ma, but actual ages range from zero to about 1000 Ma. Monte Carlo simulation indicates that the typical resurfacing 'patch' cannot exceed about 500 km in diameter without producing a crater distribution more heterogeneous than observed. Volcanic or tectonic processes within these patches must be locally intense to be able to obliterate craters completely and leave few modified. In this abstract, we describe how global geologic mapping may be used to test resurfacing hypotheses. We present preliminary evidence that the dominant mode of resurfacing on Venus is tectonism, not volcanism, and that this process must be ongoing today. Lastly, we outline a conceptual model in which to understand the relationship between global tectonics and crater distribution and preservation.
Impact craters and Venus resurfacing history
The history of resurfacing by tectonism and volcanism on Venus is reconstructed by means of an analysis of Venusian impact crater size-frequency distributions, locations, and preservation states. An atmospheric transit model for meteoroids demonstrates that for craters larger than about 30 km, the size-frequency distribution is close to the atmosphere-free case. An age of cessation of rapid resurfacing of about 500 Ma is obtained. It is inferred that a range of surface ages are recorded by the impact crater population; e.g., the Aphrodite zone is relatively young. An end-member model is developed to quantify resurfacing scenarios. It is argued that Venus has been resurfacing at an average rate of about 1 sq km/yr. Numerical simulations of resurfacing showed that there are two solution branches that satisfy the completely spatially random location restraint for Venusian craters: a is less than 0.0003 (4 deg diameter circle) and a is greater than 0.1 (74 deg diameter circle).
Anatomy of a Venusian hot spot - Geology, gravity, and mantle dynamics of Eistla Regio
Results of a study of the western and central portions of the Venusian hot spot Eistla Regio are presented. Magellan radar images were mapped to elucidate the general geologic history of the region. Radial fracture systems both on the rises and volcanoes indicate that uplift and associated faulting accompanied volcanic construction. Prominent fracture zones strike WNW to NW, parallel to the long axis of the highlands. The largest of these, Guor Linea, exhibits a progressive deformation history that may include minor clockwise rotation in addition to bulk NNE-SSW extension. Pioneer Venus line-of-sight accelerations were inverted for vertical gravity which, when combined with topography, were used to solve for mass anomalies on the crust-mantle boundary and in the upper levels of the mantle convective system.
Styles of deformation in Ishtar Terra and their implications
Styles of deformation in Ishtar Terra are examined on the basis of imaging by the Magellan radar between the start of the mission and the start of the first superior conjunction hiatus. Ishtar Terra appears to have characteristics of both plume uplifts and convergent belts, and exhibits a great variety of tectonic and volcanic activity, with large variations within distances of only a few hundred kilometers. The most prominent terrain types are the volcanic plains of Lakshmi and the mountain belts of Maxwell, Freyja, and Danu. Ishtar demonstrates three general properties of Venus. Erosional degradation is absent, leading to the preservation of patterns resulting from past activity. Many surface features are the responses of a competent layer less than 10 km thick to flows of 100 km or broader scale. These broader scale flows are controlled mainly by heterogeneities in the mantle.
Venus tectonics - An overview of Magellan observations
Magellan observations of the tectonic characteristics of highland regions on Venus are discussed with reference to competing theories for highland formation and evolution. Complex rigid terrain, or tessera, and the extent to which these elevated blocks of intensely deformed crust may be genetically related to highlands are then considered. Further, the tectonics of plains and lowland regions are examined, including deformation belts and coronae, and possible relations between such features and mantle dynamics. Implications of these observations for the global tectonics of Venus are discussed.
Martian seismicity
The design and ultimate success of network seismology experiments on Mars depends on the present level of Martian seismicity. Volcanic and tectonic landforms observed from imaging experiments show that Mars must have been a seismically active planet in the past and there is no reason to discount the notion that Mars is seismically active today but at a lower level of activity. Models are explored for present day Mars seismicity. Depending on the sensitivity and geometry of a seismic network and the attenuation and scattering properties of the interior, it appears that a reasonable number of Martian seismic events would be detected over the period of a decade. The thermoelastic cooling mechanism as estimated is surely a lower bound, and a more refined estimate would take into account specifically the regional cooling of Tharsis and lead to a higher frequency of seismic events.
Hot-spot tectonics of Eistla Regio, Venus: Results from Magellan images and Pioneer Venus gravity
Eistla Regio (ER) is a broad, low, discontinuous topographic rise striking roughly EW at low northern latitudes of Venus. Some 2000 x 7000 km in dimensions, it is the third largest rise in planform on Venus after Aphrodite Terra and Beta Phoebe Regiones. These rises are the key physiographic elements in a hot spot model of global tectonics including transient plume behavior. Since ER is the first such rise viewed by Magellan and the latitude is very favorable for Pioneer Venus gravity studies, some of the predictions of a time dependent hot spot model are tested. Western ER is defined as the rise including Gula and Sif Mons and central ER as that including Sappho Patera. Superior conjunction prevented Magellan from returning data on eastern ER (Pavlova) during the first mapping cycle. It is concluded that the western and central portions of ER, while part of the same broad topographic rise and tectonic framework, have distinctly different surface ages and gravity signatures. The western rise, including Gula and Sif Mons, is the expression of deep seated uplift with volcanism limited to the individual large shields. The eastern portion has been widely resurfaced more recently by thermal anomalies in the mantle.