Low-cost volcano surveillance from space: case studies from Etna, Krafla, Cerro Negro, Fogo, Lascar and Erebus

 Satellite data offer a means of supplementing ground-based monitoring during volcanic eruptions, especially at times or locations where ground-based monitoring is difficult. Being directly and freely available several times a day, data from the advanced very high resolution radiometer (AVHRR) offers great potential for near real-time monitoring of all volcanoes across large (3000×3000 km) areas. Herein we describe techniques to detect and locate activity; estimate lava area, thermal flux, effusion rates and cumulative volume; and distinguish types of activity. Application is demonstrated using data for active lavas at Krafla, Etna, Fogo, Cerro Negro and Erebus; a pyroclastic flow at Lascar; and open vent systems at Etna and Stromboli. Automated near real-time analysis of AVHRR data could be achieved at existing, or cheap to install, receiving stations, offering a supplement to conventional monitoring methods. Received: 21 January 1997 / Accepted: 3 April 1997     

Oxidation of pyrite during extraction of carbonate associated sulfate

The sulfur isotopic composition of carbonate associated sulfate (CAS) has been used to investigate the geochemistry of ancient seawater sulfate. However, few studies have quantified the reliability of δ³⁴S of CAS as a seawater sulfate proxy, especially with respect to later diagenetic overprinting. Pyrite, which typically has depleted δ³⁴S values due to authigenic fractionation associated with bacterial sulfate reduction, is a common constituent of marine sedimentary rocks. The oxidation of pyrite, whether during diagenesis or sample preparation, could thus adversely influence the sulfur isotopic composition of CAS. Here, we report the results of CAS extractions using HCl and acetic acid with samples spiked with varying amounts of pyrite. The results show a very strong linear relationship between the abundance of fine-grained pyrite added to the sample and the resultant abundance and δ³⁴S value of CAS. This data represents the first unequivocal evidence that pyrite is oxidized during the CAS extraction process. Our mixing models indicate that in samples with much less than 1 wt.% pyrite and relatively high δ³⁴S_pyrite values, the isotopic offset imparted by oxidation of pyrite should be much less than ? 4‰. A wealth of literature exists on the oxidation of pyrite by Fe³⁺ and we believe this mechanism drives the oxidation of pyrite during CAS extraction, during which the oxygen used to form sulfate is taken from H₂O, not O₂. Consequently, extracting CAS under anaerobic conditions would only slow, but not halt, the oxidation of pyrite. Future studies of CAS should attempt to quantify pyrite abundance and isotopic composition.     

Quantitative sinkhole hazard assessment. A case study from the Ebro Valley evaporite alluvial karst (NE Spain)

Quantitative sinkhole hazard assessments in karst areas allow calculation of the potential sinkhole risk and the performance of cost-benefit analyses. These estimations are of practical interest for planning, engineering, and insurance purposes. The sinkhole hazard assessments should include two components: the probability of occurrence of sinkholes (sinkholes/km² year) and the severity of the sinkholes, which mainly refers to the subsidence mechanisms (progressive passive bending or catastrophic collapse) and the size of the sinkholes at the time of formation; a critical engineering design parameter. This requires the compilation of an exhaustive database on recent sinkholes, including information on the: (1) location, (2) chronology (precise date or age range), (3) size, and (4) subsidence mechanisms and rate. This work presents a hazard assessment from an alluvial evaporite karst area (0.81 km²) located in the periphery of the city of Zaragoza (Ebro River valley, NE Spain). Five sinkholes and four locations with features attributable to karstic subsidence where identified in an initial investigation phase providing a preliminary probability of occurrence of 0.14 sinkholes/km² year (11.34% in annual probability). A trenching program conducted in a subsequent investigation phase allowed us to rule out the four probable sinkholes, reducing the probability of occurrence to 0.079 sinkholes/km² year (6.4% in annual probability). The information on the severity indicates that collapse sinkholes 10–15 m in diameter may occur in the area. A detailed study of the deposits and deformational structures exposed by trenching in one of the sinkholes allowed us to infer a modern collapse sinkhole approximately 12 m in diameter and with a vertical throw of 8 m. This collapse structure is superimposed on a subsidence sinkhole around 80 m across that records at least 1.7 m of synsedimentary subsidence. Trenching, in combination with dating techniques, is proposed as a useful methodology to elucidate the origin of depressions with uncertain diagnosis and to gather practical information with predictive utility about particular sinkholes in alluvial karst settings: precise location, subsidence mechanisms and magnitude, and timing and rate of the subsidence episodes.     

Horizontal spatial and temporal distribution patterns of nearshore larval fish assemblages at a temperate rocky shore

There have been no previous studies of the composition of nearshore larval fish assemblages along the coast of Portugal. We aimed to describe the composition and horizontal distribution patterns of larval fish assemblages and their temporal dynamics near a rocky reef at depths shallower than 13 m (inshore) and at two miles (3.70 km) from shore (offshore), as well as along transects perpendicular to the shoreline, from the reef to 10 miles offshore (18.52 km). Samples were taken using 5 min sub-surface trawls at the rocky shore of the Arrábida Marine Park (W Portugal). A total of 1021 larvae were collected, belonging to 61 taxa inshore and to 29 taxa offshore. Along transects, 626 larvae of 52 taxa were collected. Most larvae belonged to coastal species associated with rocky reefs. Total larval abundance and diversity were higher from May to July, which is consistent with the spawning activity of adults. Diversity and total larval abundance decreased significantly with increasing distance from shore, both in the inshore/offshore comparison and in the transects, where this decrease was evident at a very small spatial scale (within the first mile from the reef). Species assemblages differed in the pattern of distribution, with most species clearly associated to the extreme nearshore. The distribution patterns obtained were independent of the spawning mode of species. Results are discussed in the light of the possible physical mechanisms that can potentially act at the Arrábida Marine Park to facilitate larvae retention and the role of larval behaviour.     

Relationships Between Community Structure of the Intertidal Macroinfauna and Sandy Beach Characteristics Along the Chilean Coast

Abstract. Eight sandy beaches were seasonally sampled along the coast of Chile, from ca. 21 to 42° S (about 3000 km) to study the relationship between community structure of the intertidal macroinfauna and beach characteristics. Sediment samples (0.1 m², 30 cm deep) were collected (July – September 1998 and December 1998 – January 1999) with plastic cylinders at 15 equally spaced levels along three replicated transects extending from above the drift line to the swash zone. The sediment was sieved through a 1 mm mesh and the organisms collected stored in 5 % formalin. To define beach types, Dean's parameter (Ω) was calculated from wave heights and periods, and fall velocity of sand particles from the swash zone. Crustaceans (mainly peracarids) were the most diverse group with 14 species, followed by polychaetes with 5 species. The talitrid amphipod Orchestoidea tuberculata , the cirolanid isopods Excirolana braziliensis and E. hirsuticauda and the anomuran decapod Emerita analoga were the most widely distributed and common species. Regression analyses between species richness, abundance and biomass of the whole macroinfauna versus sediment characteristics, beach face slopes and morphodynamic beach states showed no significant relationships. Thus, macroinfaunal community characteristics did not increase linearly from lower intermediate to higher intermediate or dissipative beach states as had been found before in Chile or in other coasts. A comparative analysis with data from sandy beaches of other world regions showed that the number of species inhabiting Chilean sandy beaches was generally lower, whereas total population abundances were generally higher compared with values reported elsewhere.     

On the assessment of hydrodynamic coefficients of cylinders in heaving

A detailed comparison is made of the common two-dimensional methods for the evaluation of the hydrodynamic coefficients. Special attention was given to the conformal transformation method and the parameters which influence the accuracy of this method are studied. The results were compared with the Frank close-fit method with respect to accuracy and time taken. Finally the correlation between the geometric and hydrodynamic errors is examined and some relations are proposed.     

Inertial and hydrodynamic inertial loads on floating units

The analysis of the dynamic behavior of floating units usually employs a coordinate system with origin in the unit's center of gravity, which significantly simplifies the global mass matrix. Hydrodynamic coefficients are then computed considering the same coordinate system. However, to analyze other conditions of mass distribution and maintain the simplicity of a global mass matrix, it is necessary to determine again the hydrodynamic coefficients, thereby reducing the efficiency of the entire process. Another important point is that the geometries frequently used in floating units are such that the cross-terms of an added mass are relatively unimportant when compared with the main terms, and it is, therefore, common to use only some of them to analyze the unit's dynamic behavior. Recently, however, in the search for production systems suitable for water depths greater than 3000 m, other geometries have been considered in technical and economic feasibility studies. It is possible that for these new geometries all terms of the added mass matrix must be included in the analysis. This paper presents the full development used to determine the complete global mass matrix, the inertial and hydrodynamic inertial loads that make use of the added mass matrices considering any coordinate system and the six degrees of freedom, including all cross-terms.     

Dynamical conditions in planets and stars

With the available data in planets, stars and galaxies, it is studied the functions of angular momentaJ(M) and amounts of actionA _c(M) (associated to the non rotational terms in the kinetic energy). The results indicate that independently of how are these functionsJ(M),A _c(M) their ratioA _c/J remains a near invariant. It is independent also from the type of angular momenta: intrinsic spins of the bodies or the total angular (orbital) momenta of the bodies forming a system; for instance, the Solar System and the planets.The relationA _c(M) for the Solar System are analogous to these in the FGK stars of the main sequence, and the relationJ(M) (also for the Solar System) is analogous to the lower possible limit for binary stars.The different types of binary stars from the short period, detached systems to contactary systems, gives a range of functionsJ(M),A _c(M) that are the same that one can expect in stars with planetary systems. According to the detection limits given for planetary companions by Campbell, Walker and Yang (1988) (masses of less than 9 Jupiter masses and orbital periods of less than 50 years) we calculate the limits forJ(M) andA _c(M) This gives a lower limitA _c/J 1 associated to stars with planetary systems as 61 Cygni and to short period detached binaries. The upper limitA _c/J 16 correspond to planetary systems as the ours and probably to cataclysmic binaries. There are reasons to suspect that systems as the ours and in range 4 A _c/J 16 (with a lower limit analogous to contactary binaries as Algols and W Ursa Majoris) must be the most common type of planetary systems. The analogies with the functionsJ(M)A _c(M) for galaxies suggest cosmogonical conditions in the stellar formation.Independently of this, one can have boundary conditions for the Jacobi problem when applied to a collapsing cloud. Namely, from the initial stage (a molecular cloud) to the final stage (a formed stellar system: binary or planetary) the angular momenta and amounts of action decayed to 10^~4 the initial values, but in such a form thatA _c(t)/J(t) remains a near invariant.     

Coupling an ocean wave model to an atmospheric general circulation model

A coupled model, consisting of an ocean wave model and an atmospheric general circulation model (AGCM), is integrated under permanent July conditions. The wave model is forced by the AGCM wind stress, whereas the wind waves modify the AGCM surface fluxes of momentum, sensible and latent heat. We investigate the following aspects of the coupled model: how realistic are the wave fields, how strong is the coupling, and how sensitive is the atmospheric circulation to the spatially and temporally varying wave field. The wave climatology of the coupled model compares favorably with observational data. The interaction between the two models is largest (although weak) in the storm track in the Southern Hemisphere. Young windsea, which is associated with enhanced surface fluxes is generated mostly in the equatorward frontal area of an individual cyclone. However, the enhancement of the surface fluxes is too small to significantly modify the climatological mean atmospheric circulation.This paper was presented at the Second International Conference on Modelling of Global Climate Variability, held in Hamburg 7–11 September 1992 under the auspices of the Max Planck Institute for Meteorology. Guest Editor for these papers is L. Dümenil