DESDynI: Deformation, Ecosystem Structure and Dynamics of Ice
No abstract available
Engineering topics
Publications and source records attributed to Hager, Bradford H..
No abstract available
Air and space borne remote sensing have made it possible to monitor the mass and energy transport at various scales within the cryosphere-hydrosphere-atmosphere system. The recent surface mass balance (the rate of net gain of snow and ice at a geographic point) map for the Antarctic ice sheet is constructed by interpolating sparse in situ observations (about 1,800 points) with empirically calibrated satellite data of passive back emission of microwaves. The digital elevation model obtained from satellite radar altimetry is used to improve the delineation of the ice flow drainage basins. As important as these results are, the uncertainty remains up to about 2 mm/yr of eustatic sea level change with the net imbalance. In other words, we are still unable to determine even the sign of the contribution of the Antarctic ice sheet to contemporary sea level change. The problem is more likely with the discharge rather than accumulation.
In this report, we summarize what we have accomplished with five years of funding from NASA under its DOSE program, and with a comparable level of funding from NSF. We describe the development of a GPS network in the Tien Shan of Kyrgyzstan and Kazakhstan of the former Soviet Union, the analysis of data, and the main results. This discussion presents the state of the current network, which has grown significantly since the termination of our DOSE grants, with continued support both from NSF through its continental dynamics program and from NASA's SENH program. Although grants from NASA's DOSE program did not support this growth not directly, it did so indirectly by building the infrastructure that has enabled further expansion in an area where otherwise there would be only a small GPS presence. We note how the network has grown over time, but the emphasis of this discussion is on the quantity and quality of measurements that we have made.
Measurements of elastic displacements of the bedrock surrounding large ice sheets have been proposed as a means to detect mass changes in these ice sheets. However, accumulation of glacial mass on the ice sheets is a noisy process, subject to large spatial and temporal variations in precipitation. We simulated the response of the Antarctic continent to a stochastic model of interannual precipitation variations and found that interannual variations in the elastic response of the earth are large when compared to the long-term mean of displacements produced by an assumed average ice mass imbalance of 10%. If, as some scientists predict, Antarctic ice mass changes in the future become dramatic, the long-term signal should be large enough to be detected by a few years of geodetic measurements, despite climatic noise.
The singularity problem associated with a radially continuous Maxwell viscoclastic structure is investigated. A special tool called the isolation function is developed. Results calculated using the isolation function show that the discrete model assumption is no longer valid when the viscoelastic parameter becomes a continuous function of radius. Continuous variations in the upper mantle viscoelastic parameter are especially powerful in destroying the mode-like structures. The contribution to the load Love numbers of the singularities is sensitive to the convexity of the viscoelastic parameter models. The difference between the vertical response and the horizontal response found in layered viscoelastic parameter models remains with continuous models.
There were two major developments in 1994 in our collaborative GPS experiment in the Tien Shan of the Former Soviet Union (FSU). Both were motivated by our expectation that we will ultimately obtain better science at lower cost if we involve our colleagues in the FSU more deeply in (1) the collection and (2) the analysis of data. As an experimental test of the concept of having our local collaborators carry out the field work semi-autonomously, we sent 6 MIT receivers to the Tien Shan for a period of 3 months. To enable our collaborators to have the capability for data analysis, we provided computers for two data analysis centers and organized a two-week training session. This report emphasizes the rationale for deeper involvement of FSU scientists, describes the training sessions, discusses the data collection, and presents the results. We also discuss future plans. More detailed discussion of background, general scientific objectives, discussions with collaborators, and results for the campaigns in 1992 and 1993 have been given in previous reports.
Pressure load Love numbers are presented for the mantle deformation induced by the variation of the pressure field at the core mantle boundary (CNB). We find that the CMB geostrophic pressure fields, derived from 'frozen-flux' core surface flow estimates at epochs 1965 and 1975, produce a relative radial velocity (RRV) field in the range of 3mm/decade with uplift near the equator and subsidence near the poles. The contribution of this mechanism to the change in the length of day (l.o.d) is small --- about 2.3 x 10(exp -2) ms/decade. The contribution to the time variation of the ellipticity coefficient is more important --- -1.3 x 10(exp -11)/yr.
Calculating the post glacial response of a viscoelastic Earth model using the exponential decay normal mode technique leads to intrinsic singularities if viscosity varies continuously as a function of radius. We develop a numerical implementation of the Complex Real Fourier transform (CRFT) method as an accurate and stable procedure to avoid these singularities. Using CRFT, we investigate the response of a set of Maxwell Earth models to surface loading. We find that the effect of expanding a layered viscosity structure into a continuously varying structure is to destroy the modes associated with the boundary between layers. Horizontal motion is more sensitive than vertical motion to the viscosity structure just below the lithosphere. Horizontal motion is less sensitive to the viscosity of the lower mantle than the vertical motion is. When the viscosity increases at 670 km depth by a factor of about 60, the response of the lower mantle is close to its elastic limit. Any further increase of the viscosity contrast at 670 km depth or further increase of viscosity as a continuous function of depth starting from 670 km depth is unlikely to be resolved.
We formulate 2-D Cartesian finite element models that explore the fate of compositionally defined lithosphere as it encounters a viscosity increase at the boundary between the upper and lower mantle. Subducted lithosphere is represented as a cold, stiff, layered composite of denser eclogite underlain by more buoyant harzburgite. Slabs impinging on a lower mantle 30 and 100 times more viscous than the upper mantle thicken and fold strongly as they penetrate the lower mantle. Approximately a factor of two thickening occurs via pure shear just above the discontinuity, with additional enhancement due to folding by over a factor of two. No separation of the individual slab components occurs at the discontinuity, and direct comparison with models in which compositional buoyancy is explicitly ignored indicates that slab evolution is largely controlled by the thermal buoyancy. These results are at odds with hypotheses about slab evolution in which the compositional buoyancy contributions lead to component separation and the formation of slab megaliths or a compositionally layered upper mantle.
We estimate the velocity field in central and southern California using Global Positioning System (GPS) observations from 1986 to 1992 and very long baseline interferometry (VLBI) observations from 1984 to 1991. Our core network includes 12 GPS sites spaced approximately 50 km apart, mostly in the western Transverse Ranges and the coastal Borderlands. The precision and accuracy of the relative horizontal velocities estimated for these core stations are adequately described by a 95% confidence ellipse with a semiminor axis of approximately 2 mm/yr oriented roughly north-south, and a semimajor axis of approximately 3 mm/yr oriented east-west.
We have measured the deformation in the Ventura basin region, southern California, with Global Positioning System (GPS) measurements carried out over 4.6 years between 1987 and 1992. The deformation within our network is spatically variable on scales of tens of kilometers, with strain rates reaching 0.6 +/- 1 micro-rad/yr in the east-central basin. Blocklike rotations are observed south and northwest of the basin where the maximum shear strain rates are an order of magnitude lower (0.06 +/- 1 micro-rad/yr to the south). We also observed clockwise rotations of 1 deg - 7 deg/m.y. Shear strain rates determined by comparing angle changes from historical triangulation spanning several decades and GPS measurements give consistent, though less precise, results. The geodetic rates of shortening across the basin and Western Transverse Ranges are lower than those estimated from geological observations, but the patterns of deformation from the two methods agree qualitatively.
Explore the source record for details and available documents.
In this study we present the results of a geodetic experiment designed to test whether the rates and spatial distribution of deformation measured with space geodesy are compatible with the estimated geologic rates.
High-resolution radar images from the Magellan spacecraft have revealed the first details of the morphology of the Lavinia Planitia region of Venus. A number of geologic units can be distinguished, including volcanic plains units with a range of ages. Transecting these plains over much of the Lavinia region are two types of generally orthogonal features that we interpret to be compressional wrinkle ridges and extensional grooves. The dominant tectonic features of Lavinia are broad elevated belts of intense deformation that transect the plains with complex geometry. They are many tens to a few hundred kilometers wide, as much as 1000 km long, and elevated hundreds of meters above the surrounding plains. Two classes of deformation belts are seen in the Lavinia region. 'Ridge belts' are composed of parallel ridges, each a few hundred meters in elevation, that we interpret to be folds. Typical fold spacings are 5-10 km. 'Fracture belts' are dominated instead by intense faulting, with faults in some instances paired to form narrow grabens. There is also some evidence for modest amounts of horizontal shear distributed across both ridge and fracture belts. Crosscutting relationships among the belts show there to be a range in belt ages. In western Lavinia, in particular, many ridge and fracture belts appear to bear a relationship to the much smaller wrinkle ridges and grooves on the surrounding plains: ridge morphology tends to dominate belts that lie more nearly parallel to local plains wrinkle ridges, and fracture morphology tends to dominate belts that lie more nearly parallel to local plains grooves. We use simple models to explore the formation of ridge and fracture belts. We show that convective motions in the mantle can couple to the crust to cause horizontal stresses of a magnitude sufficient to induce the formation of deformation belts like those observed in Lavinia. We also use the small-scale wavelengths of deformation observed within individual ridge belts to place an approximate lower limit on the venusian thermal gradient in the Lavinia region at the time of deformation.
Geodetic observations between 1970 and 1988 indicate appreaciable strain accumulation in the Santa Barbara Channel, California. Eleven line-length changes from a six-station geodetic network spanning the eastern two-thirds of the channel were determined from electronic distance measurements in 1970/1971 and Global Positioning System observations in 1987/1988. Within this network the strains observed are spatially nonuniform. In the easternmost channel the strain is nearly uniaxial, with convergence of 6.4 +/- 0.9 mm/yr oriented N 25 deg E +/- 5 deg; this direction is consistent with the seismicity, which is dominated by thrust mechanisms with P axes directed to the northeast. In the central channel the strain is less well determined, but appears to include a significant component of shear that is left-lateral when resolved on an east-west plane.
The first measurements are reported for a major earthquake by a continuously operating GPS network, the permanent GPS Genetic ARRY (PGGA) in southern California. The Landers and Big Bear earthquakes of June 28, 1992 were monitored by daily observations. Ten weeks of measurements indicate significant coseismic motion at all PGGA sites, significant postseismic motion at one site for two weeks after the earthquakes, and no significant preseismic motion. These measurements demonstrate the potential of GPS monitoring for precise detection of precursory and aftershock seismic deformation in the near and far field.
High-resolution Magellan images and altimetry of Venus reveal a wide range of styles and scales of surface deformation that cannot readily be explained within the classical terrestrial plate tectonic paradigm. The high correlation of long-wavelength topography and gravity and the large apparent depths of compensation suggest that Venus lacks an upper-mantle low-viscosity zone. A key difference between Earth and Venus may be the degree of coupling between the convecting mantle and the overlying lithosphere. Mantle flow should then have recognizable signatures in the relationships between the observed surface topography, crustal deformation, and the gravity field. Therefore, comparison of model results with observational data can help to constrain such parameters as crustal and thermal boundary layer thicknesses as well as the character of mantle flow below different Venusian features. We explore in this paper the effects of this coupling by means of a finite element modelling technique.
The belts of the Lavinia Planitia region of Venus are discussed in detail with reference to high-resolution radar images from the Magellan spacecraft. Two types of deformation belts are identified: ridge belts and fracture belts. Ridge belts are composed of parallel ridges, each a few hundred meters in elevation, that are interpreted as folds. Typical fold spacings are 5-10 km. Fracture belts are dominated by intense faulting, with faults in some instances paired to form narrow grabens. The formation of ridge and fracture belts is explored using simple models.