The Yecheng,Xinjiang, China,earthquake of February 14, 1980

In this study, we have checked the location and focal depth of the Yecheng earthquake (m _b = 6.0, maximum intensity VII) of February 14, 1980. The result shows that this is an intermediate event with a focal depth of 90 km. The microepicenter is located at 36. 4°N, 76, 9°E, while the macroepicenter is at 37.3°N, 76.9°E, 90 km to the north of the microepicenter. This is the first destructive intermediate event in China which led to a damage as severe as of intensity VII. The focal mechanism of the event is determined to be of thrust type. Combined with the analysis of seismological and geological data in surrounding area, the possible relation between the event and plate movement has been discussed. The result in this paper indicates that in some particular place, the destructive effect of intermediate event should be considered in seismic hazard assessment.     

Age, geology, geophysics, and geochemistry of Mahukona Volcano, Hawai`i

The size, shape, and magmatic history of the most recently discovered shield volcano in the Hawaiian Islands, Mahukona, have been controversial. Mahukona corresponds to what was thought to be a gap in the paired sequence (Loa and Kea trends) of younger Hawaiian volcanoes (<4?Ma). Here, we present the results of marine expeditions to Mahukona where new bathymetry, sidescan sonar, gravity data, and lava samples were collected to address these controversies. Modeling of bathymetric and gravity data indicate that Mahukona is one of the smallest Hawaiian volcanoes (～6,000?km³) and that its magmatic system was not focused in a long-lived central reservoir like most other Hawaiian volcanoes. This lack of a long-lived magmatic reservoir is reflected by the absence of a central residual gravity high and the random distribution of cones on Mahukona Volcano. Our reconstructed subsidence history for Mahukona suggests it grew to at least ～270?m below sea level but probably did not form an island. New ⁴⁰Ar–³⁹Ar plateau ages range from 350 to 654?ka providing temporal constraints for Mahukona’s post-shield and shield stages of volcanism, which ended prematurely. Mahukona post-shield lavas have high ³He/⁴He ratios (16–21?Ra), which have not been observed in post-shield lavas from other Hawaiian volcanoes. Lava compositions range widely at Mahukona, including Pb isotopic values that straddle the boundary between Kea and Loa sequences of volcanoes. The compositional diversity of Mahukona lavas may be related to its relatively small size (less extensive melting) and the absence of a central magma reservoir where magmas would have been homogenized.     

Facets of uncertainty: epistemic uncertainty,non-stationarity,likelihood, hypothesis testing,and communication

ABSTRACT

This paper presents a discussion of some of the issues associated with the multiple sources of uncertainty and non-stationarity in the analysis and modelling of hydrological systems. Different forms of aleatory, epistemic, semantic, and ontological uncertainty are defined. The potential for epistemic uncertainties to induce disinformation in calibration data and arbitrary non-stationarities in model error characteristics, and surprises in predicting the future, are discussed in the context of other forms of non-stationarity. It is suggested that a condition tree is used to be explicit about the assumptions that underlie any assessment of uncertainty. This also provides an audit trail for providing evidence to decision makers.

Editor D. Koutsoyiannis; Associate editor S. Weijs     

Spatial distribution,source apportionment,and associated risks of trace metals (As,Pb, Cr,Cd, and Hg) from a subtropical river,Gomti, Bangladesh

Sediment samples from the Gomti River basin were investigated to determine and evaluate trace metal concentrations,their biological effect,and potential ecological and human health risks for adults and children.The mean concentrations of trace metals were organized in the descending order of mercury(Hg)(0.08 mg/kg)>chromium(Cr)(0.06 mg/kg)>lead(Pb)(0.05 mg/kg)>arsenic(As)(0.02 mg/kg)>cadmium(Cd)(0.01 mg/kg).The current study illustrated that metals were attributed to the area from natural sources and different anthropogenic sources especially from industries.However,the concentration levels were lower than the sediment quality guidelines(SQGs)based on the effect-range classifications of threshold effect concentration(TEC)and probable effect concentration(PEC).Therefore,the concentration of the elements showed no adverse biological effects on aquatic organisms.The evaluated potential ecological risk index also revealed the low toxicity to the aquatic environment.Moreover,the applied sediment quality indices,geo-accumulation index(Igeo),contamination factor(CF),and contamination severity index(CSI)suggested that the contamination levels of the elements were in the acceptable range,and the contamination had not notably impacted on the sediment quality.The risk index(HI)concerning age groups was significantly less than the threshold limit of 1 indicating that the contamination had no non-carcinogenic risk effect.The total carcinogenic risk(TCR)was less than a risk value of 1×10⁴.Hence,the current study suggests that immediate remediation is not required due to an absence of alarming conditions in the study area.Proper monitoring of the attribution of the metal elements should be done for the betterment of human and environmental health。     

Seismicity at Fuego,Pacaya, Izalco,and San Cristobal Volcanoes,Central America, 1973–1974

Seismic data collected at four volcanoes in Central America during 1973 and 1974 indicate three sources of seismicity: regional earthquakes with hypocentral distances greater than 80 km, earthquakes within 40 km of each volcano, and seismic activity originating at the volcanoes due to eruptive processes. Regional earthquakes generated by the underthrusting and subduction of the Cocos Plate beneath the Caribbean Plate are the most prominent seismic feature in Central America. Earthquakes in the vicinity of the volcanoes occur on faults that appear to be related to volcano formation. Faulting near Fuego and Pacaya volcanoes in Guatemala is more complex due to motion on a major E-W striking transform plate boundary 40 km north of the volcanoes. Volcanic activity produces different kinds of seismic signatures. Shallow tectonic or A-type events originate on nearby faults and occur both singly and in swarms. There are typically from 0 to 6 A-type events per day withb value of about 1.3. At very shallow depths beneath Pacaya, Izalco, and San Cristobal large numbers of low-frequency or B-type events are recorded with predominant frequencies between 2.5 and 4.5 Hz and withb values of 1.7 to 2.9. The relative number of B-type events appears to be related to the eruptive states of the volcanoes; the more active volcanoes have higher levels of seismicity. At Fuego Volcano, however, low-frequency events have unusually long codas and appear to be similar to tremor. High-amplitude volcanic tremor is recorded at Fuego, Pacaya, and San Cristobal during eruptive periods. Large explosion earthquakes at Fuego are well recorded at five stations and yield information on near-surface seismic wave velocities (α=3.0±0.2 km/sec.).     

Podhale,Poland, earthquake of November 30, 2004

Earthquake of November 30, 2004, in Podhale region, southern Poland, was of unexpectedly big size in this area of weak seismicity. As Poland is considered a country of low seismicity, the earthquake has caused concern about seismic hazard in Poland, especially since it took place shortly after the even more unexpected Kaliningrad Region, Russia, earthquakes of September 21, 2004, that inflicted minor damage in northern Poland. The paper presents the findings on the Podhale earthquake which reached macroseismic intensity up to 7 and magnitude 4.7 (m _b ; ISC). The event was felt up to a distance of about 100 km and inflicted slight damage to buildings in its narrow epicentral area, thus evidencing its relatively shallow depth. The quake has been located near the village of Skrzypne, about 15 km west-southwest of the district capital Nowy Targ. The source mechanism has been found to be of dip slip normal fault type, although a problem remains of association of this mechanism with known tectonic dislocations in the region. The earthquake has been followed by a long series of aftershocks. Their distribution in time is also studied and the biggest aftershocks have been located.     

Erosional unloading,hillslope geometry,and the height of the Cascade Range,Washington State,USA

Peaks in the Cascade Range in northern Washington State are on average ～800 m higher than in southern Washington. The influences of differential valley excavation and variations in hillslope length and average slope on these altitudinal trends were tested using a 3‐dimensional model for isostatic rock uplift and calculations of hillslope length and slope respectively. The magnitude of isostatic peak uplift calculated by the model is highly dependent on the flexural rigidity (D) and the related effective elastic thickness (T_e) of the crust of this region. Crustal rigidity was constrained using published estimates and by estimating the depth of the seismogenic zone in the area (D > 1 × 10²³ Nm and T_e > 24 km). With these constraints, isostatic compensation due to differential erosion added < 700 m and 300 m, or < 25% overall, of height to peaks in the northern and southern Cascades, respectively. Deeper valley incision in the northern Cascades accounts for < 300 m of the 800 m difference in peak altitudes between north and south. Similarly, variation in valley spacing and slope account for < 350 m of the difference in mean altitude between northern and southern regions. Hence, at least several hundred m difference in altitude between the northern and southern regions of the Cascades in Washington must be due to tectonic, geologic, or geophysical factors rather than surficial and geomorphic effects like isostatic response to valley incision and hillslope geometry. Copyright © 2009 John Wiley & Sons, Ltd.     

Continuing inflation at Three Sisters volcanic center,central Oregon Cascade Range,USA, from GPS,leveling, and InSAR observations

Uplift of a broad area centered ~6 km west of the summit of South Sister volcano started in September 1997 (onset estimated from model discussed in this paper) and was continuing when surveyed in August 2006. Surface displacements were measured whenever possible since August 1992 with satellite radar interferometry (InSAR), annually since August 2001 with GPS and leveling surveys, and with continuous GPS since May 2001. The average maximum displacement rate from InSAR decreased from 3–5 cm/yr during 1998–2001 to ~1.4 cm/yr during 2004–2006. The other datasets show a similar pattern, i.e., surface uplift and extension rates decreased over time but deformation continued through August 2006. Our best-fit model to the deformation data is a vertical, prolate, spheroidal point-pressure source located 4.9–5.4 km below the surface. The source inflation rate decreased exponentially during 2001–2006 with a 1/e decay time of 5.3 ± 1.1 years. The net increase in source volume from September 1997 to August 2006 was 36.5–41.9 x 10⁶ m³. A swarm of ~300 small (M _max = 1.9) earthquakes occurred beneath the deforming area in March 2004; no other unusual seismicity has been noted. Similar deformation episodes in the past probably would have gone unnoticed if, as we suspect, most are small intrusions that do not culminate in eruptions.     

Deccan plume,lithosphere rifting,and volcanism in Kutch,India

Kutch (northwest India) experienced lithospheric thinning due to rifting and tholeiitic and alkalic volcanism related to the Deccan Traps K/T boundary event. Alkalic lavas, containing mantle xenoliths, form plug-like bodies that are aligned along broadly east–west rift faults. The mantle xenoliths are dominantly spinel wehrlite with fewer spinel lherzolite. Wehrlites are inferred to have formed by reaction between transient carbonatite melts and lherzolite forming the lithosphere. The alkalic lavas are primitive (Mg# = 64–72) relative to the tholeiites (Mg# = 38–54), and are enriched in incompatible trace elements. Isotope and trace element compositions of the tholeiites are similar to what are believed to be the crustally contaminated Deccan tholeiites from elsewhere in India. In terms of Hf, Nd, Sr, and Pb isotope ratios, all except two alkalic basalts plot in a tight cluster that largely overlap the Indian Ridge basalts and only slightly overlap the field of Reunion lavas. This suggests that the alkalic magmas came largely from the asthenosphere mixed with Reunion-like source that welled up beneath the rifted lithosphere. The two alkalic outliers have an affinity toward Group I kimberlites and may have come from an old enriched (metasomatized) asthenosphere. We present a new model for the metasomatism and rifting of the Kutch lithosphere, and magma generation from a CO₂-rich lherzolite mantle. In this model the earliest melts are carbonatite, which locally metasomatized the lithosphere. Further partial melting of CO₂-rich lherzolite at about 2–2.5 GPa from a mixed source of asthenosphere and Reunion-like plume material produced the alkalic melts. Such melts ascended along deep lithospheric rift faults, while devolatilizing and exploding their way up through the lithosphere. Tholeiites may have been generated from the main plume head further south of Kutch.     

High-magnesium primary magmas from Haleakala Volcano, east Maui, Hawaii: petrography, nickel, and major-element constraints

For any given volcanic field the compositions of primary melts provide important constraints on models of magmatic processes and volcanic eruptions. In this paper, based on petrography, olivine and bulk rock compositions, two tholeiitic picrites (samples C122 and C123) from Haleakala Volcano, east Maui are evaluated as possible primary melts. Sample C122 (bulk rock MgO = 16.6%) has a high apparent Mg-Fe exchange coefficient, K_D, between olivine phenocrysts and bulk rock (0.6). However, major-elements and Ni mass-balance calculations show that the olivines in C122 are in equilibrium with the residual melt (matrix) after closed-system equilibrium fractionation of 25 wt.% olivine. Therefore, the Mg/Fe ratio, Ca content, and Ni content of C122 are consistent with the hypothesis that the bulk composition of C122 is close to a primary melt formed by partial melting of a mantle containing olivine with composition around Fo₈₉ to Fo₉₁. The uniform composition and small size (mostly 0.2–0.3 mm) of the olivine, and the glass patches in the matrix suggest fast ascent, and rapid cooling at shallow depth for C122. On the contrary, sample C123, which has an apparent K_D (between the most mafic olivine megacrysts and the bulk rock) close to the equilibrium value (0.27), the multiple planar subgrain boundaries in most of the olivine crystals indicate that it may not be a primary melt unless the deformed olivines are generated at magmatic condition as phenocrysts. If the deformed subgrain boundary texture in olivine could indeed be generated at magmatic condition, then the wide compositional range of olivine crystals in C123 (Fo₇₄ to Fo₉₁) suggests multi-stage crystallization over a wide range of cooling temperatures.The compositions of the two picrites, and a differentiated basalt which does not contain xenocrysts suggest that the Haleakala tholeiites are derived from primary melts with at least 16–17 wt.% MgO. Lavas with such high MgO content are rare in Haleakala and other Hawaiian volcanoes; therefore, most Hawaiian tholeiites must have undergone extensive fractionation histories.     

DInSAR measurements of ground deformation by sinkholes,mining subsidence,and landslides,Ebro River,Spain

Differential Interferometric Synthetic Aperture Radar (DInSAR) has been applied to detect and measure ground deformation in a stretch of the Ebro River valley (Spain) excavated in salt‐bearing evaporites. The capability of the Small Baseline Subset (SBAS) DInSAR technique to detect ground displacement is analyzed comparing the DInSAR results with the available geomorphological information. The deformation map derived from 27 European Remote Sensing (ERS) satellite images covering more than five years provides sub‐centimeter displacement measurements in zones coinciding with known active sinkholes and landslides. Moreover the map provides the first account of mining subsidence in the area. The measured deformation rates reach 1·68 cm/y for the sinkholes, 0·80 cm/y for the landslides and 1·45 cm/y for the area affected by mining subsidence. The SBAS DInSAR technique provided deformation measurements in a small proportion (5–10%) of the known active sinkholes and landslides. This limitation is mainly due to the lack of coherence in agricultural areas, the spatial resolution of the deformation map (pixel size of 90 m), and the parallelism between the ERS satellite line of sight and the linear escarpment on which most of the landslides occur. Despite this, the interferometric technique provides valuable data that complement traditional geomorphological studies including the quantification of the deformation phenomena, the identification of mining subsidence otherwise only recognizable by geodetic methods, and the detection of creep deformation which might correspond to premonitory indicators of catastrophic sinkholes and landslides capable of causing the loss of lives. Detailed DInSAR studies combined with field data would be required to improve the analysis of each deformation area. Copyright © 2009 John Wiley & Sons, Ltd.     

Longitudinal structure of pyroclastic-flow deposits,revealed by GPR survey,at Merapi Volcano,Java, Indonesia

In 2006 Merapi volcano, Indonesia, erupted for a few months, producing several block-and-ash flows reaching a maximum distance of 7.5 km from the main vent. During the eruption, we conducted a survey on those flow deposits in the Gendol Valley at Kaliadem village, about 4.5 km from the Merapi submit, using a Ground Penetrating Radar (GPR). The upper deposit was studied in its distal reaches, whereas the one below was studied in its medial reaches. The field study was carried out with a commercial RAMAC® GPR coupled with 100 MHz antennas, and the data treatment conducted with Reflex™ software. From this survey, we determined both deposits' local (1) thickness – reaching a maximum of 15 m – and (2) internal architecture. This last one is governed by long reflecting horizons extending over 20 to 30 m that delimit layers showing progradation patterns in their distal reaches. Within these layers we could also observe an internal architecture of still unknown origin. The layers are interpreted as the result of the flow pulses that progressively deposited downstream-ward by progradation. However the interpretation of those GPR profiles is a bit hazardous, because of the absence of outcrops, and we can only proceed by analogy with other studies. Nevertheless, despite numerous limitations, GPR is a helpful tool to understand pyroclastic deposits' structure when no visual observations are available.     

Numerical simulations of mantle convection: Time-dependent,three-dimensional,compressible, spherical shell

Abstract

We describe nonlinear time-dependent numerical simulations of whole mantle convection for a Newtonian, infinite Prandtl number, anelastic fluid in a three-dimensional spherical shell for conditions that approximate the Earth's mantle. Each dependent variable is expanded in a series of 4,096 spherical harmonics to resolve its horizontal structure and in 61 Chebyshev polynomials to resolve its radial structure. A semiimplicit time-integration scheme is used with a spectral transform method. In grid space there are 61 unequally-spaced Chebyshev radial levels, 96 Legendre colatitudinal levels, and 192 Fourier longitudinal levels. For this preliminary study we consider four scenarios, all having the same radially-dependent reference state and no internal heating. They differ by their radially-dependent linear viscous and thermal diffusivities and by the specified temperatures on their isothermal, impermeable, stress-free boundaries. We have found that the structure of convection changes dramatically as the Rayleigh number increases from 10⁵ to 10⁶ to 10⁷. The differences also depend on how the Rayleigh number is increased. That is, increasing the superadiabatic temperature drop, δT, across the mantle produces a greater effect than decreasing the diffusivities. The simulation with a Rayleigh number of 10⁷ is approximately 10,000 times critical, close to estimates of that for the Earth's mantle. However, although the velocity structure for this highest Rayleigh number scenario may be adequately resolved, its thermodynamic structure requires greater horizontal resolution. The velocity and thermodynamic structures of the scenarios at Rayleigh numbers of 10⁵ and 10⁶ appear to be adequately resolved. The 10⁵ Rayleigh number solution has a small number of broad regions of warm upflow embedded in a network of narrow cold downflow regions; whereas, the higher Rayleigh number solutions (with large δT) have a large number of small hot upflow plumes embedded in a broad weak background of downflow. In addition, as would be expected, these higher Rayleigh number solutions have thinner thermal boundary layers and larger convective velocities, temperatures perturbations, and heat fluxes. These differences emphasize the importance of developing even more realistic models at realistic Rayleigh numbers if one wishes to investigate by numerical simulation the type of convection that occurs in the Earth's mantle.     

Storage, migration, and eruption of magma at Kilauea volcano, Hawaii, 1971–1972

The magmatic plumbing system of Kilauea Volcano consists of a broad region of magma generation in the upper mantle, a steeply inclined zone through which magma rises to an intravolcano reservoir located about 2 to 6 km beneath the summit of the volcano, and a network of conduits that carry magma from this reservoir to sites of eruption within the caldera and along east and southwest rift zones. The functioning of most parts of this system was illustrated by activity during 1971 and 1972. When a 29-month-long eruption at Mauna Ulu on the east rift zone began to wane in 1971, the summit region of the volcano began to inflate rapidly; apparently, blockage of the feeder conduit to Mauna Ulu diverted a continuing supply of mantle-derived magma to prolonged storage in the summit reservoir. Rapid inflation of the summit area persisted at a nearly constant rate from June 1971 to February 1972, when a conduit to Mauna Ulu was reopened. The cadence of inflation was twice interrupted briefly, first by a 10-hour eruption in Kilauea Caldera on 14 August, and later by an eruption that began in the caldera and migrated 12 km down the southwest rift zone between 24 and 29 September. The 14 August and 24–29 September eruptions added about 10⁷ m³ and 8 × 10⁶ m³, respectively, of new lava to the surface of Kilauea. These volumes, combined with the volume increase represented by inflation of the volcanic edifice itself, account for an approximately 6 × 10⁶ m³/month rate of growth between June 1971 and January 1972, essentially the same rate at which mantle-derived magma was supplied to Kilauea between 1952 and the end of the Mauna Ulu eruption in 1971.The August and September 1971 lavas are tholeiitic basalts of similar major-element chemical composition. The compositions can be reproduced by mixing various proportions of chemically distinct variants of lava that erupted during the preceding activity at Mauna Ulu. Thus, part of the magma rising from the mantle to feed the Mauna Ulu eruption may have been stored within the summit reservoir from 4 to 20 months before it was erupted in the summit caldera and along the southwest rift zone in August and September.The September 1971 activity was only the fourth eruption on the southwest rift zone during Kilauea's 200 years of recorded history, in contrast to more than 20 eruptions on the east rift zone. Order-of-magnitude differences in topographic and geophysical expression indicate greatly disparate eruption rates for far more than historic time and thus suggest a considerably larger dike swarm within the east rift zone than within the southwest rift zone. Characteristics of the historic eruptions on the southwest rift zone suggest that magma may be fed directly from active lava lakes in Kilauea Caldera or from shallow cupolas at the top of the summit magma reservoir, through fissures that propagate down rift from the caldera itself at the onset of eruption. Moreover, emplacement of this magma into the southwest rift zone may be possible only when compressive stress across the rift is reduced by some unknown critical amount owing either to seaward displacement of the terrane south-southeast of the rift zone or to a deflated condition of Mauna Loa Volcano adjacent to the northwest, or both. The former condition arises when the forceful emplacement of dikes into the east rift zone wedges the south flank of Kilauea seaward. Such controls on the potential for eruption along the southwest rift zone may be related to the topographic and geophysical constrasts between the two rift zones.     

Organic and inorganic matter in Louisiana coastal waters: Vermilion, Atchafalaya, Terrebonne, Barataria, and Mississippi regions

Chromophoric dissolved organic matter (CDOM) spectral absorption, dissolved organic carbon (DOC) concentration, and the particulate fraction of inorganic (PIM) and organic matter (POM) were measured in Louisiana coastal waters at Vermilion, Atchafalaya, Terrebonne, Barataria, and Mississippi River locations, in 2007-2008. The range of CDOM was 0.092 m⁻¹ at Barataria in June 2008 to 11.225 m⁻¹ at Mississippi in February 2008. An indicator of organic matter quality was predicted by the spectral slope of absorption coefficients from 350 to 412 nm which was between 0.0087 m⁻¹ at Mississippi in May 2008 and 0.0261 m⁻¹ at Barataria in June 2008. CDOM was the dominant component of light attenuation at Terrebonne and Barataria. Detritus and CDOM were the primary components of light attenuation at Vermilion, Atchafalaya, and Mississippi. DOC ranged between 65 and 1235 μM. PIM ranged between 1.1 and 426.3 mg L⁻¹ and POM was between 0.3 and 49.6 mg L⁻¹.     

Something old,something new,something red,something blue*

Abstract

Hydrologists responsible for flood management need real-time data in order to manage imminent or ongoing floods. In this paper, innovative methods for accessing hydrological data and their spatial visualization are introduced. A multitude of relevant real-time, forecast and historical information is provided in a single, self-updating hydrological map information system. The system consists of a central database and a cartographic user interface and provides harmonized and filtered data in the form of interactive, customizable maps. Maps may also be cross-referenced with historical maps or may be animated for improved comprehension and decision making. Emphasis is placed on the development of the hydrological real-time database that manages large amounts of spatial, temporal and attributive data. The paper focuses on the cartographic user interface, its functionality and the resulting interactive hydrological maps.

Citation Lienert, C., Weingartner, R. &; Hurni, L. (2011) An interactive, web-based, real-time hydrological map information system. Hydrol. Sci. J. 56(1), 1–16     

A tale of two magmas,Fuego, Guatemala

Fuego volcano in Guatemala erupted in 1974 in a basaltic sub-Plinian event, which has been well documented and studied. In 1999, after a period of quiescence lasting 20 years, Fuego erupted again, this time less violently, but with persistent low-level activity. This study investigates the link between these episodes. Previous melt inclusion studies have shown magma erupted in 1974 to have been a volatile-rich hybrid tapped from a vertically extensive system. By contrast, magma erupted in 1999 and 2003 is similar in composition to that erupted in 1974, but melt inclusions are more evolved. Although melt inclusions from the later period are CO₂ rich (up to ∼1,500 ppm), they have low H₂O concentration (max 1.5 wt.%, compared to ∼6 wt.% in 1974). These melt inclusions have a modified H₂O concentration due to diffusive re-equilibration at shallow pressures. Despite this diffusive exchange, both eruptions show evidence of recent mingling of the same low and higher K melts, one of which was slightly cooler than the other and as a result traversed the amphibole stability field. (²¹⁰Pb/²²⁶Ra) data on selected bulk rock samples from 1974 suggest that whereas the cooler, more evolved end-member may have been degassing since the last major eruption in the 1930s, the warmer end-member intruded at most a decade prior to the 1974 eruption. The two end-members are thus batches of the same magma emplaced shallowly ∼30 years apart during which time the older batch was cooled and differentiated before mixing with the younger influx. The presence of the same two melts in the later eruptions suggests that magma in 1999 and 2003 is partly residual from 1974. The current eruptive activity is clearing the system of this residual magma prior to an expected new magma batch.     

Textural, isotopic and REE variations in spinel peridotite xenoliths, Massif Central, France

Sr and Nd isotope analyses and REE patterns are presented for a suite of well-documented mantle-derived xenoliths from the French Massif Central. The xenoliths include spinel harzburgites, spinel lherzolites and some pyroxenites. They show a wide range of textures from undeformed protogranular material through porphyroclastic to equigranular and recrystallised secondary types. Textural differences are strongly linked to trace element geochemistry and variations in radiogenic isotope ratios. Many undeformed protogranular xenoliths are Type IA LREE-depleted with MORB-type ε_Sr values between − 30.7 and − 23.6, and ε_Nd values + 13.9 to + 9.4. A second group of undeformed xenoliths are Type IB LREE-enriched with higher ε_Sr values (− 22.7 to − 10.6) and lower ε_Nd values (+ 11.9 to + 5.6). Deformed xenoliths with porphyroclastic, equigranular and secondary recrystallised textures are all Type IB (LREE-enriched, ε_Nd < 6.4, ε_Sr > 11.8). It is proposed that two separate events have given rise to the observed mixing arrays: (1) MORB-source depleted mantle was enriched by a component derived from an enriched mantle. Deformation and recrystallisation accompanied this event. (2) Subsequently, unenriched MORB-source mantle interacted with magmas chemically akin to the host basalts, and enrichment occurred with little deformation. Hypotheses of Tertiary mantle diapirism resulting in isochemical deformation and refinement of protogranular mantle to equigranular mantle are untenable because of differences in REE patterns and isotopic ratios between different textural groups.     

Deep drainage sensitivity to climate,edaphic factors,and woody encroachment,Oklahoma, USA

Because groundwater is Earth's largest pool of freshwater, understanding the sensitivity of deep drainage to climate, soils, and land cover is critical in managing water resources. To better understand controls on this critical flux in the context of woody encroachment, we determined the sensitivity of deep drainage to climate, soil texture, soil compaction, rooting depth, growing season duration, and plant–water stress response using Hydrus‐1D to simulate deep drainage. To evaluate the simulation results, we compared these results with ground measurements at two anchor sites. At both anchor sites, Hydrus‐1D predictions of deep drainage matched measured values within the errors inherent in ground measurements. Sensitivity analysis suggested greatest sensitivity of deep drainage to climate (24 mm yr^?1) and rooting depth (12 mm yr^?1), moderate sensitivity to growing season duration (5 mm yr^?1) and soil texture (4 mm yr^?1), and lowest sensitivity to topsoil compaction and plant–water stress response (3 mm yr^?1). The sensitivity analysis indicated the relative importance of the plant‐related factors considered, which, in decreasing order, were rooting depth, growing season duration, and plant–water stress response – factors that change concomitantly as a result of forestation or woody encroachment. Further ground‐truth measurements of woody encroachment effects on deep drainage are needed to confirm or refine the results of this simulation modelling study. Copyright © 2015 John Wiley & Sons, Ltd.