Gravitational effect of water circulation in the northwest Yunnan

ＧｒａｖｉｔａｔｉｏｎａｌｅｆｆｅｃｔｏｆｗａｔｅｒｃｉｒｃｕｌａｔｉｏｎｉｎｔｈｅｎｏｒｔｈｗｅｓｔＹｕｎｎａｎＭｉｎ－ＹｕＪＩＡ；（贾民育）Ｓｈａｏ－ＡｎＳＵＮ；（孙少安）Ａｉ－ＭｉｎＸＩＡＮＧ；（项爱民）ａｎｄＤｏｎｇ－ＺｈｉＬＩＵ（刘冬至）（...     

Correction to the crustal thickness of the northwestern China using the data of seis-mic tomography

Introduction The gravity anomaly is an indicator of the density distribution of the underground material. Therefore the gravity anomalies have been important data used for studying the deep crustal struc-ture for a long time. Many people have made detailed researches on the regional crustal structure inverted by Bouguer anomalies. In particular some empirical formulae and practical algorithms about the crustal thickness were brought forward, and a series of results were obtained (MENG, 1996)…     

Steric Sea-Level Change and its Impact on the Gravity Field caused by Global Climate Change

It is sometimes assumed that steric sea-level variations do not produce a gravity signal as no net mass change, thus no change of ocean bottom pressure is associated with it. Analyzing the output of two CO₂ emission scenarios over a period of 2000 years in terms of steric sea-level changes, we try to quantify the gravitational effect of steric sea-level variations. The first scenario, computed with version 2.6 of the Earth System Climate Model developed at the University of Victoria, Canada (UVic ESCM), is implemented with a linear CO₂ increase of 1% of the initial concentration of 365 ppm and shows a globally averaged steric effect of 5.2 m after 2000 years. In the second scenario, computed with UVic ESCM version 2.7, the CO₂ concentration increases quasi-exponentially to a level of 3011 ppm and is hold fixed afterwards. The corresponding globally averaged steric effect in the first 2000 years is 2.3 m. We show, due to the (vertical) redistribution of ocean water masses (expansion or contraction), the steric effect results also in a small change in the Earth’s gravity field compared to usually larger changes associated with net mass changes. Maximum effects for computation points located on the initial ocean surface can be found in scenario 1, with the effect on gravitational attraction and potential ranging from 0.0 to −0.7·10⁻⁵ m s⁻² and −3·10⁻³ to 6·10⁻³ m² s⁻², respectively. As expected, the effect is not zero but negligible for practical applications.     

Thermal structure and heat loss at the summit crater of an active lava dome

Forward-Looking Infrared (FLIR) nighttime thermal images were used to extract the thermal and morphological properties for the surface of a blocky-to-rubbley lava mass active within the summit crater of the Caliente vent at Santiaguito lava dome (Guatemala). Thermally the crater was characterized by three concentric regions: a hot outer annulus of loose fine material at 150–400°C, an inner cold annulus of blocky lava at 40–80°C, and a warm central core at 100–200°C comprising younger, hotter lava. Intermittent explosions resulted in thermal renewal of some surfaces, mostly across the outer annulus where loose, fine, fill material was ejected to expose hotter, underlying, material. Surface heat flux densities (radiative + free convection) were dominated by losses from the outer annulus (0.3–1.5 × 10⁴  s⁻¹m⁻²), followed by the hot central core (0.1–0.4 × 10⁴ J s⁻¹m⁻²) and cold annulus (0.04–0.1 × 10⁴ J s⁻¹m⁻²). Overall surface power output was also dominated by the outer annulus region (31–176 MJ s⁻¹), but the cold annulus contributed equal power (2.41–7.07 MJ s⁻¹) as the hot central core (2.68–6.92 MJ s⁻¹) due to its greater area. Cooled surfaces (i.e. the upper thermal boundary layer separating surface temperatures from underlying material at magmatic temperatures) across the central core and cold annulus had estimated thicknesses, based on simple conductive model, of 0.3–2.2 and 1.5–4.3 m. The stability of the thermal structure through time and between explosions indicates that it is linked to a deeper structural control likely comprising a central massive plug, feeding lava flow from the SW rim of the crater, surrounded by an arcuate, marginal fracture zone through which heat and mass can preferentially flow.     

Strain-geoelectric resistivity precursor and the virtual dislocation model to the Tangshan earthquake

The complete records of geoelectric resistivity before two earthquakes were analyzed, including 16 stations within 240 km around the Tangshan earthquake and 2 stations within 50–60 km from theM6.1 Datong earthquake. By eliminating various disturbances in the records and realizing the precursory anomalies to be reliable, the authors studied the distribution of the geoelectric precursor field, which proves to be physically related to the earthquake source stress field. Comparision of the sign distribution of coseismic resistivity changes with the solution of earthquake source mechanism indicates that, the coseismic resistivity changes are of opposite sign but similar spatial distribution with respect to the corresponding resistivity precursor changes. Therefore, from the resistivity observations we are of the opinion that the Tangshan earthquake is an elastic rebound process. A virtual dislocation model of geoelectric precursor for the Tangshan earthquake is proposed, in which the geoelectric precursors are supposed to be caused by the strain accumulation due to a virtual dislocation, which is opposite in sign to the actual slip taking place at the earthquake occurrence. Taking into account of the non-linear characteristics of the amplification factorK=(Δρ/ρ)/ɛ in a range of 10⁻⁷–10⁻⁵ strain changes, the theoretical distribution of geoelectric precursors for the Tangshan earthquake was calculated based on the theory of fracture mechanics and reasonably selected dislocation parameters. The results show that the semi-quantitative theoretical values are in good agreement with the observed, suggesting that the virtual dislocation model of the geoelectric precursor is appropriate to the Tangshan earthquake. Contribution No. 96A0023, Institute of Geophysics, SSB, China. This work was supported by the National Natural Science Foundation of China and the Science and Technology Department of SSB.     

Analysis of the tidal gravity observations in Three Gorges reservoir area

The comparative observations of the tidal gravity in Three Gorges reservoir area, recorded with the LCR-ET20 gravimeter from Institute of Geodesy and Geophysics, Chinese Academy of Sciences and the DZW-9 gravimeter from Institute of Seismology, China Earthquake Administration, are introduced. High-accuracy tidal gravity pa-rameters in this area are obtained, and the scale factor of the DZW-9 gravimeter is also estimated at a value of ?(756.06±0.05)×10?8 m?s?2V?1. The observed residuals of the corresponding instruments are investigated. Nu-merical results indicate that the long-term drift of the DZW-9 gravimeter appears a linear characteristic, and the observed accuracy is of the same order comparing with that of the LCR-ET20 gravimeter. The results given in the paper can provide with an effective reference model of the tidal gravity correction to the ground based and space geodesy.     

Low-field variation of magnetic susceptibility measured by the KLY-4S Kappabridge and KLF-4A magnetic susceptibility meter: Accuracy and interpretational programme

The KLY-4S Kappabridge and KLF-4A Magnetic Susceptibility Meter enable automated measurement of susceptibility variation with field in the ranges of 2–450 A/m and 5–300 A/m (in effective values), respectively. Unfortunately, the measurement accuracy decreases with decreasing field and it is not easy to decide whether the susceptibility variation at the lowest fields is natural phenomenon or results from measuring errors. To overcome this problem, the accuracies of both the above instruments were investigated experimentally using artificial specimens (mixture of pure magnetite and plaster of Paris) with variable susceptibilities ranging from 1 × 10⁻⁵ to 5 × 10⁻². The complete curve of the field variation of susceptibility of each specimen was measured 10 times and the relative error was calculated for each field. In the KLY-4S Kappabridge, in specimens with susceptibilities higher than 100 × 10⁻⁶, the relative errors are lower than 3% in all fields and lower than 1% in the fields stronger than 10 A/m. In the KLF-4A Magnetic Susceptibility Meter, in relatively strongly magnetic specimens with susceptibilities 5 × 10⁻⁴ to 5 × 10⁻², the relative error is less than 1.5% in the entire field range. While the former instrument is convenient for investigating almost all rock types, the latter instrument is convenient for measuring moderately and strongly magnetic specimens. To facilitate work with field variation of susceptibility curves, showing variable accuracies with field, the programme FieldVar was written. One of its options is plotting the measured data with corresponding field-variable error bars. In this way, a tool is offered for interpreting such susceptibility changes that are sound and reasonable from the point of view of measuring accuracy.     

Gravity change and its mechanism after the first water impoundment in Three Gorges Project

Introduction During and before the construction of the Three Gorges Project, the annual observation was continuously undertaken by the Three Gorges gravity observation network established around the head area of the reservoir (the section between Sandouping and Badong). This network has pro-vided substantial high precise data related to seismic activity. However, the annual gravity survey together monitoring of the existed observation stations can hardly effective monitor some potential sei…     

Application of <Superscript>137</Superscript>Cs fingerprinting technique to interpreting sediment production records from reservoir deposits in a small catchment of the Hilly Loess Plateau,China

According to variations of 137Cs and clay contents, 44 flood couplets were identified in a profile of reservoir deposit with a vertical length of 28.12 m in the Yuntaishan Gully. Couplet 27 at the middle of the profile had the highest average 137Cs content of 12.65 Bq kg-1, which indicated the 1963s' deposits, then 137Cs content decreased both downward and upward in the profile. The second top and bottom couplets had average 137Cs contents of 2.15 Bq kg-1 and 0.92 Bq kg-1, respectively. By integrated analysis of reservoir construction and management history, variations of 137Cs contents over the profile, sediment yields of flood couplets and rainfall data during the period of 1958-1970, individual storms related to the flood couplets were identified. 44 floods with a total sediment yield of 2.36×104 m3 occurred and flood events in a year varied between 1 and 10 times during the period of 1960-1970. 7-10 flood events occurred during the wet period of 1961-1964 with very wet autumn, while only 1-2 events during the dry period of 1965-1969. Average annual specific sediment yield was 1.29×104 t km-2 a-1 for the Yuntaishan Gully during the period of 1960-1970, which was slightly higher than 1.11 ×104 t km-2 a-1 for the Upper Yanhe River Basin above the Ganguyi Hydrological Station and slightly lower than 1.40 ×104 t km-2 a-1 for the nearby Zhifang Gully during the same period. Annual specific sediment yields for the Yuntaishan Gully were correlated to the wet season's rainfalls well.     

Interactions of river discharge and tidal modulation in a tropical estuary, NE Brazil

Observations of thermohaline properties and currents were undertaken in the Curimataú River estuary (6°18′S), Rio Grande do Norte state (RN), Brazil, during consecutive neap–spring tidal cycles in the austral autumn rainy season. Highly asymmetric neap tide along channel velocities (−0.4 to 0.9 m s⁻¹) and highly stratified conditions were generated by an increase of the buoyancy energy from the freshwater input (R _iE≈5.6). During the spring-tidal cycle the river discharge decreased and the longitudinal velocity components were higher, less asymmetrical (−0.8 to 1.1 m s⁻¹) and semidiurnal, associated with moderately stratified conditions (R _iE≈0.1) due to the increase of the kinetic tidal energy forcing mechanism. The overall salinity variation from surface to bottom during two tidal cycles was from 20.5 to 36.3 and 29 to 36.7 in the neap and spring tide experiments, respectively; in the last experiment, the tropical water (TW) mass intrusion was enhanced. The net salt transport reversed from down to up estuary during the neap and spring tide experiments, respectively, varied from 6.0 to –2.0 kg m⁻¹ s⁻¹, an indication of changes in the main forcing of the estuary dynamics. Evaluation of a classical steady analytical model, in comparison with nearly steady experimental vertical profiles of velocity, shows an agreement classifiable as reasonably fair.     

Temporal variations in gravity at Mt. Etna (Italy) associated with the 1989 and 1991 eruptions

 Results are presented from 11 microgravity surveys on Mt. Etna between 1987 and 1993, a period including the major 1989 and 1991–1993 flank eruptions and subordinate 1990 activity. Measurements were made with LaCoste and Romberg D-62 and D-157 gravity meters along a network around the volcano between 1000 and 1900 m a.s.l. and, since 1992, a N–S summit profile. Gravity changes of as much as 200 μGal were observed at scales from the size of the summit region to that of the volcano. None was associated with significant changes in ground elevation. The data show an increase in gravity for 2 years before the 1989 eruption. The increase is attributed to the accumulation of magma (0.25–1.7×10⁹m³) in an elongate zone, oriented NNW–SSE, between 2.5 and 6 km below sea level. Part of this magma was injected into the volcanic pile to supply the 1989 and 1990 eruptions. It also probably fed the start of the 1991–1993 eruption, since this event was not preceded by significant gravity changes. A large gravity increase (up to 140 μGal) detected across the volcano between June and September 1992 is consistent with the arrival in the accumulation zone of 0.32–2.2×10⁹m³ of new magma, thus favoring continued flank effusion until 1993. A large gravity decrease (200 μGal) in the summit region marked the closing stages of the 1991–1993 event and is associated with magma drainage from the upper levels of Etna's central feeding system. Received: 15 July 1995 / Accepted: 27 October 1997     

MELIMEX,an experimental heavy metal pollution study: Vertical mixing in a large limno-corral

In a limno-corral (diameter 12 m, depth to sediments 10 m), located in Baldeggersee (Switzerland), vertical mixing has been measured during more than one year and compared to the conditions in the open lake (maximum depth 65 m, surface area 5.3 km²). The temperature method by McEwen and Hutchinson yields K_z values between 5×10⁻² cm²s⁻¹ at the upper boundary of the thermocline and 2×10⁻³ cm²s⁻¹ at the bottom, a value near the molecular diffusion of heat at 4°C (1.36×10⁻³ cm²s⁻¹). K_z calculated from profiles of excess radon-222 generally agree with those from the temperature data. Compared to the open lake, the corral has a more shallow epilimnion. However, during calm meteorological conditions, vertical mixing in the upper 10 m is similar outside and inside the corral. Supported by the Swiss National Science Foundation within the framework for its National Research Program on ‘Lake Currents.’     

Rheology of the lower crust beneath the northern part of North China: Inferences from lower crustal xenoliths from Hannuoba basalts, Hebei Province, China

Lower crustal xenoliths brought up rapidly by basaltic magma onto the earth surface may provide di-rect information on the lower crust. The main purpose of this research is to gain an insight into the rheology of the lower crust through the detailed study of lower crustal xenoliths collected from the Hannuoba basalt, North China. The lower crustal xenoliths in this area consist mainly of two pyroxene granulite, garnet granulite, and light-colored granulite, with a few exception of felsic granulite. The equilibration temperature and pressure of these xenoliths are estimated by using geothermometers and geobarometers suitable for lower crustal xenoliths. The obtained results show that the equilibration temperature of these xenoliths is within the range of 785―900℃, and the equilibrium pressure is within the range of 0.8―1.2 GPa, corresponding to a depth range of 28―42 km. These results have been used to modify the previously constructed lower crust-upper mantle geotherm for the studied area. The dif-ferential stress during the deformation process of the lower crustal xenoliths is estimated by using recrystallized grain-size paleo-piezometer to be in the range of 14―20 MPa. Comparing the available steady state flow laws for lower crustal rocks, it is confirmed that the flow law proposed by Wilks et al. in 1990 is applicable to the lower crustal xenoliths studied in this paper. The strain rate of the lower crust estimated by using this flow law is within the range of 10-13―10-11 s-1, higher than the strain rate of the upper mantle estimated previously for the studied area (10-17―10-13 s-1); the equivalent viscosity is estimated to be within the range of 1017―1019Pa·s, lower than that of the upper mantle (1019―1021 Pa·s). The constructed rheological profiles of the lower crust indicate that the differential stress shows no significant linear relation with depth, while the strain rate increases with depth and equivalent vis-cosity decrease with depth. The results support the viewpoint of weak lower continental crust.     

Strombolian explosive styles and source conditions: insights from thermal (FLIR) video

Forward Looking Infrared Radiometer (FLIR) cameras offer a unique view of explosive volcanism by providing an image of calibrated temperatures. In this study, 344 eruptive events at Stromboli volcano, Italy, were imaged in 2001–2004 with a FLIR camera operating at up to 30 Hz. The FLIR was effective at revealing both ash plumes and coarse ballistic scoria, and a wide range of eruption styles was recorded. Eruptions at Stromboli can generally be classified into two groups: Type 1 eruptions, which are dominated by coarse ballistic particles, and Type 2 eruptions, which consist of an optically-thick, ash-rich plume, with (Type 2a) or without (Type 2b) large numbers of ballistic particles. Furthermore, Type 2a plumes exhibited gas thrust velocities (>15 m s⁻¹) while Type 2b plumes were limited to buoyant velocities (<15 m s⁻¹) above the crater rim. A given vent would normally maintain a particular gross eruption style (Type 1 vs. 2) for days to weeks, indicating stability of the uppermost conduit on these timescales. Velocities at the crater rim had a range of 3–101 m s⁻¹, with an overall mean value of 24 m s⁻¹. Mean crater rim velocities by eruption style were: Type 1 = 34 m s⁻¹, Type 2a = 31 m s⁻¹, Type 2b = 7 m s⁻¹. Eruption durations had a range of 6–41 s, with a mean of 15 s, similar among eruption styles. The ash in Type 2 eruptions originates from either backfilled material (crater wall slumping or ejecta rollback) or rheological changes in the uppermost magma column. Type 2a and 2b behaviors are shown to be a function of the overpressure of the bursting slug. In general, our imaging data support a broadening of the current paradigm for strombolian behavior, incorporating an uppermost conduit that can be more variable than is commonly considered.     

Observation and modeling of NPP for Pinus elliottii plantation in subtropical China

Based on the stem analysis of 59 individuals of Pinus elliottii in combination with tree biomass models, we calculated annual biomass increment of forest plots at Qianyanzhou Ecological Station, Chinese Academy of Sciences in subtropical China. In addition, canopy layer and community NPP were calculated based on 12 years’ litter fall data. NPP of the 21-year-old forest was estimated by using the BIOME BGC model; and both measured NPP and estimated NPP were compared with flux data. Community biomass was 10574 g · m⁻²; its distribution patterns in tree layer, shrub layer, herbaceous layer, tree root, herbaceous and shrub roots and fine roots were 7542, 480, 239, 1810, 230, 274 and 239 g · m⁻², respectively. From 1999 to 2004, the average annual growth rate and litter fall were 741 g · m⁻² · a⁻¹ (381.31 gC · m⁻² · a⁻¹) and 849 g · m⁻² · a⁻¹ (463 gC · m⁻² · a⁻¹), respectively. There was a significant correlation between annual litter fall and annual biomass increment; and the litter fall was 1.19 times the biomass increment of living trees. From 1985 to 2005, average NPP and GPP values based on BGC modeling were 630.88 (343.31–906.42 gC · m⁻² · a⁻¹) and 1 800 gC · m⁻² · a⁻¹ (1351.62–2318.26 gC · m⁻² · a⁻¹). Regression analysis showed a linear relationship (R ²=0.48) between the measured and simulated tree layer NPP values. NPP accounted for 30.2% (25.6%–32.9%) of GPP, while NEP accounted for 57.5% (48.1%–66.5%) of tree-layer NPP and 41.74% (37%–52%) of stand NPP. Soil respiration accounted for 77.0% of measured tree NPP and 55.9% of the measured stand NPP. NEE based on eddy covariance method was 12.97% higher than the observed NEP. Supported by the National Key Basic Research Special Foundation of China (Grant No. 2002CB4125), International Joint Research Project under Ministry of Science and Technology of China (Grant No. 2006DFB91920)     

Geochemistry and Sr-Nd isotope characteristics of tonalites in Zêtang, Tibet: New evidence for intra-Tethyan subduction

Tonalites from the island arc rock assemblage in the Zêtang segment of the Yarlung Zangbo suture zone were analyzed for major, trace elements (including REE) and Sr-Nd isotope. The experimental datademonstrate that the tonalites have the adakite-like characteristics, including high SiO2 (58%-63%),Al2O3 (18.4%-22.4%), Sr (810×10-6-940×10-6), Sr/Y (77-106), low HREE (Y=9×10-6-11×10-6, Yb=1×10-6-1.3×10-6), with LREE enrichment and faint Eu positive anomaly. Isr (0.70421-0.70487) is relatively low whereas 143Nd/144Nd (0.512896-0.512929) and εNd(t) values ( 6.7- 7.3) are high. These feainvolvement of a small amount of oceanic sediments. The identification of Z(e)tang adakites, derived from slab melting, presents new evidence for the intra-Tethyan subduction and the previous suggestion about the existence of intra-oceanic island arc within Tethys.     

Carbon budgets of three temperate forest ecosystems in Dongling Mt.,Beijing,China

There is a general agreement that forest ecosystems in the Northern Hemisphere function as signifi-cant sinks for atmospheric CO2; however, their magnitude and distribution remain large uncertainties. In this paper, we report the carbon (C) stock and its change of vegetation, forest floor detritus, and mineral soil, annual net biomass increment and litterfall production, and respiration of vegetation and soils between 1992 to 1994, for three temperate forest ecosystems, birch (Betula platyphylla) forest, oak (Quercus liaotungensis) forest and pine (Pinus tabulaeformis) plantation in Mt. Dongling, Beijing, China. We then evaluate the C budgets of these forest ecosystems. Our results indicated that total C density (organic C per hectare) of these forests ranged from 250 to 300 t C ha-1, of which 35―54 t C ha-1 from vegetation biomass C and 209―244 t C ha-1 from soil organic C (1 m depth, including forest floor detritus). Biomass C of all three forests showed a net increase, with 1.33―3.55 t C ha-1 a-1 during the study period. Litterfall production, vegetation autotrophic respiration, and soil heterotrophic respira-tion were estimated at 1.63―2.34, 2.19―6.93, and 1.81―3.49 t C ha-1 a-1, respectively. Ecosystem gross primary production fluctuated between 5.39 and 12.82 t C ha-1 a-1, about half of which (46%―59%, 3.20―5.89 t C ha-1 a-1) was converted to net primary production. Our results suggested that pine forest fixed C of 4.08 t ha-1 a-1, whereas secondary forests (birch and oak forest) were nearly in balance in CO2 exchange between the atmosphere and ecosystems.     

The dynamics of a double-cell hydrothermal system in triggering seismicity at Somma-Vesuvius: results from a high-resolution radon survey (revisited)

Data collected at Somma-Vesuvius during the 1998–1999 radon surveys have been revisited and reinterpreted in light of recent geophysical and geochemical information. The duration of selected radon anomalies, together with the decay properties of radon, have been used to estimate the permeability and porosity of rocks of the deep hydrothermal system. The current local cyclic seismicity is explained by means of a double convective-cell model. Convective cells are separated by a low-permeability horizon located at about 2–2.5 km below sea level. Fluids convecting within the upper cells show temperatures ranging 300–350°C. Rock permeabilities in this sector are estimated on the order of 10⁻¹² m², for porosities (ϕ) of about 10⁻⁵ typical of a brittle environment where fluid velocities may reach ∼800 m/day. Fluid temperatures within the lower cells may be as high as 400–450°C, consistent with supercritical regimes. The hydrodynamic parameters for these cells are lower, with permeability k ∼ 10⁻¹⁵ m², and porosity ranging from 10⁻⁶ to 10⁻⁷. Here, fluid motion toward the surface is controlled by the fracture network within a porous medium approaching brittle–ductile behaviour, and fluid velocities may reach ∼1,800 m/day. The low-permeability horizon is a layer where upper and lower convecting cells converge. In this region, fluids (convecting both at upper and lower levels) percolate through the wallrock and release their brines. Due to self-sealing processes, permeability within this horizon reaches critical values to keep the fluid pressure near lithostatic pressure (for k ∼ 10⁻¹⁸ m²). Deep fluid pressure buildups precede the onset of hydrothermally induced earthquakes. Permeability distribution and rock strength do not exclude that the next eruption at Somma-Vesuvius could be preceded by a seismic crisis, eventually leading to a precursory phreatic explosion. The coupling of these mechanisms has the potential of inducing pervasive failure within rocks of the hydrothermal shell, and may be a prelude to a magmatic eruption. It is finally emphasised that the integrated analysis of seismic and geochemical data, including radon emissions, could be successfully used in testing temperature distributions and variations of porosity and permeability in active geothermal reservoirs.     

Chemical and isotopic compositions of thermal springs,fumaroles and bubbling gases at Tacaná Volcano (Mexico–Guatemala): implications for volcanic surveillance

This study presents baseline data for future geochemical monitoring of the active Tacaná volcano–hydrothermal system (Mexico–Guatemala). Seven groups of thermal springs, related to a NW/SE-oriented fault scarp cutting the summit area (4,100m a.s.l.), discharge at the northwest foot of the volcano (1,500–2,000m a.s.l.); another one on the southern ends of Tacaná (La Calera). The near-neutral (pH from 5.8 to 6.9) thermal (T from 25.7°C to 63.0°C) HCO₃–SO₄ waters are thought to have formed by the absorption of a H₂S/SO₂–CO₂-enriched steam into a Cl-rich geothermal aquifer, afterwards mixed by Na/HCO₃-enriched meteoric waters originating from the higher elevations of the volcano as stated by the isotopic composition (δD and δ¹⁸O) of meteoric and spring waters. Boiling temperature fumaroles (89°C at ~3,600m a.s.l. NW of the summit), formed after the May 1986 phreatic explosion, emit isotopically light vapour (δD and δ¹⁸O as low as −128 and −19.9‰, respectively) resulting from steam separation from the summit aquifer. Fumarolic as well as bubbling gases at five springs are CO₂-dominated. The δ¹³C_CO2 for all gases show typical magmatic values of −3.6 ± 1.3‰ vs V-PDB. The large range in ³He/⁴He ratios for bubbling, dissolved and fumarolic gases [from 1.3 to 6.9 atmospheric ³He/⁴He ratio (R _A)] is ascribed to a different degree of near-surface boiling processes inside a heterogeneous aquifer at the contact between the volcanic edifice and the crystalline basement (⁴He source). Tacaná volcano offers a unique opportunity to give insight into shallow hydrothermal and deep magmatic processes affecting the CO₂/³He ratio of gases: bubbling springs with lower gas/water ratios show higher ³He/⁴He ratios and consequently lower CO₂/³He ratios (e.g. Zarco spring). Typical Central American CO₂/³He and ³He/⁴He ratios are found for the fumarolic Agua Caliente and Zarco gases (3.1 ± 1.6 × 10¹⁰ and 6.0 ± 0.9 R _A, respectively). The L/S (5.9 ± 0.5) and (L + S)/M ratios (9.2 ± 0.7) for the same gases are almost identical to the ones calculated for gases in El Salvador, suggesting an enhanced slab contribution as far as the northern extreme of the Central American Volcanic Arc, Tacaná.     

Echo intensity data as a directional antenna for observing processes above sloping ocean bottoms

Relative ‘echo intensity’ data (dI) from a bottom-mounted four-beam 300 kHz acoustic Doppler current profiler (ADCP) are used to infer propagation of vigorous processes above a continental slope. The 3- to 60-m horizontal beam spread and the 2-Hz sampling allow the distinction of different arrival times t _i , i = 1,..., 4, at different distances in the acoustic beams from sharp changes in dI-content associated with frontal non-linear and turbulent bores or ‘waves’. The changes in dI are partially due to variations in amounts of resuspended material carried by the near-bottom turbulence and partially due to the fast variations in density stratification (‘stratified turbulence’), as inferred from 1-Hz sampled thermistor string data above the ADCP. Such bores are observed to pass the mooring up to 80 m above the bottom, having typical propagation speeds c = 0.15–0.5 m s⁻¹, as determined from dI(t _i ). Particle speeds in the immediate environment of a bore amount to |u|_env=c ± 0.05 m s⁻¹, the equality being a necessary condition for kinematic instability, whilst the maximum particle speeds amount |u|_max = 1.2–2c. The dI-determined directions of up-, down- and alongslope processes are all to within ±10° of the ADCP’s beam-spread averaged current (particle velocity) data.