3-D Shear Wave Velocity Structure of the Crust and Upper Mantle in South China Sea and its Surrounding Regions by Surface Wave Dispersion Analysis

In this study, we construct a 3-D shear wave velocity structure of the crust and upper mantle in South China Sea and its surrounding regions by surface wave dispersion analysis. We use the multiple filter technique to calculate the group velocity dispersion curves of fundamental mode Rayleigh and Love waves with periods from 14 s to 120 s for earthquakes occurred around the Southeast Asia. We divide the study region (80° E–140° E, 16° S–32° N) into 3° × 3° blocks and use the constrained block inversion method to get the regionalized dispersion curve for each block. At some chosen periods, we put together laterally the regionalized group velocities from different blocks at the same period to get group velocity image maps. These maps show that there is significant heterogeneity in the group velocity of the study region. The dispersion curve of each block was then processed by surface wave inversion method to obtain the shear wave velocity structure. Finally, we put the shear wave velocity structures of all the blocks together to obtain the three-dimensional shear wave velocity structure of crust and upper mantle. The three-dimensional shear wave velocity structure shows that the shear wave velocity distribution in the crust and upper mantle of the South China Sea and its surrounding regions displays significant heterogeneity. There are significant differences among the crustal thickness, the lithospheric thickness and the shear wave velocity of the lid in upper mantle of different structure units. This study shows that the South China Sea Basin, southeast Sulu Sea Basin and Celebes Sea Basin have thinner crust. The thickness of crust in South China Sea Basin is 5–10 km; in Indochina is 25–40 km; in Peninsular Malaysia is 30–35 km; in Borneo is 30–35 km; in Palawan is 35 km; in the Philippine Islands is 30–35 km, in Sunda Shelf is 30–35 km, in Southeast China is 30–40 km, in West Philippine Basin is 5–10 km. The South China Sea Basin has a lithosphere with thickness of about 45–50 km, and the shear wave velocity of its lid is about 4.3–4.7 km/s; Indochina has a lithosphere with thickness of about 55–70 km, and the shear wave velocity of its lid is about 4.3–4.5 km/s; Borneo has a lithosphere with thickness of about 55–60 km, and the shear wave velocity of its lid is about 4.1–4.3 km/s; the Philippine Islands has a lithosphere with thickness of about 55–60 km, and the shear wave velocity of its lid is about 4.2–4.3 km/s, West Philippine Basin has a lithosphere with thickness of about 50–55 km, and the shear wave velocity of its lid is about 4.7–4.8 km/s, Sunda Self has a lithosphere with thickness of about 55–65 km, and the shear wave velocity of its lid is about 4.3 km/s. The Red-River Fault Zone probably penetrates to a depth of at least 200 km and is plausibly the boundary between the South China Block and the Indosinia Block.     

Lumped surface and sub‐surface runoff for erosion modeling within a small hilly watershed in northern Vietnam

Developing models to predict on‐site soil erosion and off‐site sediment transport at the agricultural watershed scale represent an on‐going challenge in research today. This study attempts to simulate the daily discharge and sediment loss using a distributed model that combines surface and sub‐surface runoffs in a small hilly watershed (< 1 km²). The semi‐quantitative model, Predict and Localize Erosion and Runoff (PLER), integrates the Manning–Strickler equation to simulate runoff and the Griffith University Erosion System Template equation to simulate soil detachment, sediment storage and soil loss based on a map resolution of 30 m × 30 m and over a daily time interval. By using a basic input data set and only two calibration coefficients based, respectively, on water velocity and soil detachment, the PLER model is easily applicable to different agricultural scenarios. The results indicate appropriate model performance and a high correlation between measured and predicted data with both Nash–Sutcliffe efficiency (Ef) and correlation coefficient (r²) having values > 0.9. With the simple input data needs, PLER model is a useful tool for daily runoff and soil erosion modeling in small hilly watersheds in humid tropical areas. Copyright © 2013 John Wiley & Sons, Ltd.     

Equation of state and pressure-induced structural changes in mirabilite (Na2SO4·10H2O) determined from ab initio density functional theory calculations

We have carried out ab initio calculations using density functional theory to determine the bulk elastic properties of mirabilite, Na₂SO₄·10H₂O, and to obtain information on structural trends caused by the application of high pressure up to ~60 GPa. We have found that there are substantial isosymmetric discontinuous structural re-organisations at ~7.7 and ~20 GPa caused by changes in the manner in which the sodium cations are coordinated by water molecules. The low-pressure and intermediate-pressure phases both have sodium in sixfold coordination but in the high-pressure phase the coordination changes from sixfold to sevenfold. These coordination changes force a re-arrangement of the hydrogen-bond network in the crystal. The trend is towards a reduction in the number of hydrogen bonds donated to the sulphate group (from twelve down to six over the range 0–60 GPa) and an increase in hydrogen bonding amongst the Na-coordinated water molecules and the two interstitial water molecules. Ultimately, we observe proton transfers from the interstitial waters (forming OH^? ions) to two of the Na-coordinated waters (forming a pair of H₃O⁺ ions). The equation of state in the athermal limit of the low-pressure phase of mirabilite, parameterised by fitting an integrated form of the third-order Birch-Murnaghan expression to the calculated energy as a function of unit-cell volume, yields the zero-pressure unit-cell volume, V ₀ = 1468.6(9) ų, the incompressibility, K ₀ = 22.21(9) GPa, and the first pressure derivative K ₀′ = (?K/?P)₀ = 5.6(1).     

Analysing the interaction between emission trading and renewable electricity support in TIMES

As the number of instruments applied in the area of energy and climate policy is rising, the issue of policy interaction needs to be explored further. This article analyses the interdependencies between the EU Emissions Trading Scheme (EU ETS) and the German feed-in tariffs (FITs) for renewable electricity in a quantitative manner using a bottom-up energy system model. Flexible modelling approaches are presented for both instruments, with which all impacts on the energy system can be evaluated endogenously. It is shown that national climate policy measures can have an effect on the supranational emissions trading system by increasing emission reduction in the German electricity sector by up to 79 MtCO₂ in 2030. As a result, emission certificate prices decline by between 1.9 €/tCO₂ and 6.1 €/tCO₂ and the burden sharing between participating countries changes, but no additional emission reduction is achieved at the European level. This also implies, however, that the cost efficiency of such a cap-and-trade system is distorted, with additional costs of the FIT system of up to €320 billion compared with lower costs for ETS emission certificates of between €44 billion and €57 billion (cumulated over the period 2013–2020).

Policy relevance

In order to fulfil ambitious emission reduction targets a large variety of climate policy instruments are being implemented in Europe. While some, like the EU ETS, directly address CO₂ emissions, others aim to promote specific low-carbon technologies. The quantitative analysis of the interactions between the EU ETS and the German FIT scheme for renewable sources in electricity generation presented in this article helps to understand the importance of such interaction effects. Even though justifications can be found for the implementation of both types of instrument, the impact of the widespread use of support mechanisms for renewable electricity in Europe needs to be taken into account when fixing the reduction targets for the EU ETS in order to ensure a credible long-term investment signal.     

Report on a recent deep-seated landslide at G��rov�� Mt., Czech Republic, triggered by a heavy rainfall: The G��rov�� Mt., Outer West Carpathians; Czech Republic

A large deep-seated landslide, triggered by a heavy rainfall, activated early in the morning on May 19, 2010 along the SW slope of Gírová Mt., NE Czech Republic. The landslide occurred within a zone of pre-existing deep-seated gravitational deformation, and it was accompanied with pronounced ground liquefying in the central and lower portion of the sliding mass. The precipitation that triggered the landslide was about 244.6 mm that fell between May 15th and 18th with an average of 61 mm/day. No properties or lives were lost. However, the landslide provided a good case study on triggering factors and process of liquefaction at mountain slopes in the E Czech Republic and adjacent areas.     

Modelling indicators of water security,water pollution and aquatic biodiversity in Europe

Abstract

The GWAVA (Global Water AVailability Assessment) model for indicating human water security has been extended with a newly developed module for calculating pollutant concentrations. This module is first described and then illustrated by being used to model nitrogen, phosphorus and organic matter concentrations. The module uses solely input variables that are likely to be available for future scenarios, making it possible to apply the module to such scenarios. The module first calculates pollutant loading from land to rivers, lakes and wetlands by considering drivers such as agriculture, industry and sewage treatment. Calculated loadings are subsequently converted to concentrations by considering aquatic processes, such as dilution, downstream transport, evaporation, human water abstraction and biophysical loss processes. Aquatic biodiversity is indicated to be at risk if modelled pollutant concentrations exceed certain water quality standards. This is indicated to be the case in about 35% of the European area, especially where lakes and wetlands are abundant. Human water security is indicated to be at risk where human water demands cannot be fulfilled during drought events. This is found to be the case in about 10% of the European area, especially in Mediterranean, arid and densely-populated areas. Modelled spatial variation in concentrations matches well with existing knowledge, and the temporal variability of concentrations is modelled reasonably well in some river basins. Therefore, we conclude that the updated GWAVA model can be used for indicating changes in human water security and aquatic biodiversity across Europe.

Editor Z.W. Kundzewicz

Citation Dumont, E., Williams, R., Keller, V., Voss, A., and Tattari, S., 2012. Modelling indicators of water security, water pollution and aquatic biodiversity in Europe. Hydrological Sciences Journal, 57 (7), 1378–1403.     

The thermal expansion and crystal structure of mirabilite (Na<Subscript>2</Subscript>SO<Subscript>4</Subscript>·10D<Subscript>2</Subscript>O) from 4.2 to 300 K,determined by time-of-flight neutron powder diffraction

We have collected high resolution neutron powder diffraction patterns from Na₂SO₄·10D₂O over the temperature range 4.2–300 K following rapid quenching in liquid nitrogen, and over a series of slow warming and cooling cycles. The crystal is monoclinic, space-group P2₁/c (Z = 4) with a = 11.44214(4) Å, b = 10.34276(4) Å, c = 12.75486(6) Å, β = 107.847(1)°, and V = 1436.794(8) Å³ at 4.2 K (slowly cooled), and a = 11.51472(6) Å, b = 10.36495(6) Å, c = 12.84651(7) Å, β = 107.7543(1)°, V = 1460.20(1) Å³ at 300 K. Structures were refined to R _P (Rietveld powder residual, \( R_{P} = {{\sum {\left| {I_{\text{obs}} - I_{\text{calc}} } \right|} } \mathord{\left/ {\vphantom {{\sum {\left| {I_{\text{obs}} - I_{\text{calc}} } \right|} } {\sum {I_{\text{obs}} } }}} \right. \kern-\nulldelimiterspace} {\sum {I_{\text{obs}} } }} \)) better than 2.5% at 4.2 K (quenched and slow cooled), 150 and 300 K. The sulfate disorder observed previously by Levy and Lisensky (Acta Cryst B34:3502–3510, 1978) was not present in our specimen, but we did observe changes with temperature in deuteron occupancies of the orientationally disordered water molecules coordinated to Na. The temperature dependence of the unit-cell volume from 4.2 to 300 K is well represented by a simple polynomial of the form V = ? 4.143(1) × 10^?7 T ³ + 0.00047(2) T² ? 0.027(2) T + 1437.0(1) Å³ (R ² = 99.98%). The coefficient of volume thermal expansion, α _V , is positive above 40 K, and displays a similar magnitude and temperature dependence to α _V in deuterated epsomite and meridianiite. The relationship between the magnitude and orientation of the principal axes of the thermal expansion tensor and the main structural elements are discussed; freezing in of deuteron disorder in the quenched specimen affects the thermal expansion, manifested most obviously as a change in the behaviour of the unit-cell parameter β.     

The effect of ferromagnetism on the equation of state of Fe3C studied by first-principles calculations

First-principles calculations have been used to determine the equation of state of Fe₃C in both its low-pressure magnetically ordered and high-pressure non-magnetically ordered states; at 0 K the ferromagnetic transition was found to occur at about 60 GPa. In the high pressure, non-magnetically ordered regime at 0 K the material may be described by a Birch-Murnaghan third-order equation of state with V₀=8.968(7) ų per atom, K₀=316.62(2) GPa and K′=4.30(2). At atmospheric pressure the ferromagnetic phase transition in Fe₃C occurs at ∼483 K; preliminary measurements of the thermal expansion by powder neutron diffraction show that this transition produces a large effect on thermoelastic properties. The volumetric thermal expansion coefficient in the paramagnetic phase was found to be 4.34×10⁻⁵ K⁻¹ at T∼550 K. By applying a thermal expansion correction to the calculated equation of state at 0 K, predicted values for the density and adiabatic incompressibility of this material at core pressures and temperatures were obtained. These results appear to be sufficiently different from seismological data so as to preclude Fe₃C as the major inner core-forming phase.