A Wavelet Analysis of Pliopleistocene Climate Indicators: A New View of Periodicity Evolution

Wavelet analysis offers an alternative to Fourier based time-series analysis, and is particularly useful when the amplitudes and periods of dominant cycles are time dependent. We analyse climatic records derived from oxygen isotopic ratios of marine sediment cores with modified Morlet wavelets. We use a normalization of the Morlet wavelets which allows direct correspondence with Fourier analysis. This provides a direct view of the oscillations at various frequencies, and illustrates the nature of the time-dependence of the dominant cycles.

Magma Extraction from the Mantle Wedge at Convergent Margins Through Dikes: A Parametric Sensitivity Analysis

We address under what conditions a magma generated by partial melting at 100 km depth in the mantle wedge above a subduction zone can reach the crust in dikes before stalling. We also address under what conditions primitive basaltic magma (Mg # >60) can be delivered from this depth to the crust. We employ linear elastic fracture mechanics with magma solidification theory and perform a parametric sensitivity analysis. All dikes are initiated at a depth of 100 km in the thermal core of the wedge, and the Moho is fixed at 35 km depth. We consider a range of melt solidus temperatures (800-1100 degrees C), viscosities (10-100 Pa s), and densities (2400-2700 kg m(-3)). We also consider a range of host rock fracture toughness values (50-300 MPa m(1/2)) and dike lengths (2-5 km) and two thermal structures for the mantle wedge (1260 and 1400 degrees C at 100 km depth and 760 and 900 degrees C at 35 km depth). For the given parameter space, many dikes can reach the Moho in less than a few hundred hours, well within the time constraints provided by U series isotope disequilibria studies. Increasing the temperature in the mantle wedge, or increasing the dike length, allows additional dikes to propagate to the Moho. We conclude that some dikes with vertical lengths near their critical lengths and relatively high solidus temperatures will stall in the mantle before reaching the Moho, and these may be returned by corner flow to depths where they can melt under hydrous conditions. Thus, a chemical signature in arc lavas suggesting partial melting of slab basalts may be partly influenced by these recycled dikes. Alternatively, dikes with lengths well above their critical lengths can easily deliver primitive magmas to the crust, particularly if the mantle wedge is relatively hot. Dike transport remains a viable primary mechanism of magma ascent in convergent tectonic settings, but the potential for less rapid mechanisms making an important contribution increases as the mantle temperature at the Moho approaches the solidus temperature of the magma.

On the Pulling Power of Ice Streams

Gravity wants to pull an ice sheet to the center of the Earth, but cannot because the Earth's crust is in the way, so ice is pushed out sideways instead. Or is it? The ice sheet "sees" nothing preventing it from spreading out except air, which is much less massive than ice. Therefore, does not ice rush forward to fill this relative vacuum; does not the relative vacuum suck ice into it, because Nature abhors a vacuum? If so, the ice sheet is not only pulled downward by gravity, it is also pulled outward by the relative vacuum. This pulling outward will be most rapid where the ice sheet encounters least resistance. The least resistance exists along the bed of ice streams, where ice-bed coupling is reduced by a basal water layer, especially if the ice stream becomes afloat and the floating part is relatively unconfined around its perimeter and unpinned to the sea floor. Ice streams are therefore fast currents of ice that develop near the margins of an ice sheet where these conditions exist. Because of these conditions, ice streams pull ice out of ice sheets and have pulling power equal to the longitudinal gravitational pulling force multiplied by the ice-stream velocity. These boundary conditions beneath and beyond ice streams can be quantified by a basal buoyancy factor that provides a life-cycle classification of ice streams into inception, growth, mature, declining and terminal stages, during which ice streams disintegrate the ice sheet. Surface profiles of ice streams are diagnostic of the stage in a life cycle and, hence, of the vitality of the ice sheet.

Stratigraphic Variation Within Polar Firn Caused by Differential Accumulation and Ice Flow: Interpretation of a 400 Mhz Short-Pulse Radar Profile from West Antarctica

We investigate causes of the stratigraphic variation revealed in a 177 km, 400 MHz short-pulse radar profile of firn from West Antarctica. The profile covers 56 m depth, and its direction was close to those of the ice flow and mean wind. The average, near-surface accumulation rates calculated from the time delays of one radar horizon consistently show minima on leeward slopes and maxima on windward slopes, confirming an earlier study based on stake observations. The stratigraphic variation includes up to 30 m depth variation in individual horizons over tens of km, fold limbs that become progressively steeper with depth, and fold-hinge loci that change direction or propagate down-ice with depth over distances far less than predicted by the ice speeds. We use an accumulation rate model to show how local rate anomalies and the effect of ice speed upon a periodic variation in accumulation rate cause these phenomena, and we reproduce two key features seen in the stratigraphic variations. We conclude that the model provides an explanation of changes in spatial stratigraphy and local measures of accumulation history given the constraints of surface topography, ice and wind velocities, and a general accumulation rate for an area.

Stable-Isotopic Composition of Precipitation Over the Northern Slope of the Central Himalaya

Stable-water-isotope data (deltaD and delta(18)O) from three groups of samples (fresh-snow and snow-pit samples collected on Qomolangma (Mount Everest) and Xixabangma during field seasons 1997,1998 and 2001, and precipitation samples collected at Tingri station during summer 2000) are presented and used to survey the isotopic composition of precipitation over the northern slope of the central Himalaya. Multi-year snow-pit samples on Qomolangma have a local meteoric water-line (slope = 8) close to the global value. Deuterium excess (d = deltaD - 8delta(18)O) values at Tingri are much lower than those in fresh snow from Qomolangma, probably due to differences in moisture source and air-mass trajectories as well as local weather conditions. There is no obvious seasonal trend for d values in the Qomolangma region. A negative relationship exists between delta(18)O and d values in both fresh snow on Qomolangma and precipitation at Tingri. Fresh-snow samples collected from different altitudes on Xixabangma allow us to investigate the altitude effect on delta(18)O values in snow. Of four storm events, only one has an obvious altitude effect on delta(18)O variation and a very low gradient of -0.1% per 100 in elevation.

Stable-Isotope and Trace Element Time Series from Fedchenko Glacier (Pamirs) Snow/Firn Cores

In summer 2005, two pilot snow/firn cores were obtained at 5365 and 5206 m a.s.l. on Fedchenko glacier, Pamirs, Tajikistan, the world's longest and deepest alpine glacier. The well-defined seasonal layering appearing in stable-isotope and trace element distribution identified the physical links controlling the climate and aerosol concentration signals. Air temperature and humidity/precipitation were the primary determinants of stable-isotope ratios. Most precipitation over the Pamirs originated in the Atlantic. In summer, water vapor was re-evaporated from semi-arid regions in central Eurasia. The semi-arid regions contribute to non-soluble aerosol loading in snow accumulated on Fedchenko glacier. In the Pamir core, concentrations of rare earth elements, major and other elements were less than those in the Tien Shan but greater than those in Antarctica, Greenland, the Alps and the Altai. The content of heavy metals in the Fedchenko cores is 2-14 times lower than in the Altai glaciers. Loess from Afghan-Tajik deposits is the predominant lithogenic material transported to the Pamirs. Trace elements generally showed that aerosol concentration tended to increase on the windward slopes during dust storms but tended to decrease with altitude under clear conditions. The trace element profile documented one of the most severe droughts in the 20th century.

Subcritical Propagation and Coalescence of Oil-Filled Cracks: Getting the Oil Out of Low-Permeability Source Rocks

We use a fracture mechanics model to study subcritical propagation and coalescence of single and collinear oil-filled cracks during conversion of kerogen to oil. The subcritical propagation distance, propagation duration, crack coalescence and excess oil pressure in the crack are determined using the fracture mechanics model together with the kinetics of kerogen-oil transformation. The propagation duration for the single crack is governed by the transformation kinetics whereas the propagation duration for the multiple collinear cracks may vary by two orders of magnitude depending on initial crack spacing. A large amount of kerogen (>90%) remains unconverted when the collinear cracks coalesce and the new, larger cracks resulting from coalescence will continue to propagate with continued kerogen-oil conversion. The excess oil pressure on the crack surfaces drops precipitously when the collinear cracks are about to coalesce, and crack propagation duration and oil pressure on the crack surfaces are strongly dependent on temperature. Citation: Jin, Z.-H., S. E. Johnson, and Z. Q. Fan (2010), Subcritical propagation and coalescence of oil-filled cracks: Getting the oil out of low-permeability source rocks, Geophys. Res. Lett., 37, L01305, doi:10.1029/2009GL041576.

Dating the Siple Dome (Antarctica) Ice Core By Manual and Computer Interpretation of Annual Layering

The Holocene portion of the Siple Dome (Antarctica) ice core was dated by interpreting the electrical, visual and chemical properties of the core. The data were interpreted manually and with a computer algorithm. The algorithm interpretation was adjusted to be consistent with atmospheric methane stratigraphic ties to the GISP2 (Greenland Ice Sheet Project 2) ice core, (BE)-B-10 stratigraphic ties to the dendrochronology C-14 record and the dated volcanic stratigraphy. The algorithm interpretation is more consistent and better quantified than the tedious and subjective manual interpretation.

Derived Bedrock Elevations, Strain Rates and Stresses from Measured Surface Elevations and Velocities - Jakobshavns-Isbrae, Greenland

Jakobshavns Isbrae (69 degrees 10'N, 49 degrees 5'W) drains about 6.5% of the Greenland ice sheet and is the fastest ice stream known. The Jakobshavns Isbrae basin of about 10 000 km(2) was mapped photogrammetrically from four sets of aerial photography, two taken in July 1985 and two in July 1986. Positions and elevations of several hundred natural features on the ice surface were determined for each epoch by photogrammetric block-aerial triangulation, and surface velocity vectors were computed from the positions. The two flights in 1985 yielded the best results and provided most common points (716) for velocity determinations and are therefore used in the modeling studies. The data from these irregularly spaced points were used to calculate ice elevations and velocity vectors at uniformly spaced grid paints 3 km apart by interpolation. The field of surface strain rates was then calculated from these gridded data and used to compute the field of surface deviatoric stresses, using the flow law of ice, for rectilinear coordinates, X, Y pointing eastward and northward. and curvilinear coordinates, L, T pointing longitudinally and transversely to the changing ice-flow direction. Ice-surface elevations and slopes were then used to calculate ice thicknesses and the fraction of the ice velocity due to basal sliding. Our calculated ice thicknesses are in fair agreement with an ice-thickness map based on seismic sounding and supplied to us by K. Echelmeyer. Ice thicknesses were subtracted from measured ice-surface elevations to map bed topography. Our calculation shows that basal sliding is significant only in the 10-15 km before Jakobshavns Isbrae becomes afloat in Jakobshavns IsfJord.

CO(2) Diffusion in Polar Ice: Observations from Naturally Formed CO(2) Spikes in the Siple Dome (Antarctica) Ice Core

One common assumption in interpreting ice-core CO(2) records is that diffusion in the ice does not affect the concentration profile. However, this assumption remains untested because the extremely small CO(2) diffusion coefficient in ice has not been accurately determined in the laboratory. In this study we take advantage of high levels of CO(2) associated with refrozen layers in an ice core from Siple Dome, Antarctica, to study CO(2) diffusion rates. We use noble gases (Xe/Ar and Kr/Ar), electrical conductivity and Ca(2+) ion concentrations to show that substantial CO(2) diffusion may occur in ice on timescales of thousands of years. We estimate the permeation coefficient for CO(2) in ice is similar to 4 x 10(-21) mol m(-1) s(-1) Pa(-1) at -23 degrees C in the top 287 m (corresponding to 2.74 kyr). Smoothing of the CO(2) record by diffusion at this depth/age is one or two orders of magnitude smaller than the smoothing in the firn. However, simulations for depths of similar to 930-950m (similar to 60-70 kyr) indicate that smoothing of the CO(2) record by diffusion in deep ice is comparable to smoothing in the firn. Other types of diffusion (e.g. via liquid in ice grain boundaries or veins) may also be important but their influence has not been quantified.