Record of enhanced late Pliocene tectono‐volcanic activity in the Lesser Antilles forearc basin: Insights into the long‐term behaviour of the Lesser Antilles subduction zone

In this study, we investigate the Lesser Antilles forearc basin, focusing on the late Pliocene to Pleistocene sedimentary archives in order to track the occurrence of extreme events triggered by enhanced subduction‐related tectono‐volcanic activity. We identify late Piacenzian deposits covering a major regional erosional surface, displaying sedimentary dykes and large marine boulders embedded in a mixed continental–marine matrix, characteristic of tsunamites. We interpret this episode of platform emersion and the successive cataclysmic deposits as resulting from enhanced tectonic activity at the interface of the subduction zone, synchronous with the initiation of the Lesser Antilles volcanic arc. We then discuss the implications in terms of the mechanical behaviour of the Lesser Antilles subduction zone.     

Numerical prediction of an intense convective system associated with the July 1987 montreal flood. part II: A trailing stratiform rainband

Abstract

In this study, the internal circulation structures of the 14 July 1987 intense mesoscale convective system (MCS) are investigated using an improved high‐resolution version of the Canadian regional finite‐element model. It is found that although the MCS is characterized by a leading convective line followed by a trailing stratiform rainband, the associated circulation structures differ substantially from those in the classical midlatitude squall system. These include the rapid propagation and separation of the leading convection from the trailing rainband, the development of a surface‐based instead of an elevated rear‐to‐front descending flow and a shallow front‐to‐rear ascending flow associated with the stratiform precipitation, the generation of low‐ and mid‐level rather than mid‐ to upper‐level stratiform cloudiness and the development of a strong anticyclonic vorticity band at the back edge of the stratiform region. It is shown that the trailing stratiform rainband is dynamically forced by frontogenetical processes, and aided by the release of conditional symmetric instability and local orographical lifting. The intense leading and trailing circulations result from latent heat released by the convective and explicit cloud schemes, respectively. Sensitivity experiments reveal that the proper coupling of these two cloud schemes is instrumental in obtaining a realistic prediction of the above‐mentioned various mesoscale components. Vorticity budget calculations show that tilting of horizontal vorticity contributes the most to the amplification of the anticyclonic vorticity band, particularly during the squall's incipient stage. The sensitivity of the simulated squall system to other model physical parameters is also examined.     

Spatial patterns of groundwater‐lake exchange – implications for acid neutralization processes in an acid mine lake

Exchange of groundwater and lake water with typically quite different chemical composition is an important driver for biogeochemical processes at the groundwater‐lake interface, which can affect the water quality of lakes. This is of particular relevance in mine lakes where anoxic and slightly acidic groundwater mixes with oxic and acidic lake water (pH < 3). To identify links between groundwater‐lake exchange rates and acid neutralization processes in the sediments, exchange rates were quantified and related to pore‐water pH, sulfate and iron concentrations as well as sulfate reduction rates within the sediment. Seepage rates measured with seepage meters (?2.5 to 5.8 L m^‐2 d^‐1) were in reasonable agreement with rates inverted from modeled chloride profiles (?1.8 to 8.1 L m^‐2 d^‐1). Large‐scale exchange patterns were defined by the (hydro)geologic setting but superimposed by smaller scale variations caused by variability in sediment texture. Sites characterized by groundwater upwelling (flow into the lake) and sites where flow alternated between upwelling and downwelling were identified. Observed chloride profiles at the alternating sites reflected the transient flow regime. Seepage direction, as well as seepage rate, were found to influence pH, sulfate and iron profiles and the associated sulfate reduction rates. Under alternating conditions proton‐consuming processes, for example, sulfate reduction, were slowed. In the uppermost layer of the sediment (max. 5 cm), sulfate reduction rates were significantly higher at upwelling (>330 nmol g^‐1 d^‐1) compared to alternating sites (<220 nmol g^‐1 d^‐1). Although differences in sulfate reduction rates could not be explained solely by different flux rates, they were clearly related to the prevailing groundwater‐lake exchange patterns and the associated pH conditions. Our findings strongly suggest that groundwater‐lake exchange has significant effects on the biogeochemical processes that are coupled to sulfate reduction such as acidity retention and precipitation of iron sulfides. Copyright © 2012 John Wiley & Sons, Ltd.     

Anticipating the effects of groundwater withdrawal on seawater intrusion and soil settlement in urban coastal areas

Intensive pumping in urban coastal areas is a common threat to water resource quality due to seawater intrusion. In those areas where subsurface water resources are not usually used for human consumption or irrigation, intensive pumping is associated with other activities like the lowering of the water table necessary to support underground structures and building foundations. This activity also increases the likelihood of soil settlement that affects building stability and the corrosion of concrete structures due to groundwater salinity. Under these circumstances, the awareness of a certain municipality (Calonge, NE Spain) of the potential effects of groundwater withdrawal upon foundations has led to an integrated approach to anticipate seawater intrusion related to urban development. Geological mapping and correlation of borehole logs, electrical resistivity tomography, and hydrochemical data provide comprehensive knowledge of the geology and hydrogeology of the area and act as screening tools necessary to discern the influence of hydrological processes in coastal areas. Developing Strack's analytical solution, new comprehensive, dimensionless expressions are herein derived to determine the critical pumping rate necessary to prevent seawater intrusion, as well as to reproduce the evolution of the wedge toe and the water table stagnation point under different withdrawal rates. Furthermore, the Dupuit–Forchheimer well discharge formula allows the estimation of the effects of the water table lowering due to such critical pumping in the surrounding building foundations. Field data from the Calonge coastal plain illustrate this approach and provide assessment criteria for future urban development and planning. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical analysis of frost effects in porous media. Benefits and limits of the finite element poroelasticity formulation

The aim of this paper is to analyse the performance of a finite element formulation usable for predicting the mechanical consequences of frost effects on porous media. It considers the characteristics of porous media and how the frost action can be assessed. The problem is then separated into two parts: thermal and poromechanical calculations. The constitutive equations developed in the framework of poromechanics are presented and the implementation in a usual finite element poroelasticity formulation based on Zuber's method is adopted. An analysis of the time‐step influence on the convergence rate is given and leads us to propose a simple method in order to obtain objectivity of the finite element response and avoid over‐long calculations. Frost effect simulations are carried out on real porous media (two fired clays) as a case study. Although the experimental behaviour of the porous media subjected to frost action is in accordance with some observations, the calculated strains appear to be overestimated compared with measurements. The problem could be largely attributable to the difficulty of assessing permeability evolution during frost development. Copyright © 2011 John Wiley & Sons, Ltd.     

Non linear waves in two dimensional keplerian flows

Abstract

Accretion discs in astrophysics are fundamental for converting gravitational binding energy into observed electromagnetic radiation. We study the behavior of waves in a two dimensional supersonic Keplerian flow inside a given gravitational potential. We present the effects of shearing and rotation on short waves, and the numerical study of the dynamical stability of such flows with respect to various perturbations. We show that a large class of dynamical effects, due to pressure and associated to short time scales, may be excited.     

TITAN2D simulations of pyroclastic flows at Cerro Machín Volcano,Colombia: Hazard implications

Cerro Machín is a dacitic tuff ring located in the central part of the Colombian Andes. It lies at the southern end of the Cerro Bravo–Cerro Machín volcanic belt. This volcano has experienced at least six major explosive eruptions during the last 5000 years. These eruptions have generated pyroclastic flows associated with Plinian activity that have traveled up to 8 km from the crater, and pyroclastic flows associated with Vulcanian activity with shorter runouts of 5 km from the source. Today, some 21,000 people live within a 8 km radius of Cerro Machín. The volcano is active with fumaroles and has shown increasing seismic activity since 2004, and therefore represents a potentially increasing threat to the local population. To evaluate the possible effects of future eruptions that may generate pyroclastic density currents controlled by granular flow dynamics we performed flow simulations with the TITAN2D code. These simulations were run in all directions around the volcano, using the input parameters of the largest eruption reported. The results show that an eruption of 0.3 km³ of pyroclastic flows from a collapsing Plinian column would travel up to 9 km from the vent, emplacing a deposit thicker than 60 m within the Toche River valley. Deposits >45 m thick can be expected in the valleys of San Juan, Santa Marta, and Azufral creeks, while 30 m thick deposits could accumulate within the drainages of the Tochecito, Bermellón, and Coello Rivers. A minimum area of 56 km² could be affected directly by this kind of eruption. In comparison, Vulcanian column-collapse pyroclastic flows of 0.1 km³ would travel up to 6 km from the vent depositing >45 m thick debris inside the Toche River valley and more than 30 m inside the valleys of San Juan, Santa Marta, and Azufral creeks. The minimum area that could be affected directly by this kind of eruption is 33 km². The distribution and thickness of the deposits obtained by these simulations are consistent with the hazard map presented by INGEOMINAS (Geological Survey of Colombia) in 2002. The composite map of the simulated flow deposits suggests that after major explosive events such as these, the generation of lahars is probable.     

Formation of intracratonic basins by lithospheric shortening and phase changes: a case study from the ultra‐deep East Barents Sea basin

The origin of large subsidence in intracratonic basins is still under debate. We propose a new and self‐consistent model for the formation of those basins, where lithospheric shortening/buckling triggers metamorphism and densification of crustal mafic heterogeneities. We use a forward thermo‐mechanical finite element technique to evaluate this mechanism for the typical example of the East Barents Sea basin (EBB) where a very large and compensated subsidence, accommodating an up to 20‐km‐thick sediment succession, is observed. The lower crust in the dynamic model is modelled with petrologic‐consistent densities for a wet mafic gabbroic composition that depend on pressure and temperature taking into account dehydration at high PT conditions. The model successfully explains the main characteristics of the EBB, notably the large anomalous and fast subsidence during the Late Permian–Early Triassic, its present‐day geometry and the absence of a significant gravity anomaly.     

Oxygen and organic carbon fluxes in sediments of the Bay of Biscay

The relationship between particulate organic carbon (POC) concentrations measured in modern sediment and fluxes of exported POC to the sediment surface needs to be understood in order to use POC content as a proxy of paleo-environmental conditions. The objective of our study was to compare POC concentrations, POC mineralization rates calculated from O₂ consumption and POC burial rates. Benthic O₂ distributions were determined in 58 fine-grained sediment cores collected at different periods at 14 stations in the southeastern part of the Bay of Biscay with depths ranging from 140 to 2800 m. Depth-dependent volume-specific oxygen consumption rates were used to assess rates of aerobic oxidation of organic matter (OM), assuming that O₂ consumption solely was related to heterotrophic activity at the sediment–water interface. Heterogeneity of benthic O₂ fluxes denoted changes in time and space of fresh organic material sedimentation. The most labile fraction of exported POC engendered a steep decrease in concentration in the upper 5 mm of vertical O₂ profiles. The rupture in the gradient of O₂ microprofile may be related to the bioturbation-induced mixing depth of fast-decaying carbon. Average diffusive O₂ fluxes showed that this fast-decaying OM flux was much higher than buried POC, although diffusive O₂ fluxes underestimated the total sediment oxygen demand, and thus the fast-decaying OM flux to the sediment surface. Sedimentary POC burial was calculated from sediment mass accumulation rate and the organic carbon content measured at the top of the sediment. The proportion of buried POC relative to total exported POC ranged at the most between 50% and 10%, depending on station location. Therefore, for a narrow geographic area like the Bay of Biscay, burial efficiency of POC was variable. A fraction of buried POC consisted of slow-decaying OM that was mineralized within the upper decimetres of sediment through oxic and anoxic processes. This fraction was deduced from the decrease with depth in POC concentration. At sites located below 500 m water depth, where the fast-decaying carbon did not reach the anoxic sediment, the slow-decaying pool may control the O₂ penetration depth. Only refractory organic material was fossilized in sedimentary records at locations where labile OM did not reach the anoxic portion of the sediment.