Scattering of surface waves modelled by the integral equation method

The integral equation method is used to model the propagation of surface waves in 3-D structures. The wavefield is represented by the Fredholm integral equation, and the scattered surface waves are calculated by solving the integral equation numerically. The integration of the Green's function elements is given analytically by treating the singularity of the Hankel function at   R = 0  , based on the proper expression of the Green's function and the addition theorem of the Hankel function. No far-field and Born approximation is made. We investigate the scattering of surface waves propagating in layered reference models imbedding a heterogeneity with different density, as well as Lamé constant contrasts, both in frequency and time domains, for incident plane waves and point sources.     

Transferring radiometric dating of the last interglacial sea level high stand to marine and ice core records

In order to derive a radiometric age marker for the end of the penultimate glacial–interglacial transition, we compiled published U-series isotope measurements on corals from the period extending from stage 6 to the middle of the last interglacial, and computed the corresponding open-system ages using Thompson et al. model (Thompson, W.G., Spiegelman, M.W., Goldstein, S.L., Speed, R.C., An open-system model for U-series age determinations of fossil corals. Earth Planet. Sci. Lett. 210 (2003) 365–381). We obtain a global mean age of 126 calendar kyr BP (ka) ± 1.7kyr (2σ) for the beginning of the last interglacial sea level high stand. After showing that the phase relationships observed between changes in sea level, North Atlantic benthic and planktonic foraminifera oxygen isotopic records, and atmospheric methane over the last deglaciation were likely also valid over the penultimate deglaciation, we derive an age of 131.2ka ± 2kyr (2σ) for the abrupt increase in atmospheric CH₄ and North Atlantic surface temperature marking the end of the penultimate glacial–interglacial transition. This age is consistent with U–Th dates of the penultimate glacial–interglacial transition recorded in speleothems from sites where speleothems isotopic records are synchronous with North Atlantic temperature records over the last deglaciation. Finally, we show that the phase obtained between the climatic response and northern hemisphere summer insolation is not constant from Termination II to Termination I, implying that northern hemisphere summer insolation alone cannot explain the timing of terminations.     

Debris-covered Glaciers and Rock Glaciers in the Nanga Parbat Himalaya, Pakistan

The origin and mobilization of the extensive debris cover associated with the glaciers of the Nanga Parbat Himalaya is complex. In this paper we propose a mechanism by which glaciers can form rock glaciers through inefficiency of sediment transfer from glacier ice to meltwater. Inefficient transfer is caused by various processes that promote plentiful sediment supply and decrease sediment transfer potential. Most debris‐covered glaciers on Nanga Parbat with higher velocities of movement and/ or efficient debris transfer mechanisms do not form rock glaciers, perhaps because debris is mobilized quickly and removed from such glacier systems. Those whose ice movement activity is lower and those where inefficient sediment transfer mechanisms allow plentiful debris to accumulate, can form classic rock glaciers. We document here with maps, satellite images, and field observations the probable evolution of part of a slow and inefficient ice glacier into a rock glacier at the margins of Sachen Glacier in c. 50 years, as well as several other examples that formed in a longer period of time. Sachen Glacier receives all of its nourishment from ice and snow avalanches from surrounding areas of high relief, but has low ice velocities and no efficient system of debris removal. Consequently it has a pronounced digitate terminus with four lobes that have moved outward from the lateral moraines as rock glaciers with prounced transverse ridges and furrows and steep fronts at the angle of repose. Raikot Glacier has a velocity five times higher than Sachen Glacier and a thick cover of rock debris at its terminus that is efficienctly removed. During the advance stage of the glacier since 1994, ice cliffs were exposed at the terminus, and an outbreak flood swept away much debris from its margins and terminus. Like the Sachen Glacier that it resembles, Shaigiri Glacier receives all its nourishment from ice and snow avalanches and has an extensive debris cover with steep margins close to the angle of repose. It has a high velocity similar to Raikot Glacier and catastrophic breakout floods have removed debris from its terminus twice in the recent past. In addition, the Shaigiri terminus blocked the Rupal River during the Little Ice Age and is presently being undercut and steepened by the river. With higher velocities and more efficient sediment transfer systems, neither the Raikot nor the Shaigiri form classic rock‐glacier morphologies.     

Complex dynamics of small-moderate volcanic events: the example of the 2011 rhyolitic Cordón Caulle eruption,Chile

After decades of repose, Puyehue-Cordón Caulle Volcano (Chile) erupted in June 2011 following a month of continuously increasing seismic activity. The eruption dispersed a large volume of rhyolitic tephra over a wide area and was characterized by complex dynamics. During the initial climactic phase of the eruption (24–30 h on 4–5 June), 11–14-km-high plumes dispersed most of the erupted tephra eastward towards Argentina, reaching as far as the Atlantic Ocean. This first eruptive phase was followed by activity of lower intensity, leading to the development of a complex stratigraphic sequence, mainly due to rapid shifts in wind direction and eruptive style. The resulting tephra deposits consist of 13 main layers grouped into four units. Each layer was characterized based on its dispersal direction, sedimentological features, and on the main characteristics of the juvenile fraction (texture, density, petrography, chemistry). The lowest part of the eruptive sequence (Unit I), corresponding to the tephra emitted between 4 and 5 June, is composed of alternating lapilli layers with a total estimated volume of ca. 0.75 km³; these layers record the highest intensity phase, during which a bent-over plume dispersed tephra towards the southeast-east, with negligible up-wind sedimentation. Products emitted during 5–6 June (Unit II) signaled an abrupt shift in wind direction towards the north, leading to the deposition of a coarse ash deposit in the northern sector (ca. 0.21 km³ in volume), followed by a resumption of easterly directed winds. A third phase (Unit III) began on 7 June and resulted in tephra deposits in the eastern sector and ballistic bombs around the vent area. A final phase (Unit IV) started after 15 June and was characterized by the emission of fine-grained white tephra from ash-charged plumes during low-level activity and the extrusion of a viscous lava flow. Timing and duration of the first eruptive phases were constrained based on comparison of the dispersal of the main tephra layers with satellite images, showing that most of the tephra was emitted during the first 72 h of the event. The analyzed juvenile material tightly clusters within the rhyolitic field, with negligible chemical variations through the eruptive sequence. Textural observations reveal that changes in eruption intensity (and consequently in magma ascent velocity within the conduit) and complex interactions between gas-rich and gas-depleted magma portions during ascent resulted in vesicular clasts with variable degrees of shear localization, and possibly in the large heterogeneity of the juvenile material.     

Interseismic Strain Localization in the San Jacinto Fault Zone

We investigate interseismic deformation across the San Jacinto fault at Anza, California where previous geodetic observations have indicated an anomalously high shear strain rate. We present an updated set of secular velocities from GPS and InSAR observations that reveal a 2–3 km wide shear zone deforming at a rate that exceeds the background strain rate by more than a factor of two. GPS occupations of an alignment array installed in 1990 across the fault trace at Anza allow us to rule out shallow creep as a possible contributor to the observed strain rate. Using a dislocation model in a heterogeneous elastic half space, we show that a reduction in shear modulus within the fault zone by a factor of 1.2–1.6 as imaged tomographically by Allam and Ben-Zion (Geophys J Int 190:1181–1196, 2012) can explain about 50 % of the observed anomalous strain rate. However, the best-fitting locking depth in this case (10.4 ± 1.3 km) is significantly less than the local depth extent of seismicity (14–18 km). We show that a deep fault zone with a shear modulus reduction of at least a factor of 2.4 would be required to explain fully the geodetic strain rate, assuming the locking depth is 15 km. Two alternative possibilities include fault creep at a substantial fraction of the long-term slip rate within the region of deep microseismicity, or a reduced yield strength within the upper fault zone leading to distributed plastic failure during the interseismic period.     

Unprecedented recent warming rate and temperature variability over the east Tibetan Plateau inferred from Alpine treeline dendrochronology

Climate Dynamics - Despite instrumental records showing recent large temperature rises on the Tibetan Plateau (TP), only a few tree-ring temperature reconstructions do capture this warming trend....     

Population Structure of Adult Blue Crabs, <Emphasis Type="Italic">Callinectes sapidus</Emphasis>, in Relation to Physical Characteristics in Barnegat Bay,New Jersey

Blue crabs (Callinectes sapidus) are an important species in coastal or lagoonal estuaries where adult population characteristics may differ as compared to drowned-river estuaries. Barnegat Bay, in southern New Jersey, is composed of two large embayments: one without and one with a salinity gradient. We tested the influence of physical characteristics on the abundance, sex ratio, and size of adult blue crabs and examined variation in measures of reproductive potential (e.g., sperm stores) in both sexes in Barnegat Bay from June to September, 2008–2009. Population structure was distinct between the embayments due to sex-specific responses to salinity: male abundance was negatively correlated with salinity whereas adult females were more abundant in high salinity because of proximity to Barnegat Inlet. This produced high sex ratios in low salinity areas and low sex ratios in high salinity areas. Summer was a growing season for adult males while in late summer-early fall, juvenile males recruited to the adult size class. The spawning season lasted from May to August and ovigerous females were concentrated near the inlets. Information on female sperm stores and ovarian development identified two cohorts of adult females: females that will spawn in the current summer and females that will not spawn until the following summer. Thus, not all adult females near the spawning grounds were members of the current spawning stock. This suggests that annual estimates of spawning stock size which overlook the proximity of females to spawning are overestimating the current spawning stock in Barnegat Bay and other estuaries.     

Past and future polar amplification of climate change: climate model intercomparisons and ice-core constraints

Climate model simulations available from the PMIP1, PMIP2 and CMIP (IPCC-AR4) intercomparison projects for past and future climate change simulations are examined in terms of polar temperature changes in comparison to global temperature changes and with respect to pre-industrial reference simulations. For the mid-Holocene (MH, 6,000 years ago), the models are forced by changes in the Earth’s orbital parameters. The MH PMIP1 atmosphere-only simulations conducted with sea surface temperatures fixed to modern conditions show no MH consistent response for the poles, whereas the new PMIP2 coupled atmosphere–ocean climate models systematically simulate a significant MH warming both for Greenland (but smaller than ice-core based estimates) and Antarctica (consistent with the range of ice-core based range). In both PMIP1 and PMIP2, the MH annual mean changes in global temperature are negligible, consistent with the MH orbital forcing. The simulated last glacial maximum (LGM, 21,000 years ago) to pre-industrial change in global mean temperature ranges between 3 and 7°C in PMIP1 and PMIP2 model runs, similar to the range of temperature change expected from a quadrupling of atmospheric CO₂ concentrations in the CMIP simulations. Both LGM and future climate simulations are associated with a polar amplification of climate change. The range of glacial polar amplification in Greenland is strongly dependent on the ice sheet elevation changes prescribed to the climate models. All PMIP2 simulations systematically underestimate the reconstructed glacial–interglacial Greenland temperature change, while some of the simulations do capture the reconstructed glacial–interglacial Antarctic temperature change. Uncertainties in the prescribed central ice cap elevation cannot account for the temperature change underestimation by climate models. The variety of climate model sensitivities enables the exploration of the relative changes in polar temperature with respect to changes in global temperatures. Simulated changes of polar temperatures are strongly related to changes in simulated global temperatures for both future and LGM climates, confirming that ice-core-based reconstructions provide quantitative insights on global climate changes. An erratum to this article can be found at     

A Protocol for the Determination of the Rare Earth Elements at Picomole Level in Rocks by ICP-MS: Results on Geological Reference Materials USGS PCC-1 and DTS-1

Rare earth element analyses are widely used in geology, environmental science and archaeology. Over the past decade inductively coupled plasma-mass spectrometry has become an important source of rare earth data on geological material. However, ICP-MS analysis of rock samples without pre-concentration can be problematic because of complex sample matrices that can generate significant molecular isobaric interferences on rare earth peaks and which need to be corrected. Such problems are exacerbated for ultramafic rocks because the low levels of rare earth elements demand more concentrated solutions in order to maintain signals above background levels. These high solid loads result in intra-run changes in instrument sensitivity which need to be monitored. Pre-concentration chemistries have been developed in order to avoid high solid loads but these are time-consuming and must offer quantitative recoveries or use a yield tracer. Here, we describe an alternative method for rare earth element analysis by ICP-MS, which involves no pre-concentration and is, therefore, able to deliver data rapidly. Our approach is to apply an external correction procedure, based on the analysis of a reference material closely matched in composition to the unknown samples, which allows correction for both interferences and variations in instrument sensitivity. Testing this method, we obtained accurate rare earth element results for basaltic rocks with a precision of about 2% (1s). We demonstrate that the method is also applicable to ultramafic rocks with abundances at ultra-trace (ng g⁻¹) level and present data for twelve separate dissolutions of the peridotite USGS PCC-1 and four separate dissolutions of the dunite DTS-1 reference materials. The repeatability of the data is between 3% and 9% (1s).