Crustal xenoliths from Calbuco Volcano,Andean Southern Volcanic Zone: implications for crustal composition and magma-crust interaction

Crustal xenoliths in the 1961 andesite flow of Calbuco Volcano, in the southern Southern Volcanic Zone (SSVZ) of the Andes, consist predominantly of pyroxene granulites and hornblende gabbronorites. The granulites contain plagioclase+pyroxene+magnetite±amphibole, and have pronounced granoblastic textures. Small amounts of relict amphibole surrounded by pyroxene-plagioclase-magnetite-glass symplectites are found in some specimens. These and similar textures in the gabbronorites are interpreted as evidence of dehydration melting. Mineral and bulk rock geochemical data indicate that the granulites are derived from an incompatible trace element depleted basaltic protolith that underwent two stages of metamorphism: a moderate pressure, high temperature stage accompanied by melting and melt extraction from some samples, followed by thermal metamorphism after entrainment in the Calbuco andesite lavas. High _Nd ^T values (+4.0 to +8.6), Nd-isotope model ages of 1.7–2.0 Ga, and trace element characteristics like chondrite normalized La/Yb< and La/Nb1 indicate that the protoliths were oceanic basalts. Similar oceanic metabasalts of greenschist to amphibolite facies are found in the Paleozoic metamorphic belt that underlies the Chilean coastal ranges. Mineral and bulk rock compositions of the gabbronorite xenoliths indicate that they are cognate, crystallizing from the basaltic andesite magma at Calbuco. Crystallization pressures for the gabbros based on total Al contents in amphibole are 6–8 kbar. These pressures point to middle to lower crustal storage of the Calbuco magma. Neither granulite nor gabbro xenoliths have the appropriate geochemical characteristics to be contaminants of Calbuco andesites, although an ancient sedimentary contaminant is indicated by the lava compositions. The presence of oceanic metabasaltic xenoliths, together with the sedimentary isotopic imprint, suggests that the lower crust beneath the volcano is analogous to the coastal metamorphic belt, which is an accretionary complex of intercalated basalts and sediments that formed along the Paleozoic Gondwanan margin. If this is the case, the geochemical composition of the lower and middle crust beneath the SSVZ is significantly different from that of most recent SSVZ volcanic rocks.     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

A decade of change in the Skidaway River estuary I. Hydrography and nutrients

The Skidaway River estuary is a tidally-dominated subtropical estuary in the southeastern USA surrounded by extensiveSpartina salt marshes. Weekly smapling at high and low tide began in 1986 for hydrography, nutrients, chlorophylla, particulate matter, and microbial and plankton biomass and composition; hydrographic and nutrient data during 1986–1996 are reported here. Salinity varied inversely with river discharge and exhibited variability at all time scales but with no long-term trend. Water temperature typically ranged over 25°C and was without apparent long-term frend. Seasonal cycles in concentrations of NO₃, NH₄, PO₄, Si(OH)₄, and DON were observed, with annual maxima generally occurring in late summer. Superimposed on seasonal cycles, all five nutrients exhibited steady increases in minimum, mean, and maximum concentrations; mean concentrations increased c. 50–150% during the decade. Nutrient concentrations were highly correlated with water temperature over the ten-year period, but weakly related to salinity and discharge. Nutrients were strongly correlated with one another, and the relative ratios among inorganic nutrients showed little long-term trend. Correlations among temperature and nutrient concentrations exhibited considerable inter-annual variability. Major spikes in organic and inorganic nutrient concentrations coincided with significant rainfall events; concentrations increased hyperbolically with rainfall. Although pristine compared to more heavily impacted waterways primarily outside the region, residential development and population density have been increasing rapidly during the past 15–20 years. Land use is apparently altering nutrient loading over the long-term (months-years), and superimposed on this are stochastic meteorological events that accelerate these changes over the short term (days-weeks).     

Kosa (yellow sand) components in precipitation collected at central Japan

In east Asia, acidic gases derived from fossil fuel combustion have increased in the past decades. On the other hand, the Asian dust, also called Kosa (yellow sand) is transported following windstorms from arid lands in the Asian continent. Many researchers have been interested in the reaction between acidic aerosols and Kosa aerosols as well as the long-range transport of these emissions. To investigate the characteristics of chemical components in precipitation on a long-term basis over Japan, precipitation was sequentially collected from April 1984 to March 1997 at Kanazawa located near the coast of the Sea of Japan. Precipitation samples were collected at 1 mm intervals for the first 5 mm rainfall and all volume of rainwater after 6 mm for all precipitation events with an automatic wet only precipitation collector. According to the analyses of precipitation including Kosa aerosols during Kosa periods, the reaction in the air between Kosa and acidic components during the long-range transport was discussed.     

Getting ‘holistic’ about land rehabilitation in hot arid regions

The paper describes Western Australian examples and causes of land degradation. It outlines shortcomings in the methodologies used to rehabilitate these areas. From this a protocol is suggested for an ‘holistic’ approach to land rehabilitation.     

Screening model for volatile pollutants in dual porosity soils

This paper develops mass fraction models for transport and fate of agricultural pollutants in structured two-region soils. Mass fraction index models, based on a semi-infinite domain solution, are derived that describe leaching at depth, vapor losses through soil surface, absorption, and degradation in the dynamic- and stagnant-water soil regions. The models predict that leaching is the result of the combined effect of the upward vapor-phase transport relative to downward advection, residence time relative to half-life, dispersion, and lateral diffusive mass transfer. Simulations show that leached fraction of volatile compounds does not always decrease monotonically with increased residence time relative to the pollutant half-life, as a result of complex interactions among the different physical and biochemical processes. The results show that leaching, volatilization, and degradation losses can be affected significantly by lateral diffusive mass transfer into immobile-water regions and advection relative to dispersion (i.e. Peclet number) in the mobile-water regions. It is shown that solute diffusion into the immobile phase and subsequent biochemical decay reduces leaching and vapor losses through soil surface. Potential use of the modified leaching index for the screening of selected pesticides is illustrated for different soil textures and infiltration rates. The analysis may be useful to the management of pesticides and the design of landfills.     

Fine structure of Pn velocity beneath Sichuan-Yunnan region

We use 23298 Pn arrival-time data from Chinese national and provincial earthquake bulletins to invert fine structure of Pn velocity and anisotropy at the top of the mantle beneath the Sichuan-Yunnan and its adjacent region. The results suggest that the Pn velocity in this region shows significant lateral variation; the Pn velocity varies from 7.7 to 8.3 km/s. The Pn-velocity variation correlates well with the tectonic activity and heat flow of the region. Low Pn velocity is observed in southwest Yunnan, Tengchong volcano area, and the Panxi tectonic area. These areas have very active seismicity and tectonic activity with high surface heat flow. On the other hand, high Pn velocity is observed in some stable regions, such as the central region of the Yangtze Platform; the most pronounced high velocity area is located in the Sichuan Basin, south of Chengdu. Pn anisotropy shows a complex pattern of regional deformation. The Pn fast direction shows a prominent clockwise rotation pattern from east of the Tibetan block to the Sichuan-Yunnan diamond block to southwest Yunnan, which may be related to southeastward escape of the Tibetan Plateau material due to the collision of the Indian Plate to the Eurasia Plate. Thus there appears to be strong correlation between the crustal deformation and the upper mantle structure in the region. The delay times of events and stations show that the crust thickness decreases from the Tibetan Plateau to eastern China, which is consistent with the results from deep seismic sounding.