Simulation of suspended sediment based on gamma test,heuristic, and regression-based techniques

In the present study, four different heuristic techniques viz. multi-layer perceptron (MLP), radial basis function (RBF), self-organizing maps (SOM), and co-active neuro-fuzzy inference system (CANFIS) with hyperbolic tangent and sigmoid transfer functions and two regression-based techniques, i.e., multiple linear regression (MLR) and sediment-rating curve (SRC), were used for suspended sediment modeling. Gamma test (GT), correlation function (CF), M test, and trail–error procedure were applied for estimation of appropriate input variables as well as training data length. The results of the GT and CF suggested the five input variables (Q_t, Q_t?1, Q_t?2, S_t?1, and S_t?2, where Q_t?1 and S_t?1 indicate the discharge and sediment values of one previous day) as the best combination. The optimal training data length (75% of total data) was estimated by M test and trail–error procedure for development of the applied models. The MLP with sigmoid transfer function (M-2) performed better than the all other models. The results of sensitivity analysis indicated that the present-day discharge (Q_t), 1-day lag discharge (Q_t?1) and 1-day lag suspended sediment (S_t?1) are the most influenced parameters in modeling current day suspended sediment (S_t).     

Effect of uncertainty in $$V_{\mathrm{S}}{-}N$$ correlations on seismic site response analysis

This paper studies the reliability of the calculated shear-wave velocity (\(V_{\mathrm{S}}\)) from different available \(V_{\mathrm{S}}\) and SPT-N correlations in terms of seismic site response analysis. In the present study, various \(V_{\mathrm{S}}\)–N correlations developed for different regions around the globe have been used to calculate the bound of \(V_{\mathrm{S}}\) variations with depth at three different sites in Kolkata city. This bound has later been used to generate the random \(V_{\mathrm{S}}\) profiles using the Monte Carlo simulation. Equivalent linear site response analysis has been performed to study the response of those generated profiles under different input motion excitations. Strong-to-weak ground motion records have been used for this purpose. The amplification spectra of the generated \(V_{\mathrm{S}}\) profiles using all soil types and specific soil-type \(V_{\mathrm{S}}\)–N correlations show significant variations. The study also shows that the \(V_{\mathrm{S}}\)–N correlation may result in quite different \(V_{\mathrm{S},30}\) values and subsequently it may lead to the different site classes according to the NEHRP 2003 classification. So, the random choice of the \(V_{\mathrm{S}}\)–N correlation, where the direct measurement of \(V_{\mathrm{S}}\) is not available, may affect the outcome of seismic hazard analysis significantly. The study points out the need for accurate estimation of the \(V_{\mathrm{S}}\) profile either from in-situ determination or using site-specific correlation.     

Pressuremeter Modulus and Limit Pressure of Clayey Soils Using <Emphasis Type="Italic">GMDH</Emphasis>-Type Neural Network and Genetic Algorithms

Pressuremeter modulus (\(E_{M}\)) and limit pressure (\(P_{L}\)) are used for the calculation of the settlement and bearing capacity of foundation respectively. As the determination of these parameters from pressuremeter test (PMT) is relatively time-consuming and expensive, various empirical correlations have been proposed to correlate the \(E_{M}\) and \(P_{L}\) to other soil parameters. For the existing equations are incapable of estimating these PMT parameters well, in present research group method of data handling type neural network is used to estimate the \(E_{M}\) and \(P_{L}\) of clayey soils. The \(E_{M}\) and \(P_{L}\) were modeled as a function of three variables including the moisture content (\(\omega\)), plasticity index and corrected SPT blow counts (\(N_{60}\)). A database containing 51 data sets have been used for training and testing of the models. The performances of proposed models are compared with those of existing empirical equations. The results demonstrate that appreciable improvement with respect to the other correlations has been achieved. At the end, sensitivity analysis of the obtained models has been performed to study the influence of input parameters on model outputs and shows that the \(N_{60}\) is the most influential parameter on the PMT parameters.     

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.     

Data-driven nonlinear modeling studies on removal of Acid Yellow 59 using Si-doped multi-walled carbon nanotubes

Data-driven modeling of removal of color index name of Acid Yellow 59 from aqueous solutions using multi-walled carbon nanotubes by multiple (non)linear regression and artificial neural networks (ANN) models based on leave-one-out cross-validation to predict the adsorbed dye amount per unit mass of adsorbent (mg g^?1) and performance evaluation of the proposed multiple (non)linear regression and ANN models is the main novel contributor of the present study. Initial dye concentration, adsorbent concentration, reaction time, and temperature were determined as explanatory variables and input neurons for multiple (non)linear regression and ANN models, respectively. The total number of experiments was determined as 1280 statistically. The results showed that multilayer perception ANN model (\(R^{2}_{\text{training}}\) = 0.9997, \(R^{2}_{\text{testing}}\) = 0.9993, RMSE = 0.7678, MAE of 0.0007) predicted q _t better than multiple (non)linear regression model (\(R^{2}_{\text{adj}}\) = 0.9645, \(R^{2}_{\text{pred}}\) = 0.9633, SE = 9.55) and MLR (R ² = 0.9543, SE = 10.87) models. The results justified the accuracy of ANN in prediction of q _t , significantly.     

Mendigite,Mn2Mn2MnCa(Si3O9)2, a new mineral species of the bustamite group from the Eifel volcanic region,Germany

A new mineral, mendigite (IMA no. 2014-007), isostructural with bustamite, has been found in the In den Dellen pumice quarry near Mendig, Laacher Lake area, Eifel Mountains, Rhineland-Palatinate (Rheinland-Pfalz), Germany. Associated minerals are sanidine, nosean, rhodonite, tephroite, magnetite, and a pyrochlore-group mineral. Mendigite occurs as clusters of long-prismatic crystals (up to 0.1 × 0.2 × 2.5 mm in size) in cavities within sanidinite. The color is dark brown with a brown streak. Perfect cleavage is parallel to (001). D _calc = 3.56 g/cm³. The IR spectrum shows the absence of H₂O and OH groups. Mendigite is biaxial (–), α = 1.722 (calc), β = 1.782(5), γ = 1.796(5), 2V _meas = 50(10)°. The chemical composition (electron microprobe, mean of 4 point analyses, the Mn²⁺/Mn³⁺ ratio determined from structural data and charge-balance constraints) is as follows (wt %): 0.36 MgO, 10.78 CaO, 37.47 MnO, 2.91 Mn₂O₃, 4.42 Fe₂O₃, 1.08 Al₂O₃, 43.80 SiO₂, total 100.82. The empirical formula is Mn_2.00(Mn_1.33Ca_0.67) (Mn_0.50 ²⁺ Mn_0.28 ³⁺ Fe_0.15 ³⁺ Mg_0.07)(Ca_0.80 (Mn_0.20 ²⁺)(Si_5.57 Fe_0.27 ³⁺ Al_0.16O₁₈). The idealized formula is Mn₂Mn₂MnCa(Si₃O₉)₂. The crystal structure has been refined for a single crystal. Mendigite is triclinic, space group \(P\bar 1\); the unit-cell parameters are a = 7.0993(4), b = 7.6370(5), c = 7.7037(4) Å, α = 79.58(1)°, β = 62.62(1)°, γ = 76.47(1)°; V = 359.29(4) ų, Z = 1. The strongest reflections on the X-ray powder diffraction pattern [d, Å (I, %) (hkl)] are: 3.72 (32) (020), 3.40 (20) (002, 021), 3.199 (25) (012), 3.000 (26), (\(01\bar 2\), \(1\bar 20\)), 2.885 (100) (221, \(2\bar 11\), \(1\bar 21\)), 2.691 (21) (222, \(2\bar 10\)), 2.397 (21) (\(02\bar 2\), \(21\bar 1\), 203, 031), 1.774 (37) (412, \(3\bar 21\)). The type specimen is deposited in the Fersman Mineralogical Museum, Russian Academy of Sciences, Moscow, registration number 4420/1.     

Soil erosion assessment on hillslope of GCE using RUSLE model

A new method for obtaining the C factor (i.e., vegetation cover and management factor) of the RUSLE model is proposed. The method focuses on the derivation of the C factor based on the vegetation density to obtain a more reliable erosion prediction. Soil erosion that occurs on the hillslope along the highway is one of the major problems in Malaysia, which is exposed to a relatively high amount of annual rainfall due to the two different monsoon seasons. As vegetation cover is one of the important factors in the RUSLE model, a new method that accounts for a vegetation density is proposed in this study. A hillslope near the Guthrie Corridor Expressway (GCE), Malaysia, is chosen as an experimental site whereby eight square plots with the size of \(8\times 8\) and \(5\times 5\) m are set up. A vegetation density available on these plots is measured by analyzing the taken image followed by linking the C factor with the measured vegetation density using several established formulas. Finally, erosion prediction is computed based on the RUSLE model in the Geographical Information System (GIS) platform. The C factor obtained by the proposed method is compared with that of the soil erosion guideline Malaysia, thereby predicted erosion is determined by both the C values. Result shows that the C value from the proposed method varies from 0.0162 to 0.125, which is lower compared to the C value from the soil erosion guideline, i.e., 0.8. Meanwhile predicted erosion computed from the proposed C value is between 0.410 and \(3.925\, \hbox {t ha}^{-1 }\,\hbox {yr}^{-1}\) compared to 9.367 to \(34.496\, \hbox {t ha}^{-1}\,\hbox {yr}^{-1 }\) range based on the C value of 0.8. It can be concluded that the proposed method of obtaining a reasonable C value is acceptable as the computed predicted erosion is found to be classified as a very low zone, i.e. less than \(10\, \hbox {t ha}^{-1 }\,\hbox {yr}^{-1}\) whereas the predicted erosion based on the guideline has classified the study area as a low zone of erosion, i.e., between 10 and \(50\, \hbox {t ha}^{-1 }\,\hbox {yr}^{-1}\).     

Comparison of earthquake source characteristics in the Kachchh Rift Basin and Saurashtra horst,Deccan Volcanic Province,western India

Seismic source parameters of small to moderate sized intraplate earthquakes that occurred during 2002–2009 in the tectonic blocks of Kachchh Rift Basin (KRB) and the Saurashtra Horst (SH), in the stable continental region of western peninsular India, are studied through spectral analysis of shear waves. The data of aftershock sequence of the 2001 Bhuj earthquake (\(M_{w}\) 7.7) in the KRB and the 2007 Talala earthquake (\(M_{w}\) 5.0) in the SH are used for this study. In the SH, the seismic moment (\(M_{o})\), corner frequency \((f_{c})\), stress drop (\(\varDelta \sigma \)) and source radius (r) vary from \(7.8\times 10^{11}\) to \(4.0\times \)10\(^{16}\) N-m, 1.0–8.9 Hz, 4.8–10.2 MPa and 195–1480 m, respectively. While in the KRB, these parameters vary from \(M_{o} \sim 1.24 \,\times \, 10^{11}\) to \(4.1 \times 10^{16}\) N-m, \(f_{c }\sim \) 1.6 to 13.1 Hz, \(\varDelta \sigma \sim 0.06\) to 16.62 MPa and \(r \sim 100\) to 840 m. The kappa (K) value in the KRB (0.025–0.03) is slightly larger than that in the SH region (0.02), probably due to thick sedimentary layers. The estimated stress drops of earthquakes in the KRB are relatively higher than those in SH, due to large crustal stress concentration associated with mafic/ultramafic rocks at the hypocentral depths. The results also suggest that the stress drop value of intraplate earthquakes is larger than the interplate earthquakes. In addition, it is observed that the strike-slip events in the SH have lower stress drops, compared to the thrust and strike-slip events.     

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.     

Associations of Mn-bearing minerals as indicators of oxygen fugacity during the metamorphism of metalliferous deposits

The paper summarizes experimental and calculation data on the effect of oxygen fugacity on the origin of mineral assemblages in Mn-bearing rocks and demonstrates the possibility of application of these data to the reconstruction of conditions under which metalliferous deposits were metamorphosed. A new variant of the T-log\(f_{O_2 } \) diagram is proposed for the Mn-Si-O system, which differs from previous ones by the location of the lines for the formation (decomposition) of braunite and tephroite. These two minerals are the most universal indicators of oxygen fugacity during the metamorphism of Mn-bearing deposits, because these minerals are widespread in nature and can be formed in diverse environments: braunite at high \(f_{O_2 } \) values in the pore solution, and tephroite at low \(f_{O_2 } \) values. The occurrence of Mn oxides and rhodonite (pyroxmangite) in a rock makes it possible to constrain the oxygen fugacity range. An original T-log\(f_{O_2 } \) diagram is constructed for the Ca-Mn-Si-O system. As follows from this diagram, a Ca admixture expands the stability field of rhodonite toward higher oxygen fugacity values. Johannsenite can be formed in these rocks at even higher \(f_{O_2 } \). The stability of both minerals is constrained in the region of low \(f_{CO_2 } \). The paper reports data on the Fe-Si-O and Mn-Fe-Si-O systems and discusses the possibility of applying the results of experiments in the Mn-Al-Si-O system to the estimation of conditions under which andalusite, spessartine, and galaxite can be formed in Mn-bearing rocks. Data on the mineralogy of numerous Mn deposits metamorphosed under various PTX parameters indicate that the origin of Mn-bearing mineral assemblages depends not so much on the temperature and pressure as on the oxygen fugacity, which is, in turn, controlled primarily by the composition of the pristine sediments (the presence or absence of organic matter in them) and host rocks and depends on the permeability of the rocks to oxygen, the P-T conditions, and the duration of the metamorphic processes.     

Linear regression models for estimating true subsurface resistivity from apparent resistivity data

Simple linear regression (SLR) models for rapid estimation of true subsurface resistivity from apparent resistivity measurements are developed and assessed in this study. The objective is to minimize the processing time and computer memory required to carry out inversion with conventional algorithms. The arrays considered are Wenner, Wenner–Schlumberger and dipole–dipole. The parameters investigated are apparent resistivity (\(\rho _a \)) and true resistivity (\(\rho _t\)) as independent and dependent variables, respectively. For the fact that subsurface resistivity is nonlinear, the datasets were first transformed into logarithmic scale to satisfy the basic regression assumptions. Three models, one each for the three array types, are thus developed based on simple linear relationships between the dependent and independent variables. The generated SLR coefficients were used to estimate \(\rho _t\) for different \(\rho _a\) datasets for validation. Accuracy of the models was assessed using coefficient of determination (\(R^{2})\), F-test, standard error (SE) and weighted mean absolute percentage error (wMAPE). The model calibration \(R^{2}\) and F-value are obtained as 0.75 and 2286, 0.63 and 1097, and 0.47 and 446 for the Wenner, Wenner–Schlumberger and dipole–dipole array models, respectively. The SE for calibration and validation are obtained as 0.12 and 0.13, 0.16 and 0.25, and 0.21 and 0.24 for the Wenner, Wenner–Schlumberger and dipole–dipole array models, respectively. Similarly, the wMAPE for calibration and validation are estimated as 3.27 and 3.49%, 3.88 and 5.72%, and 5.35 and 6.07% for the three array models, respectively. When compared with standard constraint least-squares (SCLS) inversion and Incomplete Gauss–Newton (IGN) algorithms, the SLR models were found to reduce about 80–96.5% of the processing time and memory space required to carry out the inversion with the SCLS algorithm. It is concluded that the SLR models can rapidly estimate \(\rho _t\) for the various arrays accurately.     

Rheological properties of clayey soils originating from flow-like landslides

Flow-like landslides in clayey soils represent serious threats for populations and infrastructures and have been the subject of numerous studies in the past decade. However, despite the rising need for landslide mitigation with growing urbanization, the transient mechanisms involved in the solid-fluid transition are still poorly understood. One way of characterizing the solid-fluid transition is to carry out rheometrical tests on clayey soil samples to assess the evolution of viscosity with the shear stress. In this study, we carried out geotechnical and rheometrical tests on clayey samples collected from six flow-like landslides in order to assess if these clayey soils exhibit similar characteristics when they fluidize (solid-fluid transition). The results show that (1) all tested soils except one exhibit a yield-stress fluid behavior that can be associated with a bifurcation in viscosity (described by the critical shear rate \( \dot{\gamma_c} \)) and in shear modulus G; (2) the larger the amplitude of the viscosity bifurcation, the larger the associated drop in G; and (3) the water content (w) deviation from the Atterberg liquid limit (LL) seem a key parameter controlling a common mechanical behavior of these soils at the solid-fluid transition. We propose exponential laws describing the evolution of the critical shear stress τ_c, the critical shear rate \( \dot{\gamma_c} \), and the shear modulus G as a function of the deviation w-LL.     

Shear heating by translational brittle reverse faulting along a single,sharp and straight fault plane

Shear heating by reverse faulting on a sharp straight fault plane is modelled. Increase in temperature (T _i ) of faulted hangingwall and footwall blocks by frictional/shear heating for planar rough reverse faults is proportional to the coefficient of friction (μ), density and thickness of the hangingwall block (ρ). T _i increases as movement progresses with time. Thermal conductivity (K _i ) and thermal diffusivity (\(k_{\mathrm {i}}^{\prime }\)) of faulted blocks govern T _i but they do not bear simple relation. T _i is significant only near the fault plane. If the lithology is dry and faulting brings adjacent hangingwall and footwall blocks of the same lithology in contact, those blocks undergo the same rate of increase in shear heating per unit area per unit time.     

Equation of state of synthetic qandilite Mg<Subscript>2</Subscript>TiO<Subscript>4</Subscript> at ambient temperature

Using a diamond-anvil cell and synchrotron X-ray diffraction, the compressional behavior of a synthetic qandilite Mg_2.00(1)Ti_1.00(1)O₄ has been investigated up to about 14.9 GPa at 300 K. The pressure–volume data fitted to the third-order Birch–Murnaghan equation of state yield an isothermal bulk modulus (K _T0) of 175(5) GPa, with its first derivative \(K_{T0}^{{\prime }}\) attaining 3.5(7). If \(K_{T0}^{{\prime }}\) is fixed as 4, the K _T0 value is 172(1) GPa. This value is substantially larger than the value of the adiabatic bulk modulus (K _S0) previously determined by an ultrasonic pulse echo method (152(7) GPa; Liebermann et al. in Geophys J Int 50:553–586, 1977), but in general agreement with the K _T0 empirically estimated on the basis of crystal chemical systematics (169 GPa; Hazen and Yang in Am Miner 84:1956–1960, 1999). Compared to the K _T0 values of the ulvöspinel (Fe₂TiO₄; ~148(4) GPa with \(K_{T0}^{{\prime }} = 4\)) and the ringwoodite solid solutions along the Mg₂SiO₄–Fe₂SiO₄ join, our finding suggests that the substitution of Mg²⁺ for Fe²⁺ on the T sites of the 4–2 spinels can have more significant effect on the K _T0 than that on the M sites.     

A mathematical model on depth-averaged <Emphasis Type="Italic">β</Emphasis>-factor in open-channel turbulent flow

The well-known Rouse equation is the most widely used equation to determine the vertical distribution of suspended sediment concentration in an open-channel flow. The exponent of Rouse equation, known as Rouse number, contains the parameter β defined by the ratio of sediment diffusion coefficient to turbulent diffusion coefficient. As such to measure sediment concentration accurately, an appropriate expression for β is essentially needed. The present study, therefore, focuses on the derivation of depth-averaged β through modified expressions of sediment and turbulent diffusion coefficients. A regression analysis is done to establish the relation between β and normalized settling velocity, and the relation is used to determine suspension concentration.     

Rethinking communication in risk interpretation and action

We have studied the attenuation characteristics of eastern Himalaya and southern Tibet by using local earthquake data set that consists of 123 well-located events, recorded by the Himalayan Nepal Tibet Seismic Experiment operated during 2001–2003. We have used single backscattering model to calculate frequency-dependent values of coda Q (\(Q_\mathrm{c}\)). The estimation of \(Q_\mathrm{c}\) is made at central frequencies 2, 4, 8 and 12 Hz through five lapse time windows from 10 to 50 s starting at double the travel time of the S-wave. The observed \(Q_\mathrm{c}\) is found to be strongly frequency-dependent and follows a similar trend as observed in other tectonically active parts of the Himalaya. The trend of variation of \(Q_\mathrm{c}\) with lapse time and the corresponding apparent depths is also studied. Increase in \(Q_\mathrm{c}\) values with the lapse time suggests that the deeper part of the study region is less heterogeneous than the shallower part. The observed values of \(Q_0\) (\(Q_\mathrm{c}\) at 1 Hz) and frequency parameter n indicate that the medium beneath the study area is highly heterogeneous and tectonically very active. A regionalization of the estimated \(Q_0\) is carried out, and a contour map is prepared for the whole region. Some segments of Lesser Himalaya and Sub-Himalaya exhibit very low \(Q_0\) , while the whole Tethyan Himalaya and some parts of Greater Himalaya are characterized by low \(Q_0\) values. Our results are comparable with those obtained from tectonically active regions in the world.     

Functional Performance of Three Invasive <Emphasis Type="Italic">Marenzelleria</Emphasis> Species Under Contrasting Ecological Conditions Within the Baltic Sea

A 4-week laboratory experiment investigated the behaviour (survival and bioirrigation) and impact of the invasive polychaetes Marenzelleria viridis, M. neglecta and M. arctia on sediment-water solutes exchange, porewater chemistry, and Fe and P interactions in high-salinity sandy sediment (HSS) and low-salinity muddy sediment (LSM) from the Baltic Sea. M. viridis showed deep burrowing with efficient bioirrigation (11 L m^?2 day^?1) and high survival (71%) in HSS, while M. arctia exhibited shallow burrowing with high bioirrigation (12 L m^?2 day^?1) and survival (88%) in LSM. M. neglecta behaved poorly in both ecological settings (bioirrigation, 5–6 L m^?2 day^?1; survival, 21–44%). The deep M. viridis bioirrigation enhanced total microbial CO₂ (TCO₂) production in HSS by 175% with a net efflux of NH₄⁺ and PO₄^3?, at rates 3- to 27-fold higher than for the other species. Although the shallow and intense bioirrigation of M. arctia in LSM stimulated microbial TCO₂ production to some extent (61% enhancement), the nutrient fluxes close to zero indicate that it effectively prevented the P release. Porewater Fe:PO₄^3? ratios revealed that the oxidizing effect of M. arctia bioirrigation increased the PO₄^3? adsorption capacity of LSM twofold relative to defaunated controls while no buffering of PO₄^3? was detected in M. viridis HSS treatment. Therefore, the different behaviour of the three species in various environments and the sharp contrast between M. viridis and M. arctia effects on C, N and P cycling must be considered carefully when the ecological role of Marenzelleria species in the Baltic Sea is evaluated.     

Experimental determination of liquidus H<Subscript>2</Subscript>O contents of haplogranite at deep-crustal conditions

The liquidus water content of a haplogranite melt at high pressure (P) and temperature (T) is important, because it is a key parameter for constraining the volume of granite that could be produced by melting of the deep crust. Previous estimates based on melting experiments at low P (≤0.5 GPa) show substantial scatter when extrapolated to deep crustal P and T (700–1000 °C, 0.6–1.5 GPa). To improve the high-P constraints on H₂O concentration at the granite liquidus, we performed experiments in a piston–cylinder apparatus at 1.0 GPa using a range of haplogranite compositions in the albite (Ab: NaAlSi₃O₈)—orthoclase (Or: KAlSi₃O₈)—quartz (Qz: SiO₂)—H₂O system. We used equal weight fractions of the feldspar components and varied the Qz between 20 and 30 wt%. In each experiment, synthetic granitic composition glass + H₂O was homogenized well above the liquidus T, and T was lowered by increments until quartz and alkali feldspar crystalized from the liquid. To establish reversed equilibrium, we crystallized the homogenized melt at the lower T and then raised T until we found that the crystalline phases were completely resorbed into the liquid. The reversed liquidus minimum temperatures at 3.0, 4.1, 5.8, 8.0, and 12.0 wt% H₂O are 935–985, 875–900, 775–800, 725–775, and 650–675 °C, respectively. Quenched charges were analyzed by petrographic microscope, scanning electron microscope (SEM), X-ray diffraction (XRD), and electron microprobe analysis (EMPA). The equation for the reversed haplogranite liquidus minimum curve for Ab_36.25Or_36.25Qz_27.5 (wt% basis) at 1.0 GPa is \(T = - 0.0995 w_{{{\text{H}}_{ 2} {\text{O}}}}^{ 3} + 5.0242w_{{{\text{H}}_{ 2} {\text{O}}}}^{ 2} - 88.183 w_{{{\text{H}}_{ 2} {\text{O}}}} + 1171.0\) for \(0 \le w_{{{\text{H}}_{ 2} {\text{O}}}} \le 17\) wt% and \(T\) is in °C. We present a revised \(P - T\) diagram of liquidus minimum H₂O isopleths which integrates data from previous determinations of vapor-saturated melting and the lower pressure vapor-undersaturated melting studies conducted by other workers on the haplogranite system. For lower H₂O (<5.8 wt%) and higher temperature, our results plot on the high end of the extrapolated water contents at liquidus minima when compared to the previous estimates. As a consequence, amounts of metaluminous granites that can be produced from lower crustal biotite–amphibole gneisses by dehydration melting are more restricted than previously thought.