An elastic-perfectly plastic analysis of the bending of the lithosphere at a trench

A number of authors have modelled the flexure of the lithosphere at an oceanic trench using a thin elastic plate with a hydrostatic restoring force. In some cases good agreement with observed topography is obtained but in other cases the slope of the lithosphere within the trench is greater than that predicted by the elastic theory. In this paper the bending of a thin plate is considered using an elastic-perfectly plastic rheology. It is found that the lithosphere behaves elastically seaward of the trench, but that plasticity decreases the radius of curvature within the trench. The results are compared with a number of observed trench profiles. The elastic-perfectly plastic profiles are in excellent agreement with those profiles that deviate from elastic behavior.     

Sedimentary Deposits from the 17 July 2006 Western Java Tsunami,Indonesia: Use of Grain Size Analyses to Assess Tsunami Flow Depth,Speed, and Traction Carpet Characteristics

The 2006 western Java tsunami deposited a discontinuous sheet of sand up to 20 cm thick, flooded coastal southern Java to a depth of at least 8 m and inundated up to 1 km inland. In most places the primarily heavy mineral sand sheet is normally graded, and in some it contains complex internal stratigraphy. Structures within the sand sheet probably record the passage of up to two individual waves, a point noted in eyewitness accounts. We studied the 2006 tsunami deposits in detail along a flow parallel transect about 750 m long, 15 km east of Cilacap. The tsunami deposit first becomes discernable from the underlying sediment 70 m from the shoreline. From 75 to 300 m inland the deposit has been laid down in rice paddies, and maintains a thickness of 10–20 cm. Landward of 300 m the deposit thins dramatically, reaching 1 mm by 450 m inland. From 450 m to the edge of deposition (around 700 m inland) the deposit remains <1 mm thick. Deposition generally attended inundation—along the transect, the tsunami deposited sand to within about 40 m of the inundation limit. The thicker part of the deposit contains primarily sand indistinguishable from that found on the beach 3 weeks after the event, but after about 450 m (and roughly coinciding with the decrease in thickness) the tsunami sediment shifts to become more like the underlying paddy soil than the beach sand. Grain sizes within the deposit tend to fine upward and landward, although overall upward fining takes place in two discrete pulses, with an initial section of inverse grading followed by a section of normal grading. The two inversely graded sections are also density graded, with denser grains at the base, and less dense grains at the top. The two normally graded sections show no trends in density. The inversely graded sections show high density sediment to the base and become less dense upward and represents traction carpet flows at the base of the tsunami. These are suggestive of high shear rates in the flow. Because of the grain sorting in the traction carpet, the landward-fining trends usually seen in tsunami deposits are masked, although lateral changes of mean sediment grain size along the transect do show overall landward fining, with more variation as the deposit tapers off. The deposit is also thicker in the more seaward portions than would be produced by tsunamis lacking traction carpets.     

Variations in the accuracy of gravity recovery due to ground track variability: GRACE,CHAMP, and GOCE

Following an earlier recognition of degraded monthly geopotential recovery from GRACE (Gravity Recovery And Climate Experiment) due to prolonged passage through a short repeat (low order resonant) orbit, we extend these insights also to CHAMP (CHAllenging Minisatellite Payload) and GOCE (Gravity field and steady state Ocean Circulation Explorer). We show wide track-density variations over time for these orbits in both latitude and longitude, and estimate that geopotential recovery will be as widely affected as well within all these regimes, with lesser track density leading to poorer recoveries. We then use recent models of atmospheric density to estimate the future orbit of GRACE and warn of degraded performance as other low order resonances are encountered in GRACE’s free fall. Finally implications for the GOCE orbit are discussed.     

Headless submarine canyons and fluid flow on the toe of the Cascadia accretionary complex

Headless submarine canyons with steep headwalls and shallowly sloping floors occur on both the second and third landward vergent anticlines on the toe of the Cascadia accretionary complex off central Oregon (45 °N, 125° 30′W). In September 1993, we carried out a series of nine deep tow camera sled runs and nine ALVIN dives to examine the relationship between fluid venting, structure and canyon formation. We studied four canyons on the second and third landward vergent anticlines, as well as the apparently unfailed intercanyon regions along strike. All evidence of fluid expulsion is associated with the canyons; we found no evidence of fluid flow between canyons. Even though all fluid seeps are related to canyons, we did not find seeps in all canyons, and the location of the seeps within the canyons differed. On the landward facing limb of the second landward vergent anticline a robust cold seep community occurs at the canyon’s inflection point. This seep is characterized by chemosynthetic vent clams, tube worms and extensive authigenic carbonate. Fluids for this seep may utilize high-permeability flow paths either parallel to bedding within the second thrust ridge or along the underlying thrust fault before leaking into the overriding section. Two seaward facing canyons on the third anticlinal ridge have vent clam communities near the canyon mouths at approximately the intersection between the anticlinal ridge and the adjacent forearc basin. No seeps were found along strike at the intersection of the slope basin and anticlinal ridge. We infer that the lack of seepage along strike and the presence of seeps in canyons may be related to fluid flow below the forearc basin/slope unconformity (overpressured by the impinging thrust fault to the west?) directed toward canyons at the surface.     

A study of the chemistry of pore fluids and authigenic carbonates in methane seep environments: Kodiak Trench, Hydrate Ridge, Monterey Bay, and Eel River Basin

Analyses of the chemical and isotopic composition of carbonates rocks recovered from methane seepage areas of the Kodiak Trench, Hydrate Ridge, Monterey Bay Clam Flats, and the Eel River Basin, coupled with the studies of the chemistry of the pore fluids, have shown that these carbonates have grown within the sediment column. Geochemical profiles of pore fluids show that, in deep water seeps (Kodiak Trench—4450 m; Monterey Bay—1000 m; Hydrate Ridge—650 m), δ¹³C (DIC) values are low (isotopically light), whereas in the Eel River area ( 350–500 m), δ¹³C (DIC) values are much higher (isotopically heavier). In all cases, the δ¹³C values indicate that processes of methane oxidation, associated with sulfate reduction, are dominant in the shallow sediments. Data on the isotopic composition of authigenic carbonates found at sites in Kodiak Trench, Eel River Basin South, and Eel River Basin North indicate a variable composition and origin in different geochemical environments. Some of the authigenic carbonates from the study sites show a trend in their δ¹³C values similar to those of the pore fluids obtained in their vicinity, suggesting formation at relatively shallow depths, but others indicate formation at greater sediment depths. The latter usually consist of high magnesium calcite or dolomite, which, from their high values of δ¹³C (up to 23‰;) and δ¹⁸O (up to 7.5‰), suggest formation in the deeper horizons of the sediments, in the zone of methanogenesis. These observations are in agreement with observations by other workers at Hydrate Ridge, in Monterey Bay, and in the Eel River Basin.     

Approximate axial alignment times for spinning bodies

Using a different approximate scheme for each of two separate body shapes, we compute for torque-free rotating quasi-rigid bodies the rate of damping of the nutation angle due to internal energy dissipation. The first method calculates the alignment rate as well as the free precession period of a Maxwell body; it modifies, and thereby makes clearer, a similar calculation by Gerstenkorn (1967b). The second method deals with an elongated body and models it by symmetrically placing masses along a massless rod, where all energy dissipation occurs through the hysteresis loss of the strain energy stored during the precession motion. For both cases, the damping time is proportional to the material strength of the body and to its damping properties, and inversely proportional to the square of the body's radius and the cube of the rotation frequency.The results agree well with those of another approximate solution (Burns and Safronov, 1973) and numerically check that solution for the alignment times of asteroids. They show that essentially all observed asteroids align quite rapidly and thus explain why asteroids are observed in pure spin despite numerous misaligning collisions.     

Field Survey of the March 28, 2005 Nias-Simeulue Earthquake and Tsunami

On the evening of March 28, 2005 at 11:09?p.m. local time (16:09 UTC), a large earthquake occurred offshore of West Sumatra, Indonesia. With a moment magnitude (M _w) of 8.6, the event caused substantial shaking damage and land level changes between Simeulue Island in the north and the Batu Islands in the south. The earthquake also generated a tsunami, which was observed throughout the source region as well as on distant tide gauges. While the tsunami was not as extreme as the tsunami of December 26th, 2004, it did cause significant flooding and damage at some locations. The spatial and temporal proximity of the two events led to a unique set of observational data from the earthquake and tsunami as well as insights relevant to tsunami hazard planning and education efforts.     

Slope basins, headless canyons, and submarine palaeoseismology of the Cascadia accretionary complex

A combination of geomorphological, seismic reflection and geotechnical data constrains this study of sediment erosion and deposition at the toe of the Cascadia accretionary prism. We conducted a series of ALVIN dives in a region south of Astoria Canyon to examine the interrelationship of fluid flow and slope failure in a series of headless submarine canyons. Elevated head gradients at the inflection point of canyons have been inferred to assist in localized failures that feed sediment into a closed slope basin. Measured head gradients are an order of magnitude too low to cause seepage-induced slope failure alone; we therefore propose transient slope failure mechanisms. Intercanyon slopes are uniformly unscarred and smooth, although consolidation tests indicate that up to several metres of material may have been removed. A sheet-like failure would remove sediment uniformly, preserving the observed smooth intercanyon slope. Earthquake-induced liquefaction is a likely trigger for this type of sheet failure as the slope is too steep and short for sediment flow to organize itself into channels. Bathymetric and seismic reflection data suggest sediment in a trench slope basin between the second and third ridges from the prism's deformation is derived locally. A comparison of the amounts of material removed from the slopes and that in the basin shows that the amount of material removed from the slopes may slightly exceed the amount of material in the basin, implying that a small amount of sediment has escaped the basin, perhaps when the second ridge was too low to form a sufficient dam, or through a gap in the second ridge to the south. Regardless, almost 80% of the material shed off the slopes around the basin is deposited locally, whereas the remaining 20% is redeposited on the incoming section and will be re-accreted.
     

On the elastic-perfectly plastic bending of the lithosphere under generalized loading with application to the Kurd Trench

Summary. The bending of the lithosphere as it approaches an oceanic trench is modelled using an elastic-perfectly plastic rheology. The rock within the elastic lithosphere is assumed to behave elastically until a yield stress is reached; for larger strains the stress is assumed to remain at the yield value (no strain hardening). The deflection of the lithosphere is obtained numerically for a range of applied moments, vertical forces and horizontal forces at the trench axis. The results of the analysis are compared with a series of topographic and gravity profiles across the Kuril Trench. In order to model the large observed curvatures of the lithosphere on the outer trench slope plasticity is required. The location of the region of high curvatures gives a value for the compressive axial load. One conclusion is that the lithosphere is under compression seaward of the Kuril Trench.     

Error calibration of geopotential harmonics in recent and past gravitational fields

The geopotential coefficients of a number of recent and past models have been examined with the aim of calibrating their formal (or published) errors, principally by direct comparison with the same coefficients of a precise and wholly independent satellite-only reference field determined (in 2010) predominantly from over 7 years of μ/s low–low GRACE inter-satellite range rate observations (ITG-GRACE2010S). In all of these comparisons the reference field used, over specified spectral ranges, has much smaller reported errors than the ones to be calibrated. In particular we find that a recently published field (in 2010) using gravity gradient and position data from ESA’s GOCE satellite, GO_CONS_GCF_2_TIM_R1, has formal errors which are significantly optimistic for the lowest degrees (n) but with increasing realism where the gradiometer gains influence over the position information (15 < n < 120). Other GOCE models do not afford an unbiased error calibration. Validating error calibration with the independent reference model, of two past comprehensive fields (containing both satellite and surface data), confirms that one of them (JGM3, published in 1994, complete to n = 70) reported generally realistic formal errors while another (EGM96, published in 1998, complete to n = 360, tested only for n < 150) had significantly optimistic ones for most n < 100 but with better realism when affected only by the surface information.