Reconnaissance geophysics and geology of the Pinacate craters,Sonora, Mexico

The Pinacate volcanic field of northwest Sonora, Mexico contains eleven tuff-rimmed explosion craters as well as hundreds of cinder cones. A multiphase eruptive and collapse history for each of the craters is evidenced by: cinder cones exposed in cross-section in crater rims, rim-breaking cinder cones, major eccentric vents, erosion hiatuses between tuff deposition, and coalescing craters exhibiting differeng degrees of erosion. Basaltic intrusions occurred along ring fractures and in the center of at least one crater. Magnetic surveys indicate that explosion craters have tuffaceous feeder vents whereas cinder cones have basaltic, highly magnetic plugs. Computer models illustrate that magnetic terrain effects can contribute to anomalies observed associated with craters, and need to be considered where magnetic surveys include irregularities such as rims or peaks.     

Two‐dimensional nonlinear diffusive numerical simulation of geomorphic modifications to cinder cones

The temporal evolution of simple landforms such as cinder cones by nonlinear diffusive processes is studied through the use of a new 2D numerical model using well‐established and accurate numerical mathematics and high‐resolution digital elevation models (DEMs). Extending 1D (profile) nonlinear diffusion analyses used in cinder cone, hillslope and fault scarp evolution studies, we have implemented a 2D numerical model with a spatially and temporally varying sediment transport rate coefficient scaled nonlinearly by the ratio of local slope to critical slope. The high accuracy and efficient numerical implementation are documented in the paper and the MATLAB toolkit developed is used to solve for the developmentof an initial 2D cone form. First, we examine the nonlinear transport rule and suggest a refinement that accounts explicitly for flux at threshold slopes. We find that the maximum diffusion (necessarily introduced in the numerical model to avoid infinite rates) at the critical slope controls the final morphology, especially approaching steady state. Secondly, solving the landscape evolution problem in 2D enables a natural accounting for sediment flux convergence or divergence in the profile. Thirdly, the boundary behavior of a given landscape element controls much of what happens in that domain and so we allow for arbitrary flux magnitude or elevation boundary conditions. Fourthly, landscapes are heterogeneous in their surface cover and so we allow for spatially and temporally varying transport rate k and we permit an arbitrary vertical displacement field within the model domain. To test the new formulation for the nonlinear term, the effect of variable diffusivity k and the numerical schemes implemented, we apply the model to cinder cones built on the flanks of Mount Etna in 2001 and 2002–2003. We explore the effects of DEM resolution with data from the 2001 cone and the utility of spatially variable diffusivity to explain the variation in erosion measured by differencing repeat light detection and ranging (LIDAR) surveys gathered in 2004 and 2007 over the 2002–2003 cone complex. Copyright © 2013 John Wiley & Sons, Ltd.     

The influence of climate and microclimate (aspect) on soil creep efficiency: Cinder cone morphology and evolution along the eastern Mediterranean Golan Heights

Although hillslope evolution has been subject to much investigation for more than a century, the effect of climate on the morphology of soil‐mantled hillslopes remains poorly constrained. In this study, we perform numerical simulations of volcanic cinder cones in the Golan Heights (eastern Mediterranean) to estimate soil transport efficiency across a significant north–south gradient in mean annual precipitation (1100 to 500 mm). We use the initial cinder cone morphology (constrained by stratigraphy), the modern hillslope form (surveyed with sub‐meter accuracy) and the eruption age (based on ⁴⁰Ar–³⁹Ar chronology) to predict the best‐fit value of the soil transport coefficient (‘diffusivity’) based on a nonlinear transport model. Our results indicate that the best‐fit diffusivity (K ) varies from 1 to 6 m² ka^?1 among the five cinder cones in our field area. Diffusivity (K ) values vary systematically with precipitation and hillslope aspect; specifically, K is higher on south‐facing (drier) hillslopes and decreases with mean annual precipitation. We interpret this climate dependency to reflect vegetation‐driven variations in apparent soil cohesion, which increases with root network density, and attenuation of rain splash and overland flow erosion, which increases with vegetative ground cover. To assess how vegetative root mass and ground cover vary with precipitation and aspect, we quantified the spatial distribution of NDVI (normalized difference vegetation index) from ASTER satellite images and observed spatial variations that correlate with our calibrated values of K . Analysis of previously studied cinder cones in the USA can be used to extend our framework to arid domains. This endeavor suggests a humped relationship between K and precipitation with maximum diffusivity at mean annual precipitation of 400–600 mm. Copyright © 2017 John Wiley & Sons, Ltd.     

A transition from strombolian to phreatomagmatic activity induced by a lava flow damming water in a valley

The Golan Heights is a Plio-Pleistocene volcanic plateau. Cinder cones of Late Pleistocene age are very common in the eastern and northern Golan, while phreatomagmatic deposits are relatively rare and occur just in two structures — the maar of Birket Ram and the tuff ring of Mt. Avital. The complex of Mt. Avital includes two large cinder cones, a tuff ring with an elongated central depression and several basaltic flows, some of them breach the cinder cones. The (exposed) eruptive history of the complex includes (1) an early stage of basaltic lava flows, (2) strombolian activity and the buildup of the southern cinder cone, (3) a second stage of basaltic flows and the buildup of the northern cinder cone, and then a transition to (4) phreatomagmatic explosions. The phreatomagmatic deposits include surges, lapilli fallout deposits and coarse-grained lithic tuff breccias, which were found up to 200 m above the central depression. Basaltic and scoriaceous clasts are the main component of all deposits, while juvenile material is usually a minor component, almost absent in the lapilli deposits.It is suggested that the phreatomagmatic events in Mt. Avital were induced by the infiltration of water from a lake that existed in a nearby topographic low (Quneitra Valley). The lake was formed or significantly expanded at about 300 ka due to a lava flow that blocked the drainage of the valley to the west. The interlayering of tuff and scoria at the top of the northern cinder cone and the good preservation of a lava flow top breccia under the surges imply that the phreatomagmatic activity immediately followed and even coincided with the last stages of strombolian activity. It is suggested that the dry–wet transition was triggered by the effusion of the second stage lavas and the buildup of the northern cinder cone, which probably caused a reduction of pressure in the magmatic system and allowed the lake water an access to the magmatic system. The minimum age of the phreatomagmatic events is determined by a 54 ka Musterian site which lies directly on top of the tuff in the Quneitra Valley.     

Spatial soil loss risk and reservoir siltation in semi-arid Tunisia

Abstract

Soil erosion vulnerability and extreme rainfall characteristics over the Mediterranean semi-arid region of Tunisia are crucial input for estimation of siltation rate in artificial reservoirs. A comprehensive high-resolution database on erosive rainfall, together with siltation records for 28 small reservoirs, were analysed for this region, the Tunisian Dorsal (the easternmost part of the Atlas Mountains). The general life-span of these reservoirs is only about 14 years. Depending on the soil degradation in the different catchments, the corresponding reservoirs display a wide range of soil erosion rates. The average soil loss was 14.5 t ha^?1 year^?1 but some catchments display values of up to 36.4 t ha^?1 year^?1. The maximum 15-min duration rainfall intensity was used to determine the spatial distribution of rainfall erosivity. The northwestern parts of the Tunisian Dorsal display the most extreme rainfall erosivity. Spatial erosion patterns are to some extent similar; however, they vary greatly according to their location in the “soil degradation cycle”. This cycle determines the soil particle delivery potential of the catchment. In general, the northwestern parts of the Dorsal display modest soil erosion patterns due to the already severely degraded soil structure. Here, the soil surface is often the original bedrock. However, the greatest soil erosion occurs in the mid-eastern parts of the Dorsal, which represents the “degradation front”. The latter corresponds to the area with highest erosion, which is continuously progressing westward in the Dorsal. The large variation between the erosive rainfall events and the annual soil loss rates was explained by two important factors. The first relates to the soil degradation cycle. The second factor corresponds to the degradation front with the highest soil loss rates. At present this front is located at 300 m altitude and appears to be moving along an 80-km westward path starting from the east coast. A better understanding of the above can be used to better manage soils and soil covers in the Tunisian Dorsal area and, eventually, to decrease the soil erosion and reservoir siltation risk.

Citation Jebari, S., Berndtsson, R., Bahri, A. & Boufaroua, M. (2010) Spatial soil loss risk and reservoir siltation in semi-arid Tunisia. Hydrol. Sci. J. 55(1), 121–137.     

Subaqueous,rhyolitic dome-top tuff cones: a model based on the Devonian Bunga Beds,southeastern Australia and a modern analogue

The relics of two Late Devonian subaqueous rhyolitic dome-top tuff and pumice cone successions are preserved in the Bunga Beds outlier of the Boyd Volcanic Complex, southeastern Australia. These cone successions and other rhyolitic volcanics of the Bunga Beds are associated with turbidite and other deep-water massflow sedimentary rocks. The two cone successions have a generally similar stratigraphy. At the base, flow-banded, variably autobrecciated and quench-fragmented rhyolite, representing an intrusive to extrusive dome, is overlain by rhyolitesediment breccia, representing extrusion of the dome through the deep-water sediment pile and resedimentation down its flanks. In the northern cone succession an overlying, succession of bedded pumiceous crystal-rich to crystal-poor tuffs represents the onset of pyroclastic activity and growth of a tuff cone. An overyling debris flow deposit represents degradation of part of the cone. The topmost unit, a stratified pumice succession, is thought to represent another cone-building eruptive phase, and is separated from the underlying strata by a major slide surface. The southern cone succession contains less tuff and abundant pumice, and is also terminated by a debris-flow deposit, indicating cone degradation. A modern analogue for the inferred eruptive style and sequence is the 1953–1957 rhyolite eruption that formed the Tuluman Island lava-tuff cone complex in the Bismarck Sea. The eruptions were often cyclical consisting of an initial inferred submarine-lava-forming stage, passing into a pumicecone-forming stage, in some cases a subaeriallava-forming stage, and a final stage, following the cessation of volcanism, during which the cones collapsed gravitationally or were destroyed by wave erosion. Using observations from both the Tuluman Island eruptions and the preserved stratigraphies of the Devonian tuff cones, a dynamic model is proposed for the formation of subaqueous rhyolitic dome-top tuff and pumice cones.     

Relations between rainfall–runoff‐induced erosion and aeolian deposition at archaeological sites in a semi‐arid dam‐controlled river corridor

Process dynamics in fluvial‐based dryland environments are highly complex with fluvial, aeolian, and alluvial processes all contributing to landscape change. When anthropogenic activities such as dam‐building affect fluvial processes, the complexity in local response can be further increased by flood‐ and sediment‐limiting flows. Understanding these complexities is key to predicting landscape behavior in drylands and has important scientific and management implications, including for studies related to paleoclimatology, landscape ecology evolution, and archaeological site context and preservation. Here we use multi‐temporal LiDAR surveys, local weather data, and geomorphological observations to identify trends in site change throughout the 446‐km‐long semi‐arid Colorado River corridor in Grand Canyon, Arizona, USA, where archaeological site degradation related to the effects of upstream dam operation is a concern. Using several site case studies, we show the range of landscape responses that might be expected from concomitant occurrence of dam‐controlled fluvial sand bar deposition, aeolian sand transport, and rainfall‐induced erosion. Empirical rainfall‐erosion threshold analyses coupled with a numerical rainfall–runoff–soil erosion model indicate that infiltration‐excess overland flow and gullying govern large‐scale (centimeter‐ to decimeter‐scale) landscape changes, but that aeolian deposition can in some cases mitigate gully erosion. Whereas threshold analyses identify the normalized rainfall intensity (defined as the ratio of rainfall intensity to hydraulic conductivity) as the primary factor governing hydrologic‐driven erosion, assessment of false positives and false negatives in the dataset highlight topographic slope as the next most important parameter governing site response. Analysis of 4+ years of high resolution (four‐minute) weather data and 75+ years of low resolution (daily) climate records indicates that dryland erosion is dependent on short‐term, storm‐driven rainfall intensity rather than cumulative rainfall, and that erosion can occur outside of wet seasons and even wet years. These results can apply to other similar semi‐arid landscapes where process complexity may not be fully understood. Published 2015. This article is a U.S. Government work and is in the public domain in the USA     

大兴安岭焰山、高山火山——一种新的火山喷发型式

在野外地质资料基础上,利用火山形态学方法,探讨了大兴安岭焰山、高山火山的喷发型式。结果表明,大兴安岭哈拉哈河-绰尔河火山群中的焰山和高山火山不同于斯通博利式喷发形成的火山,其早期爆破喷发的火山碎屑形成火山渣锥、空降火山碎屑席和小型火山碎屑流,晚期溢出大量熔岩。两火山具有较高大的锥体(标高200～300m以上),在结构上,松散火山砾、火山弹等构成下部的降落锥,熔结集块岩构成上部的溅落锥。由火山砾和火山灰组成的空降火山碎屑席分布在火山锥体周围。两火山溢出的熔岩经历了从结壳熔岩→翻花石→渣状熔岩的演变。根据喷发产物可推断焰山和高山火山具有以下喷发特征:爆破喷发形成持续的喷发柱→斯通博利式喷发→熔岩喷泉喷溢,其中以持续时间较长的喷发柱区别于典型的斯通博利式喷发。类似焰山、高山火山的喷发特征,在龙岗第四纪火山群、镜泊湖全新世火山群中也都有个例,这是中国大陆火山作用中一种新的喷发型式。     

Relations of cinder cones to the magmatic evolution of Mount Mazama,Crater Lake National Park,Oregon

Cinder cones at Crater Lake are composed of high-alumina basaltic to andesitic scoria and lavas. The Williams Crater Complex, a basaltic cinder cone with andesitic to dacitic lava flows, stands on the western edge of the caldera, against an andesite flow from Mount Mazama. Bombs erupted from Williams Crater contain cores of banded andesite and dacite, similar to those erupted during the climatic eruption of Mount Mazama.Major- and trace-element variations exhibit an increase in incompatible elements and a decrease in compatible elements, consistent with crystal fractionation of olivine, plagioclase, clinopyroxene, orthopyroxene, and magnetite. LREE patterns in the rocks are irregular; each successive basalt is enriched in LREE relative to the preceding andesite.Compositional variations in the magmas of the cinder cones suggest that three magmatic processes were involved, partial melting, fractional crystallization, and magma mixing. Partial melting of more than one source produced primary basaltic magma(s). Subsequent mixing and fractional crystallization produced the more differentiated basaltic to andesitic magmas.     

Computer-simulation models of scoria cone degradation

Long-term erosional modifications of the relatively simple morphology of scoria (‘cinder') cones are ideally suited for study by field and computer-simulation methods. A series of temporally-distinct cones in the San Francisco and Springerville volcanic fields of Arizona provides the foundation for documenting the degradational evolution of scoria cones in a semi-arid climate. Progressive changes due to erosion are illustrated by the systematic decrease with increasing age of various morphometric parameters, including scoria cone height, cone height/width ratio (H_co/W_co), crater depth/width ratio, and slope angle. For example, Holocene–latest Pleistocene cones in the San Francisco field have a mean H_co/W_co value of 0.178±0.041, a mean maximum slope angle of 29.7±4.2°, and a mean average slope angle of 26.4±7.3°, whereas the group of Pliocene cones have values of 0.077±0.024, 20.5±5.8°, and 8.7±2.7°, respectively. Comparative morphology of scoria cones is a potentially useful dating tool for mapping volcanic fields.In order to better understand the degradational modifications of these volcanic landforms, we have developed a numerical approach to simulate the surficial processes responsible for the erosion of a typical scoria cone. The simulation algorithm can apply either a linear diffusion-equation model or a model with a nonlinear transport law. Using a finite-difference formulation, the simulation operates upon a three-dimensional scoria cone input as a matrix of elevation values. Utilizing both field and model results, the correlation between changing H_co/W_co value, cone age, and computer time step was expressed graphically to derive comprehensive values of the transport or diffusion coefficient (D_f) for both volcanic fields. For the San Francisco volcanic field, D_f had a calculated value of 21.4 m²/kyr for the linear model and 5.3 m/kyr for the nonlinear model, while for the Springerville volcanic field D_f had a calculated value of 24.4 m²/kyr for the linear model and 6.3 m/kyr for the nonlinear model.     

Physical modelling of rainfall‐ and snowmelt‐induced erosion of stony slope underlain by permafrost

Physical modelling experiments have been carried out in a cold room to test on a small scale, the effects of water supply during the thaw of an experimental slope with permafrost. Permafrost was maintained at depth and a thin active layer was frozen and thawed from the surface. Data from the experiments relate to two different conditions, first with moderate rainfall, and second with heavy rainfall during the thaw period. When moderate rainfall is applied during thaw phases, the experimental slope is slightly degraded. At the scale of the experiment, erosion processes involve frost jacking of the coarse blocks, frost creep and gelifluction that induce slow and gradual down slope displacements of the active layer, but also small landslides leading to large but slow mass movements with short displacements. Changes in experimental slope morphology are marked by the initiation of a small‐scale drainage network and the development of a little crest line which shows a progressive upslope migration. With such boundary conditions, there is not enough water supply to evacuate downslope the whole of the eroded material and a topographic smoothing is observed. When heavy rainfall is applied during thaw periods, rapid mass wasting (small mud‐flows and debris flows) become prominent. Slope failures are largely controlled by the water saturation of the active layer and by the occurrence of steeper slopes. At the scale of the experiment, rates of erosion and maximum incision increase by about 100% leading to significant slope degradation with marked and specific scars comparable to gullying. These morphological changes are dependant on both the size and the frequency of catastrophic events. These experiments provide detailed data that could improve the knowledge of the physical parameters that control the initiation, at a small‐scale, of erosion processes on periglacial slopes with a thin active layer and/or with thin cover of mobilizable slope deposits. Copyright © 2010 John Wiley & Sons, Ltd.     

Hydrology of stockpiled industrial coal exposed to rainfall

The movement of water and particulate matter through a simulated coal pile using a rainfall–coal‐pile simulator is investigated. In conjunction with coal‐pile management practices, this study indicates that coal‐pile systems change their internal and surficial structure with successive rainstorm events. This change in coal‐pile structure alters the proportions, rates and suspended solid concentrations within surface runoff and infiltration during subsequent storm events. Results demonstrate the need to determine the optimal coal‐pile slope that minimizes water contact time, but at the same time minimizes the potential for mass wasting and rill erosion. Promotion of a compacted surface is beneficial for the minimization of solids removal, and a stormwater detention pond with potential flocculant additive is suggested as a viable remedial measure for the containment and treatment of coal‐pile runoff. Copyright © 2002 John Wiley & Sons, Ltd.     

Monogenetic vulcanism in sierra Chichinautzin,Mexico

The variation in the activity patterns of the Chichinautzin volcanic rocks is discussed. This sequence of lavas and pyroclastic deposits is located in the central part of the Mexican Volcanic Belt, directly south of Mexico City, and is typical of its Quaternary monogenetic vulcanism. One-hundred and fourty-six volcanoes and their deposits covering 952 km² were mapped. Cone density is 0.15 km² with heights ranging from to 315 m and crater diameters from 50 to 750 m. Ratios of cone height/diameter decreased from 0.20 to 0.12 with age. Basal diameters varied from 0.1 km to 2 km. Lavas are mainly blocky andesites but some dacites and basalts were found. Lengths of flows range from 1.0 to 21.5 km with heights of 0.5 to 300 m and aspect rations of 21.4 to 350. Three types of volcanic structures are found in the area: scoria cones, lavas cones and thick flows lacking a cone. Pyroclastic deposits are basically Strombolian although some deposits were produced by more violent activity and lava cones seem to have formed by activity transitional to Hawaiian-type vulcanism. Therre is a dominant E-W trend shown mainly by the orientation of cone clusters. The Chichinautzin volcanic centers are compared to the monogenetic volcanoes of the Toluca and Paricutin areas which are similar.     

Relative scales of time and effectiveness of climate in watershed geomorphology

Peak rainfalls and peak runoff rates per unit area are comparable over a worldwide spectrum of climates. However, while the magnitude of the external contribution of energy or force in diverse regions is similar, the impact on the landscape varies markedly between regions. Absolute magnitudes of climatic events and absolute time intervals between such events do not provide satisfactory measures of the geomorphic effectiveness of events of different magnitudes and recurrence intervals. Although geomorphic processes are driven by complex sets of interrelated climatic, topographic, lithologic, and biologic factors, the work done by individual extreme events can be scaled as a ratio to mean annual erosion and the effectiveness of such events in forming landscape features can be related to the rate of recovery of channel form or mass wasting scars following alteration by the extreme event. Thus, a time scale for effectiveness may relate the recurrence interval of an event to the time required for a landform to recover the form existing prior to the event. River channels in temperate regions widened by floods of recurrence intervals from 50 to more than 200 years may regain their original width in matters of months or years. In semi-arid regions, recovery of channel form depends not only upon flows but upon climatic determinants of the growth of bottomland vegetation resulting in variable rates of recovery, on the order of decades, depending upon coincidence of average flows and strengthened vegetation. In truly arid regions the absence of vegetation and flow precludes recovery and the width of channels increases in drainage areas up to 100 km² but remains relatively constant at larger drainage areas. Area as well as time controls the effectiveness of specific events inasmuch as the likelihood of simultaneous peak discharges or rainfalls and large areas is less, particularly in arid regions where events spanning areas of more than several thousand km² are extremely rare if experienced at all. To some extent a decrease in area in a humid region is comparable with a regional change from humid toward more arid climate reflected in the increase in importance of episodic as contrasted with more continuous processes. Exceedingly rare floods of extreme magnitudes, estimated recurrence intervals of 500 years or longer, may exceed thresholds of competence otherwise unattainable in the ‘normal’ record resulting in ‘irreparable’ transformations of valley landforms. Denudation of hillslopes by mass wasting during relatively rare events can also be related to mean rates of denudation and to recovery of hillslope surfaces after scarring by different kinds of landslides. Measured recovery times described in the literature vary from less than a decade for some tropical regions to decades or more in temperate regions. Recurrence intervals of high magnitude storms which trigger mass wasting range from 1 to 2 years in some tropical areas, to 3 or 4 per hundred years in some areas of seasonal rainfall and to 100 or more years in some temperate regions. The effectiveness of climatic events on both hillslopes and rivers is not separable from gradient, lithology or other variables which control both thresholds of activity and recovery rate.     

Calibration of an oscillating nozzle‐type rainfall simulator

Nozzle‐type rainfall simulators are commonly used in hydrologic and soil erosion research. Simulated rainfall intensity, originating from the nozzle, increases as the distance between the point of measurement and the source is decreased. Hence, rainfall measured using rain gauges would systematically overestimate the rainfall received at the ground level. A simple model was developed to adjust rainfall measured anywhere under the simulator to plot‐wide average rainfall at the ground level. Nozzle height, plot width, gauge diameter and height, and gauge location are required to compute this adjustment factor. Results from 15 runs at different rain intensities and durations, and with different rain gauge layouts, showed that a simple average of measured rain would overestimate the plot‐wide rain by about 20 per cent. Using the adjustment factor to convert measured rainfall for individual gauges before averaging improved the estimate of plot‐wide rainfall considerably. For the 15 runs considered, overall discrepancy between actual and measured rain is reduced to less than 1 per cent with a standard error of 0·97 mm. This model can be easily tested in the ?eld by comparing rainfall depths of different sized gauges. With the adjustment factor they should all give very similar values. Copyright © 2003 John Wiley & Sons, Ltd.     

Computer program for stability analysis of steep,cohesive riverbanks

The ability to predict the stability of eroding riverbanks is a prerequisite for modelling alluvial channel width adjustment and a requirement for predicting bank erosion rates and sediment yield associated with bank erosion. Mass‐wasting of bank materials under gravity occurs through a variety of specific mechanisms, with a separate analysis required for each type of failure. This paper presents a computer program for the analysis of the stability of steep, cohesive riverbanks with respect to planar‐type failures. Planar‐type failures are common along stream channels destabilized by severe bed degradation. Existing stability analyses for planar‐type failures have a number of limitations that affect their physical basis and predictive ability. The computer program presented here is based on an analysis developed by Darby and Thorne. The software takes account of the geotechnical characteristics of the bank materials, the shape of the bank profile, and the relative elevations of the groundwater and surface water to estimate stability with respect to mass failure along a planar‐type failure surface. Results can be displayed either in terms of a factor of safety (ratio of resisting to driving forces), or probability of failure. The computer analysis is able to determine the relative amounts of bed degradation and bank‐toe erosion required to destabilize an initially stable bank. Data for the analysis are supplied in the form of either HEC‐2 hydrographic survey data files or user‐supplied bank profile data, in conjunction with user‐supplied geotechnical parameter values. Some examples, using data from the Upper Missouri River in Montana, are used to demonstrate potential applications of the software. Copyright © 2000 John Wiley & Sons, Ltd.     

Geology and K-Ar dating of the Tuxtla Volcanic Field,Veracruz, Mexico

The Tuxtla Volcanic Field (TVF) is located on the coast of the Gulf of Mexico in the southern part of the state of Veracruz, Mexico. Volcanism began about 7 my ago, in the Late Miocene, and continued to recent times with historical eruptions in ad 1664 and 1793. The oldest rocks occur as highly eroded remnants of lava flows in the area surrounding the historically active cone of San Martín Tuxtla. Between about 3 and 1 my ago, four large composite volcanoes were built in the eastern part of the area. Rocks from these structures are hydrothermally altered and covered with lateritic soils, and their northern slopes show extensive erosional dissection that has widened preexisting craters to form erosional calderas. The eastern volcanoes are composed of alkali basalts, hawaiites, mugearites, and benmoreites, with less common calc-alkaline basaltic andesites and andesites. In the western part of the area, San Martín Tuxtla Volcano and its over 250 satellite cinder cones and maars produced about 120 km³ of lava over the last 0.8 my. A ridge of flank cinder cones blocked drainage to the north to form Laguna Catemaco. Lavas erupted from San Martín and its flank vents are restricted to compositions between basanite and alkali basalt. The alignment of major volcanoes and flank vents along a N55°W trend suggests an extensional stress field in the crust with a minimum compressional stress orientation of N35° E. In total, about 800 km³ of lava has been erupted in the TVF in the last 7 my. This gives a magma output rate of about 0.1 km³/1000 year, a value smaller than most composite cones, but similar to cinder cone fields that occur in central Mexico. Individual eruptions over the last 5000 years had volumes on the order of 0.1km³, with average recurrence intervals of 600 years. The alkaline compositions of the TVF lavas contrast markedly with the calc-alkaline compositions erupted in the subduction-related Mexican Volcanic Belt to the west, leading previous workers to suggest that the TVF is not related to subduction. Trace-element signatures of TVF lavas indicate, however, that they are probably related to subduction. We suggest that the alkaline character of the TVF lavas is the result of low degrees of melting of a mantle source coupled with a stress regime that allows these small-volume melts to reach the surface in the TVF.     

Detecting short-term evolution of Etnean scoria cones: a LIDAR-based approach

The 2001 and 2002–2003 flank eruptions on Mount Etna (Italy) were characterized by intense explosive activity which led to the formation of two large monogenetic scoria cones (one from each eruption) on the upper southern flank of the volcano. Continuous monitoring of Etna, especially during flank eruptions, has provided detailed information on the growth of these cones. They differ in genesis, shape, and size. A set of high resolution (1 m) digital elevation models (DEMs) derived from light detection and ranging (LIDAR) data collected during four different surveys (2004, 2005, 2006, and 2007) has been used to map morphology and to extract the morphometric parameters of the scoria cones. By comparing LIDAR-derived DEMs with a pre-eruption (1998) 10 m DEM, the volume of the two scoria cones was calculated for the first time. Comparison of the LIDAR-derived DEMs revealed in unprecedented detail morphological changes during scoria cone degradation. In particular, the morphologically more exposed and structurally weaker 2002–2003 cone was eroded rapidly during the first few years after its emplacement mainly due to gravitational instability of slopes and wind erosion.     

Vulnerability of water resources,vegetation productivity and soil erosion to climate change in Mediterranean watersheds

Climate change is expected to increase temperatures and lower rainfall in Mediterranean regions; however, there is a great degree of uncertainty as to the amount of change. This limits the prediction capacity of models to quantify impacts on water resources, vegetation productivity and erosion. This work circumvents this problem by analysing the sensitivity of these variables to varying degrees of temperature change (increased by up to 6·4 °C), rainfall (reduced by up to 40%) and atmospheric CO₂ concentrations (increased by up to 100%). The SWAT watershed model was applied to 18 large watersheds in two contrasting regions of Portugal, one humid and one semi‐arid; incremental changes to climate variables were simulated using a stochastic weather generator. The main results indicate that water runoff, particularly subsurface runoff, is highly sensitive to these climate change trends (down by 80%). The biomass growth of most species showed a declining trend (wheat down by 40%), due to the negative impacts of increasing temperatures, dampened by higher CO₂ concentrations. Mediterranean species, however, showed a positive response to milder degrees of climate change. Changes to erosion depended on the interactions between the decline in surface runoff (driving erosion rates downward) and biomass growth (driving erosion rates upward). For the milder rainfall changes, soil erosion showed a significant increasing trend in wheat fields (up to 150% in the humid watersheds), well above the recovery capacity of the soil. Overall, the results indicate a shift of the humid watersheds to acquire semi‐arid characteristics, such as more irregular river flows and increasingly marginal conditions for agricultural production. Copyright © 2007 John Wiley & Sons, Ltd.     

Morphology and mechanism of eruption of postglacial shield volcanoes in Iceland

Postglacial Icelandic shield volcanoes were formed in monogenetic eruptions mainly in the early Holocene epoch. Shield volcanoes vary in their cone morphology and in the areal extent of the associated lava flows. This paper presents the results of a study of 24 olivine tholeiite and 7 picrite basaltic shield volcanoes. For the olivine tholeiitic shields the median slope is 2.7°, the median height 60 m, the median diameter 3.6 km, the median aspect ratio (height against diameter) 0.019, and the median cone volume 0.2 km³. The picritic shield volcanoes are considerably steeper and smaller. A shield-volcano cone forms from successive lava lake overflows which are of shelly-type pahoehoe. A widespread apron surrounding the cone forms from tube-fed P-type pahoehoe. The slopes of the cones have (a) a planar or slightly convex form, (b) a concave form, or (c) a convex-concave form. A successive stage of a shield volcano is determined on the basis of cone morphology and lava assemblages. A shield-producing eruption has alternating episodes of lava lake overflows and tube-fed delivery to the distal parts of the flow field. In the late stages of eruption, the cone volume increases in response to the increased amount of rootless outpouring on the cone flanks. Normally, only a small percentage of the total erupted volume of a shield volcano, sometimes as little as 1–3%, is in the shield volcano cone itself, the main volume being in the apron of the shield.