Hydrothermal calderas

The standard model of caldera formation is related to the emptying of a magma chamber and ensuing roof collapse during large eruptions or subsurface withdrawal. Although this model works well for numerous volcanoes, it is inappropriate for many basaltic volcanoes (with the notable exception of Hawaii), as these have eruptions that involve volumes of magma that are small compared to the collapse. Many arc volcanoes also have similar oversized depressions, such as Poas (Costa Rica) and Aoba (Vanuatu). In this article, we propose an alternative caldera model based on deep hydrothermal alteration of volcanic rocks in the central part of the edifice. Under certain conditions, the clay-rich altered and pressurized core may flow under its own weight, spread laterally, and trigger very large caldera-like collapse. Several specific mechanisms can generate the formation of such hydrothermal calderas. Among them, we identify two principal modes: mode 1: ripening with summit loading and flank spreading and mode II: unbuttressing with flank subsidence and flank sliding. Processes such as summit loading or flank subsidence may act simultaneously in hybrid mechanisms. Natural examples are shown to illustrate the different modes of formation. For ripening, we give Aoba (Vanuatu) as an example of probable summit loading, while Casita (Nicaragua) is the type example of flank spreading. For unbuttressing, Nuku Hiva Island (Marquesas) is our example for flank subsidence and Piton de la Fournaise (La Réunion) is our example of flank sliding. The whole process is slow and probably needs (a) at least a few tens of thousands of years to deeply alter the edifice and reach conditions suitable for ductile flow and (b) a few hundred years to achieve the caldera collapse. The size and the shape of the caldera strictly mimic that of the underlying weak core. Thus, the size of the caldera is not controlled by the dimensions of the underlying magma reservoir. A collapsing hydrothermal caldera could generate significant phreatic activity and trigger major eruptions from a coexisting magmatic complex. As the buildup to collapse is slow, such caldera-forming events could be detected long before their onset.     

THE IMPACT-FLOOD CONNECTION: DOES IT EXIST?*

In a recent article in Terra Nova, Kristan-Tollmann and Tollmann (1994) suggested that the Biblical Flood can be explained by seven fragments of a comet that impacted the ocean at seven locations on Earth at 03.00h (C.E.T.) on 23 September, 9545 yr BP. We demonstrate that all the ‘geological proofs’ that allegedly support their conclusions are not supported by the available data on impact cratering. Their hypothesis is based on insufficient and ambiguous data, selective citation, and incomplete comprehension of previous research.     

A statistical method for varve verification using seasonal pollen deposition

We developed a fast, inexpensive, statistically rigorous method of varve verification that uses prospective varve layer-splitting and seasonal pollen deposition. This method can be used on any sediment that contains seasonally deposited pollen, and avoids the need for radio-isotopic, optical, or thermoluminescence dating. The method uses a χ² test and non-parametric regression, together with recorded plant bloom times from pollen traps, to assess the differences in pollen abundances in light/dark (i.e. summer/winter-spring) sediment layers. A statistical test is required because such seasonally deposited data are inherently noisy, with a low signal-to-noise ratio. To illustrate this approach, visible laminations of Lake Mina, Minnesota, USA, were assessed in two separate regions of a core spanning 900 years. Results show that the laminations were deposited annually.     

Development of a Darcy-Brinkman model to simulate water flow and tracer transport in a heterogeneous karstic aquifer (Val d’Orléans, France)

Darcy’s law is the equation of reference widely used to model aquifer flows. However, its use to model karstic aquifers functioning with large pores is problematic. The physics occurring within the karstic conduits requires the use of a more representative macroscopic equation. A hydrodynamic model is presented which is adapted to the karstic aquifer of the Val d’Orléans (France) using two flow equations: (1) Darcy’s law, used to describe water flow within the massive limestone, and (2) the Brinkman equation, used to model water flow within the conduits. The flow equations coupled with the transport equation allow the prediction of the karst transfer properties. The model was tested by using six dye tracer tests and compared to a model that uses Darcy’s law to describe the flow in karstic conduits. The simulations show that the conduit permeability ranges from 5?×?10^?6 to 5.5?×?10^?5?m² and the limestone permeability ranges from 8?×?10^?11 to 6?×?10^?10?m². The dispersivity coefficient ranges from 23 to 53 m in the conduits and from 1 to 5 m in the limestone. The results of the simulations carried out using Darcy’s law in the conduits show that the dispersion towards the fractures is underestimated.     

First results from the North Iceland experiment

Between June 2004 and September 2004 a temporary seismic network was installed on the northern insular shelf of Iceland and onshore in north Iceland. The seismic setup aimed at resolving the subsurface structure and, thus, the geodynamical transition from Icelandic crust to typical oceanic crust along Kolbeinsey Ridge. The experiment recorded about 1,000 earthquakes. The region encloses the Tjörnes Fracture Zone containing the Husavik–Flatey strike-slip fault and the extensional seismic Grimsey Lineament. Most of the seismicity occurs in swarms offshore. Preliminary results reveal typical mid-ocean crust north of Grimsey and a heterogeneous structure with major velocity anomalies along the seismic lineaments and north–south trending subsurface features. Complementary bathymetric mapping highlight numerous extrusion features along the Grimsey Lineament and Kolbeinsey Ridge. The seismic dataset promises to deliver new insights into the tectonic framework for earthquakes in an extensional transform zone along the global mid-ocean ridge system.     

Structure of the upper mantle in the Northeastern Atlantic close to the azores-gibraltar ridge from surface-wave and body-wave observations

Phase velocities of teleseismic Rayleigh waves have been measured in the central North Atlantic on both sides of the Azores-Gibraltar Ridge (AGR) by means of a specially designed long-period station network. The dispersion data obtained were regionalized and then subjected to a “hedgehog” inversion, which gives a set of upper mantle models compatible with the observational data within specified error bounds.Reasonable model solutions were selected by using regional body-wave observations, such as P_n- and S_n-wave velocities determined from earthquakes along the AGR. The S_(itn) velocities measured indicate that the shear-wave velocity in the mantle part of the lithosphere is much higher on the northern side of the AGR. Strongly negative P-wave residuals in this area indicate faster seismic propagation than implied by the Jeffreys-Bullen travel-time tables, while propagation is much slower in the Gulf of Cadiz area. Furthermore the residuals show a clear difference for paths through oceanic and continental domains and suggest that the transition between these two domains extends much further into the ocean on the southern side of the AGR than on the northern side.The proposed model for the structure of the upper mantle in that region shows that there exists a pronounced velocity contrast across the AGR. Thickening of the lithospheric plate with increasing plate age is indicated to the south of the ridge. The greatest thickness is reached close to the continental margin within a zone about 500 km wide, whose velocity close to the Canary Islands and Madeira is significantly lower, probably due to the well-known volcanic activity there. These observations together with the travel time residuals reveal that this zone seems to be of a transitional nature somewhere between a continental and oceanic structure.     

The Relationship of Green-Line Transients to White-Light Coronal Mass Ejections

We report observations by the Large Angle Spectrometric Coronagraph (LASCO) on the SOHO spacecraft of three coronal green-line transients that could be clearly associated with coronal mass ejections (CMEs) detected in Thomson-scattered white light. Two of these events, with speeds >25 km s^-1, may be classified as ‘whip-like’ transients. They are associated with the core of the white-light CMEs, identified with erupting prominence material, rather than with the leading edge of the CMEs. The third green-line transient has a markedly different appearance and is more gradual than the other two, with a projected outward speed <10 km s^-1. This event corresponds to the leading edge of a ‘streamer blowout’ type of CME. A dark void is left behind in the emission-line corona following each of the fast eruptions. Both fast emission-line transients start off as a loop structure rising up from close to the solar surface. We suggest that the driving mechanism for these events may be the emergence of new bipolar magnetic regions on the surface of the Sun, which destabilize the ambient corona and cause an eruption. The possible relationship of these events to recent X-ray observations of CMEs is briefly discussed. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1004981125702     

Geomorphological models as a tool for environmental studies

Relief is understood as a fundamental component of the environmental system. The structure of geomorphodynamic processes and of the relief sphere as the main energy transformation surface is a determining factor for environmental conditions of both sites and areas. Modelling is the instrument of the applied and systems analytic approach to geomorphology. Geomorphological tools are appropriate means of determining, structuring, assessing and predicting the geoecological potential of the environment. This is demonstrated by quantitative and multivariate models to describe spatial patterns of the environment (fluvial network analysis) and to evaluate and present environmentally relevant geomorphological disposition (slope stability analysis).     

Application of the Canadian regional climate model to the Laurentian great lakes region: Implementation of a lake model

Abstract

This study reports on the implementation of an interactive mixed‐layer/thermodynamic‐ice lake model coupled with the Canadian Regional Climate Model (CRCM). For this application the CRCM, which uses a grid mesh of 45 km on a polar stereographic projection, 10 vertical levels, and a timestep of 15 min, is nested with the second generation Canadian General Circulation Model (GCM) simulated output. A numerical simulation of the climate of eastern North America, including the Laurentian Great Lakes, is then performed in order to evaluate the coupled model. The lakes are represented by a “mixed layer” model to simulate the evolution of the surface water temperature, and a thermodynamic ice model to simulate evolution of the ice cover. The mixed‐layer depth is allowed to vary spatially. Lake‐ice leads are parametrized as a function of ice thickness based on observations. Results from a 5‐year integration show that the coupled CRCM/lake model is capable of simulating the seasonal evolution of surface temperature and ice cover in the Great Lakes. When compared with lake climatology, the simulated mean surface water temperature agrees within 0.12°C on average. The seasonal evolution of the lake‐ice cover is realistic but the model tends to underestimate the monthly mean ice concentration on average. The simulated winter lake‐induced precipitation is also shown, and snow accumulation patterns on downwind shores of the lakes are found to be realistic when compared with observations.     

New brown dwarf candidates in the Pleiades

We have performed deep, wide‐field imaging on a ∼0.4 deg² field in the Pleiades (Melotte 22). The selected field was not yet target of a deep search for low mass stars and brown dwarfs. Our limiting magnitudes are R ∼ 22 mag and I ∼ 20 mag, sufficient to detect brown dwarf candidates down to 40 M_J. We found 197 objects, whose location in the (I, R – I) color magnitude diagram is consistent with the age and the distance of the Pleiades. Using CTK R and I as well as JHK photometry from our data and the 2MASS survey we were able to identify 7 new brown dwarf candidates. We present our data reduction technique, which enables us to resample, calibrate, and co‐add many images by just two steps. We estimate the interstellar extinction and the spectral type from our optical and the NIR data using a two‐dimensional χ² fitting (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)