The origin,typology, spatial distribution and detrimental effects of the sinkholes developed in the alluvial evaporite karst of the Ebro River valley downstream of Zaragoza city (NE Spain)

Three types of sinkhole have been mapped in a 50 km² stretch of the Ebro River valley downstream of Zaragoza: large collapse sinkholes, large shallow subsidence depressions and small cover-collapse sinkholes. The sinkholes relate to the karstification of evaporitic bedrock that wedges out abruptly downstream, giving way to a shale substratum. Twenty-three collapse sinkholes, up to 50 m in diameter by 6 m deep, and commonly hosting saline ponds, have been identified in the floodplain. They have been attributed to the upward stoping of dissolutional cavities formed within the evaporitic bedrock by rising groundwater flows. Twenty-four large shallow subsidence depressions were mapped in the floodplain. These may reach 850 m in length and were formed by structurally controlled interstratal karstification of soluble beds (halite or glauberite? and gypsum) by rising groundwater flow and the progressive settlement of the overlying bedrock and overburden sediments. A total of 447 small cover-collapse, or dropout, sinkholes have been recognized in a perched alluvial level along the southern margin of the valley. These sinkholes result from the upward propagation of voids through the alluvial mantle caused by the downward migration of detrital sediments into dissolutional voids. The majority of these sinkholes, commonly 1·5–2 m in diameter, are induced by human activities. Over the karstic bedrock, there is a significant increase in sinkhole density downstream. This is interpreted as being a result of the evaporitic bedrock wedging out and the convergence of the groundwater flow lines in the karstic aquifer. The collapse sinkholes in this area, locally with a probability of occurrence higher than 140 sinkholes/km²/year, cause substantial damage to the linear infrastructures, buildings and agriculture, and they might eventually cause the loss of human lives. Copyright © 2006 John Wiley & Sons, Ltd.     

A first estimate of ground water ages for the deep aquifer of the Kathmandu Basin, Nepal, using the radioisotope chlorine-36

The Kathmandu Basin in Nepal contains up to 550 m of Pliocene-Quaternary fluvio-lacustrine sediments which have formed a dual aquifer system. The unconfined sand and gravel aquifer is separated by a clay aquitard, up to 200 m thick, from the deeper, confined aquifer, comprised of Pliocene sand and gravel beds, intercalated with clay, peat, and lignite. The confined aquifer currently provides an important water supply to the central urban area but there are increasing concerns about its sustainability due to overexploitation. A limited number of determinations of the radioisotope 36Cl have been made on bore waters in the basin, allowing us to postulate on the age of ground water in the deeper, confined aquifer. Ground water evolution scenarios based on radioisotope decay, gradual dissolution of formational salts as the ground waters move downgradient, and flow velocity estimations produce comparable ground water ages for the deep waters, ranging from 200,000 to 400,000 years. From these ages, we deduce a mean ground water flow velocity of only 45 mm/year from recharge in the northeast to the main extraction region 15 km to the southwest. We thus estimate current recharge at about 5 to 15 mm/year, contributing 40,000 to 1.2 million m3/year to the ground water system. Current ground water extraction is estimated to be 20 times this amount. The low specific discharge confirms that the resource is being mined, and, based on current projections, reserves will be used up within 100 years.     

Rippled scour depressions on continental shelf bank slopes off Nordland and Troms,Northern Norway

Multibeam bathymetry acquired under the MAREANO programme from the continental shelf off Nordland and Troms, northern Norway, show bedforms that we have interpreted as rippled scour depressions. They occur in three areas offshore on bank slopes facing southeast, more than 15 km from land. They are generally found where the slope gradient is low, in water depths of 70–160 m. Individual depressions are up to 3 km long, 1 m deep and up to 300 m wide. They occur in areas where sediments evolve quickly from glacial deposits on the banks to post-glacial muddy sediments on the glacial troughs. Multibeam backscatter and underwater video data show that depression floors are covered by rippled, gravelly, shelly sand. Ripple crests are parallel or slightly oblique to the depression axis orientation. Sand without bedforms is observed between the depressions. TOPAS seismic lines show that the uppermost seismic unit consists of the sand between the depressions. The base of this unit may be the last transgressive/tidal/wave ravinement surface. Physical oceanographic modelling indicates that maximum current velocities are up to 0.6 m/s in the rippled scour depression areas. Stronger currents appear to inhibit the building of these features. Tidal currents play an important role as they trend parallel to the southeast banks slopes and are likely responsible of the gravelly ripples formation inside the depressions as well as the persistence of these depressions which are not covered by finer sediments. On Malangsgrunnen bank, some of the rippled scour depressions are in the extension of NW–SE furrows located on the bank. Simulated bottom currents indicate currents mainly perpendicular to these furrows, as for the rippled scour depressions on the bank slopes. Nevertheless, these features could also highlight currents coming from the northwest which reach the bank margin and continue down to the areas of the rippled scour depressions. These currents could be responsible for the formation of some of the bedforms, together with tidal currents.     

Solution-collapse depressions and suspensates in the limnocrenic lake of banyoles (ne Spain)

A high resolution seismic reflection survey in the Banyoles limnocrenic solution lake allowed penetration of dense suspensates occupying cone-like bottom depressions of different size. The depressions result from the dissolution and collapse of underlying Eocene calcareous and gypsiferous materials over which lacustrine sediments of varying thicknesses have accumulated. The suspensates occupying the depressions present three main types of seismic signatures: stratified, semistratified, and transparent. The densities of the suspensates and the water depths of their tops, which fluctuate continuously, vary from one depression to another. The maximum seismically recorded suspensate thickness is 44 m. Morphological and structural features, seismic characters, and variable degrees of hydraulic activity, point to the existence of different stages of maturity in the lake bottom depressions. This work brings new insight on the dynamics and evolution of limnocrenic solution lakes.     

Detection of preferential infiltration pathways in sinkholes using joint inversion of self-potential and EM-34 conductivity data

The percolation of water in the ground is responsible for measurable electric potentials called self‐potentials. These potentials are influenced by the distribution of the electrical conductivity of the ground. Because sinkholes are associated both with self‐potential and electrical conductivity anomalies, a joint inversion of EM‐34 conductivity and self‐potential data is proposed as a way of delineating the location of these features. Self‐potential and EM conductivity data were obtained at a test site in Normandy (France) where sinkholes and crypto‐sinkholes are present over a karstic area in a chalk substratum overlain by clay‐with‐flint and loess covers. The presence of sinkholes and crypto‐sinkholes is associated with negative self‐potential anomalies with respect to a reference electrode located outside the area where the sinkholes are clustered. The sinkholes also have a conductivity signature identified by the EM‐34 conductivity data. We used the simulated‐annealing method, which is a global optimization technique, to invert jointly EM‐34 conductivity and self‐potential data. Self‐potential and electrical conductivity provide clear complementary information to determine the interface between the loess and clay‐with‐flint formations. The sinkholes and crypto‐sinkholes are marked by depressions in this interface, focusing the groundwater flow towards the aquifer contained in the chalk substratum.     

TEM study of the geoelectrical structure and groundwater salinity of the Nahal Hever sinkhole site, Dead Sea shore, Israel

Since the 1990s a large number of sinkholes have appeared in the Dead Sea (DS) coastal area. Sinkhole development was triggered by the lowering of the DS level. In the literature the relationship between the sinkholes and the DS level is explained by intrusion of relatively fresh water into the aquifer thereby dramatically accelerating the salt dissolution with creation of subsurface caverns, which in turn cause sinkholes. The main goal of our project was detection and localization of relatively fresh groundwater. During our study we used the transient electromagnetic method (TEM) to measure the electrical resistivity of the subsurface. As a test site we selected Nahal Hever South which is typical for the DS coast. Our results show that resistivity of the shallow subsurface reflects its vertical and lateral structure, e.g., its main hydrogeological elements explain the inter-relations between geology, hydrogeology, and sinkholes. The TEM method has allowed detailed differentiation of layers (clay, salt, etc.) in the subsurface based on their bulk resistivity. The 10 m-thick salt layer composed of idiomorphic crystals of halite deposited during the earlier Holocene period was extrapolated from borehole HS-2 through the study area. It was found that in Nahal Hever the typical value of the bulk resistivity of clay saturated with the DS brine varies between 0.2 and 0.3 Ωm, whereas saturated gravel and sandy sediments are characterized by resistivity between 0.4 and 0.6 Ωm. The high water salinity of the aquifer (enveloping the salt layer) expressed in terms of resistivity is also an important characterization of the sinkhole development mechanism. The electrical resistivity of the aquifer in the vicinity of the salt unit and its western border did not exceed 1 Ωm (in most cases aquifer resistivity was 0.2-0.6 Ωm) proving that, in accordance with existing criteria, the pores of the alluvial sediments are filled with highly mineralized DS brine. However, we suggest that the criterion of the aquifer resistivity responsible for the salt dissolution should be decreased from 1 Ωm to 0.6 Ωm corresponding to the chloride concentration of approximately 100 g/l (the chloride saturation condition reaches 224 g/l in the northern DS basin and 280 g/l in the southern one).Based on TEM results we can reliably conclude that in 2005, when most of sinkholes had appeared at the surface, salt was located within a very low resistivity environment inside sediments saturated with DS brine. Intrusion of relatively fresh groundwater into the aquifer through the 600 × 600 m²area affected by sinkholes has not been observed.     

Imaging a shallow salt diapir using ambient seismic vibrations beneath the densely built-up city area of Hamburg,Northern Germany

Salt diapirs are common features of sedimentary basins. If close to the surface, they can bear a significant hazard due to possible dissolution sinkholes, karst formation and collapse dolines or their influence on ground water chemistry. We investigate the potential of ambient vibration techniques to map the 3-D roof morphology of shallow salt diapirs. Horizontal-to-vertical (H/V) spectral peaks are derived at more than 900 positions above a shallow diapir beneath the city area of Hamburg, Germany, and are used to infer the depth of the first strong impedance contrast. In addition, 15 small-scale array measurements are conducted at different positions in order to compute frequency-dependent phase velocities of Rayleigh waves between 0.5 and 25 Hz. The dispersion curves are inverted together with the H/V peak frequency to obtain shear-wave velocity profiles. Additionally, we compare the morphology derived from H/V and array measurements to borehole lithology and a gravity-based 3-D model of the salt diapir. Both methods give consistent results in agreement with major features indicated by the independent data. An important result is that H/V and array measurements are better suited to identify weathered gypsum caprocks or gypsum floaters, while gravity-derived models better sample the interface between sediments and homogeneous salt. We further investigate qualitatively the influence of the 3-D subsurface topography of the salt diapir on the validity of local 1-D inversion results from ambient vibration dispersion curve inversion.     

Probabilistic sinkhole modelling for hazard assessment

A method for quantitatively assessing sinkhole susceptibility (spatial probability) and hazard (spatio‐temporal probability) has been developed and independently tested in a 50 km² sector of the Ebro Valley evaporite karst. Three genetic types of sinkholes have been mapped in the floodplain and a terrace surface: 947 small cover‐collapse sinkholes (type 1, terrace), large collapse sinkholes (type 2, floodplain) and large subsidence depressions (type 3, floodplain). The type 1 sinkhole inventory includes two temporal populations: 447 sinkholes formed before 24 November 2005, and 500 between that date and 2 November 2006. Sinkhole susceptibility models have been elaborated analysing the statistical relationships between the sinkholes of the 2005 inventory and a set of potential conditioning factors. The independent evaluation (validation) of the susceptibility models by means of several strategies (random, sequentially excluded, and temporal) has allowed us to select the most significant variables for each sinkhole type and assess quantitatively the quality of models; which are reasonable for the three sinkhole types. Validation has also provided information on the contribution of specific variables and the effect of changing their accuracy to the prediction capability of models. Susceptibility models for type 3 sinkholes have been validated satisfactorily with the 2006 sinkhole inventory (temporal validation). The best susceptibility model has been transformed into a hazard map considering the frequency of sinkholes that occurred in each susceptibility class between 2005 and 2006, as well as their average size. The susceptibility and hazard models obtained could be used as an objective basis for the application of mitigation measures, either of preventive or corrective nature. Copyright © 2008 John Wiley & Sons, Ltd.     

Geoelectric structure of the Ein Gedi sinkhole occurrence site at the Dead Sea shore in Israel

Hazardous sinkholes started to appear in alluvial fans and unconsolidated sediments along the western Dead Sea coast in 1990. Since then hundreds of sinkholes have appeared from north to south along the shoreline. The Electrical Resistivity Tomography (ERT) method was used to achieve a better understanding of the subsurface geoelectric structure at the sinkhole development sites, taking into account that electric parameters (such as resistivity or conductivity) are very sensitive to formation properties and their variations in time. Fifteen image lines were surveyed at the Ein Gedi area during a period of active sinkhole development (in 2001–2002) over an area of 300 × 550 m². Resistivity cross-sections and maps were constructed from 2-D linear surveys. The process of sinkhole formation in the surveyed area is located in a strip 50–70 m wide and 300–500 m long, extending approximately in a north–south direction. The sinkholes are arranged along a tortuous line within this strip. On resistivity maps and sections this U-shaped zone appears as an alternation of high resistivity anomalies of 350–1000 Ωm (at sinkhole group locations) with narrow background resistivity zones of 50–100 Ωm. The large size of resistivity anomalies (250 × 300 m²), which are considerably greater than those of the sinkholes, form one of the features of the sinkhole sites in the Ein Gedi area. The anomalies continue down to the water table or even deeper (maximum of 25–35 m depth). A low resistivity layer of 1–8 Ωm underlies them. The combined analysis of the image results and other geophysical data shows that high resistivity anomalies are associated with the decompaction of the soil mass at the sinkhole development sites and surrounding areas. Recent studies have shown that sinkholes in the Ein Gedi area are developing along the salt western edge located at a depth of 50 m. The subsurface high resistivity anomaly conforms to the sinkhole line (and salt boundary). They are presumably located above the great dissolution caverns at the salt edge. The heterogeneity of the resistivity structure within the high resistivity anomaly (seen in both lateral and vertical planes) confirms that a disintegration of internal formation structure takes place. Away from the sinkhole sites the subsurface resistivity distribution is homogeneous.     

Exploring landscape and geologic controls on spatial patterning of streambank groundwater discharge in a mixed land use watershed

Preferential groundwater discharge features along stream corridors are ecologically important at local and stream network scales, yet we lack quantification of the multiscale controls on the spatial patterning of groundwater discharge. Here we identify physical attributes that best explain variation in the presence and lateral extent of preferential groundwater discharges along two 5th order streams, the Housatonic and Farmington Rivers, and 32 1st to 4th order reaches across the Farmington River network. We mapped locations of preferential groundwater discharge exposed along streambanks using handheld thermal infrared cameras paired with high-resolution topographic and land use land cover datasets, surficial soil characteristic maps, and depth-to-bedrock geophysical measurements. The unconfined Housatonic River, MA, USA (12 km) had fewer discharge locations and less lateral extent (41 discharge locations with 38 m of active discharge/km of river) compared to the partially confined Farmington River, CT, USA (26 km; 169 discharge locations with 129 m of active discharge/km of river). Using a moving window analysis, we found along both rivers that discharge was more likely to occur where bank slopes were steeper, floodplain extent was narrower, and degree of confinement was higher. Along the Farmington River, groundwater discharge was more likely to occur where saturated hydraulic conductivity was higher and depth-to-bedrock was shallower. Among the 32 stream reaches surveyed (33.2 km of total stream length) within the Farmington River watershed, preferential discharge was observed in all but two stream reaches, varied from 0 to 25% of lateral extent along stream banks (mean = 6%), and was more likely to occur where stream reach slopes were steep, saturated hydraulic conductivity was high, and watershed urbanization was low. Our results show that, though both surface (e.g., topographic, land use land cover) and subsurface (e.g., soil characteristics, bedrock depth) factors control the prevalence of streambank preferential groundwater discharge, the dominant controls vary across valley settings and stream sizes.     

Using δ18O and δ2H to Detect Hydraulic Connection Between a Sinkhole Lake and a First-Magnitude Spring

Oxygen and hydrogen isotopes were used in this study to detect a hydraulic connection between a sinkhole lake and a karst spring. In karst areas, surface water that flows to a lake can drain through sinkholes in the lakebed to the underlying aquifer, and then flows in karst conduits and through aquifer matrix. At the study site located in northwest Florida, USA, Lake Miccosukee immediately drains into two sinkholes. Results from a dye tracing experiment indicate that lake water discharges at Natural Bridge Spring, a first-magnitude spring 32 km downgradient from the lake. By collecting weekly water samples from the lake, the spring, and a groundwater well 10 m away from the lake during the dry period between October 2019 and January 2020, it was found that, when rainfall effects on isotopic signature in spring water are removed, increased isotope ratios of spring water can be explained by mixing of heavy-isotope-enriched lake water into groundwater, indicating hydraulic connection between the lake and the spring. Such a detection of hydraulic connection at the scale of tens of kilometers and for a first-magnitude spring has not been previously reported in the literature. Based on the isotope ratio data, it was estimated that, during the study period, about 8.5% the spring discharge was the lake water that drained into the lake sinkholes.     

Role of ground water in geomorphology,geology, and paleoclimate of the Southern High Plains,USA

Study of ground water in the Southern High Plains is central to an understanding of the geomorphology, deposition of economic minerals, and climate change record in the area. Ground water has controlled the course of the Canadian and Pecos rivers that isolated the Southern High Plains from the Great Plains and has contributed significantly to the continuing retreat of the westward escarpment. Evaporative and dissolution processes are responsible for current plateau topography and the development of the signature 20,000 small playa basins and 40 to 50 large saline lake basins in the area. In conjunction with eolian processes, ground water transport controls the mineralogy of commercially valuable mineral deposits and sets up the distribution of fine efflorescent salts that adversely affect water quality. As the water table rises and retreats, lunette and tufa formation provides valuable paleoclimate data for the Southern High Plains. In all these cases, an understanding of ground water processes contributes valuable information to a broad range of geological topics, well beyond traditional interest in water supply and environmental issues.     

Salt karst and tectonics: sinkholes development along tension cracks between parallel strike‐slip faults,Dead Sea,Jordan

This work deals with the tectonic interpretation of an alignment of more than 300 sinkholes stretching along the Jordanian coast of the Dead Sea, Ghor Al Haditha area. Its dimensions are 6 km long with a width of 600 m. Sinkholes appeared during the last decades as a consequence of the very rapid lowering of the lake level. The linear shape was inferred from ground collapse inventories carried out between 1991 and 2008. The lineament is replaced and analyzed in its structural setting at regional and local scales. Its direction (N 24° E) is sub‐parallel to the ones displayed by many focal mechanisms, especially the one associated with the earthquake of the 23 April, 1979 (mb = 5·1; N 20° E ± 5°), which is representative of all focal mechanisms calculated on a fault plane compatible with the general direction of the Jordan‐Dead Sea Transform fault system for the east coast of the Dead Sea area. The alignment of sinkholes is constituted by 13 minor linear segments separated by as many empty spaces. Four minor linear units present an en‐echelon arrangement from which one can deduce the presence of a local extensional stress field. In this context, the sinkhole locations provide information of subsurface discontinuities interpreted as hidden fractures. In a close future, such results could support the work of decision‐makers and engineers in the projected development of the area. Copyright © 2009 John Wiley & Sons, Ltd.     

Dilute gravity current and rain-flushed ash deposits in the 1.8 ka Hatepe Plinian deposit,Taupo, New Zealand

Two groups of poorly sorted ash-rich beds, previously interpreted as rain-flushed ashes, occur in the ca. AD 180 Hatepe Plinian pumice fall deposit at Taupo volcano, New Zealand. Two ash beds with similar dispersal patterns and an aggregate thickness of up to 13 cm make up the lowermost group (A). Group A beds extend 45 km north-east of the vent and cover 290 km². In the southern part of the group A distribution area, a coarse ash to lapilli-size Plinian pumice bed (deposit B) separates the two group A beds. The scarcity of lapilli (material seen elsewhere from the still-depositing pumice fall) in group A beds indicates that they were rapidly transported and deposited. However, this rapid transportation and deposition did not produce cross-bedding, nor did it erode the underlying deposits. It is proposed that thick (>600 m) but dilute gravity currents generated from the collapsing outer margin of the otherwise buoyant Hatepe Plinian eruption column deposited the group A beds. The upper ash beds (group C) consist of one to seven layers, attain an aggregate thickness of 35 cm, and vary considerably in thickness and number of beds with respect to distance from vent. Group C beds contain variable amounts of ash mixed with angular Plinian pumices and are genuine rain-flushed ashes. Several recent eruptions at other volcanoes (Ukinrek Maars, Vulcan, Rabaul, La Soufrère de Guadeloupe and Soufrière, St Vincent) have produced gravity currents similar in style, but much smaller than those envisaged for group A deposits. The overloaded margins of otherwise buoyant eruption plumes generated these gravity currents. Laboratory studies have produced experimental gravity current analogues. Hazards from dilute gravity currents are considerable but often overlooked, thus the recognition of gravity current deposits will contribute to more thorough volcanic hazard assessment of prehistoric eruption sequences.     

Morphology and distribution of dolines on ultramafic rocks from airborne LiDAR data: the case of southern Grande Terre in New Caledonia (SW Pacific)

Dolines are closed geomorphological depressions which are surface manifestations of karstic systems. Usually developed on limestones, they also typify the morphology of the New Caledonian landscape, particularly on the southern massif of the main island (known as Massif du Sud). The specificity of dolines here lies in their development on ultramafic rocks. They are evidences of subsidence, suffosion and collapse phenomena resulting from dissolution weathering of peridotites. However, extensive underground drainage systems are still not yet recognized. Semi‐automatic mapping of dolines is carried out on a 148 km² area of the Massif du Sud from a high accuracy LiDAR digital elevation model. In this area 8601 dolines ranging from 1 m² to 2 km² are identified and morphologically characterized with precision. Most are small, shallow and round‐shaped, yet more complex shapes are locally observed. Size distribution analysis allows the setting of a threshold of 20 000 m² above which surface processes rather than chemical weathering control doline evolution. Doline density analysis reveals high concentrations on flat areas where ferricrete overlies the complete weathering profile, especially in the case of elevated rainy watersheds. Dolines are aligned and elongated along a north 135° ± 5° major fracture direction, which is inherited from the obduction of the Pacific Plate upper mantle in the Late Eocene. Finally, we propose a pioneering morphometric typology of dolines that provides important clues as to pseudokarstic activity. We define collapse, bowl‐shaped and flat bottom dolines. Collapse and bowl‐shaped dolines are assumed to denote active pseudokarst. They may widen and deepen, or eventually be filled by sediments. They are distinguished from flat bottom dolines that are partially to completely filled, which suggests that they are associated with paleo‐pseudokarsts. However the groundwater flow paths associated with the genesis and evolution of dolines must be clarified, thus collapse and bowl‐shaped dolines should be hydrologically monitored. Copyright © 2016 John Wiley & Sons, Ltd.     

Caldera formation at Volcán Colima, Mexico, by a large large holocene volcanic debris avalanche

About 4,300 years ago, 10 km³ of the upper cone of ancestral Volcán Colima collapsed to the southwest leaving a horseshoe-shaped caldera 4 km in diameter. The collapse produced a massive volcanic debris avalanche deposit covering over 1550 km² on the southern flanks of the volcano and extending at least 70 km from the former summit. The avalanche followed a steep topographic gradient unobstructed by barriers, resulting in an unusually high area/volume ratio for the Colima deposit. The apparent coefficient of friction (fall height/distance traveled) for the Colima avalanche is 0.06, a low value similar to those of other large-volume deposits. The debris avalanche deposit contains 40–75% angular volcanic clasts from the ancestral cone, a small proportion of vesicular blocks that may be juvenile, and in distal exposures, rare carbonate clasts plucked from the underlying surface by the moving avalanche. Clasts range in size to over 20 m in diameter and are brecciated to different degrees, pulverized, and surrounded by a rock-flour matrix. The upper surface of the deposit shows prominent hummocky topography with closed depressions and surface boulders. A thick, coarse-grained, compositionally zoned scoria-fall layer on the upper northeastern slope of the volcano may have erupted at the time of collapse. A fine-grained surge layer is present beneath the avalanche deposit at one locality, apparently representing an initial blast event. Most of the missing volume of the ancestral volcano has since been restored at an average rate of 0.002 km³/yr through repeated eruptions from the post-caldera cone. As a result, the southern slope of Volcán Colima may again be susceptible to collapse. Over 200,000 people are now living on primary or secondary deposits of the debris avalanche, and a repetition of this event would constitute a volcanic disaster of great magnitude.Ancestral Volcán Colima grew on the southern, trenchward flank of the earlier and larger volcano Nevado de Colima. Trenchward collapse was favored by the buttressing effect of Nevado, the rapid elevation drop to the south, and the intrusion of magma into the southern flank of the ancestral volcano. Other such trenchward-younging, paired volcanoes are known from Mexico, Guatemala, El Salvador, Chile, and Japan. The trenchward slopes of the younger cones are common sites for cone collapse to form avalanche deposits, as occurred at Colima and Popocatepetl in Mexico and at San Pedro Volcano in Chile.     

Segmentation of the southern Mariana back-arc spreading center

SeaBeam multibeam bathymetry obtained during cruise SO-69 of research vessel (R/V) Sonne defines the segmentation and structure of ∼ 300 km of the Mariana back-arc spreading center south of the Pagan fracture zone at 17°33'N. Eight ridge segments, ranging from 14 to 64 km in length, are displaced as much as 2.7–14.5 km by both right- (predominantly) and left-lateral offsets and transform faults. An axial ridge commonly occupies the middle portion of the rift valley and rises from 200 to 700 m above the adjacent sea floor, in places shoaling to a water depth of 3200 m. An exception is the 60-km-long segment between 16°58' and 17°33'N where single peaks only a few tens of meters high punctuate the rift axis. Photographic evidence and rock samples reveal the presence of mostly pillow lavas outcropping on the axial ridges or peaks whereas the deeper parts of the rift valley floor (max. depth 4900 m) are heavily to totally sedimented. Abundant talus ramps along fault scarps testify to ongoing disruption of the crust. Lozenge-shaped collapse structures are covered by layers of sediment up to tens of centimeters thick on the rift valley floor. The presence of discrete volcanic ridges in the southern Mariana back-arc spreading region suggests that emplacement of oceanic crust at this slow spreading center occurs by `multi-site' injection of magma. Along-axis variations in length, crestal depth, and size of the axial ridges can be best explained by different stages in the cyclicity of magma supply along-axis.     

沁水盆地大地热流与地温场特征

分析了20口井的温度数据和39个岩石样品的热导率数据,计算了20个大地热流值.结果表明：沁水盆地具有偏高的大地热流背景,现今热流变化于448～1018 mW/m2之间,平均627±152 mW/m2,现今地温梯度介于209～476 ℃/km之间,平均地温梯度282±103 ℃/km. 热流的分布与本区上、下主煤层含气量的分布特征相似,在一定程度上反映了现代热流继承了煤层气生气期的古地热背景.     

On the morphological characteristics of overdeepenings in high‐mountain glacier beds

Overdeepenings, i.e. closed topographic depressions with adverse slopes in the direction of flow, are characteristic for glacier beds and glacially sculpted landscapes. Quantitative information about their morphological characteristics, however, has so far hardly been available. The present study provides such information by combining the analysis of (a) numerous bed overdeepenings below still existing glaciers of the Swiss Alps and the Himalaya‐Karakoram region modelled with a robust shear stress approximation and (b) detailed bathymetries from recently exposed lakes in the Peruvian Andes. The investigated overdeepenings exist where glacier surface slopes are low (< 5°–10°), occur in bedrock or morainic material and are most commonly a fraction of a kilometre squared in surface area, hundreds of metres long, about half the length in width and tens of metres deep. They form under conditions of low to high basal shear stresses, at cirque, confluence, trunk valley and terminus positions. The most striking phenomenon, however, is the high variability of their geometries: Depths, surface areas, lengths and widths of the overdeepenings vary over orders of magnitude and are only weakly – if at all – interrelated. Inclinations of adverse slopes do not differ significantly from those of forward slopes and are in many cases higher than so far assumed theoretical limits for supercooling of ascending water and corresponding closure of sub‐glacial channels. Such steep adverse slopes are a robust observation and in support of recently developed new concepts concerning the question about where supercooling of sub‐glacial water and closure of ice channels can or must occur. However, the question of when and under what climatic, topographic and ice conditions the overdeepenings had formed remains unanswered. This open question constitutes a key problem concerning the interpretation of observed overdeepenings, the understanding of the involved glacio‐hydraulic processes and the possibility of realistic predictive modelling of overdeepening formation. Copyright © 2016 John Wiley & Sons, Ltd.     

Heat flow and ground water flow in the Great Plains of the United States

Regional groundwater flow in deep aquifers adds advective components to the surface heat flow over extensive areas within the Great Plains province. The regional groundwater flow is driven by topographically controlled piezometric surfaces for confined aquifers that recharge either at high elevations on the western edge of the province or from subcrop contacts. The aquifers discharge at lower elevations to the east. The assymetrical geometry for the Denver and Kennedy Basins is such that the surface areas of aquifer recharge are small compared to the areas of discharge. Consequently, positive advective heat flow occurs over most of the province. The advective component of heat flow in the Denver Basin is on the order of 15 mW m⁻² along a zone about 50 km wide that parallels the structure contours of the Dakota aquifer on the eastern margin of the Basin. The advective component of heat flow in the Kennedy Basin is on the order of 20 mW m⁻² and occurs over an extensive area that coincides with the discharge areas of the Madison (Mississippian) and Dakota (Cretaceous) aquifers. Groundwater flow in Paleozoic and Mesozoic aquifers in the Williston Basin causes thermal anomalies that are seen in geothermal gradient data and in oil well temperature data. The pervasive nature of advective heat flow components in the Great Plains tends to mask the heat flow structure of the crust, and only heat flow data from holes drilled into the crystalline basement can be used for tectonic heat flow studies.