Measuring and interpreting heat fluxes from shallow volcanic bodies using vertical temperature profiles: a preliminary test

 To test the potential of heat flux prospecting in active volcanic areas using shallow temperature data taken along vertical profiles, we carried out two thermal profile surveys, one not far from Yasur cone on Tanna Island, and another inside the caldera of Ambrym (New Hebrides arc, southwestern Pacific). The basic steady heat flux of internal volcanic origin was determined, taking into account both conductive and convective heat transfers. At both locations there exists, over small distances, significant differences in the heat flux. These differences correspond to shallow sources of heat. The use of a network of vertical profiles allowed: (a) heat flux mapping; (b) location of shallow volcanic heat sources; and (c) observation of the detailed structure of the heat release at quiescent but active volcanoes. Received: 18 July 1997 / Accepted: 13 May 1998     

Mapping landslide susceptibility from small datasets: A case study in the Pays de Herve (E Belgium)

A landslide susceptibility map is proposed for the Pays de Herve (E Belgium), where large landslides affect Cretaceous clay outcrop areas. Based on a Bayesian approach, this GIS-supported probabilistic map identifies the areas most susceptible to deep landslides. The database is comprised of the source areas of ten pre-existing landslides (i.e. a sample of 154 grid cells) and of six environmental data layers, namely lithology, proximity to active faults, slope angle and aspect, elevation and distance to the nearest valley-floor. A 30-m-resolution DEM from the Belgian National Geographical Institute is used for the analysis. Owing to the small size of the sample, a special cross-validation procedure of the susceptibility map is performed, which uses in an iterative way each of the landslides to test the predictive power of the map derived from the other landslides. Four different sets of variables are used to produce four susceptibility maps, whose prediction curves are compared. While the prediction rates associated with the models not involving the “proximity to active fault” criterion are comparable to those of the models considering this variable, strong weaknesses inherent in the fault data on which the latter rely suggest that the final susceptibility map should be based on a model that excludes any reference to fault. This highlights the difference between a triggering factor and determining factors, and in the same time broadens the scope of the produced map. A single reactivated slide is also used to test the possibility of predicting future reactivation of existing landslides in the area. Finally, the need for geomorphological control over the mathematical treatment is underlined in order to obtain realistic prediction maps.     

Twin structure in coesite studied by high resolution electron microscopy

A synthetic crystal of coesite, SiO₂, has been found to contain microtwins. The (100) reflection twin has been made evident on incoherent (010) planes. The twin structure is directly imaged by high resolution electron microscopy at 2,5 Å resolution, allowing an atomic model to be proposed. This twin could be incoherent on large surfaces and the dislocation density is predicted. The X-ray crystallographic data are analyzed critically in view of the presence of such defects in single crystals. The mutual orientation between a silicon matrix and coesite precipitates in Czochralski-silicon annealed specimens is proposed to be non unique as a consequence of the easy twinning of coesite.     

Multi-scale aquifer characterization and groundwater flow model parameterization using direct push technologies

Direct push (DP) technologies are typically used for cost-effective geotechnical characterization of unconsolidated soils and sediments. In more recent developments, DP technologies have been used for efficient hydraulic conductivity (K) characterization along vertical profiles with sampling resolutions of up to a few centimetres. Until date, however, only a limited number of studies document high-resolution in situ DP data for three-dimensional conceptual hydrogeological model development and groundwater flow model parameterization. This study demonstrates how DP technologies improve building of a conceptual hydrogeological model. We further evaluate the degree to which the DP-derived hydrogeological parameter K, measured across different spatial scales, improves performance of a regional groundwater flow model. The study area covers an area of ~60 km² with two overlying, mainly unconsolidated sand aquifers separated by a 5–7 m thick highly heterogeneous clay layer (in north-eastern Belgium). The hydrostratigraphy was obtained from an analysis of cored boreholes and about 265 cone penetration tests (CPTs). The hydrogeological parameter K was derived from a combined analysis of core and CPT data and also from hydraulic direct push tests. A total of 50 three-dimensional realizations of K were generated using a non-stationary multivariate geostatistical approach. To preserve the measured K values in the stochastic realizations, the groundwater model K realizations were conditioned on the borehole and direct push data. Optimization was performed to select the best performing model parameterization out of the 50 realizations. This model outperformed a previously developed reference model with homogeneous K fields for all hydrogeological layers. Comparison of particle tracking simulations, based either on the optimal heterogeneous or reference homogeneous groundwater model flow fields, demonstrate the impact DP-derived subsurface heterogeneity in K can have on groundwater flow and solute transport. We demonstrated that DP technologies, especially when calibrated with site-specific data, provide high-resolution 3D subsurface data for building more reliable conceptual models and increasing groundwater flow model performance.     

Late Quaternary kinematics of the Pallatanga strike-slip fault (Central Ecuador) from topographic measurements of displaced morphological features

The northeast-trending Pallatanga right-lateral strike-slip fault runs across the Western Cordillera connecting N50E-N70E trending normal faults in the Gulf of Guayaquil with N-S reverse faults in the Interandean Depression. Over most of its length, the fault trace has been partly obscured by erosional processes and can be inferred in the topography only at the large scale. Only the northern fault segment, which follows the upper Rio Pangor valley at elevations above 3600 m, is prominent in the morphology. Valleys and ridges cut and offset by the fault provide an outstanding record of right-lateral cumulative fault displacement. The fault geometry and kinematics of this particular fault segment can be determined from detailed topographic levellings. The fault strikes N30E and dips 75 to the NW. Depending on their size and nature, transverse morphological features such as tributaries of the Rio Pangor and intervening ridges, reveal right-lateral offsets which cluster around 27 ± 11m, 41.5 ± 4 m, 590 ± 65 m and 960 ± 70 m. The slip vector deduced from the short-term offsets shows a slight reverse component with a pitch of about 11.5 SW. The 41.5 ± 4 m displacements are assumed to be coeval with the last glacial termination, yielding a mean Holocene slip-rate of 2.9- 4.6 mm yr⁻¹. Assuming a uniform slip rate on the fault in the long term, the 27 m offset appears to correlate with an identified middle Holocene morphoclimatic event, and the long term offsets of 590 m and 960 m coincide with the glacial terminations at the beginning of the last two interglacial periods.     

A new method for sampling and analyzing volcanic sublimates — application to Merapi volcano, Java

A new method for the sampling of sublimates from high-temperature volcanic gases has been used at Merapi volcano, Java, in 1978. The sublimates were collected on the inner walls of silica tubes introduced into fumarolic vents. Volcanic gases were allowed to move freely through the tubes and as they cooled, a fraction of the volatile components condensed on the inner walls of the tubes along the temperature gradient. The sublimates were then analyzed by a combination of light microscopy, scanning electron microscopy, electron microprobe and X-ray diffraction.Six successive zones of different compositions and mineralogical associations have been identified along the covered range of temperatures (900° to around 400°C). From the high to the low temperatures, these zones are composed of: (1) cristobalite, magnetite, hercynite; (2) molybdenite; (3) acmite; (4) halite, sylvite; (5) sphalerite, pyrite; and (6) galena. Equilibrium calculations show that these crystalline phases are stable for pS₂, pC₁, and pO₂, values typical of magma-buffered gases that have not been contaminated by atmospheric oxygen.The deposits observed in the tubes may be useful in aiding the understanding of the mechanisms acting during the cooling of the gaseous phase on its way to the surface and before its emission into the atmosphere.     

Tar ball frequency data and analytical results from a long-term beach monitoring program

Following the spill of fuel oils from the New Carissa in February 1999, approximately 300 km of beaches on the Pacific coast of North America were surveyed. A long-term observation program focused on the documentation of stranded tar balls in the vicinity of the spill site. Systematic beach surveys which were conducted over the period March 1999 to April 2001 and semi-logarithmic scale, time-series plots proved the most useful format for identifying trends. Beach monitoring continued through to August 2001. by which time 212 tar balls had been analyzed by GC/MS for their chemical characteristics. The samples of tar balls collected between February 1999 and August 2001 were qualitatively compared with New Carissa source oils (NCSO) and 101 (48%) were not consistent with NSCO. The presence of tar balls that are not related to an incident can confound attempts to define cleanup or endpoint criteria and to assess possible injury to natural resources.     

Post-collision neogene volcanism of the Eastern Rif (Morocco): magmatic evolution through time

Neogene volcanism in the Eastern Rif (Morocco) comprises a series of calc-alkaline, potassic calc-alkaline, shoshonitic and alkaline volcanic rocks. According to new stratigraphical, along with new and previous chronological and geochemical data, the orogenic volcanism was successively (1) calc-alkaline (basaltic andesites and andesites: 13.1 to 12.5 Ma, rhyolites: 9.8 Ma), (2) K-calc-alkaline (basaltic andesitic to rhyolitic lavas and granodiorites: 9.0 to 6.6 Ma), and (3) shoshonitic (absarokites, shoshonites, latites, trachytes: 7.0 to 5.4 Ma). The later Pliocene volcanism was basaltic and alkaline (5.6 to 1.5 Ma). The calc-alkaline and K-calc-alkaline series exhibit lower K₂O (0.7–5.3 wt.%), Nb (8–19 ppm) contents and higher ⁸⁷Sr/⁸⁶Sr (0.70773–0.71016) than the shoshonitic series (K₂O: 2.4–7.2 wt.%, Nb: 21–38 ppm, ⁸⁷Sr/⁸⁶Sr: 0.70404–0.70778). Pliocene alkaline basalts have a sodic tendency (Na₂O/K₂O: 1.7–3.5), high Nb content (up to 52 ppm), and low ⁸⁷Sr/⁸⁶Sr ratio (0.70360–0.70413). The variations through time of K₂O, Nb and Sr isotopic ratio reflect different mantle sources: (i) calc-alkaline, potassic calc-alkaline and shoshonitic series are derived from a mantle source modified by older subduction, (ii) alkaline basalts are derived mainly from an enriched mantle source. Through time, incompatible elements such as Nb increased while ⁸⁷Sr/⁸⁶Sr decreased, suggesting a decreasing influence of metasomatized mantle (inherited subduction). Such evolution is related to the post-collision regimes operating in this area, and could be linked to the succession of extensional, compressional and strike-slip fault tectonics.     

Cratering and modification of wet‐target craters: Projectile impact experiments and field observations of the Lockne marine‐target crater (Sweden)

Abstract— Marine impacts are one category of crater formation in volatile targets. At target water depths exceeding the diameter of the impactor, the zones of vaporization, melting, and excavation of the standard land‐target cratering model develop partially or entirely in the water column. The part of the crater that has a potential of being preserved (seafloor crater) may to a great extent be formed by material emplacement and excavation processes that are very different from land‐target craters. These processes include a high‐energy, water‐jet‐driven excavation flow. At greater water depths, the difference in strength of the target layers causes a concentric crater to evolve. The crater consists of a wide water cavity with a shallow excavation flow along the seabed surrounding a nested, deeper crater in the basement. The modification of the crater is likewise influenced by the water through its forceful resurge to fill the cavity in the water mass and the seafloor. The resurge flow is strongly erosive and incorporates both ejecta and rip‐up material from the seabed surrounding the excavated crater. A combination of field observations and impact experiments has helped us analyze the processes affecting the zone between the basement crater and the maximum extent of the water cavity. The resurge erosion is facilitated by fragmentation of the upper parts of the solid target caused by a) spallation and b) vibrations from the shallow excavation flow and, subsequently, c) the vertical collapse of the water cavity rim wall. In addition, poorly consolidated and saturated sediments may collapse extensively, possibly aided by a violent expansion of the pore water volume when it turns into a spray during passage of the rarefaction wave. This process may also occur at impacts into water‐saturated targets without an upper layer of seawater present. Our results have implications for impacts on both Earth and Mars, and possibly anywhere in the solar system where volatiles exist/have existed in the upper part of the target.     

The 15 m.y. geomagnetic reversal periodicity: a quantitative test

We have analyzed the three most commonly used geomagnetic polarity time scales that cover the past 100 m.y.. In the three cases, different spectral analyses show a 13–16 m.y. periodicity in the rate of reversal occurrence. To test whether this periodicity is real or simply arises from a random generator we have compared these polarity time scales with a large number of synthetic sequences produced by a random process, characterized by a linear time variation of its mean activity. Geomagnetic and generated sequences were regularly sampled by using sliding windows, anf then the Fourier spectra of the obtained frequency signals were compared. This test shows that the detected periodicity is presumably not a simple statistical fluctuation of an aperiodic generator, and consequently that a long-term periodicity in the geodynamo must be seriously considered.