古尔班通古特沙漠季节性冻土入渗特性试验研究

以古尔班通古特沙漠南缘融雪期冻结土壤入渗试验为依据,分析了沙漠地区季节性冻土水分入渗特性及主要影响因素.结果表明:沙漠地区季节性冻土具有较高的入渗能力,沙漠冻结风沙土的稳渗率为0.26～0.30 mm/min,是田间冻结壤土的10～20倍,可以保证融雪水及时入渗进入土壤,为融雪水的高效储存创造了有利条件.入渗能力随着土壤含水率的升高而减小.沙地土壤初始的低含水率、土壤大孔隙结构特征是沙地冻结土壤具有较高入渗能力的主要原因.研究结果为进一步研究春季融雪水在沙地的再分配过程奠定了基础.     

Quaternary landscape evolution and erosion rates for an intramontane Neogene basin (Guadix–Baza basin, SE Spain)

The landscape evolution in Neogene intramontane basins is a result of the interaction of climatic, lithologic, and tectonic factors. When sedimentation ceases and a basin enters an erosional stage, estimating erosion rates across the entire basin can offer a good view of landscape evolution. In this work, the erosion rates in the Guadix–Baza basin have been calculated based on a volumetric estimate of sediment loss by river erosion since the Late Pleistocene. To do so, the distribution of a glacis surface at ca. 43 kyr, characterised by a calcrete layer that caps the basin infilling, has been reconstructed. To support this age, new radiometric data of the glacis are presented. The volume of sediment loss by water erosion has been calculated for the entire basin by comparing the reconstructed geomorphic surface and the present-day topography. The resulting erosion rates vary between 4.28 and 6.57 m³ ha^− 1 yr^− 1, and are the consequence of the interaction of climatic, lithologic, topographic, and tectonic factors. Individual erosion rates for the Guadix and Baza sub-basins (11.80 m³ ha^− 1 yr^− 1 and 1.77 m³ ha^− 1 yr^− 1 respectively) suggest different stages of drainage pattern evolution in the two sub-basins. We attribute the lower values obtained in the Baza sub-basin to the down-throw of this sub-basin caused by very recent activity along the Baza fault.     

Comparison of long-term and short-term uplift rates along an active blind reverse fault zone (Bajo Segura,Se Spain)

The seismic activity of the Guardamar-Torrevieja zone (Eastern Betic Cordillera, SE Spain) can be associated with the Bajo Segura fault zone, an E-W reverse blind fault with secondary NW-SE dextral faults. A high-precision levelling profile 30 kilometers long was set up and levelled in 1997 to monitor the vertical displacement of this active fault zone. This profile runs parallel to an older high-precision levelling line included in the Spanish first order levelling network measured by the Instituto Geográfico Nacional (IGN) in two different campaigns (1934 and 1976). The 1997 line was relevelled in 2003 and 10 new benchmarks were set up, both to increase benchmark density and to restore some of the 1997 benchmarks that had been destroyed. We have used the 1976 IGN and the 2003 measurement to construct a recent vertical movements profile with a significant time difference (27 years). This recent vertical movements profile shows that the vertical movements are very small, nearly equal to the error bars, with a 0.2 mm/year rise in the town of Guardamar, and a 0.2 mm/year subsidence of the southern part of the profile (Punta Prima) respect to the town of Torrevieja. These movements could be related to the activity of the Bajo Segura and the San Miguel faults respectively. Using geological markers we have deduced uplift rates of 0.1 mm/year during the last 3 million years. Therefore, these preliminary results indicate that geodetically (short-term deformation) determined uplift rates are similar to those estimated from geological markers (long term deformation).     

Perchlorate accumulation and migration in the deep vadose zone in a semiarid region

For 25 years, a plant in Israel manufacturing ammonium perchlorate disposed of untreated wastewater in four unlined ponds. This study explores the transport mechanisms of perchlorate infiltrated from 1965 to 1990 from one of these active storage ponds into a deep (40 m) layered vadose zone and the underlying Israeli coastal aquifer. Perchlorate migration from 1990, when wastewater disposal ceased, until today, with infiltration due only to natural rain (500 mm y⁻¹), was also studied. Several indirect methods were used, including: mass balance in the unsaturated zone profile, δ¹⁸O and δ²H profiles below the pond, and a comparison of the same sediment profiles in 2005 and 2007. The isotopic composition of the pore water could be divided into two separate groups: lighter (depleted) and heavier (enriched) samples. All samples in the lighter group were from the shallow vadose zone, above two clayey layers, and represent natural infiltration of rainwater. The enriched samples were from the deeper section of the unsaturated zone (20–40 m) and represent water used for perchlorate manufacturing 14 years prior to drilling. Consequently, the overall maximum infiltration rate was estimated to be 1.4 m y⁻¹. Below the clayey layer almost identical perchlorate concentrations were found along the sediment profile in 2005 and 2007 (two boreholes, 3 m apart). Very different perchlorate profiles were observed above the clayey layers. This suggests that perchlorate below the clay layers (20–40 m) is practically stagnant under the current natural conditions. The reduction in perchlorate concentration in groundwater below the ponds vs. its increased concentration further downgradient supports the contention that the current migration of perchlorate from the vadose zone to the groundwater is very small. We estimate that perchlorate concentration in the groundwater under the infiltration pond, which was 187 mg l⁻¹ in 2004, will reach 10 μg l⁻¹ within about 14 years. The existence of a clayey layer crossing the thick vadose zone was thus found to significantly change the infiltration rate when ponded conditions were replaced with natural precipitation.     

Construction dynamics of a lava channel

We use a kinematic GPS and laser range finder survey of a 200 m-long section of the Muliwai a Pele lava channel (Mauna Ulu, Kilauea) to examine the construction processes and flow dynamics responsible for the channel–levee structure. The levees comprise three packages. The basal package comprises an 80–150 m wide ′a′a flow in which a ∼2 m deep and ∼11 m wide channel became centred. This is capped by a second package of thin (<45 cm thick) sheets of pahoehoe extending no more than 50 m from the channel. The upper-most package comprises localised ′a′a overflows. The channel itself contains two blockages located 130 m apart and composed of levee chunks veneered with overflow lava. The channel was emplaced over 50 h, spanning 30 May–2 June, 1974, with the flow front arriving at our section (4.4 km from the vent) 8 h after the eruption began. The basal ′a′a flow thickness yields effusion rates of 35 m³ s⁻¹ for the opening phase, with the initial flow advancing across the mapped section at ∼10 m/min. Short-lived overflows of fluid pahoehoe then built the levee cap, increasing the apparent channel depth to 4.8 m. There were at least six pulses at 90–420 m³ s⁻¹, causing overflow of limited extent lasting no more than 5 min. Brim-full flow conditions were thus extremely short-lived. During a dominant period of below-bank flow, flow depth was ∼2 m with an effusion rate of ∼35 m³ s⁻¹, consistent with the mean output rate (obtained from the total flow bulk volume) of 23–54 m³ s⁻¹. During pulses, levee chunks were plucked and floated down channel to form blockages. In a final low effusion rate phase, lava ponded behind the lower blockage to form a syn-channel pond that fed ′a′a overflow. After the end of the eruption the roofed-over pond continued to drain through the lower blockage, causing the roof to founder. Drainage emplaced inflated flows on the channel floor below the lower blockage for a further ∼10 h. The complex processes involved in levee–channel construction of this short-lived case show that care must be taken when using channel dimensions to infer flow dynamics. In our case, the full channel depth is not exposed. Instead the channel floor morphology reflects late stage pond filling and drainage rather than true channel-contained flow. Components of the compound levee relate to different flow regimes operating at different times during the eruption and associated with different effusion rates, flow dynamics and time scales. For example, although high effusion rate, brim-full flow was maintained for a small fraction of the channel lifetime, it emplaced a pile of pahoehoe overflow units that account for 60% of the total levee height. We show how time-varying volume flux is an important parameter in controlling channel construction dynamics. Because the complex history of lava delivery to a channel system is recorded by the final channel morphology, time-varying flow dynamics can be determined from the channel morphology. Developing methods for quantifying detailed flux histories for effusive events from the evidence in outcrop is therefore highly valuable. We here achieve this by using high-resolution spatial data for a channel system at Kilauea. This study not only indicates those physical and dynamic characteristics that are typical for basaltic lava flows on Hawaiian volcanoes, but also a methodology that can be widely applied to effusive basaltic eruptions.     

Discontinuous Galerkin methods for advective transport in single-continuum models of fractured media

Accurate simulation of flow and transport processes in fractured rocks requires that flow in fractures and shear zones to be coupled with flow in the porous rock matrix. To this end, we will herein consider a single-continuum approach in which both fractures and the porous rock are represented as volumetric objects, i.e., as cells in an unstructured triangular grid with a permeability and a porosity value associated with each cell. Hence, from a numerical point of view, there is no distinction between flow in the fractures and the rock matrix. This enables modelling of realistic cases with very complex structures. To compute single-phase advective transport in such a model, we propose to use a family of higher-order discontinuous Galerkin methods. Single-phase transport equations are hyperbolic and have an inherent causality in the sense that information propagates along streamlines. This causality is preserved in our discontinuous Galerkin discretization. We can therefore use a simple topological sort of the graph of discrete fluxes to reorder the degrees-of-freedom such that the discretized linear system gets a lower block-triangular form, from which the solution can be computed very efficiently using a single-pass forward block substitution. The accuracy and utility of the resulting transport solver is illustrated through several numerical experiments.     

High-resolution dating of Co-rich crusts: A comparative study using the methods of orbital pacing and <Superscript>230</Superscript>Th<Subscript>ex</Subscript>/<Superscript>232</Superscript>Th dating

Due to their slow growth rates, seamount Co-rich crusts are very difficult to date with high resolution and precision. This paper is to test the use of orbital pacing on the growth profile of crusts to determine high-resolution age and growth rate. Crust CB14 from the central Pacific Ocean was selected for this study. We first examined the growth pattern in detail under a reflected-light microscope and ascertained that the growth environment was stable for the sub-layer 1 (0–3 mm). We then used electron microprobe line-scanning to obtain elemental profiles. The pattern of the power spectrum analysis of the Al-profile revealed that there are significant cycles of 113.9, 87.8, 51.5, 42.2 and 25.8 μm. These cycles correspond to the Milankovitch cycles of 53.1, 41, 24, 19.7 and 12 ka, respectively, and yield the growth rate of about 2.14 mm/Ma and an age of about 1.40 Ma for the boundary between the sub-layer 1 and sub-layer 2. We also used a drilling machine with a numerically controlled drive to obtain high-resolution samples at 0.1mm intervals, and used the ²³⁰Th_ex/²³²Th method to date the samples. For the uppermost 1.3 mm, the growth rate was about 2.15 mm/Ma, and the age for the layer at the depth of 3 mm was about 1.40 Ma, which coincides perfectly with the results obtained from orbital pacing. Thus, it is considered that orbital pacing is a new and effective method to determine the growth rate of the seamount Co-rich crust. This method is applicable for establishing a high-resolution age frame for the crusts of the world’s oceans. Supported by China Ocean Mineral Resources R & P Association (Grant No. DY105-01-01-08) and National Natural Science Foundation of China (Grant Nos. 40106005, 40476050)     

High-resolution dating of Co-rich crusts: A comparative study using the methods of orbital pacing and ~(230)Th_(ex)/~(232)Th dating

Due to their slow growth rates, seamount Co-rich crusts are very difficult to date with high resolution and precision. This paper is to test the use of orbital pacing on the growth profile of crusts to determine high-resolution age and growth rate. Crust CB14 from the central Pacific Ocean was selected for this study. We first examined the growth pattern in detail under a reflected-light microscope and ascertained that the growth environment was stable for the sub-layer 1 (0-3 mm). We then used electron mic...     

Probabilistic seismic hazard assessment of Switzerland: best estimates and uncertainties

We present the results of a new genera tion of probabilistic seismic hazard assessment for Switzerland. This study replaces the previous intensity-based generation of national hazard maps of 1978. Based on a revised moment-magnitude earthquake catalog for Switzerland and the surrounding regions, covering the period 1300–2003, sets of recurrence parameters (a and b values, M _max ) are estimated. Information on active faulting in Switzerland is too sparse to be used as source model. We develop instead two models of areal sources. The first oriented towards capturing historical and instrumental seismicity, the second guided largely by tectonic principles and express ing the alterative view that seismicity is less stationary and thus future activity may occur in previously quiet regions. To estimate three alterna tive a and b value sets and their relative weighting, we introduce a novel approach based on the modified Akaike information criterion, which allows us to decide when the data in a zone deserves to be fitted with a zone-specific b value. From these input parameters, we simulate synthetic earthquake catalogs of one-million-year duration down to magnitude 4.0, which also reflect the difference in depth distribution between the Alpine Foreland and the Alps. Using a specific predictive spectral ground motion model for Switzerland, we estimate expected ground motions in units of the 5% damped acceleration response spectrum at frequencies of 0.5–10 Hz for all of Switzerland, referenced to rock sites with an estimated shear wave velocity of 1,500 m/s² in the upper 30 m. The highest hazard is found in the Wallis, in the Basel region, in Graubünden and along the Alpine front, with maximum spectral accelerations at 5 Hz frequency reaching 150 cm/s² for a return period of 475 years and 720 cm/s² for 10,000 years. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users.     

Volcanic evolution of Volcán Aucanquilcha: a long-lived dacite volcano in the Central Andes of northern Chile

Volcán Aucanquilcha, northern Chile, has produced ∼37 km³ of dacite (63–66 wt% silica), mainly as lavas with ubiquitous magmatic inclusions (59–62 wt% silica) over the last ∼1 million years. A pyroclastic flow deposit related to dome collapse occurs on the western side of the edifice and a debris avalanche deposit occurs on the eastern side. The >6,000-m high edifice defines a 9-km E–W ridge and lies at the center of a cluster of more than 15 volcanoes, the Aucanquilcha Volcanic Cluster, that has been active for at least the past 11 million years. The E–W alignment of vents is nearly orthogonal to the arc axis. A majority of Volcán Aucanquilcha was constructed during the first 200,000 years of eruption, whereas the last 800,000 years have added little additional volume. The peak eruptive rate during the edifice-building phases was ∼0.16 km³/ka and the later eruptive rate was ∼0.02 km³/ka. Comparable dacite volcanoes elsewhere show a similar pattern of high volcanic productivity during the early stages and punctuated rather than continuous activity. Volcán Aucanquilcha lavas are dominated by phenocrysts of plagioclase, accompanied by two populations of amphibole, biotite, clinopyroxene, Fe–Ti oxides and (or) orthopyroxene. Accessory phases include zircon, apatite and rare quartz and sanidine. One amphibole population is pargasite and the other is hornblende. The homogeneity of dacite lava from Volcán Aucanquilcha contrasts with the heterogeneity (52–66 wt% silica) at nearby Volcán Ollagüe, which has been active over roughly the same period of time. We attribute this homogeneity at Aucanquilcha to the thermal development of the crust underneath the volcano resulting from protracted magmatism there, whereas Volcán Ollagüe lacks this magmatic legacy.