Trend analysis of evapotranspiration over India: Observed from long-term satellite measurements

Owing to the lack of consistent spatial time series data on actual evapotranspiration (ET), very few studies have been conducted on the long-term trend and variability in ET at a national scale over the Indian subcontinent. The present study uses biome specific ET data derived from NOAA satellite’s advanced very high resolution radiometer to investigate the trends and variability in ET over India from 1983 to 2006. Trend analysis using the non-parametric Mann–Kendall test showed that the domain average ET decreased during the period at a rate of \(0.22\,\hbox {mm year}^{-1}\). A strong decreasing trend (\(m = -1.75\, \hbox {mm year}^{-1}\), \(F = 17.41\), \(P\) 0.01) was observed in forest regions. Seasonal analyses indicated a decreasing trend during southwest summer monsoon (\(m= -0.320\, \hbox {mm season}^{-1}\,\hbox {year}^{-1})\) and post-monsoon period (\(m= -0.188\, \hbox {mm season}^{-1 }\,\hbox {year}^{-1})\). In contrast, an increasing trend was observed during northeast winter monsoon (\(m = 0.156 \,\hbox {mm season}^{-1 }\,\hbox {year}^{-1})\) and pre-monsoon (\(m = 0.068\, \hbox {mm season}^{-1 }\,\hbox {year}^{-1})\) periods. Despite an overall net decline in the country, a considerable increase ( \(4 \,\hbox {mm year}^{-1}\)) was observed over arid and semi-arid regions. Grid level correlation with various climatic parameters exhibited a strong positive correlation (\(r \!>\!0.5\)) of ET with soil moisture and precipitation over semi-arid and arid regions, whereas a negative correlation (\(r\) \(-0.5\)) occurred with temperature and insolation in dry regions of western India. The results of this analysis are useful for understanding regional ET dynamics and its relationship with various climatic parameters over India. Future studies on the effects of ET changes on the hydrological cycle, carbon cycle, and energy partitioning are needed to account for the feedbacks to the climate.     

Timescales for mechanisms for the dynamical evolution of open star clusters

We have obtained the stellar velocity dispersion in three mutually perpendicular directions in the halos and cores of clusters as a function of time for several non-stationary open-cluster models. During the dynamical evolution of the open-cluster models, the velocity dispersions undergo oscillations that do not decay during 5–10 violent-relaxation timescales, τ _vr . We estimated the time for synchronization of the rotation of the open-cluster models and their motion around the center of the Galaxy, t _s , which, depending on the model parameters, is t _s ? (5–27)τ _vr . Synchronization mechanisms for the models are discussed. The disruption of the systems in the force field of the Galaxy is strongly affected by tidal friction. We have also estimated the time for the formation of a spherical stellar-velocity distribution in the cluster models, t _σ ? (6 ? 25)τ _vr . The impact of instability in the stellar motions in a cluster on the formation of a spherical velocity distribution in the open-cluster models is discussed. We have noted a tendency for a weakening of the dependence of the coarse phase density of the cluster on small initial perturbations of the stellar phase coordinates in the model cluster cores for times about five times longer than the violent-relaxation time.     

H-bonding scheme and cation partitioning in axinite: a single-crystal neutron diffraction and Mössbauer spectroscopic study

The crystal chemistry of a ferroaxinite from Colebrook Hill, Rosebery district, Tasmania, Australia, was investigated by electron microprobe analysis in wavelength-dispersive mode, inductively coupled plasma–atomic emission spectroscopy (ICP–AES), ⁵⁷Fe Mössbauer spectroscopy and single-crystal neutron diffraction at 293 K. The chemical formula obtained on the basis of the ICP–AES data is the following: \( ^{X1,X2} {\text{Ca}}_{4.03} \,^{Y} \left( {{\text{Mn}}_{0.42} {\text{Mg}}_{0.23} {\text{Fe}}^{2 + }_{1.39} } \right)_{\varSigma 2.04} \,^{Z1,Z2} \left( {{\text{Fe}}^{3 + }_{0.15} {\text{Al}}_{3.55} {\text{Ti}}_{0.12} } \right)_{\varSigma 3.82} \,^{T1,T2,T3,T4} \left( {{\text{Ti}}_{0.03} {\text{Si}}_{7.97} } \right)_{\varSigma 8} \,^{T5} {\text{B}}_{1.96} {\text{O}}_{30} \left( {\text{OH}} \right)_{2.18} \). The ⁵⁷Fe Mössbauer spectrum shows unambiguously the occurrence of Fe²⁺ and Fe³⁺ in octahedral coordination only, with Fe²⁺/Fe³⁺ = 9:1. The neutron structure refinement provides a structure model in general agreement with the previous experimental findings: the tetrahedral T1, T2, T3 and T4 sites are fully occupied by Si, whereas the T5 site is fully occupied by B, with no evidence of Si at the T5, or Al or Fe³⁺ at the T1–T5 sites. The structural and chemical data of this study suggest that the amount of B in ferroaxinite is that expected from the ideal stoichiometry: 2 a.p.f.u. (for 32 O). The atomic distribution among the X1, X2, Y, Z1 and Z2 sites obtained by neutron structure refinement is in good agreement with that based on the ICP–AES data. For the first time, an unambiguous localization of the H site is obtained, which forms a hydroxyl group with the oxygen atom at the O16 site as donor. The H-bonding scheme in axinite structure is now fully described: the O16–H distance (corrected for riding motion effect) is 0.991(1) Å and an asymmetric bifurcated bonding configuration occurs, with O5 and O13 as acceptors [i.e. with O16···O5 = 3.096(1) Å, H···O5 = 2.450(1) Å and O16–H···O5 = 123.9(1)°; O16···O13 = 2.777(1) Å, H···O13 = 1.914(1) Å and O16–H···O13 = 146.9(1)°].     

Stable carbon and oxygen isotope study on benthic foraminifera: Implication for microhabitat preferences and interspecies correlation

Stable isotopes of benthic foraminifera have widely been applied in micropalaeontological research to understand vital effects in foraminifera. Isotopic fractionations are mainly controlled by ontogeny, bottom/pore water chemistry, habitat preference, kinetic effect and respiration. Discontinuous abundance of a species for isotopic analysis has forced us to select multiple species from down-core samples. Thus standardisation factors are required to convert isotopic values of one species with respect to other species. The present study is pursued on isotopic values of different pairs of benthic foraminifera from the Krishna–Godavari basin and Peru offshore to understand habitat-wise isotopic variation and estimation of isotopic correction factors for the paired species (Cibicides wuellerstorfi–Bulimina marginata, Ammonia spp.–Loxostomum amygdalaeformis and Bolivina seminuda–Nonionella auris). Infaunal species (B. marginata, Ammonia spp. and N. auris) show a lighter carbon isotopic excursion with respect to the epifaunal to shallow infaunal forms (C. wuellerstorfi, L. amygdalaeformis and B. seminuda). These lighter \(\updelta ^{13}\) \(\hbox {C}\) values are related to utilisation of \(\hbox {CO}_{2}\) produced by anaerobic remineralisation of organic matter. However, enrichment of \(\updelta ^{18}\) \(\hbox {O}\) for the deeper microhabitat (bearing lower pH and decreased \({\hbox {CO}_{3}}^{2-})\) is only recorded in case of B. marginata. It is reverse in case of N. auris and related to utilisation of respiratory \(\hbox {CO}_{2}\) and internal dissolve inorganic carbon pool. Estimation of interspecies isotopic correction factors for the species pairs (\(\updelta ^{13}\) \(\hbox {C}\) of C. wuellerstorfi–B. marginata, L. amygdalaeformis–Ammonia spp., N. auris–B. seminuda) and \(\updelta ^{18}\) \(\hbox {O}\) of C. wuellerstorfi–B. marginata are statistically reliable and may be used in palaeoecological studies.     

Fractal dimension of pore space in carbonate samples from Tushka area (Egypt)

To investigate inhomogeneous and porous structures in nature, the concept of fractal dimension was established. This paper briefly introduces the definition and measurement methods of fractal dimension. Three different methods including mercury injection capillary pressure (MICP), nuclear magnetic resonance (NMR), and nitrogen adsorption (BET) were applied to determine the fractal dimensions of the pore space of eight carbonate rock samples taken from West Tushka area, Egypt. In the case of fractal behavior, the capillary pressure P _c and cumulative fraction V _c resulting from MICP are linearly related with a slope of D-3 in a double logarithmic plot with D being the value of fractal dimension. For NMR, the cumulative intensity fraction V _c and relaxation time T ₂ show a linear relation with a slope of 3-D in a double logarithmic plot. Fractal dimension can also be determined by the specific surface area S _por derived from nitrogen adsorption measurements and the effective hydraulic radius. The fractal dimension D shows a linear relation with the logarithm of S _por . The fractal dimension is also used in models of permeability prediction. To consider a more comprehensive data set, another 34 carbonate samples taken from the same study area were integrated in the discussion on BET method and permeability prediction. Most of the 42 rock samples show a good agreement between measured permeability and predicted permeability if the mean surface fractal dimension for each facies is used.     

Diffusive fractionation of volatiles and their isotopes during bubble growth in magmas

Bubbles grow in decompressing magmas by simple expansion and by diffusive supply of volatiles to the bubble/melt interface. The latter phenomenon is of significant geochemical interest because diffusion can fractionate elements and isotopes (or isotopologues) of dissolved components. This raises the possibility that the character of volatile components in bubbles may not reflect that of volatiles dissolved in the host melt over the lifetime of a bubble—even in the absence of equilibrium vapor/melt isotopic fractionation. Recent experiments have confirmed the existence of an isotope mass effect on diffusion of the volatile element Cl in silicate melt [Fortin et al. (Isotopic fractionation of chlorine during chemical diffusion in a dacitic melt and its implications for isotope behavior during bubble growth (abstract), 2016 Fall AGU Meeting, 2016)], so there is a clear need to understand the efficacy of diffusive fractionation during bubble growth. In this study, numerical models of diffusion and mass redistribution during bubble growth were implemented for both “passive” volatiles—those whose concentrations are generally well below saturation levels—and “active” volatiles such as CO₂ and H₂O, whose elevated concentrations and limited solubilities are the cause of bubble nucleation and growth. Both diffusive and convective bubble-growth scenarios were explored. The magnitude of the isotope mass effect on passive volatiles partitioned into bubbles growing at a constant rate R in a static system depends upon R/D _L, K _d and D _H/D _L (K _d = bubble/melt partition coefficient; D _H/D _L = diffusivity ratio of the heavy and light isotopes). During convective bubble growth, the presence of a discrete (physical) melt boundary layer against the growing bubble (of width x _BL) simplifies outcomes because it leads to the quick onset of steady-state fractionation during growth, the magnitude of which depends mainly upon R?x _BL/D _L and D _H/D _L (bubble/melt fractionation is maximized at R?x _BL/D _L ≈0.1). Constant R is unrealistic for most real systems, so other scenarios were explored by including the solubility and EOS of an “active” volatile (e.g., CO₂) in the numerical simulations. For plausible decompression paths, R increases exponentially with time—leading, potentially, to larger isotopic fractionation of species partitioned into the growing bubble. For volatile species whose isotope mass effects on diffusion have been measured (Cl, Li), predicted isotope fractionation in the exsolved vapor can be as large as ?4‰ for Cl and ?25‰ for Li.     

Spots and the activity of M dwarfs from observations with the Kepler Space Telescope

A method for estimating the spottedness parameter S (the spotted area as a fraction of the surface of an active star) proposed earlier is applied to an analysis of activity in 1570 M dwarf stars. The analysis is based on observational material obtained with the Kepler Space Telescope, as well as data on the fluxes of the studied objects in the near and far ultraviolet (NUV and FUV) based on data from the GALEX space telescope. The variations of S with the ages of the stars are studied (four groups with different ages are distinguished), as well as variations of S with their rotational periods. A diagram characterizing the relationship between S and the Rossby number Ro resembles the classical dependence of the X-ray luminosities of active stars on Ro, and a saturation regime is attained at the same value, Ro = 0.13. Moreover, objects with ages of more than 100 million years do not form a single sequence (and stars older than 900 million years possess surface spottednesses of order 1%). The S?Ro dependence obtained could expand possibilities for analyzing the dependence of the X-ray luminosities of active stars on their Rossby numbers, and could also be applied to refine parameters characterizing the action of dynamo mechanisms, such as the dynamo number N _D . A comparison of the GALEX NUV and FUV brightness estimates with the activity parameters of the stars suggests that younger, more rapidly rotating active stars are brighter in the NUV, and that the FUV flux grows and the difference of the FUV and NUV brightnesses decreases with increasing spottedness S.     

Water diffusion in Mount Changbai peralkaline rhyolitic melt

Diffusion couple experiments with wet half (up to 4.6 wt%) and dry half were carried out at 789–1,516 K and 0.47–1.42 GPa to investigate water diffusion in a peralkaline rhyolitic melt with major oxide concentrations matching Mount Changbai rhyolite. Combining data from this work and a related study, total water diffusivity in peralkaline rhyolitic melt can be expressed as:

$ D_{{{\text{H}}_{ 2} {\text{O}}_{\text{t}} }} = D_{{{\text{H}}_{ 2} {\text{O}}_{\text{m}} }} \left( {1 - \frac{0.5 - X}{{\sqrt {[4\exp (3110/T - 1.876) - 1](X - X^{2} ) + 0.25} }}} \right), $

$ {\text{with}}\;D_{{{\text{H}}_{ 2} {\text{O}}_{\text{m}} }} = \exp \left[ { - 1 2. 7 8 9- \frac{13939}{T} - 1229.6\frac{P}{T} + ( - 27.867 + \frac{60559}{T})X} \right], $

where D is in m² s^?1, T is the temperature in K, P is the pressure in GPa, and X is the mole fraction of water and calculated as X = (C/18.015)/(C/18.015 + (100 ? C)/33.14), where C is water content in wt%. We recommend this equation in modeling bubble growth and volcanic eruption dynamics in peralkaline rhyolitic eruptions, such as the ~1,000-ad eruption of Mount Changbai in North East China. Water diffusivities in peralkaline and metaluminous rhyolitic melts are comparable within a factor of 2, in contrast with the 1.0–2.6 orders of magnitude difference in viscosities. The decoupling of diffusivity of neutral molecular species from melt viscosity, i.e., the deviation from the inversely proportional relationship predicted by the Stokes–Einstein equation, might be attributed to the small size of H₂O molecules. With distinct viscosities but similar diffusivity, bubble growth controlled by diffusion in peralkaline and metaluminous rhyolitic melts follows similar parabolic curves. However, at low confining pressure or low water content, viscosity plays a larger role and bubble growth rate in peralkaline rhyolitic melt is much faster than that in metaluminous rhyolite.

     

Shear wave velocity as function of cone penetration resistance and grain size for Holocene-age uncemented soils: a new perspective

For feasibility studies and preliminary design estimates, field measurements of shear wave velocity, V _s, may not be economically adequate and empirical correlations between V _s and more available penetration measurements such as cone penetration test, CPT, data turn out to be potentially valuable at least for initial evaluation of the small-strain stiffness of soils. These types of correlations between geophysical (Vs) and geotechnical (N-SPT, q _c-CPT) measurements are also of utmost importance where a great precision in the calculation of the deposit response is required such as in liquefaction evaluation or earthquake ground response analyses. In this study, the stress-normalized shear wave velocity V _s1 (in m/s) is defined as statistical functions of the normalized dimensionless resistance, Q _tn-CPT, and the mean effective diameter, D ₅₀ (in mm), using a data set of different uncemented soils of Holocene age accumulated at various sites in North America, Europe, and Asia. The V _s1–Q _tn data exhibit different trends with respect to grain sizes. For soils with mean grain size (D ₅₀) < 0.2 mm, the V _s1/Q _tn ^0.25 ratio undergoes a significant reduction with the increase in D ₅₀ of the soil. This trend is completely reversed with further increase in D ₅₀ (D ₅₀ > 0.2 mm). These results corroborate earlier results that stressed the use of different CPT-based correlations with different soil types, and those emphasized the need to impose particle-size limits on the validity of the majority of available correlations.     

The morphological characteristics of gully systems and watersheds in Dry-Hot Valley,SW China

Gully systems and watersheds are geomorphic units with clear boundaries that are relatively independent of basin landscapes and play an important role in natural geography. In order to explore the morphological characteristics of gully systems and watersheds in the Dry-Hot Valley [South West (SW) China], gullies are interpreted from online Google images with high resolution and watersheds are extracted from digital elevation model at a scale of 1:50,000. The results show that: (1) There are 17,382 gullies (with a total area of 1141.66 km²) and 42 watersheds in the study area. (2) The average gully density of the study area (D) is 4.29 km/km², gully frequency (F) is 14.39 gullies/km², the branching ratio (B) is 5.13, the length ratio (L) is 3.12, and the coefficient of the main and tributary gullies (M) is 0.06. The degree of gully erosion is strong to extremely strong, the main development intensity of gully erosion ranges from intense to moderate, and the type of gully system is tributary. (3) The watershed areas (A) are between 0.39 and 96.43 km², the relief ratio (R) is from 0.10 to 0.19, the circularity ratio (C) is from 0.30 to 0.83, the texture ratio (T) is from 0.82 to 39.35, and the dominant geomorphological texture type is fine. (4) There is a quantitative relationship between F and D:F?=?0.624D² (R?=0.84) and T is closely related to D, F, M (R²?>?0.7). A, R and C are related to M (R²?>?0.5). The development of gully systems is the result of coupling effects between multiple factors. In this area, the degree of erosion and the condition of the main and tributary gullies can be controlled by the degree of topographic breakage in the watershed, which provides some theoretical basis for the evaluation of gully erosion by the latter. In addition, the scale, relief, and shape have a significant impact on the locations of the main and tributary gullies. For tributary gullies, attention should be paid to the interception and control of runoff and sediment in the small confluence branches in order to prevent gully expansion and head advance. These features can inform the development of targeted measures for the control of soil erosion.     

Experimental calibration of Forsterite–Anorthite–Ca-Tschermak–Enstatite (FACE) geobarometer for mantle peridotites

The crystallization of plagioclase-bearing assemblages in mantle rocks is witness of mantle exhumation at shallow depth. Previous experimental works on peridotites have found systematic compositional variations in coexisting minerals at decreasing pressure within the plagioclase stability field. In this experimental study we present new constraints on the stability of plagioclase as a function of different Na₂O/CaO bulk ratios, and we present a new geobarometer for mantle rocks. Experiments have been performed in a single-stage piston cylinder at 5–10 kbar, 1050–1150?°C at nominally anhydrous conditions using seeded gels of peridotite compositions (Na₂O/CaO?=?0.08–0.13; X _Cr = Cr/(Cr?+?Al)?=?0.07–0.10) as starting materials. As expected, the increase of the bulk Na₂O/CaO ratio extends the plagioclase stability to higher pressure; in the studied high-Na fertile lherzolite (HNa-FLZ), the plagioclase-spinel transition occurs at 1100?°C between 9 and 10 kbar; in a fertile lherzolite (FLZ) with Na₂O/CaO?=?0.08, it occurs between 8 and 9 kbar at 1100?°C. This study provides, together with previous experimental results, a consistent database, covering a wide range of P–T conditions (3–9 kbar, 1000–1150?°C) and variable bulk compositions to be used to define and calibrate a geobarometer for plagioclase-bearing mantle rocks. The pressure sensitive equilibrium:

$$\mathop {{\text{M}}{{\text{g}}_{\text{2}}}{\text{Si}}{{\text{O}}_{\text{4}}}^{{\text{Ol}}}}\limits_{{\text{Forsterite}}} +\mathop {{\text{CaA}}{{\text{l}}_{\text{2}}}{\text{S}}{{\text{i}}_{\text{2}}}{{\text{O}}_{\text{8}}}^{{\text{Pl}}}}\limits_{{\text{Anorthite}}~} =\mathop {{\text{CaA}}{{\text{l}}_{\text{2}}}{\text{Si}}{{\text{O}}_{\text{6}}}^{{\text{Cpx}}}}\limits_{{\text{Ca-Tschermak}}} +{\text{ }}\mathop {{\text{M}}{{\text{g}}_{\text{2}}}{\text{S}}{{\text{i}}_{\text{2}}}{{\text{O}}_{\text{6}}}^{{\text{Opx}}}}\limits_{{\text{Enstatite}}} ,$$

has been empirically calibrated by least squares regression analysis of experimental data combined with Monte Carlo simulation. The result of the fit gives the following equation:

$$P=7.2( \pm 2.9)+0.0078( \pm 0.0021)T{\text{ }}+0.0022( \pm 0.0001)T{\text{ }}\ln K,$$

$${R^2}=0.93,$$

where P is expressed in kbar and T in kelvin. K is the equilibrium constant K?=?a _CaTs × a _en/a _an × a _fo, where a _CaTs, a _en, a _an and a _fo are the activities of Ca-Tschermak in clinopyroxene, enstatite in orthopyroxene, anorthite in plagioclase and forsterite in olivine. The proposed geobarometer for plagioclase peridotites, coupled to detailed microstructural and mineral chemistry investigations, represents a valuable tool to track the exhumation of the lithospheric mantle at extensional environments.

     

The influence of materials on the breaching process of natural dams

In this study, a series of natural dam overtopping laboratory tests are reported. In these tests, the effect of seven different sediment mixtures on the breaching process was investigated. According to the test results, three stages of the breaching process of natural dams made of different materials were observed. Backward erosion was the primary cause for the incising slopes. The effects of backward erosion became stronger with the larger fines contents of the materials. With an increase in the median diameter (d ₅₀) of particles, the breaching time became longer. However, the peak discharge became smaller. With an increase in the fines contents (p), the median diameter of the particles and the void ratio were changed, which resulted in a decrease in the breaching time and an increase in the peak discharge. The breaching time and peak discharge were more sensitive to the median diameter than to the fines contents. The relation between breach width and depth was found to follow a logistic function \( W\kern0.5em =\kern0.5em \frac{\zeta }{1\kern0.5em +\kern0.5em {e}^{\left(-k\left(D\kern0.5em -\kern0.5em {D}_0\right)\right)}} \). The parameters ζ, k, and D ₀ are defined by a linear relationship with the median diameter and fines content. A breach of the side slope occurred as a tensile failure when the fines contents of the materials were large; otherwise, shear failure occurred. Furthermore, when the materials had fewer fines contents, the volume of the collapsed breach side slope became larger.     

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 use of ultrasonic frequencies to control the bloom formation,regrowth, and eco-toxicity in <Emphasis Type="Italic">Microcystis aeruginosa</Emphasis>

This study is focused on the use of ultrasound to disrupt Microcystis aeruginosa growth, with consideration for the gap between laboratory-scale experiment and field application. Laboratory-scale sonication systems with different frequencies (i.e., 20, 584, 869, and 1137 kHz) for 30 min at 10-min intervals were conducted to investigate their effectiveness at disrupting M. aeruginosa growth. The toxicological effect of sonicated M. aeruginosa including microcystin-LR was also evaluated using a Daphnia magna acute toxicity test. High frequencies, 869 and 1137 kHz, resulted in more than a 60% reduction of cells after 30 min of sonication. Low to middle frequencies, 20 and 584 kHz, only showed a 30% reduction of cells after 30 min of sonication. High frequencies also led to the inactivation of cell proliferation during M. aeruginosa regrowth, due to cellular destruction and finally cell death. However, the concentration of microcystin-LR and the potential adverse effects of M. aeruginosa on D. magna could not be controlled using ultrasonic frequencies. Therefore, our results suggest that ultrasonic frequencies between 869 and 1137 kHz are effective at controlling bloom formation in M. aeruginosa and the regrowth of M. aeruginosa after sonication, but not at controlling microcystin-LR concentrations and its adverse effects on D. magna. Consequently, a combined purification technology to reduce the cyanotoxins such as microcystin-LR, rather than ultrasonic frequency alone, will be needed to control M. aeruginosa growth and its toxicity levels in the aquatic environment.     

Laboratory parallel-beam transmission X-ray powder diffraction investigation of the thermal behavior of nitratine NaNO<Subscript>3</Subscript>: spontaneous strain and structure evolution

Present work provides in-situ structural data at a fine temperature scale from RT to the melting point of nitratine, NaNO₃. From the analysis of log e ₃₃ versus log t plots, it is possible to prove that an univocal indication on the R \( \overline{3} \) c (low temperature, LT) → R \( \overline{3} \) m (high temperature, HT) transition mechanism cannot be obtained because of the relevant role played by the arbitrary assumptions required for defining the c ₀ dependence from temperature of the HT phase. This is due to the occurrence of excess thermal expansion for the HT phase. A significantly better fit for an Ising-spin structural model over a non-Ising rigid-body one has been obtained for the LT phase. Moreover, the Ising model led to a smooth variation of the oxygen site x fractional coordinate throughout the transition. The structure of the HT polymorph has been successfully refined considering an oxygen site at x, 0, ½, with 50% occupancy. Such model was the only acceptable one from the crystal chemical point of view as the alternative model (oxygen site at x, y, z with 25% occupancy) led to unrealistically aplanar \( {\text{NO}}_{3}^{ - } \) groups.     

Distribution of REE,LILE, and HFSE between biotite,feldspar, and the melt in the granulite facies migmatite,nimnyr block,Aldan shield

The behavior of trace elements under conditions of partial melting of granitoid rocks has been studied. The element’s partition coefficients between minerals and the melt D_i^min/melt depends, in the first place, on the composition of the primary melt. In biotite the HREE D_i are a little below 1, while those of LREE, especially D_i for Ce, are 1–3 orders of magnitude less. This leads to an efficient differentiation of REEs in anatexic melts especially when biotite is the main mineral phase of restite. On the contrary, there are feldspars, the D_i of which cannot provide such a magnitude of differentiation. Unlike garnets and pyroxenes, whose stability in restite permits enrichment of anatexic melts produced in migmatization zones with Nb, Ti, and Cr, the presence of biotite in restite causes depletion of melts with those elements as well as with Rb. Feldspars, under conditions of their fractional crystallization or during differentiation of an anatexic melt, deplete the latter with Sr, Ba, and Rb, but enrich it with Nb, Ti, Cr, Y, Zr, and V.     

Heat capacity of hydrous trachybasalt from Mt Etna: comparison with CaAl<Subscript>2</Subscript>Si<Subscript>2</Subscript>O<Subscript>8</Subscript> (An)–CaMgSi<Subscript>2</Subscript>O<Subscript>6</Subscript> (Di) as basaltic proxy compositions

The specific heat capacity (C _p) of six variably hydrated (~3.5 wt% H₂O) iron-bearing Etna trachybasaltic glasses and liquids has been measured using differential scanning calorimetry from room temperature across the glass transition region. These data are compared to heat capacity measurements on thirteen melt compositions in the iron-free anorthite (An)–diopside (Di) system over a similar range of H₂O contents. These data extend considerably the published C _p measurements for hydrous melts and glasses. The results for the Etna trachybasalts show nonlinear variations in, both, the heat capacity of the glass at the onset of the glass transition (i.e., C _p ^g ) and the fully relaxed liquid (i.e., C _p ^l ) with increasing H₂O content. Similarly, the “configurational heat capacity” (i.e., C _p ^c  = C _p ^l  ? C _p ^g ) varies nonlinearly with H₂O content. The An–Di hydrous compositions investigated show similar trends, with C _p values varying as a function of melt composition and H₂O content. The results show that values in hydrous C _p ^g , C _p ^l and C _p ^c in the depolymerized glasses and liquids are substantially different from those observed for more polymerized hydrous albitic, leucogranitic, trachytic and phonolitic multicomponent compositions previously investigated. Polymerized melts have lower C _p ^l and C _p ^c and higher C _p ^g with respect to more depolymerized compositions. The covariation between C _p values and the degree of polymerization in glasses and melts is well described in terms of SM_hydrous and NBO/T _hydrous. Values of C _p ^c increase sharply with increasing depolymerization up to SM_hydrous ~ 30–35 mol% (NBO/T _hydrous ~ 0.5) and then stabilize to an almost constant value. The partial molar heat capacity of H₂O for both glasses (\( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{g}} \)) and liquids (\( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{l}} \)) appears to be independent of composition and, assuming ideal mixing, we obtain a value for \( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{l}} \) of 79 J mol^?1 K^?1. However, we note that a range of values for \( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{l}} \) (i.e., ~78–87 J mol^?1 K^?1) proposed by previous workers will reproduce the extended data to within experimental uncertainty. Our analysis suggests that more data are required in order to ascribe a compositional dependence (i.e., nonideal mixing) to \( C_{{{\text{p}}\;{\text{H}}_{2} {\text{O}}}}^{\text{l}} \).     

Towards a better understanding of the evolution of the flood risk in Mediterranean urban areas: the case of Barcelona

Dam failure constitutes a grave threat to human life. However, there is still a lack of systematic and comprehensive research on the loss of life (L) caused by dam break in China. From the perspective of protecting human life, a new calculation method for L occurred in dam break floods is put forward. Fourteen dam failure cases in China are selected as the basic data by three-dimensional stratified sampling, balancing spatial, vertical elevation and temporal representations, as well as considering various conditions of the dam collapse. The method includes three progressive steps: Firstly, some impact factors of loss of life (IFL) are selected by literature survey, i.e., severity of dam break flood (S _F), population at risk (P _R), understanding of dam break (U _B), warning time (T _W) and evacuation condition (E _C). And the other IFL of weather during dam break (W _B), dam break mode (M _B), water storage (S _W), building vulnerability (V _B), dam break time (T _B) and average distance from affected area to dam (D _D) are also taken into account to get a more comprehensive consideration. According to disaster system and disaster risk, these eleven IFL are divided into four categories. Through the improved entropy method, eight key IFL are further selected out of the eleven. Secondly, four L modules are built based on four categories, which are L-causing factor module (M ₁), L-prone environment module (M ₂), affected body module (M ₃) and rescue condition module (M ₄). Eventually, by using two methods of multivariate nonlinear regression and leave-one-out cross-validation in combination with coupled four modules, the calculation method for L is established. Compared with the results of Graham method and D&M method, the result of the proposed one is much closer to the actual value and performs better in fitting effect and regional applicability. In the application, L calculation and consequence assessment are carried out in the example of Hengjiang reservoir that has already broken down. At the same time, L calculation and risk prediction are used in the analysis of Yunshan reservoir, which is under planning. The proposed method can not only be applied to estimate L and its rate (f _L ) under various types of dam break conditions in China, but also provide a reliable consequence assessment and prediction approach to reduce the risk of L.