Stability of sandy slopes under seepage conditions

Stability against shallow mass sliding in saturated sandy slopes under seepage depends on the flow direction and hydraulic gradient, particularly near the ground surface. Two modes of instability i.e., Coulomb sliding and liquefaction have been studied and the critical flow directions discussed. The utility of the numerical approach in solving complex flow problems with irregular boundaries and surface topography is demonstrated by means of two slope examples with different internal drainage conditions. The numerical results for the seepage gradients at different points are compared with those predicted by the simple expression derived in this study, and the corresponding effects on the stability are evaluated.     

Comparison of numerical models of two debris flows in the Cortina d’ Ampezzo area,Dolomites, Italy

The accurate prediction of runout distances, velocities and the knowledge of flow rheology can reduce the casualties and property damage produced by debris flows, providing a means to delineate hazard areas, to estimate hazard intensities for input into risk studies and to provide parameters for the design of protective measures. The application of most of models that describe the propagation and deposition of debris flow requires detailed topography, rheological and hydrological data that are not always available for the debris-flow hazard delineation and estimation. In the Cortina d’Ampezzo area, Eastern Dolomites, Italy, most of the slope instabilities are represented by debris flows; 325 debris-flow prone watersheds have been mapped in the geomorphological hazard map of this area. We compared the results of simulations of two well-documented debris flows in the Cortina d’Ampezzo area, carried on with two different single-phase, non-Newtonian models, the one-dimensional DAN-W and the two-dimensional FLO-2D, to test the possibility to simulate the dynamic behaviour of a debris flow with a model using a limited range of input parameters. FLO-2D model creates a more accurate representation of the hazard area in terms of flooded area, but the results in terms of runout distances and deposits thickness are similar to DAN-W results. Using DAN-W, the most appropriate rheology to describe the debris-flow behaviour is the Voellmy model. When detailed topographical, rheological and hydrological data are not available, DAN-W, which requires less detailed data, is a valuable tool to predict debris-flow hazard. Parameters obtained through back-analysis with both models can be applied to predict hazard in other areas characterized by similar geology, morphology and climate.     

Provenance of the Bosnian Flysch

Sandwiched between the Adriatic Carbonate Platform and the Dinaride Ophiolite Zone, the Bosnian Flysch forms a c. 3000 m thick, intensely folded stack of Upper Jurassic to Cretaceous mixed carbonate and siliciclastic sediments in the Dinarides. New petrographic, heavy mineral, zircon U/Pb and fission-track data as well as biostratigraphic evidence allow us to reconstruct the palaeogeology of the source areas of the Bosnian Flysch basin in late Mesozoic times. Middle Jurassic intraoceanic subduction of the Neotethys was shortly followed by exhumation of the overriding oceanic plate. Trench sedimentation was controlled by a dual sediment supply from the sub-ophiolitic high-grade metamorphic soles and from the distal continental margin of the Adriatic plate. Following obduction onto Adria, from the Jurassic–Cretaceous transition onwards a vast clastic wedge (Vranduk Formation) was developed in front of the leading edge, fed by continental basement units of Adria that experienced Early Cretaceous synsedimentary cooling, by the overlying ophiolitic thrust sheets and by redeposited elements of coeval Urgonian facies reefs grown on the thrust wedge complex. Following mid-Cretaceous deformation and thermal overprint of the Vranduk Formation, the depozone migrated further towards SW and received increasing amounts of redeposited carbonate detritus released from the Adriatic Carbonate Platform margin (Ugar Formation). Subordinate siliciclastic source components indicate changing source rocks on the upper plate, with ophiolites becoming subordinate. The zone of the continental basement previously affected by the Late Jurassic–Early Cretaceous thermal imprint has been removed; instead, the basement mostly supplied detritus with a wide range of pre-Jurassic cooling ages. However, a c. 80 Ma, largely synsedimentary cooling event is also recorded by the Ugar Formation, that contrasts the predominantly Early Cretaceous cooling of the Adriatic basement and suggests, at least locally, a fast exhumation.     

Long-Term Changes in Water Quality and Productivity in the Patuxent River Estuary: 1985 to 2003

We conducted a quantitative assessment of estuarine ecosystem responses to reduced phosphorus and nitrogen loading from sewage treatment facilities and to variability in freshwater flow and nonpoint nutrient inputs to the Patuxent River estuary. We analyzed a 19-year dataset of water quality conditions, nutrient loading, and climatic forcing for three estuarine regions and also computed monthly rates of net production of dissolved O₂ and physical transport of dissolved inorganic nitrogen (DIN) and phosphorus (DIP) using a salt- and water-balance model. Point-source loading of DIN and DIP to the estuary declined by 40–60% following upgrades to sewage treatment plants and correlated with parallel decreases in DIN and DIP concentrations throughout the Patuxent. Reduced point-source nutrient loading and concentration resulted in declines in phytoplankton chlorophyll-a (chl-a) and light-saturated carbon fixation, as well as in bottom-layer O₂ consumption for upper regions of the estuary. Despite significant reductions in seaward N transport from the middle to lower estuary, chl-a, turbidity, and surface-layer net O₂ production increased in the lower estuary, especially during summer. This degradation of water quality in the lower estuary appears to be linked to a trend of increasing net inputs of DIN into the estuary from Chesapeake Bay and to above-average river flow during the mid-1990s. In addition, increased abundance of Mnemiopsis leidyi significantly reduced copepod abundance during summer from 1990 to 2002, which favored increases in chl-a and allowed a shift in total N partitioning from DIN to particulate organic nitrogen. These analyses illustrate (1) the value of long-term monitoring data, (2) the need for regional scale nutrient management that includes integrated estuarine systems, and (3) the potential water quality impacts of altered coastal food webs.     

Tracing the exhumation of the Eclogite Zone (Tauern Window,Eastern Alps) by <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar dating of white mica in eclogites

New radiometric ages from the Subpenninic nappes (Eclogite Zone and Rote Wand – Modereck Nappe, Tauern Window) show that phengites formed under eclogite-facies metamorphic conditions retain their initial isotopic signature, even when associated lithologies were overprinted by greenschist- to amphibolite-facies metamorphism. Different stages of the eclogite-facies evolution can be dated provided ⁴⁰Ar/³⁹Ar dating is combined with micro-structural analyses. An age of 39 Ma from the Rote Wand – Modereck Nappe is interpreted to be close to the burial age of this unit. Eclogite deformation within the Eclogite Zone started at the pressure peak along distinct shear zones, and prevailed along the exhumation path. An age of ca. 38 Ma is only observed for eclogites not affected by subsequent deformation and is interpreted as maximum age due to the possible influence of homogenously distributed excess argon. During exhumation deformation was localised along distinct mylonitic shear zones. This stage is mainly characterised by the formation of dynamically recrystallized omphacite2 and phengite. Deformation resulted in the resetting of the Ar isotopic system within the recrystallized white mica. Flat argon release spectra showing ages of 32 Ma within mylonites record the timing of cooling along the exhumation path, and the emplacement onto the Venediger Nappe. Ar-release patterns and ³⁶Ar/⁴⁰Ar vs.³⁹Ar/⁴⁰Ar isotope correlation analyses indicate no significant ⁴⁰Ar-loss after initial closure, and only a negligible incorporation of excess argon. From the pressure peak onwards, eclogitic conditions prevailed for almost 8–10 Ma.     

Probabilistic modeling of seismically triggered landslides using Monte Carlo simulations

The 2004 Chuetsu earthquakes of Niigata (Japan) triggered numerous landslides, and the most widespread types of landslides were highly disrupted, relatively shallow slides and soil (debris) flows. This paper presented a method to evaluate slope instability using Newmark displacement on a pixel-by-pixel basis in a given area. The proposed method was able to integrate Newmark displacement modeling and Monte Carlo simulations within geographical information systems. In the modeling, an empirical attenuation relationship was utilized to calculate Arias intensity over this study area, and the variability of geotechnical parameters was taken into account to calculate coseismic landslide displacement. Before deriving the displacement from related inputs, the Monte Carlo simulations ran 1,500 times and generated 1,500 displacement values for each grid cell, and then means and standard deviations of displacement were calculated and probabilistic distributions can be obtained. Finally, given 10 cm as a threshold value of displacement, estimated probabilities of displacement exceeding 10 cm were shown as a map of seismic landslide hazards. The resulting hazard map was classified into four categories from very low to high level.     

Slope inclinometers for landslides

Slope inclinometers/indicators are used to determine the magnitude, rate, direction, depth, and type of landslide movement. This information is usually vitally important for understanding the cause, behavior, and remediation of a landslide. However, many inclinometer measurements fail to achieve these intended aims because of lack of appreciation of the many factors that need to be correctly implemented during installation, monitoring, and data reduction to yield useful data. This paper presents some guidelines for understanding, installing, and interpreting slope inclinometers and presents three case histories that illustrate some of the pitfalls that can develop if these guidelines are not followed.     

印度板块和亚洲大陆在何时何地碰撞

印度板块和亚洲大陆的初始碰撞时间是所有相关的喜马拉雅-西藏造山体系演化模式的主控条件,并严重影响到对众多与青藏高原隆升和东亚大陆挤出相关的地质过程速率的解释,以及对新生代全球气候变化的理解。尽管印度板块和亚洲大陆汇聚的速率在55Ma突然减缓被广泛地认为是初始碰撞的标志,但这次碰撞所造成的主要构造效应直到20多个百万年以后才显现出来。对印度板块和亚洲大陆相对位置的重新估算,表明它们在55Ma时并没有达到可以彼此发生碰撞的距离。基于来自西藏新的野外证据和对已有数据的重新评估,认为初始碰撞发生在始新世—渐新世之交（约34Ma）,并对55Ma时发生的地质事件提出了另一种解释     

Automated Image Registration for Hydrologic Change Detection in the Lake-Rich Arctic

Multitemporal remote sensing provides a unique tool to track lake dynamics at the pan-Arctic scale but requires precise registration of thousands of satellite images. This is a challenging task owing to a dearth of stable features to be used as tie points [(TPs), i.e., control points] in the dynamic landscapes. This letter develops an automated method to precisely register images in the lake-rich Arctic. The core premise of the method is that the centers of lakes are generally stable even if their shorelines are not. The proposed procedures first extract lakes in multitemporal satellite images, derive lake centroids and match them between images, and then use the centroids of stable lakes as TPs for image registration. The results show that this approach can achieve subpixel registration accuracy, outcompeting the conventional manual methods in both efficiency and accuracy. The proposed method is fully automated and represents a feasible way to register images for lake change detection at the pan-Arctic scale.      

Using Sequential Self-Calibration Method to Identify Conductivity Distribution: Conditioning on Tracer Test data

An iterative inverse method, the sequential self-calibration method, is developed for mapping spatial distribution of a hydraulic conductivity field by conditioning on nonreactive tracer breakthrough curves. A streamline-based, semi-analytical simulator is adopted to simulate solute transport in a heterogeneous aquifer. The simulation is used as the forward modeling step. In this study, the hydraulic conductivity is assumed to be a deterministic or random variable. Within the framework of the streamline-based simulator, the efficient semi-analytical method is used to calculate sensitivity coefficients of the solute concentration with respect to the hydraulic conductivity variation. The calculated sensitivities account for spatial correlations between the solute concentration and parameters. The performance of the inverse method is assessed by two synthetic tracer tests conducted in an aquifer with a distinct spatial pattern of heterogeneity. The study results indicate that the developed iterative inverse method is able to identify and reproduce the large-scale heterogeneity pattern of the aquifer given appropriate observation wells in these synthetic cases.