中国新疆北天山和台湾南部陆地泥火山研究进展

在系统介绍中国新疆北天山地区和台湾南部地区泥火山研究进展的基础上,对其地质特征进行了对比分析。结果显示,北天山和台湾南部地区的泥火山均沿着断裂带分布,主要位于背斜轴部,泥火山分布区地层多出露为含泥岩层。对两个地区泥火山喷出物物理特征进行了对比分析,固体喷发物的矿物成分相似,如石英、蒙脱石等;液体喷出物的泥浆温度与冒泡频率相近,但最大气泡直径与气体流量有很大差别。又分别对两地区液体、气体喷出物的化学特征进行了对比分析,液体喷出物均盐度高;甲烷是大多数泥火山喷发气体的主要成分,一些泥火山喷发的气体主要是二氧化碳。区域构造地质和气候条件不同,导致两地泥火山喷出物存在差异。从现有研究来看,两地泥火山的喷发都是岩层的孔隙压力增大造成的。两个地区泥火山与当地地震活动之间表现出良好的对应关系。泥火山的地球化学参数可能是地震活动的潜在指标。     

Automatic extraction of geometrical characteristics hidden in ground‐penetrating radar sectional images using simultaneous perturbation artificial bee colony algorithm

Ground‐penetrating radar is one of the most effective methods of detecting shallow buried objects. Ground‐penetrating radar radargram is a vertical map of the radar pulse reflections that returns from subsurface objects, and in the case of cylindrical objects, it would be a hyperbola. In order to get clear and accurate information on the presence, location, and geometry of the buried objects, the radargrams need to be interpreted. Interpretation of the results is a time‐consuming task and needs an expert with vast knowledge. Development of an automatic interpretation method of B‐scan ground‐penetrating radar images would be an effective and efficient solution to this problem. A novel automatic interpretation method of ground‐penetrating radar images, based on simultaneous perturbation artificial bee colony algorithm using tournament selection strategy, simultaneous perturbation stochastic approximation method, and new search equations, is introduced in this paper. The proposed algorithm is used to extract geometrical parameters, i.e. depth, location, and radius, of buried cylindrical objects in order to assess its accuracy. Synthetic data, simulated using GprMax2D forward modelling program, and real data, surveyed in the campus of Isfahan University of Technology, are used in the assessment. The performance of the proposed method in detecting synthetic hyperbolas is compared with that of the original artificial bee colony algorithm, genetic algorithm, and modified Hough transform. The results show superiority of the proposed algorithm, in detecting synthetic hyperbolas. Furthermore, the performance of the proposed method in estimating depth and radius of pipes in real ground‐penetrating radar images is compared with that of the modified Hough transform. The results indicate higher accuracy of the proposed method in estimating geometrical parameters of the buried cylindrical objects.     

Chemical evolution and volcanic activity of the active crater lake of Poás volcano, Costa Rica, 1993–1997

Concentrations of chloride and sulfate and pH in the hot crater lake (Laguna Caliente) at Poás volcano and in acid rain varied over the period 1993–1997. These parameters are related to changes in lake volume and temperature, and changes in summit seismicity and fumarole activity beneath the active crater. During this period, lake level increased from near zero to its highest level since 1953, lake temperature declined from a maximum value of 70°C to a minimum value of 25°C, and pH of the lake water increased from near zero to 1.8. In May 1993 when the lake was nearly dry, chloride and sulfate concentrations in the lake water reached 85,400 and 91,000 mg l⁻¹, respectively. Minimum concentrations of chloride and sulfate after the lake refilled to its maximum volume were 2630 and 4060 mg l⁻¹, respectively. Between January 1993 and May 1995, most fumarolic activity was focused through the bottom of the lake. After May 1995, fumarolic discharge through the bottom of the lake declined and reappeared outside the lake within the main crater area. The appearance of new fumaroles on the composite pyroclastic cone coincided with a dramatic decrease in type B seismicity after January 1996. Between May 1995 and December 1997, enhanced periods of type A seismicity and episodes of harmonic tremor were associated with an increase in the number of fumaroles and the intensity of degassing on the composite pyroclastic cone adjacent to the crater lake. Increases in summit seismic activity (type A, B and harmonic tremor) and in the height of eruption plumes through the lake bottom are associated with a period of enhanced volcanic activity during April–September 1994. At this time, visual observations and remote fumarole temperature measurements suggest an increase in the flux of heat and gases discharged through the bottom of the crater lake, possibly related to renewed magma ascent beneath the active crater. A similar period of enhanced seismic activity that occurred between August 1995 and January 1996, apparently caused fracturing of sealed fumarole conduits beneath the composite pyroclastic cone allowing the focus of fumarolic degassing to migrate from beneath the lake back to the 1953–1955 cone. Changes in the chemistry of summit acid rain are correlated changes in volcanic activity regardless of whether fumaroles are discharging into the lake or are discharging directly into the atmosphere.     

Origin, characteristics, and behaviour of lahars following the 1990 eruption of Kelud volcano, eastern Java (Indonesia)

 In contrast to most twentieth-century eruptions of Kelud volcano (eastern Java), the 10 February 1990 plinian eruption was not accompanied by lake-outburst lahars. However, at least 33 post-eruption lahars occurred between 15 February and 28 March 1990. They swept down 11 drainage systems and travelled as far as 24 km at an estimated mean peak velocity in the range of 4–11 m s^–1. The deposits (volume ≥30 000 000 m³) were approximately 7 m thick 2 km from vent, and 3 m thick 10 km from vent, on the volcaniclastic apron surrounding the volcano. Subtle but significant sedimentological differences in the deposits relate to four flow types: (a) Early, massive deposits are coarse, poorly sorted, slightly cohesive, and commonly inversely graded. They are inferred to record hot lahars that incorporated pumice and scoria from pyroclastic-flow deposits, probably by rapid remobilization of hot proximal pyroclastic flow deposits by rainfall runoff. Sedimentary features, such as clasts subparallel to bedding and thick, reversely to ungraded beds, suggest that these flows were laminar. (b) Abundant, very poorly sorted deposits include non-cohesive, clast-supported, inversely graded beds and ungraded, finer-grained, and cohesive matrix-supported beds. These beds display layering and vertical segregation/density stratification, suggesting unsteady properties of pulsing debris flows. They are interpreted as deposited from segments of flow waves at a middle distance downstream that incorporated pre-eruption sediments. Sedimentological evidence suggests unsteady flow properties during progressive aggradation. (c) Fine-grained, poorly sorted and ungraded deposits are interpreted as recording late hyperconcentrated streamflows that formed in the waning stage of an overflow and transformed downcurrent into streamflows. (d) Ungraded, crudely stratified deposits were emplaced by flows transitional between hyperconcentrated flows and streamflows that traveled farther downvalley (as far as 27 km from the vent). At Kelud, the transformation of flow and behavior occurs within only 10 km of the source, at the apex of the alluvial fans. The rapid change of flow behavior is attributed to the low fines content and to the unsteady flow regime, which may be due to: (a) the rapid deposition of bedload, owing to the break in channel gradient close to the vent and to changes in channel cross-section and roughness; and (b) the very low silt+clay content in the non-cohesive deposits. These deposits mix with water to produce streamflows. Received: 27 June 1997 / Accepted: 5 January 1998     

Architecture of the deep critical zone in the Río Icacos watershed (Luquillo Critical Zone Observatory,Puerto Rico) inferred from drilling and ground penetrating radar (GPR)

In the critical zone, surficial bedrock interactions result in the formation of a mantle of chemically‐ and physically‐altered material defined here as regolith. In the watershed of the Río Icacos, an upland river draining the Luquillo Mountains in tropical Puerto Rico, we explored the influence of lithology (quartz diorite versus hornfels‐facies volcaniclastic rock) on weathering. Regolith profiles were studied by drilling boreholes and imaging the subsurface using ground penetrating radar (GPR). Overall, the regolith structure is not laterally continuous but rather is punctuated by zones of deep fractures that host in situ weathering, corestones, and colluvial material. GPR images of these vertical zones show reflectors at 15–20 m depth. Thus, the architecture of the critical zone in the upper Luquillo Mountains is highly dependent on lithology and its influence on fracture development. At the highest elevations where hornfels overlies quartz diorite, positive feedbacks occur when the water table drops so that oxidative weathering of biotite in the more felsic rock creates microfractures and allows deeper infiltration of meteoric waters. Such exposure results in some of the fastest weathering rocks in the world and may contribute to formation of the knickpoint in the Río Icacos watershed. This work represents the first study combining GPR and drilling to look at the structure of the deep critical zone and demonstrates: (1) the importance of combining direct methods (such as drilling) with indirect methods (such as GPR) to understand the architecture of the critical zone in tropical systems; (2) the interplay of the surficial stress regime, lithology and climate in dictating the architecture of weathering. Copyright © 2016 John Wiley & Sons, Ltd.     

Hydraulic,physical and rheological characteristics of rain‐triggered lahars at Semeru volcano,Indonesia

Lahars (volcanic debris flows) have been responsible for 40% of all volcanic fatalities over the past century. Mount Semeru (East Java, Indonesia) is a persistently active composite volcano that threatens approximately one million people with its lahars and pyroclastic flows. Despite their regularity, the behaviour and the propagation of these rain‐triggered lahars are poorly understood. In situ samples were taken from lahars in motion at two sites in the Curah Lengkong River, on the southeast flank of Semeru, providing estimates of the particle concentration, grain size spectrum, grain density and composition. This enables us to identify flow sediment from three categories of lahars: (a) hyperconcentrated flow, (b) non‐cohesive, clast‐ and matrix‐supported debris flow, and (c) muddy flood. To understand hyperconcentrated flow sediment transport processes, it is more appropriate to sample the active flows than the post‐event lahar deposits because in situ sampling retains the full spectrum of the grain‐size distribution. Rheometrical tests on materials sampled from moving hyperconcentrated flows were carried out using a laboratory vane rheometer. Despite technical difficulties, results obtained on the <63, <180, and <400 µm fractions of the sampled sediment, suggest a purely frictional behaviour. Importantly, and contrary to previous experiments conducted with monodisperse suspensions, our results do not show any transition towards a viscous behaviour for high shear rates. These data provide important constraints for future physical and numerical modelling of lahar flows. Copyright © 2010 John Wiley & Sons, Ltd.     

Major and trace element geochemistry of neutral and acidic thermal springs at El Chichón volcano,Mexico: Implications for monitoring of the volcanic activity

Four groups of thermal springs with temperatures from 50 to 80 °C are located on the S–SW–W slopes of El Chichón volcano, a composite dome-tephra edifice, which exploded in 1982 with a 1 km wide, 160 m deep crater left. Very dynamic thermal activity inside the crater (variations in chemistry and migration of pools and fumaroles, drastic changes in the crater lake volume and chemistry) contrasts with the stable behavior of the flank hot springs during the time of observations (1974–2005). All known groups of hot springs are located on the contact of the basement and volcanic edifice, and only on the W–SW–S slopes of the volcano at almost same elevations 600–650 m asl and less than 3 km of direct distance from the crater. Three groups of near-neutral (pH ≈ 6) springs at SW–S slopes have the total thermal water outflow rate higher than 300 l/s and are similar in composition. The fourth and farthest group on the western slope discharges acidic (pH ≈ 2) saline (10 g/kg of Cl) water with a much lower outflow rate (< 10 l/s).     

Preliminary sensitivity study of eruption source parameters for operational volcanic ash cloud transport and dispersion models — A case study of the August 1992 eruption of the Crater Peak vent, Mount Spurr, Alaska

Ash clouds are one of the major hazards that result from volcanic eruptions. Once an eruption is reported, volcanic ash transport and dispersion (VATD) models are used to forecast the location of the ash cloud. These models require source parameters to describe the ash column for initialization. These parameters include: eruption cloud height and vertical distribution, particle size distribution, and start and end time of the eruption. Further, if downwind concentrations are needed, the eruption mass rate and/or volume of ash need to be known. Upon notification of an eruption, few constraints are typically available on many of these source parameters. Recently, scientists have defined classes of eruption types, each with a set of pre-defined eruption source parameters (ESP). We analyze the August 18, 1992 eruption of the Crater Peak vent at Mount Spurr, Alaska, which is the example case for the Medium Silicic eruption type. We have evaluated the sensitivity of two of the ESP – the grain size distribution (GSD) and the vertical distribution of ash – on the modeled ash cloud. HYSPLIT and Puff VATD models are used to simulate the ash clouds from the different sets of source parameters. We use satellite data, processed through the reverse absorption method, as reference for computing statistics that describe the modeled-to-observed comparison. With the grain size distribution, the three options chosen, (1) an estimated distribution based on past eruption studies, (2) a distribution with finer particles and (3) the National Oceanic and Atmospheric Administration HYSPLIT GSD, have little effect on the modeled ash cloud. For the initial vertical distribution, both linear (uniform concentration throughout the vertical column) and umbrella shapes were chosen. For HYSPLIT, the defined umbrella distribution (no ash below the umbrella), apparently underestimates the lower altitude portions of the ash cloud and as a result has a worse agreement with the satellite detected ash cloud compared to that with the linear vertical distribution for this particular eruption. The Puff model, with a Poisson function to represent the umbrella cloud, gave similar results as for a linear distribution, both having reasonable agreement with the satellite detected cloud. Further sensitivity studies of this eruption, as well as studies using the other source parameters, are needed.