Organic carbon deposition and phosphorus accumulation during Oceanic Anoxic Event 2 in Tarfaya, Morocco

With a multi-proxy approach, an attempt was made to constrain productivity and bottom-water redox conditions and their effects on the phosphorus accumulation rate at the Mohammed Plage section on the Tarfaya coast, Morocco, during the Cenomanian-Turonian Anoxic Event (OAE 2). A distinct δ¹³C_org isotope excursion of +2.5‰ occurs close to the top of the section. The unusually abrupt shift of the isotope excursion and disappearance of several planktonic foraminiferal species (e.g. Rotalipora cushmani and Rotalipora greenhornensis) in this level suggests a hiatus of between 40–60 kyrs at the excursion onset. Nevertheless, it was possible to determine both the long-term environmental history as well as the processes that took place immediately prior to and during OAE 2. TOC% values increase gradually from the base of the section to the top (from 2.5% to 10%). This is interpreted as the consequence of a long-term eustatic sea-level rise and subsidence causing the encroachment of less oxic waters into the Tarfaya Basin. Similarly a reduction in the mineralogically constructed ‘detrital index’ can be explained by the decrease in the continental flux of terrigenous material due to a relative sea-level rise. A speciation of phosphorus in the upper part of the section, which spans the start and mid-stages of OAE 2, shows overall higher abundances of P_reactive mass accumulation rates before the isotope excursion onset and lower values during the plateau. Due to the probable short hiatus, the onset of the decrease in phosphorus content relative to the isotope excursion is uncertain, although the excursion plateau already contains lower concentrations. The C_org/P_total and V/Al ratios suggest that this reduction was mostly likely caused by a decrease in the available bottom oxygen content (probably as a result of higher productivity) and a corresponding fall in the phosphorus retention ability of the sediment. Productivity appears to have remained high during the isotope plateau possibly due to a combination of ocean-surface fertilisation via increased aridity (increased K/Al and Ti/Al ratios) and/or higher dissolved inorganic phosphorus content in the water column as a result of the decrease in sediment P retention. The evidence for decreased P-burial has been observed in many other palaeoenvironments during OAE 2. Tarfaya's unique upwelling paleosituation provides strong evidence that the nutrient recycling was a global phenomenon and therefore a critical factor in starting and sustaining OAE 2.     

Late quaternary deformation and palaeoseismicity: Insight into geomorphotectonic evolution of Kachchh, western Indian Subcontinent

Kachchh in western Indian Shield, according to the Bureau of Indian Standard (IS:1893:2002), falls in Seismic Zone V. This is intriguing considering that the region is far away from active Plate margin. Apart from the recent incidences of earthquakes, there are several pre-historic/archaeological records of earthquakes in the region. Beyond these, the geological evidence of earth-movements (causing earthquakes) is provided by the occurrence of several’ active’ faults, which are considered geological markers of palaeoseismicity. There are records of innumerable incidences of faulting in the region in not so distant geological past. Study of fault features especially the scarp faces marking abrupt change in physical relief proves that the different levels of topography in the entire terrain are fault-bound features. Studies also confirm that the topographic difference between the high and ’sunken’ features have formed due to uplift and relative down-sagging during the geomorphotectonic evolution of the terrain. Features that make the region unique are: (i) restriction of fault-related deformation zone to a narrow strip between the southern margin of Thar Desert and the south coast line of the Kachchh Peninsula; (ii) overall sub-horizontality of bedding and other topographic and planation surfaces over the entire region; (iii) evidence of fault-controlled geomorphology indicating vertical movement along fault planes; (iv) evidence constraining the time of geomorphological evolution of the terrain only during the Late Quaternary, making it the youngest neotectonically evolved terrain in the Precambrian Indian Shield.     

Influence of seismic motions on behavior of piles in liquefied soils

In the present study an analytical procedure based on finite element technique is proposed to investigate the influence of vertical load on deflection and bending moment of a laterally loaded pile embedded in liquefiable soil, subjected to permanent ground displacement. The degradation of subgrade modulus due to soil liquefaction and effect of nonlinearity are also considered. A free headed vertical concrete elastic nonyielding pile with a floating tip subjected to vertical compressive loading, lateral load, and permanent ground displacement due to earthquake motions, in liquefiable soil underlain by nonliquefiable stratum, is considered. The input seismic motions, having varying range of ground motion parameters, considered here include 1989 Loma Gilroy, 1995 Kobe, 2001 Bhuj, and 2011 Sikkim motions. It is calculated that maximum bending moment occurred at the interface of liquefiable and nonliquefiable soil layers and when thickness of liquefiable soil layer is around 60% of total pile length. Maximum bending moment of 1210 kNm and pile head deflection of 110 cm is observed because of 1995 Kobe motion, while 2001 Bhuj and 2011 Sikkim motions amplify the pile head deflection by 14.2 and 14.4 times and bending moment approximately by 4 times, when compared to nonliquefiable soil. Further, the presence of inertial load at the pile head increases bending moment and deflection by approximately 52% when subjected to 1995 Kobe motion. Thus, it is necessary to have a proper assessment of both kinematic and inertial interactions due to free field seismic motions and vertical loads for evaluating pile response in liquefiable soil.     

Variations in shear wave velocity and soil site class in Kolkata city using regression and sensitivity analysis

The detrimental effects of an earthquake are strongly influenced by the response of soils subjected to dynamic loading. The behavior of soils under dynamic loading is governed by the dynamic soil properties such as shear wave velocity, damping characteristics and shear modulus. Worldwide, it is a common practice to obtain shear wave velocity (V _s in m/s) using the correlation with field standard penetration test (SPT) N values in the absence of sophisticated dynamic field test data. In this paper, a similar but modified advanced approach has been proposed for a major metro city of eastern India, i.e., Kolkata city (latitudes 22°20′N–23°00′N and longitudes 88°04′E–88°33′E), to obtain shear wave velocity profile and soil site classification using regression and sensitivity analyses. Extensive geotechnical borehole data from 434 boreholes located across 75 sites in the city area of 185 km² and laboratory test data providing information on the thickness of subsoil strata, SPT N values, consistency indices and percentage of fines are collected and analyzed thoroughly. A correlation between shear wave velocity (V _s) and SPT N value for various soil profiles of Kolkata city has been established by using power model of nonlinear regression analysis and compared with existing correlations for other Indian cities. The present correlations, having regression coefficients (R ²) in excess of 0.96, indicated good prediction capability. Sensitivity analysis predicts that significant influence of soil type exists in determining V _s values, for example, typical silty sand shows 30.4 % increase in magnitude of V _s as compared to silt of Kolkata city. Moreover, the soil site classification shows Class D and Class E category of soil that exists typically in Kolkata city as per NEHRP (Recommended provisions for seismic regulations for new buildings and other structures—Part 1: Provisions. Prepared by the Building Seismic Safety Council for the Federal Emergency Management Agency (Report FEMA 450), Washington, DC, 2003) guidelines and thereby highlighting the seismic vulnerability of the city. The results presented in this study can be utilized for seismic microzonation, ground response analysis and hazard assessment for Kolkata city.

     

Lithium, boron and chlorine as tracers for metasomatism in high-pressure metamorphic rocks: a case study from Syros (Greece)

High-pressure metamorphic (HPM) rocks (derived from igneous protoliths) and their metasomatised rinds from the island of Syros (Greece) were analysed for their B and Cl whole-rock abundances and their H₂O content by prompt-gamma neutron-activation analysis (PGNAA) and for their Li and Be whole-rock abundances by ICP-OES. In the HPM rocks, B?/Be and Cl?/Be ratios correlate with H₂O contents and appear to be controlled by extraction of B and Cl during dehydration and prograde metamorphism. In contrast, samples of the metasomatised rinds show no such correlation. B?/Be ratios in the rinds are solely governed by the presence or absence of tourmaline, and Cl?/Be ratios vary significantly, possibly related to fluid inclusions. Li/Be ratios do not correlate with H₂O contents in the HPM rocks, which may in part be explained by a conservative behaviour of Li during dehydration. However, Li abundances exceed the vast majority of published values for Li abundances in fresh, altered, or differentiated oceanic igneous rocks and presumably result from metasomatic enrichment of Li. High Li concentrations and highly elevated Li/Be ratios in most metasomatised samples demonstrate an enrichment of Li in the Syros HP mélange during fluid infiltration. This study suggests that B and Cl abundances of HPM meta-igneous rocks can be used to trace prograde dehydration, while Li concentrations seem to be more sensitive for retrograde metasomatic processes in such lithologies.     

Bayes estimation: A novel approach to derivation of internally consistent thermodynamic data for minerals,their uncertainties,and correlations. Part II: Application

Internally consistent thermodynamic data, including their uncertainties and correlations, are reported for 22 phases of the quaternary system CaO-Al₂O₃-SiO₂-H₂O. These data have been derived by simultaneous evaluation of the appropriate phase properties (PP) and reaction properties (RP) by the novel technique of Bayes estimation (BE). The thermodynamic model used and the theory of BE was expounded in Part I of this paper. Part II is the follow-up study illustrating an application of BE. The input for BE comprised, among others, the a priori values for standard enthalpy of formation of the i-th phase, Δ_f H _i ⁰ , and its standard entropy, S _i ⁰ , in addition to the reaction reversal constraints for 33 equilibria involving the relevant phases. A total of 269 RP restrictions have been processed, of which 107 turned out to be non-redundant. The refined values for Δ_f H _i ⁰ and S _i ⁰ obtained by BE, including their 2σ-uncertainties, appear in Table 4; the Appendix reproduces the corresponding correlation matrix. These data permit generation of computed phase diagrams with 2σ-uncertainty envelopes based on conventional error propagation; Fig. 3 depicts such a phase diagram for the system CaO-Al₂O₃-SiO₂. It shows that the refined dataset is capable of yielding phase diagrams with uncertainty envelopes narrow enough to be geologically useful. The results in Table 4 demonstrate that the uncertainties of the prior values for Δ_f H _i ^Emphasis>0 , given in Table 1, have decreased by up to an order of magnitude, while those for S _i ⁰ improved by a factor of up to two. For comparison, Table 4 also lists the refined Δ_f H _i ⁰ and S _i ⁰ data obtained by mathematical programming (MAP), minimizing a quadratic objective function used earlier by Berman (1988). Examples of calculated phase diagrams are given to demonstrate the advantages of BE for deriving internally consistent thermodynamic data. Although P-T curves generated from both MAP and BE databases will pass through the reversal restrictions, BE datasets appear to be better suited for extrapolations beyond the P-T range explored experimentally and for predicting equilibria not constrained by reversals.     

An internally consistent set of thermodynamic data for twentyone CaO-Al2O3-SiO2- H2O phases by linear parametric programming

The technique of linear parametric programming has been applied to derive sets of internally consistent thermodynamic data for 21 condensed phases of the quaternary system CaO-Al₂O₃-SiO₂-H₂O (CASH) (Table 4). This was achieved by simultaneously processing:

1. calorimetric data for 16 of these phases (Table 1), and
2. experimental phase equilibria reversal brackets for 27 reactions (Table 3) involving these phases.

Calculation of equilibrium P-T curves of several arbitrarily picked reactions employing the preferred set of internally consistent thermodynamic data from Table 4 shows that the input brackets are invariably satisfied by the calculations (Fig. 2a). By contrast, the same equilibria calculated on the basis of a set of thermodynamic data derived by applying statistical methods to a large body of comparable input data (Haas et al. 1981; Hemingway et al. 1982) do not necessarily agree with the experimental reversal brackets. Prediction of some experimentally investigated phase relations not included into the linear programming input database also appears to be remarkably successful. Indications are, therefore, that the thermodynamic data listed in Table 4 may be used with confidence to predict geologic phase relations in the CASH system with considerable accuracy. For such calculated phase diagrams and their petrological implications, the reader's attention is drawn to the paper by Chatterjee et al. (1984).     

Mildly alkaline basalts from Pavagadh Hill, India: Deccan flood basalts with an asthenospheric origin

Summary Of twelve flows at Pavagadh Hill, the two three-phenocryst-basalt flows with Mg#0.70 and Ni/MgO33 are the most primitive and perhaps as primitive as any basalts in the Deccan province. Scatter on variation diagrams and the occurrence of primitive flows at two different levels in the volcanic sequence implies that most rocks are probably not, strictly speaking, comagmatic. Nevertheless, mass balance calculations indicate a generalized differentiation scheme from primitive basalt to hawaiite that involved removal of olivine, augite, plagioclase and Fe-Ti oxides in the proportions 40:33:22:5 with 50% of the magma remaining. Crustal assimilation had a minimal effect on evolution of the basalts but rhyolites at the top of the volcanic sequence may have been produced by crustal melting following prolonged heat release from alkali basalt pooled along fault zones in the continental crust. Major element based calculations indicate that the most primitive basalts were generated by 7 to 10% melting of mantle peridotite. These low percentages of melting, typical of alkali basalts, are consistent with the steep slopes on chondrite-normalized REE diagrams. Low heavy REE concentrations point to residual garnet in the source region. Incompatible element concentrations (e.g. Rb, Ba, Zr, La) in Pavagadh basalts exceed those in Deccan tholeiitic basalts but are substantially lower than those reported for some other Deccan alkali basalts. Obviously Pavagadh basalts do not reflect the lowest percentages of melting and greatest amount of source region metasomatic enrichment attained in the Deccan province. Deccan tholeiitic and alkali basalts are largely characterized by low La/Nb ratios and high La/Ba ratios similar to those in oceanic island basalts. This indicates minimal involvement of the subcontinental lithospheric mantle in their petrogenesis. Comparison with continental mafic magma provinces where a subcontinental lithospheric mantle imprint is common indicates long periods of extension and/or melting of mantle lithosphere still hot from pre-extension subduction are more likely to produce magmas bearing the lithospheric imprint.

---

Alkalische Basalte von Pavagadh Hill, Indien: Deccan-Flutbasalte von Astenosphärischer Herkunft

Zusammenfassung Im Gebiet von Pavagadh Hill, Indien, treten 12 Spät-Deccan und rhyolithische alkalibasaltische Ergüsse und Intrusiva auf. Variationsdiagramme zeigen, daß die Abfolge nicht komagmatisch ist. Zusammen mit Berechnungen der Massenbilanz unterstützen sie vielmehr ein Zwei-Stadienmodell für die Entstehung von Hawaiiten aus sehr primitiven (i.e. Mg#=Mg/(Mg+.(0.9*Fe_totat)) at.%0.70) Basalten. Olivin und Augit dominierten die frühe Fraktionierung während Augit vorherrschte als der Magmaanteil von 65% auf 50% sank. Die Entfernung von Plagioklas spielte bei der Differentiation nur eine geringe Rolle. Niedrige Th/Nb (0,2), Rb/Sr(<0,12) und K/NbVerhältnisse geben keine Hinweise auf signifikante Assimilation von Krustenmaterial. Die Seltene-Erd-Verteilungsmuster (SEE), niedrige Gehalte an schweren SEE sowie die Hauptelementspektren der Alkalibasalte weisen auf eine granatführende Ursprungsregion und auf einen Aufschmelzungsgrad von nur 7% bis 10% hin. Es gibt jedoch auch stärker alkalische (höhere Rb, Zr etc.) Deccanbasalte (i.e. Rajpipla). Die Assoziation von Deccanalkalibasalten, Rhyolithen und Störungszonen zeigt, daß letztere die Extraktion von Magma aus dem Mantel erleichterten und dazu führten, daß Magma aus Magmenkammern Krustenschmelzen (Rhyolithe) produzierte. Deccanbasalte tendieren zu hohen La/Ba und niedrigen La/Nb-Verhältnissen; dies weist auf eine asthenosphärische Herkunft hin, selbst wenn die Gesteine verhältnismäßig spät gebildet wurden (i.e. Pavagadh). Längere Perioden von Krustenextension oder von Subduktion, die der Extension vorhergeht, führt offensichtlich zur Entstehung von Magmen mit einer lithosphärischen Komponente.

     

A statistical procedure for the analysis of seismotectonically induced hydrochemical signals: A case study from the Eastern Carpathians, Romania

Changes in the stress field of an aquifer system induced by seismotectonic activity may change the mixing ratio of groundwaters with different compositions in a well, leading to hydrochemical signals which in principle could be related to discrete earthquake events. Due to the complexity of the interactions and the multitude of involved factors the identification of such relationships is a difficult task. In this study we present an empiric statistical approach suitable to analyse if there is an interdependency between changes in the chemical composition of monitoring wells and the regional seismotectonic activity of a considered area. To allow a rigorous comparison with hydrochemistry the regional earthquake time series was aggregated into an univariate time series. This was realized by expressing each earthquake in form of a parameter “e”, taking into consideration both energetic (magnitude of a seismic event) and spatial parameters (position of epi/hypocentrum relative to the monitoring site). The earthquake and the hydrochemical time-series were synchronised aggregating the e-parameters into “earthquake activity” functions E, which takes into account the time of sampling relative to the earthquakes which occurred in the considered area. For the definition of the aggregation functions a variety of different “e” parameters were considered. The set of earthquake functions E was grouped by means of factor analysis to select a limited number of significant and representative earthquake functions E to be used further on in the relation analysis with the multivariate hydrochemical data set. From the hydrochemical data a restricted number of hydrochemical factors were extracted. Factor scores allow to represent and analyse the variation of the hydrochemical factors as a function of time. Finally, regression analysis was used to detect those hydrochemical factors which significantly correlate with the aggregated earthquake functions.This methodological approach was tested with a hydrochemical data set collected from a deep well monitored for two years in the seismically active Vrancea region, Romania. Three of the hydrochemical factors were found to correlate significantly with the considered earthquake activities. A screening with different time combinations revealed that correlations are strongest when the cumulative seismicity over several weeks was considered. The case study also showed that the character of the interdependency depends sometimes on the geometrical distribution of the earthquake foci. By using aggregated earthquake information it was possible to detect interrelationships which couldn't have been identified by analysing only relations between single geochemical signals and single earthquake events. Further on, the approach allows to determine the influence of different seismotectonic patterns on the hydrochemical composition of the sampled well. The method is suitable to be used as a decision instrument in assessing if a monitoring site is suitable or not to be included in a monitoring net within a complex earthquake prediction strategy.     

Riddle of arsenic in groundwater of Bengal Delta Plain—role of non-inland source and redox traps

In the Bengal Delta Plain (BDP) the primary arsenic sourcing appears to be different from the global scenario. Here, the Terminal Pleistocene–Holocene depositional platform, the interactive early Holocene depositional morphology with fluvio-estuarine and marine incursions played a crucial role for arsenic sourcing and enrichment. The lenticular silt-fine sand layer between anoxic clay beds favoured entrapment of dissolved organic carbon with decayed phyto-planktons debris. The Terminal Pleistocene–Holocene transgression and regression processes may have acted as major events in the BDP. Interestingly, at the end of the last glacial maxima, the Pleistocene delta had undergone block movements, wherein some parts of the platform were raised above the level of Holocene deposition. Those blocks were found to be free from arsenic in the groundwater. The sea, during re-emerging inundation (10–7 ka BP), has witnessed a monsoon-induced environment in the BDP with the resultant oscillation of sea level leading to higher upsurge towards the north. This might have resulted in the marine incursion and inundation in pre-existing land depressions. Meanwhile arsenic entrapments through marine incursion as well as enrichment in the presence of organic carbon/DOC and/or Fe/Mn/Al catalytic agents could have developed into localised redox traps. It may be of relevance that due to the repetitive transgressive–regressive phases in Holocene, resulting in periodic exposure and weathering of iron-bearing minerals and consequent iron enrichment in the aquifer system. The iron, thus present, had free charge to host arsenic as a sink. It appears that arsenic, wherever found, would likely be of atypical localised exhaustible phenomenon, both in horizontal and vertical context. It also rationalises the cause of the absence of arsenic in the other nearby Pleistocene platform, which has not come across Holocene interaction and marine incursion, as to the likely limiting condition for the search for arsenic in the BDP or beyond.