Considering large earthquake clustering in seismic risk analysis

Earthquakes cluster in space and time resulting in nonlinear damage effects. We compute earthquake interactions using the Coulomb stress transfer theory and dynamic vulnerability from the concept of ductility capacity reduction. We combine both processes in the generic multi-risk framework where risk scenarios are simulated using a variant of the Markov chain Monte Carlo method. We apply the proposed approach to the thrust fault system of northern Italy, considering earthquakes with characteristic magnitudes in the range ~[6, 6.5], different levels of tectonic loading \(\dot{\tau }\) = {10⁻⁴, 10⁻³, 10⁻²} bar/year and a generic stock of fictitious low-rise buildings with different ductility capacities μ _Δ = {2, 4, 6}. We describe the process’ stochasticity by non-stationary Poisson earthquake probabilities and by binomial damage state probabilities. We find that earthquake clustering yields a tail fattening of the seismic risk curve, the effect of which is amplified by damage-dependent fragility due to clustering. The impact of clustering alone is in average more important than dynamic vulnerability, the spatial extent of the former phenomenon being greater than of the latter one.

     

Thermal equation of state of natural tourmaline at high pressure and temperature

Synchrotron-based in situ angle-dispersive X-ray diffraction experiments were conducted on a natural uvite-dominated tourmaline sample by using an external-heating diamond anvil cell at simultaneously high pressures and temperatures up to 18 GPa and 723 K, respectively. The angle-dispersive X-ray diffraction data reveal no indication of a structural phase transition over the P–T range of the current experiment in this study. The pressure–volume–temperature data were fitted by the high-temperature Birch–Murnaghan equation of state. Isothermal bulk modulus of K ₀ = 96.6 (9) GPa, pressure derivative of the bulk modulus of \(K_{0}^{\prime } = 12.5 \;(4)\), thermal expansion coefficient of α ₀ = 4.39 (27) × 10^?5 K^?1 and temperature derivative of the bulk modulus (?K/?T) _P  = ?0.009 (6) GPa K^?1 were obtained. The axial thermoelastic properties were also obtained with K _a0 = 139 (2) GPa, \(K_{a0}^{\prime }\) = 11.5 (7) and α _a0 = 1.00 (11) × 10^?5 K^?1 for the a-axis, and K _c0 = 59 (1) GPa, \(K_{c0}^{\prime }\) = 11.4 (5) and α _c0 = 2.41 (24) × 10^?5 K^?1 for the c-axis. Both of axial compression and thermal expansion exhibit large anisotropic behavior. Thermoelastic parameters of tourmaline in this study were also compared with that of the other two ring silicates of beryl and cordierite.     

The effects of undercooling and deformation rates on the crystallization kinetics of Stromboli and Etna basalts

We have investigated the effect of undercooling and deformation on the evolution of the texture and the crystallization kinetics of remelted basaltic material from Stromboli (pumice from the March 15, 2007 paroxysmal eruption) and Etna (1992 lava flow). Isothermal crystallization experiments were conducted at different degrees of undercooling and different applied strain rate (T = 1,157–1,187 °C and $ \dot{\gamma }_{i} $ γ · i  = 4.26 s^?1 for Stromboli; T = 1,131–1,182 °C and $ \dot{\gamma }_{i} $ γ · i  = 0.53 s^?1 for Etna). Melt viscosity increased due to the decrease in temperature and the increase in crystal content. The mineralogical assemblage comprises Sp + Plg (dominant) ± Cpx with an overall crystal fraction (?) between 0.06 and 0.27, increasing with undercooling and flow conditions. Both degree of undercooling and deformation rate deeply affect the kinetics of the crystallization process. Plagioclase nucleation incubation time strongly decreases with increasing ΔT and flow, while slow diffusion-limited growth characterizes low ΔT—low deformation rate experiments. Both Stromboli (high strain rate) and Etna (low strain rate) plagioclase growth rates (G) display relative small variations with Stromboli showing higher values (4.8 ± 1.9 × 10^?9 m s^?1) compared to Etna (2.1 ± 1.6 × 10^?9 m s^?1). Plagioclase average nucleation rates J continuously increase with undercooling from 1.4 × 10⁶ to 6.7 × 10⁶ m^?3 s^?1 for Stromboli and from 3.6 × 10⁴ to 4.0 × 10⁶ m^?3 s^?1 for Etna. The extremely low value of 3.6 × 10⁴ m^?3 s^?1 recorded at the lowest undercooling experiment for Etna (ΔT = 20 °C) indicates that the crystallization process is growth-dominated and that possible effects of textural coarsening occur. G values obtained in this paper are generally one or two orders of magnitude higher compared to those obtained in the literature for equivalent undercooling conditions. Stirring of the melt, simulating magma flow or convective conditions, facilitates nucleation and growth of crystals via mechanical transportation of matter, resulting in the higher J and G observed. Any modeling pertaining to magma dynamics in the conduit (e.g., ascent rate) and lava flow emplacement (e.g., flow rate, pāhoehoe–‘a‘ā transition) should therefore take the effects of dynamic crystallization into account.     

Behavior of 10-Å phase at low temperatures

The thermal evolution of 10-Å phase Mg₃Si₄O₁₀(OH)₂·H₂O, a phyllosilicate which may have an important role in the storage/release of water in subducting slabs, was studied by X-ray single-crystal diffraction in the temperature range 116–293 K. The lattice parameters were measured at several intervals both on cooling and heating. The structural model was refined with intensity data collected at 116 K and compared to the model refined at room temperature. As expected for a layer silicate on cooling in this temperature range, the a and b lattice parameters undergo a small linear decrease, α _a  = 1.7(4) 10^?6 K^?1 and α _b  = 1.9(4) 10^?6 K^?1, where α is the linear thermal expansion coefficient. The greater variation is along the c axis and can be modeled with the second order polynomial c _T  = c ₂₉₃(1 + 6.7(4)10^?5 K^?1ΔT + 9.5(2.5)10^?8 K^?2(ΔT)²) where ΔT = T ? 293 K; the monoclinic angle β slightly increased. The cell volume thermal expansion can be modeled with the polynomial V _T  = V ₂₉₃ (1 + 8.0 10^?5 K^?1 ΔT + 1.4 10^?7 K^?2 (ΔT)²) where ΔT = T ? 293 is in K and V in ų. These variations were similar to those expected for a pressure increase, indicating that T and P effects are approximately inverse. The least-squares refinement with intensity data measured at 116 K shows that the volume of the SiO₄ tetrahedra does not change significantly, whereas the volume of the Mg octahedra slightly decreases. To adjust for the increased misfit between the tetrahedral and octahedral sheets, the tetrahedral rotation angle α changes from 0.58° to 1.38°, increasing the ditrigonalization of the silicate sheet. This deformation has implications on the H-bonds between the water molecule and the basal oxygen atoms. Furthermore, the highly anisotropic thermal ellipsoid of the H₂O oxygen indicates positional disorder, similar to the disorder observed at room temperature. The low-temperature results support the hypothesis that the disorder is static. It can be modeled with a splitting of the interlayer oxygen site with a statistical distribution of the H₂O molecules into two positions, 0.6 Å apart. The resulting shortest O_bas–O_W distances are 2.97 Å, with a significant shortening with respect to the value at room temperature. The low-temperature behavior of the H-bond system is consistent with that hypothesized at high pressure on the basis of the Raman spectra evolution with P.     

Dewatering Behaviour of a Uranium Ore Slurry Containing Clays

The main objective of this paper was to investigate the dewatering behaviour of a clayey uranium ore slurry. The slurry (containing 28% clay size) exhibited moderate water adsorption (w _l  = 83% and w _p  = 30%). Primarily composed of muscovite (46%) and quartz (30%), the clay minerals included illite (8%), chlorite (5%) and kaolinite (2%) alongside a CEC of 41 (cmol(+)/kg) with Ca²⁺ and Mg²⁺ as the dominant cations. Likewise, the high EC (17,600 μS/cm) and ionic strength (1.15 mol/L) indicated a flocculated microstructure due to the presence of SO₄ ^2? (22,600 mg/L) and Mg²⁺ (1340 mg/L) in the slurry water. Settling included sedimentation and consolidation at low initial solids condition (25–35%) whereas only consolidation was observed at high initial solids contents (40–50%). The average k reduced from 1.2 × 10^?6 m/s (initial s = 25%) to 5.3 × 10^?8 m/s (initial s = 50%) along with a void ratio reduction from 7.4 to 2.6. Due to thixotropic strength, volume compressibility during consolidation showed apparent pre-consolidation at low effective stress (0.3–2 kPa) with a reduction in void ratio from 2.6 to 2.5. The e _s was found to be 2.46 at σ′ = 2 kPa and was followed by a steeper slope with the void ratio reducing to 2.1 at σ′ = 31 kPa. Likewise, the hydraulic conductivity during consolidation decreased from 2.6 × 10^?9 m/s (at e = 2.6) to 2.0 × 10^?10 m/s (at e = 2.1).     

The P–V–T equation of state of CaPtO3 post-perovskite

Orthorhombic post-perovskite CaPtO₃ is isostructural with post-perovskite MgSiO₃, a deep-Earth phase stable only above 100 GPa. Energy-dispersive X-ray diffraction data (to 9.4 GPa and 1,024 K) for CaPtO₃ have been combined with published isothermal and isobaric measurements to determine its P–V–T equation of state (EoS). A third-order Birch–Murnaghan EoS was used, with the volumetric thermal expansion coefficient (at atmospheric pressure) represented by α(T) = α₀ + α₁(T). The fitted parameters had values: isothermal incompressibility, $ K_{{T_{0} }} $  = 168.4(3) GPa; $ K_{{T_{0} }}^{\prime } $  = 4.48(3) (both at 298 K); $ \partial K_{{T_{0} }} /\partial T $  = ?0.032(3) GPa K^?1; α₀ = 2.32(2) × 10^?5 K^?1; α₁ = 5.7(4) × 10^?9 K^?2. The volumetric isothermal Anderson–Grüneisen parameter, δ _T , is 7.6(7) at 298 K. $ \partial K_{{T_{0} }} /\partial T $ for CaPtO₃ is similar to that recently reported for CaIrO₃, differing significantly from values found at high pressure for MgSiO₃ post-perovskite (?0.0085(11) to ?0.024 GPa K^?1). We also report axial P–V–T EoS of similar form, the first for any post-perovskite. Fitted to the cubes of the axes, these gave $ \partial K_{{aT_{0} }} /\partial T $  = ?0.038(4) GPa K^?1; $ \partial K_{{bT_{0} }} /\partial T $  = ?0.021(2) GPa K^?1; $ \partial K_{{cT_{0} }} /\partial T $  = ?0.026(5) GPa K^?1, with δ _T  = 8.9(9), 7.4(7) and 4.6(9) for a, b and c, respectively. Although $ K_{{T_{0} }} $ is lowest for the b-axis, its incompressibility is the least temperature dependent.     

Thermal equation of state of CaIrO<Subscript>3</Subscript> post-perovskite

The pressure–volume–temperature (P–V–T) relation of CaIrO₃ post-perovskite (ppv) was measured at pressures and temperatures up to 8.6 GPa and 1,273 K, respectively, with energy-dispersive synchrotron X-ray diffraction using a DIA-type, cubic-anvil apparatus (SAM85). Unit-cell dimensions were derived from the Le Bail full profile refinement technique, and the results were fitted using the third-order Birth-Murnaghan equation of state. The derived bulk modulus \( K_{T0} \) at ambient pressure and temperature is 168.3 ± 7.1 GPa with a pressure derivative \( K_{T0}^{\prime } \) = 5.4 ± 0.7. All of the high temperature data, combined with previous experimental data, are fitted using the high-temperature Birch-Murnaghan equation of state, the thermal pressure approach, and the Mie-Grüneisen-Debye formalism. The refined thermoelastic parameters for CaIrO₃ ppv are: temperature derivative of bulk modulus \( (\partial K_{T} /\partial T)_{P} \) = ?0.038 ± 0.011 GPa K^?1, \( \alpha K_{T} \) = 0.0039 ± 0.0001 GPa K^?1, \( \left( {\partial K_{T} /\partial T} \right)_{V} \) = ?0.012 ± 0.002 GPa K^?1, and \( \left( {\partial^{2} P/\partial T^{2} } \right)_{V} \) = 1.9 ± 0.3 × 10^?6 GPa² K^?2. Using the Mie-Grüneisen-Debye formalism, we obtain Grüneisen parameter \( \gamma_{0} \) = 0.92 ± 0.01 and its volume dependence q = 3.4 ± 0.6. The systematic variation of bulk moduli for several oxide post-perovskites can be described approximately by the relationship K _T0  = 5406.0/V(molar) + 5.9 GPa.     

Evaluating of the earthquake hazard parameters with Bayesian method for the different seismic source regions of the North Anatolian Fault Zone

The earthquake hazard parameters and earthquake occurrence probabilities are computed for the different regions of the North Anatolia Fault Zone (NAFZ) using Bayesian method. A homogenous earthquake catalog for M _S magnitude which is equal or larger than 4.0 is used for a time period between 1900 and 2015. Only two historical earthquakes (1766, M _S = 7. 3 and 1897, M _S = 7. 0) are included in Region 2 (Marmara Region) where a large earthquake is expected in the near future since no large earthquake has been observed for the instrumental period. In order to evaluate earthquake hazard parameters for next 5, 10, 20, 50, 100 years, M _max (maximum regional magnitude), β value, λ (seismic activity or density) are computed for the different regions of NAFZ. The computed M _max values are changed between 7.11 and 7.89. While the highest magnitude value is calculated in the Region 9 related to Tokat-Erzincan, the lowest value in the Region 10 including the eastern of Erzincan. The “quantiles” of “apparent” and “true” magnitudes of future time intervals of 5, 10, 20, 50, and 100 years are calculated for confidence limits of probability levels of 50, 70 and 90 % of the 10 different seismic source regions. The region between Tokat and Erzincan has earthquake hazard level according to the determined parameters. In this region the expected maximum earthquake size is 7.8 with 90 % occurrence probability in next 100 years. While the regional M _max value of Marmara Region is computed as 7.61, expected maximum earthquake size is 7.37 with 90 % occurrence probability in next 100 years.     

Thermoelastic properties of chromium oxide Cr<Subscript>2</Subscript>O<Subscript>3</Subscript> (eskolaite) at high pressures and temperatures

A new synchrotron X-ray diffraction study of chromium oxide Cr₂O₃ (eskolaite) with the corundum-type structure has been carried out in a Kawai-type multi-anvil apparatus to pressure of 15 GPa and temperatures of 1873 K. Fitting the Birch–Murnaghan equation of state (EoS) with the present data up to 15 GPa yielded: bulk modulus (K _0,T0), 206 ± 4 GPa; its pressure derivative K′_0,T , 4.4 ± 0.8; (?K _0,T /?T) = ?0.037 ± 0.006 GPa K^?1; a = 2.98 ± 0.14 × 10^?5 K^?1 and b = 0.47 ± 0.28 × 10^?8 K^?2, where α _0,T  = a + bT is the volumetric thermal expansion coefficient. The thermal expansion of Cr₂O₃ was additionally measured at the high-temperature powder diffraction experiment at ambient pressure and α _0,T0 was determined to be 2.95 × 10^?5 K^?1. The results indicate that coefficient of the thermal expansion calculated from the EoS appeared to be high-precision because it is consistent with the data obtained at 1 atm. However, our results contradict α ₀ value suggested by Rigby et al. (Brit Ceram Trans J 45:137–148, 1946) widely used in many physical and geological databases. Fitting the Mie–Grüneisen–Debye EoS with the present ambient and high-pressure data yielded the following parameters: K _0,T0 = 205 ± 3 GPa, K′_0,T  = 4.0, Grüneisen parameter (γ ₀) = 1.42 ± 0.80, q = 1.82 ± 0.56. The thermoelastic parameters indicate that Cr₂O₃ undergoes near isotropic compression at room and high temperatures up to 15 GPa. Cr₂O₃ is shown to be stable in this pressure range and adopts the corundum-type structure. Using obtained thermoelastic parameters, we calculated the reaction boundary of knorringite formation from enstatite and eskolaite. The Clapeyron slope (with \({\text{d}}P/{\text{d}}T = - 0.014\) GPa/K) was found to be consistent with experimental data.     

http://dx.doi.org/10.1016/j.gsf.2016.07.001

In the tandem planet formation regime,planets form at two distinct sites where solid particles are densely accumulated due to the on/off state of the magnetorotational instability(MRI).We found that tandem planet formation can reproduce the solid component distribution of the Solar System and tends to produce a smaller number of large planets through continuous pebble flow into the planet formation sites.In the present paper,we investigate the dependence of tandem planet formation on the vertical magnetic field of the protoplanetary disk.We calculated two cases of B_Z 3.4 × 10~(-3) G and B_Z = 3.4 × 10~(-5) G at 100 AU as well as the canonical case of B_Z = 3.4 × 10~(-4) G.We found that tandem planet formation holds up well in the case of the strong magnetic field(B_Z 3.4 × 10~(-3) G).On the other hand,in the case of a weak magnetic field(B_Z= 3.4 × 10~(-5) G) at 100 AU,a new regime of planetary growth is realized:the planets grow independently at different places in the dispersed area of the MRl-suppressed region of r-8-30 AU at a lower accretion rate of M 10~(-7.4)M_⊙yr~(-1).We call this the "dispersed planet formation" regime.This may lead to a system with a larger number of smaller planets that gain high eccentricity through mutual collisions.     

<Emphasis Type="Italic">P</Emphasis>–<Emphasis Type="Italic">V</Emphasis>–<Emphasis Type="Italic">T</Emphasis> equation of state of CaAl<Subscript>4</Subscript>Si<Subscript>2</Subscript>O<Subscript>11</Subscript> CAS phase

The thermoelastic parameters of the CAS phase (CaAl₄Si₂O₁₁) were examined by in situ high-pressure (up to 23.7 GPa) and high-temperature (up to 2,100 K) synchrotron X-ray diffraction, using a Kawai-type multi-anvil press. P–V data at room temperature fitted to a third-order Birch–Murnaghan equation of state (BM EOS) yielded: V _0,300 = 324.2 ± 0.2 Å³ and K _0,300 = 164 ± 6 GPa for K′ _0,300 = 6.2 ± 0.8. With K′ _0,300 fixed to 4.0, we obtained: V _0,300 = 324.0 ± 0.1 Å³ and K _0,300 = 180 ± 1 GPa. Fitting our P–V–T data with a modified high-temperature BM EOS, we obtained: V _0,300 = 324.2 ± 0.1 Å³, K _0,300 = 171 ± 5 GPa, K′ _0,300 = 5.1 ± 0.6 (?K _0,T /?T) _P  = ?0.023 ± 0.006 GPa K^?1, and α_0,T  = 3.09 ± 0.25 × 10^?5 K^?1. Using the equation of state parameters of the CAS phase determined in the present study, we calculated a density profile of a hypothetical continental crust that would contain ~10 vol% of CaAl₄Si₂O₁₁. Because of the higher density compared with the coexisting minerals, the CAS phase is expected to be a plunging agent for continental crust subducted in the transition zone. On the other hand, because of the lower density compared with lower mantle minerals, the CAS phase is expected to remain buoyant in the lowermost part of the transition zone.     

Toxinogenicity and cytotoxicity of <Emphasis Type="Italic">Alternaria,Aspergillus</Emphasis> and <Emphasis Type="Italic">Penicillium</Emphasis> moulds isolated from working environments

There is currently limited research available on the secondary metabolites of moulds in workplaces. The aim of this study was to determine the mould contamination in museums (N = 4), composting plants (N = 4) and tanneries (N = 4) and the secondary metabolite profiles of Alternaria, Aspergillus and Penicillium isolates from these workplaces. Alternaria, Aspergillus and Penicillium species were identified using the ITS1/2 sequence of the rDNA region. Mould metabolites were quantitatively analysed on standard laboratory medium and mineral medium containing materials specific to each workplace using liquid chromatography-mass spectrometry. We also examined the cytotoxicity of the moulds using MTT (3-(4,5-dimethylthiazolyl-2)-2,5-diphenyltetrazolium bromide) assays. Air microbiological contamination analyses showed a number of microorganisms, ranging from 2.4 × 10³ CFU m^?3 (composting plants) to 6.8 × 10⁴ CFU m^?3 (tanneries). We identified high percentages of Alternaria, Aspergillus and Penicillium moulds (air 57–59%, surfaces 10–65%) in all workplaces. The following moulds were the most cytotoxic (>90%): Alternaria alternata, A. limoniasperae, Aspergillus flavus, Penicillium biourgeianum, P. commune and P. spinulosum. The same mould species isolated from different working environments exhibited varying toxigenic and cytotoxic properties. Modifying the culture medium to simulate environmental conditions most often resulted in the inhibition of secondary metabolite production. Moulds isolated from the working environments produced the following mycotoxins (ng g^?1): chanoclavines (0.28–204), cyclopiazonic acid (27.1–1045), fumigaclavines (0.33–10,640,000), meleagrin (0.57–13,393), roquefortins (0.01–16,660), rugulovasines (112–220), viridicatin (0.12–957), viridicatol (4.23–2753) and quinocitrinines (0.07–1104), which may have a negative impact on human health.     

Seasonal disparity in the co-occurrence of arsenic and fluoride in the aquifers of the Brahmaputra flood plains,Northeast India

Arsenic (As) and fluoride (F^?) in groundwater are increasing global water quality and public health concerns. The present study provides a deeper understanding of the impact of seasonal change on the co-occurrence of As and F^?, as both contaminants vary with climatic patterns. Groundwater samples were collected in pre- and post-monsoon seasons (n = 40 in each season) from the Brahmaputra flood plains (BFP) in northeast India to study the effect of season on As and F^? levels. Weathering is a key hydrogeochemical process in the BFP and both silicate and carbonate weathering are enhanced in the post-monsoon season. The increase in carbonate weathering is linked to an elevation in pH during the post-monsoon season. A Piper diagram revealed that bicarbonate-type water, with Na⁺, K⁺, Ca²⁺, and Mg²⁺ cations, is common in both seasons. Correlation between Cl^? and NO₃ ^? (r = 0.74, p = 0.01) in the post-monsoon indicates mobilization of anthropogenic deposits during the rainy season. As was within the 10 µg L^?1 WHO limit for drinking water and F^? was under the 1.5 mg L^?1 limit. A negative correlation between oxidation reduction potential and groundwater As in both seasons (r = ?0.26 and ?0.49, respectively, for pre-monsoon and post-monsoon, p = 0.05) indicates enhanced As levels due to prevailing reducing conditions. Reductive hydrolysis of Fe (hydr)oxides appears to be the predominant process of As release, consistent with a positive correlation between As and Fe in both seasons (r = 0.75 and 0.73 for pre- and post-monsoon seasons, respectively, at p = 0.01). Principal component analysis and hierarchical cluster analysis revealed grouping of Fe and As in both seasons. F^? and sulfate were also clustered during the pre-monsoon season, which could be due to their similar interactions with Fe (hydr)oxides. Higher As levels in the post-monsoon appears driven by the influx of water into the aquifer, which drives out oxygen and creates a more reducing condition suitable for reductive dissolution of Fe (hydr)oxides. An increase in pH promotes desorption of As oxyanions AsO₄ ^3? (arsenate) and AsO₃ ^3? (arsenite) from Fe (hydr)oxide surfaces. Fluoride appears mainly released from F^?-bearing minerals, but Fe (hydr)oxides can be a secondary source of F^?, as suggested by the positive correlation between As and F^? in the pre-monsoon season.     

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 β.     

Sulfate Reduction and Sulfur Cycles at Two Seagrass Beds Inhabited by Cold Affinity <Emphasis Type="Italic">Zostera marina</Emphasis> and Warm Affinity <Emphasis Type="Italic">Halophila nipponica</Emphasis> in Temperate Coastal Waters

To evaluate the impact of invading seagrass on biogeochemical processes associated with sulfur cycles, we investigated the geochemical properties and sulfate reduction rates (SRRs) in sediments inhabited by invasive warm affinity Halophila nipponica and indigenous cold affinity Zostera marina. A more positive relationship between SRR and below-ground biomass (BGB) was observed at the H. nipponica bed (SRR = 0.6809 × BGB ? 4.3162, r ² = 0.9878, p = 0.0006) than at the Z. marina bed (SRR = 0.3470 × BGB ? 4.0341, r ² = 0.7082, p = 0.0357). These results suggested that SR was more stimulated by the dissolved organic carbon (DOC) exuded from the roots of H. nipponica than by the DOC released from the roots of Z. marina. Despite the enhanced SR in spring-summer, the relatively lower proportion (average, 20%) of acid-volatile sulfur (AVS) in total reduced sulfur and the strong correlation between total oxalate-extractable Fe (Fe_(oxal)) and chromium-reducible sulfur (CRS = 0.2321 × total Fe_(oxal) + 1.8180, r ² = 0.3344, p = 0.0076) in the sediments suggested the rapid re-oxidation of sulfide and precipitation of sulfide with Fe. The turnover rate of the AVS at the H. nipponica bed (0.13 day^?1) was 2.5 times lower than that at the Z. marina bed (0.33 day^?1). Together with lower AVS turnover, the stronger correlation of SRR to BGB in the H. nipponica bed suggests that the extension of H. nipponica resulting from the warming of seawater might provoke more sulfide accumulation in coastal sediments.     

Evaluating the susceptibility of landslide landforms in Japan using slope stability analysis: a case study of the 2016 Kumamoto earthquake

This study analyzed 267 landslide landforms (LLs) in the Kumamoto area of Japan from the database of about 0.4 million LLs for the whole of Japan identified from aerial photos by the National Research Institute for Earth Science and Disaster Resilience of Japan (NIED). Each LL in the inventory is composed of a scarp and a moving mass. Since landslides are prone to reactivation, it is important to evaluate the sliding-recurrence susceptibility of LLs. One possible approach to evaluate the susceptibility of LLs is slope stability analysis. A previous study found a good correlation (R ² = 0.99) between the safety factor (F _s ) and slope angle (α) of F _s  = 17.3α ^?0.843. We applied the equation to the analysis of F _s for 267 LLs in the area affected by the 2016 Kumamoto earthquake (M _j  = 7.3). The F _s was calculated for the following three cases of failure: scarps only, moving mass only, and scarps and moving mass together. Verification with the 2016 Kumamoto earthquake event shows that the most appropriate method for the evaluation of LLs is to consider the failure of scarps and moving mass together. In addition, by analyzing the relationship between the factors of slope of entire landslide and slope of scarp for LLs and comparing the results with the Aso-ohashi landslide, the largest landslide caused by the 2016 Kumamoto earthquake, we also found that morphometric analysis of LLs is useful for forecasting the travel distance of future landslides.     

Numerical analysis under seismic loads of soils improvement with floating stone columns

Geotechnical Engineering has developed many methods for soil improvement so far. One of these methods is the stone column method. The structure of a stone column generally refers to partial change of suitable subsurface ground through a vertical column, poor stone layers which are completely pressed. In general terms, to improve bearing capacity of problematic soft and loose soil is implemented for the resolution of many problems such as consolidation and grounding problems, to ensure filling and splitting slope stability and liquefaction that results from a dynamic load such as earthquake. In this study, stone columns method is preferred as an improvement method, and especially load transfer mechanisms and bearing capacity of floating stone column are focused. The soil model, 32 m in width and 8 m in depth, used in this study is made through Plaxis 2D finite element program. The clay having 5° internal friction angle with different cohesion coefficients (c 10, c 15, c 20 kN/m²) are used in models. In addition, stone columns used for soil improvement are modeled at different internal friction angles (? 35°, ? 40°, ? 45°) and in different s/D ranges (s/D 2, s/D 3), stone column depths (B, 2B, 3B) and diameters (D 600 mm, D 800 mm, D 1000 mm). In the study, maximum acceleration (a _max = 1.785 m/s²) was used in order to determine the seismic coefficient used. In these soil models, as maximum acceleration, maximum east–west directional acceleration value of Van Muradiye earthquake that took place in October 23, 2011 was used. As a result, it was determined that the stone column increased the bearing capacity of the soil. In addition, it is observed that the bearing capacity of soft clay soil which has been improved through stone column with both static and earthquake load effect increases as a result of increase in the diameter and depth of the stone column and decreases as a result of the increase in the ranges of stone column. In the conducted study, the bearing capacity of the soil models, which were improved with stone column without earthquake force effect, was calculated as 1.01–3.5 times more on the average, compared to the bearing capacity of the soil models without stone column. On the other hand, the bearing capacity of the soil models with stone columns, which are under the effect of earthquake force, was calculated as 1.02–3.7 times more compared to the bearing capacity of the soil models without stone column.     

Static elasticity of cordierite I: Effect of heavy ion irradiation on the compressibility of hydrous cordierite

The effect of ion beam irradiations on the elastic properties of hydrous cordierite was investigated by means of Raman and X-ray diffraction experiments. Oriented single crystals were exposed to swift heavy ions (Au, Bi) of various specific energies (10.0–11.1 MeV/u and 80 MeV/u), applying fluences up to 5 × 10¹³ ions/cm². The determination of unit-cell constants yields a volume strain of 3.4 × 10^?3 up to the maximum fluence, which corresponds to a compression of non-irradiated cordierite at ~480 ± 10 MPa. The unit-cell contraction is anisotropic (e ₁ = 1.4 ± 0.1 × 10^?3, e ₂ = 1.5 ± 0.1 × 10^?3, and e ₃ = 7 ± 1 × 10^?4) with the c-axis to shrink only half as much as the axes within the ab-plane. The lattice elasticity for irradiated cordierite (? = 1 × 10¹² ions/cm²) was determined from single-crystal XRD measurements in the diamond anvil cell. The fitted third-order Birch–Murnaghan equation-of-state parameters of irradiated cordierite (V ₀ = 1548.41 ± 0.16 Å³, K ₀ = 117.1 ± 1.1 GPa, ?K/?P = ?0.6 ± 0.3) reveal a 10–11 % higher compressibility compared to non-irradiated cordierite. While the higher compressibility is attributed to the previously reported irradiation-induced loss of extra-framework H₂O, the anomalous elasticity as expressed by elastic softening (β _a ^?1 , β _b ^?1 , β _c ^?1  = 397 ± 9, 395 ± 28, 308 ± 11 GPa, ?(β ^?1)/?P = ?4.5 ± 2.7, ?6.6 ± 8.4, ?5.4 ± 3.0) appears to be related to the framework stability and to be independent of the water content in the channels and thus of the ion beam exposure.     

Biosorption optimization and equilibrium isotherm of industrial dye compounds in novel strains of microscopic fungi

The aim of this study was to evaluate the biosorption capacity of selected strains of microscopic fungi. We optimized the biosorption process and used the Freundlich isotherm for three strains: H. haematococca BwIII43, K37 and T. harzianum BsIII33 to describe the biosorption equilibrium of anthraquinone dye, Alizarin Blue Black B (ABBB) and alkali lignin (AL). In optimal conditions (1 g of mycelium biomass, pH = 7.0, 28 °C) for ABBB and AL sorption, the live biomass of H. haematococca BwIII43 was characterized by a higher sorption capacity, amounting to 247.47 and 161.00 mg g^?1, respectively. The highest sorption properties toward anthraquinone dye (K _F = 19.96 mg g^?1) were shown for the biomass of H. haematococca K37. In the presence of alkali lignin, the highest sorption capacity and bond strength exhibited the biomass of H. haematococca BwIII43 (K _F = 28.20 mg g^?1, n = 3.46). Effective decolorization of ABBB and AL by the selected strains of microscopic fungi indicated that the biosorption process additionally enhanced the removal of color compounds from the solution.     

Destructive near-fault strong ground motion from the 2016 Kumamoto prefecture,Japan, M7.3 earthquake

The sequence of the 2016 Kumamoto, Japan, earthquake, which included an initial M6.5 foreshock on April 14, followed by a larger M7.3 mainshock on April 16, and subsequently occurred high aftershock activity, caused significant damage in Kumamoto and neighboring regions. The near-field strong motion record by strong motion network (K-NET and KiK-net) and the intensity meter network demonstrated clearly the characteristics of the strong ground motion developed by the shallow (H = 12 km), inland earthquake comprising short-time duration (<15–20 s) but large (>1G) ground accelerations. The velocity response spectra of the near-fault motion at Mashiki and Nishihara showed large levels (>300–550 cm/s) in the short-period range (T = 1–2 s), several times larger than that of the near-field record of the destructive 1995 Kobe earthquake (M7.3) and that of the 2004 Mid-Niigata earthquake (M6.8). This period corresponds to the collapse vulnerability of Japanese wooden-frame houses, and is the major cause of severe damage during the Kumamoto earthquake. The response spectra also showed extremely large levels (>240–340 cm/s) in the long-period (T > 3 s) band, which is potentially disastrous for high-rise buildings, large oil storage tanks, etc. to have longer resonant period. Such long-period motion was, for the most parts, developed by the static displacement of the fault movement rather than by the seismic waves radiating from the source fault. Thus, the extreme near-fault long-period motion was hazardous only close to the fault but it attenuated very rapidly away from the fault.