Radiative forcing and climate impact resulting from SO2 injections based on a 200,000-year record of Plinian eruptions along the Central American Volcanic Arc

We present for the first time a self-consistent methodology connecting volcanological field data to global climate model estimates for a regional time series of explosive volcanic events. Using the petrologic method, we estimated SO₂ emissions from 36 detected Plinian volcanic eruptions occurring at the Central American Volcanic Arc (CAVA) during the past 200,000 years. Together with simple parametrized relationships collected from past studies, we derive estimates of global maximum volcanic aerosol optical depth (AOD) and radiative forcing (RF) describing the effect of each eruption on radiation reaching the Earth’s surface. In parallel, AOD and RF time series for selected CAVA eruptions are simulated with the global aerosol model MAECHAM5-HAM, which shows a relationship between stratospheric SO₂ injection and maximum global mean AOD that is linear for smaller volcanic eruptions (<5 Mt SO₂) and nonlinear for larger ones (≥5 Mt SO₂) and is qualitatively and quantitatively consistent with the relationship used in the simple parametrized approximation. Potential climate impacts of the selected CAVA eruptions are estimated using an earth system model of intermediate complexity by RF time series derived by (1) directly from the global aerosol model and (2) from the simple parametrized approximation assuming a 12-month exponential decay of global AOD. We find that while the maximum AOD and RF values are consistent between the two methods, their temporal evolutions are significantly different. As a result, simulated global maximum temperature anomalies and the duration of the temperature response depend on which RF time series is used, varying between 2 and 3 K and 60 and 90 years for the largest eruption of the CAVA dataset. Comparing the recurrence time of eruptions, based on the CAVA dataset, with the duration of climate impacts, based on the model results, we conclude that cumulative impacts due to successive eruptions are unlikely. The methodology and results presented here can be used to calculate approximate volcanic forcings and potential climate impacts from sulfur emissions, sulfate aerosol or AOD data for any eruption that injects sulfur into the tropical stratosphere.     

Anthropogenic influence of trace metals in sediments of the Al-Qilt catchment,West Bank,Palestine: 1. Contamination factor and bonding forms

Al-Qilt catchment, located east of both the Ramallah and Jerusalem districts is strongly influenced by waste waters discharged from Israeli settlements and Palestinian urban centers. In this research, our aim was to verify the chemical status of the surface sediments from the Al-Qilt catchment during 2008/2009. For this purpose, 36 surface sediment samples were collected and analyzed with ICP/MS within the fraction of below 63 μm. The extent of contamination of trace metals from Al-Qilt sediments sites were measured by evaluating the contamination factors (CF) from digestion by aqua regia and additionally by sequential extraction steps using the BCR-method. On the basis of the calculated CF, sediments appear particularly contaminated with Zn, Cu, Ag, Sn, Cd, Hg, Bi, and B. Three regions, Wadi Sweanit, Ras Al-Qilt, and Qalandiah have been identified as critical points of contamination. Trace metal inputs to the Al-Qilt catchment need to be kept under strict control in the future since Ras Al-Qilt is considered as one of the important springs in the area and is used for domestic purposes. Moreover, it will be susceptible to pollution if no action is taken to decrease the pollution at the upstream of Al-Qilt.     

Spatial heterogeneous response of land use and landscape functions to ecological restoration: the case of the Chinese loess hilly region

Ecological restorations over time may have profound effects on ecological and socio-economic systems. However, land-use changes and landscape functions that accompany ecological restorations can have spatial differentiations due to varied biophysical and socio-economic contexts. Therefore, these spatial differentiations caused by ecological restoration must be understood for better planning and management of restoration activities. The Baota District, with 576 villages in the center of the Chinese Loess Plateau, was selected as the study area because of its dramatic transition from cropland to grassland and shrubland from 1990 to 2010. Using the ArcGIS software and a k-means clustering analysis, an approach to identify types of land-use change patterns (TLCPs) at the village level was developed, and four TLCPs were delineated. The analysis indicated a general pattern of cropland decline by 21.6 %, but revealed significant spatial variations between villages in different TLCPs. Vegetation cover and soil retention, which are key proxies for landscape functions, increased by 22.70 and 108 %, respectively, from 2000 to 2010 with significant spatial heterogeneity. The Universal Soil Loss Equation was employed for the assessment of soil retention. The analysis of landscape metrics revealed a major trend of fragmentation and regularity on the county and village scale; however, spatial variations remained. Physical attributes were used to characterize different TLCPs, and notable differences were found. The spatial heterogeneous change in land use and landscape functions on the village scale may be useful for land use and ecological restoration management policy makers.     

Can a Training Image Be a Substitute for a Random Field Model?

In most multiple-point simulation algorithms, all statistical features are provided by one or several training images (TI) that serve as a substitute for a random field model. However, because in practice the TI is always of finite size, the stochastic nature of multiple-point simulation is questionable. This issue is addressed by considering the case of a sequential simulation algorithm applied to a binary TI that is a genuine realization of an underlying random field. At each step, the algorithm uses templates containing the current target point as well as all previously simulated points. The simulation is validated by checking that all statistical features of the random field (supported by the simulation domain) are retrieved as an average over a large number of outcomes. The results are as follows. It is demonstrated that multiple-point simulation performs well whenever the TI is a complete (infinitely large) realization of a stationary, ergodic random field. As soon as the TI is restricted to a limited domain, the statistical features cannot be obtained exactly, but integral range techniques make it possible to predict how much the TI should be extended to approximate them up to a prespecified precision. Moreover, one can take advantage of extending the TI to reduce the number of disruptions in the execution of the algorithm, which arise when no conditioning template can be found in the TI.     

Responses of microorganisms and enzymes to soil contamination with metazachlor

A greenhouse pot experiment was conducted at the University of Warmia and Mazury in Olsztyn, with the aim of describing the influence of metazachlor on counts and biodiversity of soil microorganisms, soil enzymatic activity, physicochemical properties of soil and yield of spring oilseed rape. The first experimental factor was soil contamination with increasing rates of metazachlor: 0 (soil without the herbicide), 0.333 (recommended by the manufacturer), 6.666, 13.332, 26.665, 53.328, 106.656 and 213.312 mg kg^?1 dm of soil, while the second factor comprised two dates of determination: 30 and 60 days after starting the experiment. The tested herbicide had an adverse effect on reproduction of all analyzed microorganisms: oligotrophic bacteria and their endospore forms, Azotobacter spp. bacteria, organotrophic bacteria, actinomycetes and fungi. The values of the colony development and eco-physiological diversity indices decreased under the influence of excessive doses of the herbicide. Also, significant fluctuation in the enzymatic activity of soil was observed in response to the higher doses of metazachlor and depending on the date of determinations. The herbicide had an invariably negative influence on the activity of soil enzymes, causing the inhibition of dehydrogenases, catalase, urease, acid phosphatase, alkaline phosphatase, arylsulfatase and β-glucosidase. The physicochemical status of soil depended significantly on the degree of soil contamination with the herbicide, same as the yield of spring oilseed rape.     

Streamflow response to regional climate model output in the mountainous watershed: a case study from the Swiss Alps

Regional climate model (RCM) outputs are often used in hydrological modeling, in particular for streamflow forecasting. The heterogeneity of the meteorological variables such as precipitation, temperature, wind speed and solar radiation often limits the ability of the hydrological model performance. This paper assessed the sensitivity of RCM outputs from the PRUDENCE project and their performance in reproducing the streamflow. The soil and water assessment tool was used to simulate the streamflow of the Rhone River watershed located in the southwestern part of Switzerland, with the climate variables obtained from four RCMs. We analyzed the difference in magnitude of precipitation, maximum and minimum air temperature, and wind speed with respect to the observed values from the meteorological stations. In addition, we also focused on the impact of the grid resolution on model performance, by analyzing grids with resolutions of 50 × 50 and 25 × 25 km². The variability of the meteorological inputs from various RCMs is quite severe in the studied watershed. Among the four different RCMs, the Danish Meteorological Institute provided the best performance when simulating runoff. We found that temperature lapse rate is significantly important in the mountainous snow and glacier dominated watershed as compared to other variables like precipitation, and wind speed for hydrological performance. Therefore, emphasis should be given to minimum and maximum temperature in the bias correction studies for downscaling climatic data for impact modeling in the mountainous snow and glacier dominated complex watersheds.     

Sedimentological and stratigraphic evolution of northern Lebanon since the Late Cretaceous: implications for the Levant margin and basin

This paper presents an updated review of the Upper Mesozoic and Cenozoic sedimentological and stratigraphic evolution of the Levant margin with a focus on the northern Lebanon. Facies and microfacies analysis of outcrop sections and onshore well cores (i.e., Kousba and Chekka) supported by nannofossil and planktonic foraminifers biostratigraphy, allowed to constrain the depositional environments prevailing in the Turonian to Late Miocene. The “Senonian” (a historical term used to define the Coniacian to Maastrichtian) source rock interval was subdivided into four sub-units with related outer-shelfal facies: (1) Upper Santonian, (2) Lower, (3) Upper Campanian, and (4) Lower Maastrichtian. This Upper Cretaceous rock unit marks the major drowning of the former Turonian rudist platform. This paper confirms the Late Lutetian to Late Burdigalian hiatus, which appears to be a direct consequence of major geodynamic events affecting the Levant region (i.e., the continued collision of Afro-Arabia with Eurasia), potentially enhanced by regressional cycles (e.g., Rupelian lowstand). The distribution of Late Burdigalian–Serravallian rhodalgal banks identified in northern Lebanon was controlled by pre-existing structures inherited from the pulsating onshore deformation. Reef barriers facies occur around the Qalhat anticline, separating an eastern, restricted back-reef setting from a western, coastal to open marine one. The acme of Mount Lebanon’s uplift and exposure is dated back to the Middle–Late Miocene; it led to important erosion of carbonates that were subsequently deposited in paleo-topographic lows. The Late Cretaceous to Cenozoic facies variations and hiatuses show that the northern Lebanon was in a higher structural position compared to the south since at least the Late Cretaceous.     

Behaviour of radiogenic Pb in zircon during ultrahigh-temperature metamorphism: an ion imaging and ion tomography case study from the Kerala Khondalite Belt,southern India

Zircon crystals from a locally charnockitized Paleoproterozoic high-K metagranite from the Kerala Khondalite Belt (KKB) of southern India have been investigated by high-spatial resolution secondary ion mass spectrometry analysis of U–Th–Pb and rare earth elements (REE), together with scanning ion imaging and scanning ion tomography (depth-profiled ion imaging). The spot analyses constrain the magmatic crystallization age of the metagranite to ca. 1,850 Ma, with ultrahigh-temperature (UHT) metamorphism occurring at ca. 570 Ma and superimposed charnockite formation at ca. 520–510 Ma, while the ion imaging reveals a patchy distribution of radiogenic Pb throughout the zircon cores. Middle- to heavy-REE depletion in ca. 570 Ma zircon rims suggests that these grew in equilibrium with garnet and therefore date the UHT metamorphism in the KKB. The maximum apparent ²⁰⁷Pb/²⁰⁶Pb age obtained from the unsupported radiogenic Pb concentrations is also consistent with formation of the Pb patches during this event. The superimposed charnockitization event appears to have caused additional Pb-loss in the cores and recrystallization of the rims. The results of depth-profiling of the scanning ion tomography image stack show that the Pb-rich domains range in size from <5 nm to several 10 nm (diameter if assumed to be spherical). The occurrence of such patchy Pb has previously been documented only from UHT metamorphic zircon, where it likely results from annealing of radiation-damaged zircon. The formation of a discrete, heterogeneously distributed and subsequently immobile Pb phase effectively arrests the normal Pb-loss process seen at lower grades of metamorphism.     

Petrogenesis of a voluminous Quaternary adakitic volcano: the case of Baru volcano

The origin of adakite magmas remains controversial because initially the term adakite had petrogenetic significance implying an origin by direct melting of the eclogitized subducting oceanic crust. Many models have been produced for their origin, and until now there has not been a straightforward method to discriminate between these models in a given adakite suite. Here, we use detailed chronological and geochemical studies of selected adakitic edifices that allows for the determination of the magmatic output rate parameter (Qe), which has been correlated with the rates of magma generation deep within subduction zones. By providing temporal and eruption rate estimates, we provide constraints on the possible petrogenetic processes involved in the generation of adakite-like signatures. Adakite magmas derived from the melting of the subducting slab should be volumetrically insignificant when compared to the adakite-like magmas produced by typical arc magma generation processes. In this study, we use this observation and the extraordinary stratigraphic exposure from Miocene to present in an adakitic volcano in Panama and to study the temporal and chemical variation in erupted magmas to estimate rates of magma generation. Detailed chemical and geochronological analyses of Baru volcano indicate that the volcanic edifice was constructed in its entirety during the Quaternary and magmas display adakite-like features such as steep rare earth elements patterns, pronounced depletions in the heavy rare earth elements, low Y, high Sr, and high Sr/Y. The magmatic output rates (Qe) that we have calculated show that compared to other typical adakitic volcanoes, most of the volcanic edifice of Baru volcano was constructed extremely rapidly (<~213 k.a.) and in time frames that are similar to typical arc volcanoes. The observed chemical and mineralogical variation, coupled with the high magma production rates, indicate that Baru volcano is more representative of a typical arc volcano than a small-volume melt of the subducting oceanic crust. The technique we outline may have broader application in determining the petrogenetic conditions of other adakite suites.     

Origin and significance of poikilitic and mosaic peridotite xenoliths in the western Pannonian Basin: geochemical and petrological evidences

Peridotite xenoliths erupted by alkali basaltic volcanoes in the western Pannonian Basin can be divided into two fundamentally contrasting groups. Geochemical characteristics of the abundant protogranular, porphyroclastic and equigranular nodules suggest that these samples originate from an old consolidated and moderately depleted lithospheric mantle domain. In contrast, the geochemical features of the worldwide rare, but in the Pannonian Basin relatively abundant, poikilitic xenoliths attest to a more complex evolution. It has been argued that the origin of the peculiar texture and chemistry may be intimately linked to melt/rock reactions at successively decreasing liquid volumes in a porous melt flow system. The most likely site where such reactions can take place is the asthenosphere–lithosphere boundary. In this context, poikilitic xenoliths may provide petrological and geochemical evidence for reactions between magmatic liquids issued from the uprising asthenosphere and the solid mantle rocks of the lithosphere. These reactions are important agents of the thermal erosion of the lithosphere; thus, they could have considerably contributed to the thinning of the lithosphere in the Pannonian region. We suggest that in the Pannonian Basin, there could be a strong relation between the unusual abundance of poikilitic mantle xenoliths and the strongly eroded lithosphere.