Parameterization of historical earthquakes in Switzerland

Macroseismic earthquake parameters of historical events have been reassessed in the framework of the update of the Earthquake Catalogue of Switzerland ECOS-09. The Bakun and Wentworth method (Bakun and Wentworth 1997) has been used to assess location, magnitude, and, when possible, focal depth. We apply a two-step procedure. Intensity attenuation is assessed first by fitting a model with a logarithmic and a linear term, using a set of 111 earthquakes. The magnitude range is 3 and 5.8. Then, intensity to magnitude relation is developed. A subset of the 111 events, all having an instrumental moment magnitude, was used to perform this intensity to magnitude calibration. Five final calibration strategies were developed based on different intensity calibration datasets, regionalized or non-regionalized models, and fixed or variable source depth. The final assessment of the macroseismic earthquake parameters is based on an expert judgment procedure, using the results derived from all five strategies, and taking into consideration the historical knowledge available for the particular earthquake. A bootstrap procedure has been applied to assess the uncertainty of parameters. Indicative lower and upper bounds of uncertainty are derived from distributions of location and magnitude for a number of events, obtained through bootstrap sampling of the intensity field and of the single intensity values. The final uncertainties are given in terms of parameter uncertainty classes already used in previous versions of the earthquake catalogue of Switzerland.     

Age modeling of young non-varved lake sediments: methods and limits. Examples from two lakes in Central Chile

High-resolution and highly precise age models for recent lake sediments (last 100–150 years) are essential for quantitative paleoclimate research. These are particularly important for sedimentological and geochemical proxies, where transfer functions cannot be established and calibration must be based upon the relation of sedimentary records to instrumental data. High-precision dating for the calibration period is most critical as it determines directly the quality of the calibration statistics. Here, as an example, we compare radionuclide age models obtained on two high-elevation glacial lakes in the Central Chilean Andes (Laguna Negra: 33°38′S/70°08′W, 2,680 m a.s.l. and Laguna El Ocho: 34°02′S/70°19′W, 3,250 m a.s.l.). We show the different numerical models that produce accurate age-depth chronologies based on ²¹⁰Pb profiles, and we explain how to obtain reduced age-error bars at the bottom part of the profiles, i.e., typically around the end of the 19th century. In order to constrain the age models, we propose a method with five steps: (i) sampling at irregularly-spaced intervals for ²²⁶Ra, ²¹⁰Pb and ¹³⁷Cs depending on the stratigraphy and microfacies, (ii) a systematic comparison of numerical models for the calculation of ²¹⁰Pb-based age models: constant flux constant sedimentation (CFCS), constant initial concentration (CIC), constant rate of supply (CRS) and sediment isotope tomography (SIT), (iii) numerical constraining of the CRS and SIT models with the ¹³⁷Cs chronomarker of AD 1964 and, (iv) step-wise cross-validation with independent diagnostic environmental stratigraphic markers of known age (e.g., volcanic ash layer, historical flood and earthquakes). In both examples, we also use airborne pollutants such as spheroidal carbonaceous particles (reflecting the history of fossil fuel emissions), excess atmospheric Cu deposition (reflecting the production history of a large local Cu mine), and turbidites related to historical earthquakes. Our results show that the SIT model constrained with the ¹³⁷Cs AD 1964 peak performs best over the entire chronological profile (last 100–150 years) and yields the smallest standard deviations for the sediment ages. Such precision is critical for the calibration statistics, and ultimately, for the quality of the quantitative paleoclimate reconstruction. The systematic comparison of CRS and SIT models also helps to validate the robustness of the chronologies in different sections of the profile. Although surprisingly poorly known and under-explored in paleolimnological research, the SIT model has a great potential in paleoclimatological reconstructions based on lake sediments.     

Frequency–Magnitude Distribution of −3.7 ≤ M W ≤ 1 Mining-Induced Earthquakes Around a Mining Front and b Value Invariance with Post-Blast Time

We investigated frequency-magnitude distribution (FMD) of acoustic emissions (AE) occurring near an active mining front in a South African gold mine, using a catalog developed from an AE network, which is capable of detecting AEs down to M _W  ?5. When records of blasts were removed, FMDs of AEs obeyed a Gutenberg?Richter law with similar b values, not depending on post-blasting time from the initial 1-min interval through more than 30 h. This result denies a suggestion in a previous study (Richardson and Jordan Bull Seismol Soc Am, 92:1766–1782, 2002) that new fractures generated by blasting disturb the size distribution of background events, which they interpreted as slip events on existing weak planes. Our AE catalog showed that the GR law with b ～ 1.2 was valid between M _W  ?3.7 and 0 for AEs around the mining front. Further, using the mine’s seismic catalog, which covers a longer time period of the same area, we could extend the validity range of the GR law with the same b value up to M _W 1.     

Primitive andesites from the Taupo Volcanic Zone formed by magma mixing

Andesites with Mg# >45 erupted at subduction zones form either by partial melting of metasomatized mantle or by mixing and assimilation processes during melt ascent. Primitive whole rock basaltic andesites from the Pukeonake vent in the Tongariro Volcanic Centre in New Zealand’s Taupo Volcanic Zone contain olivine, clino- and orthopyroxene, and plagioclase xeno- and antecrysts in a partly glassy matrix. Glass pools interstitial between minerals and glass inclusions in clinopyroxene, orthopyroxene and plagioclase as well as matrix glasses are rhyolitic to dacitic indicating that the melts were more evolved than their andesitic bulk host rock analyses indicate. Olivine xenocrysts have high Fo contents up to 94%, δ¹⁸O_(SMOW) of +5.1‰, and contain Cr-spinel inclusions, all of which imply an origin in equilibrium with primitive mantle-derived melts. Mineral zoning in olivine, clinopyroxene and plagioclase suggest that fractional crystallization occurred. Elevated O isotope ratios in clinopyroxene and glass indicate that the lavas assimilated sedimentary rocks during stagnation in the crust. Thus, the Pukeonake andesites formed by a combination of fractional crystallization, assimilation of crustal rocks, and mixing of dacite liquid with mantle-derived minerals in a complex crustal magma system. The disequilibrium textures and O isotope compositions of the minerals indicate mixing processes on timescales of less than a year prior to eruption. Similar processes may occur in other subduction zones and require careful study of the lavas to determine the origin of andesite magmas in arc volcanoes situated on continental crust.     

A new occurrence of extremely rare Ferro-Chloro-Pargasite in Tudor,Ontario

International Journal of Earth Sciences -     

A mesocosm experiment investigating the effects of substratum quality and wave exposure on the survival of fish eggs

In a mesocosm experiment, the attachment of bream (Abramis brama) eggs to spawning substrata with and without periphytic biofilm coverage and their subsequent survival with and without low-intensity wave exposure were investigated. Egg attachment was reduced by 73% on spawning substrata with a natural periphytic biofilm, compared to clean substrata. Overall, this initial difference in egg numbers persisted until hatching. The difference in egg numbers was even increased in the wave treatment, while it was reduced in the no-wave control treatment. Exposure to a low-intensity wave regime affected egg development between the two biofilm treatments differently. Waves enhanced egg survival on substrata without a biofilm but reduced the survival of eggs on substrata with biofilm coverage. In the treatment combining biofilm-covered substrata and waves, no attached eggs survived until hatching. In all treatments, more than 75% of the eggs became detached from the spawning substrata during the egg incubation period, and <1% of these detached eggs survived within the substratum interstices. Hence, detached eggs contributed little to the reproductive success we observed. The implications of these results on the spawning success of fish species that use freshly inundated spawning substrata in shallow waters are discussed.     

Earthquakes in Switzerland and surrounding regions during 2008

This report of the Swiss Seismological Service summarizes the seismic activity in Switzerland and surrounding regions during 2008. During this period, 451 earthquakes and 75 quarry blasts were detected and located in the region under consideration. The three strongest events occurred in the Valais, near Lac des Toules (M_L 3.6), and in Graubünden, near Ilanz (M_L 3.7) and Paspels (M_L 4.0). Although felt by the population, they were not reported to have caused any damage. However, with a total of only 15 events with M_L ≥ 2.5, the seismic activity in the year 2008 was far below the average over the previous 33 years.     

Weathering and the mobility of phosphorus in the catchments and forefields of the Rhône and Oberaar glaciers, central Switzerland: Implications for the global phosphorus cycle on glacial-interglacial timescales

In this study we evaluate the dynamics of the biophile element phosphorus (P) in the catchment and proglacial areas of the Rhône and Oberaar glaciers (central Switzerland). We analysed erosion and dissolution rates of P-containing minerals in the subglacial environment by sampling water and suspended sediment in glacier outlets during three ablation and two accumulation seasons. We also quantified biogeochemical weathering rates of detrital P in proglacial sedimentary deposits using two chronosequences of samples of fresh, suspended, material obtained from the Oberaar and Rhône water outlets, Little-Ice-Age (LIA) moraines and Younger Dryas (YD) tills in each catchment. Subglacial P weathering is mainly a physical process and detrital P represents more than 99% of the precipitation-corrected total P denudation flux (234 and 540 kg km⁻² yr⁻¹ for the Rhône and Oberaar catchments, respectively). The calculated detrital P flux rates are three to almost five times higher than the world average flux. The precipitation-corrected soluble reactive P (SRP) flux corresponds to 1.88-1.99 kg km⁻² yr⁻¹ (Rhône) and 2.12-2.44 kg km⁻² yr⁻¹ (Oberaar), respectively. These fluxes are comparable to those of tropical rivers draining transport-limited, tectonically inactive weathering areas.In order to evaluate the efficiency of detrital P weathering in the Rhône and Oberaar proglacial areas, we systematically graded apatite grains extracted from the chronosequence in each catchment relative to weathering-induced changes in their surface morphologies (grades 1-4). Fresh apatite grains are heavily indented and dissolution rounded (grade 1). LIA grains from two 0-10 cm deep moraine samples show extensive dissolution etching, similar to surface grains from the YD profile (mean grades 2.7, 3.5 and 3.5, respectively). In these proglacial deposits, the weathering front deepens progressively as a function of time due to biocorrosion in the evolving acidic pedosphere, with mechanical indentations on grains acting as sites of preferential dissolution. We also measured iron-bound, organic and detrital P concentrations in the chronosequence and show that organic and iron-bound P has almost completely replaced detrital P in the top layers of the YD profiles. Detrital P weathering rates are calculated as 310 and 280 kg km⁻² yr⁻¹ for LIA moraines and 10 kg km⁻² yr⁻¹ for YD tills. During the first 300 years of glacial sediment exposure P dissolution rates are shown to be approximately 70 times higher than the mean global dissolved P flux from ice-free continents. After 11.6 kyr the flux is 2.5 times the global mean. These data strengthen the argument for substantial changes in the global dissolved P flux on glacial-interglacial timescales. A crude extrapolation from the data described here suggests that the global dissolved P flux may increase by 40-45% during the first few hundred years of a deglaciation phase.     

Possible Effects of Altered Growth Behaviour of Norway Spruce (Picea abies) on Carbon Accounting

Implementing the Kyoto Protocol necessitates precise and practical inventory methods for estimating the carbon reservoir in forests. An analysis of the German state North Rhine-Westphalia forest inventory data showed enormous deviations from the yield tables of Wiedemann (1936, 1942). Specifically the ratio of height/diameter at breast height was up to 30% lower than in the yield tables, which indicated wider annual rings and a lower basic density. Such differences most likely pose consequences for the calculation of the amount of carbon in these forests. In this study the aboveground green weight, moisture and the carbon content of fifteen 46-year old Norway spruces (Picea abies) were measured, and two approaches for calculating the aboveground carbon content of trees were discussed: the common methodology of applying expansion and conversion factors as well as using optimized regression equations. The average expansion factor deviated approximately 1.2% from the yield tables of Grundner and Schwappach (1906). The wood density was 368 kg/m³ and within the European scope of 320–420 kg/m³ (Bosshard, 1984; Hakkila, 1989; Knigge and Schulz, 1966; Trendelenburg and Mayer-Wegelin, 1955). The woody carbon content was found to be 50.42%, which corresponded well with reference values of 50.3% (Bosshard, 1984) and 51.4% (Knigge and Schulz, 1966). However, the average percentage of crown material carbon content was 0.75% higher than in the tree boles. This study demonstrated that the most reliable methodology for calculating aboveground tree carbon was a power regression function, which calculated the growing stock volume, converting it directly into tree carbon content. The results deviated by less than 1.3% (R² = 0.99) from the measured tree carbon. Using volume expansion factors and factors of basic density and carbon percentage, as in Burschel et al. (1993), the calculated results diverged –3.7% from the measured outcome. The adapted form according to Fang et al. (2001) deviated –1.2% (R² = 0.63) and the form following Brown and Schroeder (1999) about –2.5% (R² = 0.61) from the measured carbon of the sample trees. However, coefficients of variation from 12% and 13% for the presented algorithm illustrated that in reality the calculated carbon of single trees could diverge considerably from the model results.