Application of a multi‐temporal,LiDAR‐derived,digital terrain model in a landslide‐volume estimation

Sediments produced by landslides are crucial in the sediment yield of a catchment, debris flow forecasting, and related hazard assessment. On a regional scale, however, it is difficult and time consuming to measure the volumes of such sediment. This paper uses a LiDAR‐derived digital terrain model (DTM) taken in 2005 and 2010 (at 2 m resolution) to accurately obtain landslide‐induced sediment volumes that resulted from a single catastrophic typhoon event in a heavily forested mountainous area of Taiwan. The landslides induced by Typhoon Morakot are mapped by comparison of 25 cm resolution aerial photographs taken before and after the typhoon in an 83.6 km² study area. Each landslide volume is calculated by subtraction of the 2005 DTM from the 2010 DTM, and the scaling relationship between landslide area and its volume are further regressed. The relationship between volume and area are also determined for all the disturbed areas (V_L = 0.452A_L^1.242) and for the crown areas of the landslides (V_L = 2.510A_L^1.206). The uncertainty in estimated volume caused by use of the LiDAR DTMs is discussed, and the error in absolute volume estimation for landslides with an area >10⁵ m² is within 20%. The volume–area relationship obtained in this study is also validated in 11 small to medium‐sized catchments located outside the study area, and there is good agreement between the calculation from DTMs and the regression formula. By comparison of debris volumes estimated in this study with previous work, it is found that a wider volume variation exists that is directly proportional to the landslide area, especially under a higher scaling exponent. Copyright © 2013 John Wiley & Sons, Ltd.     

Southern high latitude climate variability in the Late Cretaceous greenhouse world

A palynological study of oil exploration wells in the Gippsland Basin southeastern Australia has provided a record of southern high latitude climate variability for the last 12 million years of the Cretaceous greenhouse world. During this time, the vegetation was dominated by a cool to temperate flora of Podocarpaceae, Proteaceae and Nothofagidites spp. at a latitude of 60°S. Milankovitch forced cyclic alternations from drier to wetter climatic periods caused vegetation variability from 72 to 77 Ma. This climate change was probably related to the waxing and waning of ephemeral (100 ky) small ice sheets in Antarctica during times of insolation minima and maxima. Drying and cooling after 72 Ma culminated from 68 to 66 Ma, mirroring trends in global δ¹⁸O data. Quantitative palynofloral analyses have the potential to provide realistic proxies for small-scale climate variability in the predominantly ice-free Late Cretaceous.     

Mass density of the upper atmosphere derived from Starlette’s Precise Orbit Determination with Satellite Laser Ranging

The atmospheric mass density of the upper atmosphere from the spherical Starlette satellite’s Precise Orbit Determination is first derived with Satellite Laser Ranging measurements at 815 to 1115 km during strong solar and geomagnetic activities. Starlette’s orbit is determined using the improved orbit determination techniques combining optimum parameters with a precise empirical drag application to a gravity field. MSIS-86 and NRLMSISE-00 atmospheric density models are compared with the Starlette drag-derived atmospheric density of the upper atmosphere. It is found that the variation in the Starlette’s drag coefficient above 800 km corresponds well with the level of geomagnetic activity. This represents that the satellite orbit is mainly perturbed by the Joule heating from geomagnetic activity at the upper atmosphere. This result concludes that MSIS empirical models strongly underestimate the mass density of the upper atmosphere as compared to the Starlette drag-derived atmospheric density during the geomagnetic storms. We suggest that the atmospheric density models should be analyzed with higher altitude acceleration data for a better understanding of long-term solar and geomagnetic effects.     

Chemical and textural equilibration of garnet during amphibolite facies metamorphism: The influence of coupled dissolution–reprecipitation

Metamorphic equilibration requires chemical communication between minerals and may be inhibited through sluggish volume diffusion and or slow rates of dissolution in a fluid phase. Relatively slow diffusion and the perceived robust nature of chemical growth zoning may preclude garnet porphyroblasts from readily participating in low‐temperature amphibolite facies metamorphic reactions. Garnet is widely assumed to be a reactant in staurolite‐isograd reactions, and the evidence for this has been assessed in the Late Proterozoic Dalradian pelitic schists of the Scottish Highlands. The 3D imaging of garnet porphyroblasts in staurolite‐bearing schists reveals a good crystal shape and little evidence of marginal dissolution; however, there is also lack of evidence for the involvement of either chlorite or chloritoid in the reaction. Staurolite forms directly adjacent to the garnet, and its nucleation is strongly associated with deformation of the muscovite‐rich fabrics around the porphyroblasts. “Cloudy” fluid inclusion‐rich garnet forms in both marginal and internal parts of the garnet porphyroblast and is linked both to the production of staurolite and to the introduction of abundant quartz inclusions within the garnet. Such cloudy garnet typically has a Mg‐rich, Mn‐poor composition and is interpreted to have formed during a coupled dissolution–reprecipitation process, triggered by a local influx of fluid. All garnet in the muscovite‐bearing schists present in this area is potentially reactive, irrespective of the garnet composition, but very few of the schists contain staurolite. The staurolite‐producing reaction appears to be substantially overstepped during the relatively high‐pressure Barrovian regional metamorphism reflecting the limited permeability of the schists in peak metamorphic conditions. Fluid influx and hence reaction progress appear to be strongly controlled by subtle differences in deformation history. The remaining garnet fails to achieve chemical equilibrium during the reaction creating distinctive patchy compositional zoning. Such zoning in metamorphic garnet created during coupled dissolution–reprecipitation reactions may be difficult to recognize in higher grade pelites due to subsequent diffusive re‐equilibration. Fundamental assumptions about metamorphic processes are questioned by the lack of chemical equilibrium during this reaction and the restricted permeability of the regional metamorphic pelitic schists. In addition, the partial loss of prograde chemical and textural information from the garnet porphyroblasts cautions against their routine use as a reliable monitor of metamorphic history. However, the partial re‐equilibration of the porphyroblasts during coupled dissolution–reprecipitation opens possibilities of mapping reaction progress in garnet as a means of assessing fluid access during peak metamorphic conditions.     

Mid-Miocene (post 12 Ma) displacement along the central Karakoram fault zone in the Nubra Valley,Ladakh, India from spot LA-ICPMS U/Pb zircon ages of granites

The Karakoram fault zone is a prominent right lateral fault that connects the frontal thrust of the North Pamir with the Indus suture zone near Mount Kailas. Its nature and age of initiation is controversial. In the Nubra valley, Ladakh, India, a Karakoram range granite is thrust over Cretaceous magmatic arc rocks and this thrust is cut by a western strand of the Karakoram fault zone. Three different lithologies from this granite gave weighted mean zircon U/Pb ages of 12.92±0.77 Ma, 12.41±0.43 Ma, and 11.72±0.31 Ma. The ages indicate a relatively short intrusive history of about 1 Ma for the phases: the geochemistry is practically identical to the Pangong leucogranites in the same tectonic block. The Karakoram fault zone in this area is thus less than ~12 Ma old which supports a post middle Miocene (Serravallian) age of Karakoram fault initiation in this area.     

Experimental and numerical study of soil-reinforcement effects on the low-strain stiffness and bearing capacity of shallow foundations

This paper presents results of meticulous laboratory testing and numerical simulations on the effect of reinforcement on the low-strain stiffness and bearing capacity of shallow foundations on dry sand. The effect of the location and the number of reinforcement layers is studied in the laboratory, whereas numerical simulations are used to study the reinforcement-foundation interaction. Laboratory tests show an increase of 100, 200, and 275% not only in bearing capacity but also in low-strain stiffness (linear load–displacement behaviour) of a square foundation when one, two, and three layers of reinforcement are used, respectively. The specimen preparation technique is found to be crucial for the repeatability and reliability of the laboratory results (less than 5% variability). Numerical simulations demonstrate that if reinforcements are placed up to a depth of one footing width (B) below the foundation, better re-distribution of the load to deeper layers is achieved, thus reducing the stresses and strains underneath the foundation. Numerical simulations and experimental results clearly identify a critical zone between 0.3 and 0.5B, where maximum benefits not only on the bearing capacity but also on the low-strain stiffness of the foundation are obtained. Therefore, soil reinforcement can also be used to reduce low-strain vibrations of foundations.     

Seasonal changes in imposex and tissue burden of butyltin compounds in Thais clavigera populations along the coastal area of Mirs Bay, China

We assessed the current status of tributyltin (TBT) contamination of Thais clavigera (Gastropoda) along the coastal area of Mirs Bay, China for the first time. The snail samples were collected from 10 different sites in the summer (June) and winter (December) of 2006, respectively. They were analyzed for imposex status, i.e. relative penis size index (RPSI) and vas deferens sequence index (VDSI), followed by quantification of butyltins in their tissues. Most of the collected females suffered from imposex, but the level of butyltin contamination varied with the distance from Yantian Port, which is currently the fourth busiest container port in the world. The tissue concentration of TBT varied with season. For a particular site, the winter samples in general contained much higher concentration of TBT than the summer samples. RPSI, VDSI and organotin concentrations were higher in T. clavigera collected from sites closer to the Port, such as the Seafood Street and Kat O. Both RPSI and VDSI were positively correlated with the tissue burden of tributyltin. VDSI exhibited little seasonal variability, whereas RPSI showed marked seasonal variability, with lower values in the summer samples. The current results will serve as an important reference for long-term monitoring of butyltin contamination in this area.     

210Pb in the Pacific: the GEOSECS measurements of particulate and dissolved concentrations

We report here on particulate and dissolved²¹⁰Pb profiles at 16 stations, and on total²¹⁰Pb profiles at 3 stations, all occupied during the Pacific GEOSECS expedition. Comparison with measurements at Yale on GEOSECS library samples indicates that during separation of particulate lead from dissolved lead, our filtered water samples suffered some loss of²¹⁰Pb in the filtration system; this effect appears to have reduced the dissolved²¹⁰Pb activities by ～ 20% in stations where the water was filtered. However, for these first Pacific data on the²¹⁰Pb distribution between the two phases, this effect does not significantly interfere with our recognition of the major features of both particulate and dissolved²¹⁰Pb distributions.The dissolved²¹⁰Pb profiles in general vary geographically, following the²²⁶Ra profiles. In deep water,²²⁶Ra increases northward and eastward from the southwest Pacific, from ～ 22dpm/100kg, to over 40 dpm/100 kg in the northeast Pacific. Our dissolved²¹⁰Pb profiles show a similar increase in deep water, varying from about 10 to 20 dpm/100 kg along this line, and are commonly characterized by a mid-depth maximum. This²¹⁰Pb maximum reflects the mid-depth²²⁶Ra maximum of the Pacific Deep Water observed along the western boundary current.In surface water at low latitudes there is a significant²¹⁰Pb flux from the atmosphere, which produces a²¹⁰Pb/²²⁶Ra activity ratio generally greater than unity. This flux penetrates as deep as 600 m, as indicated by an “induced”²¹⁰Pb minimum caused by the surface maximum. The surface water²¹⁰Pb excess decreases toward high southern latitudes and vanishes in the Circumpolar region.The particulate²¹⁰Pb profiles show a general increase with depth, from ～ 0.3dpm/100kg in subsurface water to ～ 1.5dpm/100kg in bottom water, with or without a mid-depth maximum that reflects the²²⁶Ra or dissolved²¹⁰Pb maximum. The particulate²¹⁰Pb normally comprises about 2% of the total²¹⁰Pb in subsurface water, and this fraction increases to about 10% near the bottom. As the filtration loss is not taken into account, the fraction of particulate²¹⁰Pb quoted here is an upper limit. Since the particulate matter concentrations are quite uniform in the water column below a few hundred meters, the²¹⁰Pb activity of the particulate matter also increases with depth. The particulate matter has a²¹⁰Pb concentration of ～ 100dpm/g in subsurface water, but the concentration increases to ～ 500dpm/g or more toward the bottom. This indicates that there is a cumulative adsorption of Pb onto the suspended particles as they are sinking through the water column.     

Application of bivariate frequency analysis to the derivation of rainfall–frequency curves

Bivariate distributions have been recently employed in hydrologic frequency analysis to analyze the joint probabilistic characteristics of multivariate storm events. This study aims to derive practical solutions of application for the bivariate distribution to estimate design rainfalls corresponding to the desired return periods. Using the Gumbel mixed model, this study constructed rainfall–frequency curves at sample stations in Korea which provide joint relationships between amount, duration, and frequency of storm events. Based on comparisons and analyses of the rainfall–frequency curves derived from univariate and bivariate storm frequency analyses, this study found that conditional frequency analysis provides more appropriate estimates of design rainfalls as it more accurately represents the natural relationship between storm properties than the conventional univariate storm frequency analysis.