The recording of floods and earthquakes in Lake Chichój,Guatemala during the twentieth century

Laguna Chichój (Lake Chichój) is the only deep permanent lake in the central highlands of Guatemala. The lake is located in the boundary zone between the North American and Caribbean plates. The lake has been struck by devastating earthquakes and tropical cyclones in historical times. We investigated the imprint of twentieth century extreme events on the sedimentary record of this tropical lake using a bathymetric survey of the lake, coring the lake floor, and providing a chronology of sediment accumulation. The lake occupies a series of circular depressions likely formed by the rapid dissolution of a buried body of gypsum. ²¹⁰Pb and ¹³⁷Cs inventories and varve counting indicate high rates of sedimentation (1–2 cm year^?1). The annually layered sediment is interrupted by turbidites of two types: a darker-colored turbidite, enriched in lake-derived biogenic constituents, and interpreted as a seismite, and a lighter-colored type, enriched in catchment-derived constituents, interpreted as a flood layer. Comparison of our ¹³⁷Cs-determined layer ages with a catalog of twentieth century earthquakes shows that an earthquake on the Motagua fault in 1976 generated a conspicuous darker-colored turbidite and slumped deposits in separate parts of the lake. The entire earthquake inventory further reveals that mass movements in the lake are triggered at Modified Mercalli Intensities higher than V. Tropical cyclonic depressions known to have affected the lake area had limited effect on the lake, including Hurricane Mitch in 1998. One storm however produced a significantly thicker flood layer in the 1940s. This storm is reportedly the only event to have generated widespread slope failures in the lake catchment. It is thus inferred that abundant landsliding provided large amounts of concentrated sediment to the lake, through hyperpycnal flows.     

Numerical modelling of strain in lava tubes

Received: 5 November 1998 / Accepted: 17 November 1999     

Numerical simulation of the propagation of P waves in fractured media

We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q ⁻¹ of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q ⁻¹ of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) ⁻¹ in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.     

Offshore experiments on styrene spillage in marine waters for risk assessment

Within the context of risk evaluation of chemical spillages into the marine environment, this paper reports on an offshore experiment to study the behaviour of styrene spilled into sea under natural conditions and discusses theoretical approaches. Floating structures were used to enclose the spillage and the gaseous cloud formation, and dissolution processes were in situ monitored. The identification of spill risks for man and marine environment through GESAMP's hazard profile is described for styrene: Styrene is rated as a chemical with a significant health hazard that will float but also evaporate. However, monitoring of the water column in the experiments showed that the concentration of styrene in water during the first hour represents 50% of the product spilled. For the potentially exposed public, the GESAMP hazard rating recommends the closure of beaches and evacuation. The risk assessment developed from experimental data confirms this safety advice.     

Effect of the processing methodology on satellite altimetry-based global mean sea level rise over the Jason-1 operating period

Determining how the global mean sea level (GMSL) evolves with time is of primary importance to understand one of the main consequences of global warming and its potential impact on populations living near coasts or in low-lying islands. Five groups are routinely providing satellite altimetry-based estimates of the GMSL over the altimetry era (since late 1992). Because each group developed its own approach to compute the GMSL time series, this leads to some differences in the GMSL interannual variability and linear trend. While over the whole high-precision altimetry time span (1993–2012), good agreement is noticed for the computed GMSL linear trend (of $3.1\pm 0.4$  mm/year), on shorter time spans (e.g., ${<}10~\hbox {years}$ ), trend differences are significantly larger than the 0.4 mm/year uncertainty. Here we investigate the sources of the trend differences, focusing on the averaging methods used to generate the GMSL. For that purpose, we consider outputs from two different groups: the Colorado University (CU) and Archiving, Validation and Interpretation of Satellite Oceanographic Data (AVISO) because associated processing of each group is largely representative of all other groups. For this investigation, we use the high-resolution MERCATOR ocean circulation model with data assimilation (version Glorys2-v1) and compute synthetic sea surface height (SSH) data by interpolating the model grids at the time and location of “true” along-track satellite altimetry measurements, focusing on the Jason-1 operating period (i.e., 2002–2009). These synthetic SSH data are then treated as “real” altimetry measurements, allowing us to test the different averaging methods used by the two processing groups for computing the GMSL: (1) averaging along-track altimetry data (as done by CU) or (2) gridding the along-track data into $2^{\circ }\times 2^{\circ }$ meshes and then geographical averaging of the gridded data (as done by AVISO). We also investigate the effect of considering or not SSH data at shallow depths $({<}120~\hbox {m})$ as well as the editing procedure. We find that the main difference comes from the averaging method with significant differences depending on latitude. In the tropics, the $2^{\circ }\times 2^{\circ }$ gridding method used by AVISO overestimates by 11 % the GMSL trend. At high latitudes (above $60^{\circ }\hbox {N}/\hbox {S}$ ), both methods underestimate the GMSL trend. Our calculation shows that the CU method (along-track averaging) and AVISO gridding process underestimate the trend in high latitudes of the northern hemisphere by 0.9 and 1.2 mm/year, respectively. While we were able to attribute the AVISO trend overestimation in the tropics to grid cells with too few data, the cause of underestimation at high latitudes remains unclear and needs further investigation.     

Dating carbonate rocks with in-situ produced cosmogenic Be: Why it often fails

Surface exposure dating of carbonate rocks using cosmogenic ¹⁰Be is problematic. We have performed step-wise leaching of calcite-rich samples in order to investigate the reasons for this. Results on different grain size fractions clearly indicate the source of atmospheric ¹⁰Be is clay. We demonstrate that partial-leaching procedures, which result in moderate pH levels will not release ¹⁰Be (in-situ produced or atmospheric) due to the instant re-absorption on grain surfaces. By contrast, under strongly acidic conditions, all absorbed ¹⁰Be is leached from aluminosilicates giving abnormally high ¹⁰Be concentrations and consequently exposure ages that are too old. Dating is only possible if samples do not contain any clay minerals or if they can be removed prior to carbonate dissolution.     

Rapid and cyclic aeolian deposition during the Last Glacial in European loess: a high-resolution record from Nussloch, Germany

This paper reports the results of an investigation of the Weichselian Upper Pleniglacial loess sequences of Nussloch (Rhine Valley, Germany) based on stratigraphy, palaeopedology, sedimentology, palynology, malacology and geochemistry (δ¹³C), supported by radiocarbon, TL and OSL dating. Grain-size and magnetic susceptibility records are taken at 5 cm intervals from the Upper Pleniglacial (UPG) loess. The data indicate cyclic variations in loess deposition between ca 34 and 17 ka, when the sedimentation rate is especially high (1.0–1.2 m per ka for more than 10 m). The grain-size index (GSI: ratio of coarse silt versus fine silt and clay) shows variations, which are assumed to be an indirect measurement of wind intensity. The sedimentation rate, interpreted from the profiles, indicates high values in loess (Loess events LE-1 to LE-7) and low or negligible values in tundra gley horizons G1 to G8. OSL ages from the loess and ¹⁴C dates from organic matter in the loess show that loess deposition was rapid but was interrupted by shorter periods of reduced aeolian sedimentation. Comparison between the data from Nussloch and other European sequences demonstrates a progressive coarsening of the loess deposits between ca 30 and 22 ka. This coarsening trend ends with a short but major grain-size decrease and is followed by an increase to a new maximum at 20 ± 2 ka (“W” shape). Correlation between the loess GSI and the Greenland ice-core dust records, suggests a global connection between North Atlantic and Western European global atmospheric circulation and wind regimes. In addition, the typical Upper Pleniglacial loess deposition begins at ca 30–31 ka, close to Heinrich event (HE) 3, and the main period of loess sedimentation at about 25 ± 2 ka is coeval to HE 2. Correlation of magnetic susceptibility and grain-size records shows that the periods, characterised by high GSI, coincide with an increase in the amount of ferromagnetic minerals reworked from the Rhine alluvial plain. They suggest enhancement in the frequency of the storms from N–NW. These results are integrated within a palaeogeographical model of dust transport and deposition in Western Europe for the Weichselian Upper Pleniglacial (or Late Pleniglacial).     

Experimental survey of the hydrodynamic performance of a small spar buoy

The objective in this experimental work is to evaluate the capability of several cylindrical buoys to follow and measure waves. Eleven configurations have been studied, eight among which were furnished with disks of different diameters at the waterline. The experiments took place in the test basin at ISITV. On the one hand, we have obtained a set of transfer coefficients for each frequency and each configuration in regular waves. On the other hand, we have determined the corresponding transfer functions, then we have used the latter to measure irregular waves. The time signals as measured by the buoys with and without corrections have been compared with wave gage measurements and subjected to a wave-by-wave analysis. The results allow the conclusion that the buoy with a medium size disk constitutes the best compromise. In fact, it permits the buoy to follow the free surface while minimizing parasitic motions.