Semi-automatic shallow landslide detection by the integration of airborne imagery and laser scanning data

Landslide mapping is essential for effective watershed management. In Taiwan, a typhoon or earthquake event can trigger hundreds, even thousands, of shallow landslides in mountainous watersheds. Thus, improving the efficiency of landslide mapping by means of remote sensing techniques is an important issue. This study proposes a new method that uses concurrent aerial laser scanning (ALS) data and color ortho-imagery as input data: the topographic indices of slope, surface roughness, and object height model can be derived from the ALS data and the Green–Red Vegetation Index from the ortho-images. The method first uses these topographic and spectral indices in a global, semi-automatic algorithm to separate landslide from non-landslide pixels. It then offers a region growing tool and a 3D Eraser/Painter to edit detected landslides locally. These global and local operations are designed with a user interface, which is intuitive and user-friendly. Results from four test sites in a mountainous watershed prove that the method is easy, accurate, and suitable for landslide mapping in Taiwan.     

Does sedimentary organic delta 13C record variations in quaternary ocean [CO2(aq)]?

Ocean surface water [CO2(aq)] variations based on glacial/interglacial changes in sediment delta 13Corg are shown to compare favorably with reconstructions based on ice core [CO2]. In particular, an approximate 80 microatmospheres increase in atmospheric pCO2 during the last glacial-interglacial transition is calculated to correspond to a 3-4 micromolar increase in ocean surface water [CO2(aq)] at atmospheric equilibrium. A widespread marine delta 13Corg decrease of 1-2% accompanied this event and was not preceded by an equivalent isotopic change in surface water total dissolved inorganic carbon. These observations support the hypothesis that [CO2(aq)] influences photosynthetic isotope fractionation between marine inorganic and organic carbon pools, and therefore that plankton/sediment delta 13Corg may serve as a proxy for surface water [CO2(aq)].     

Long-term spatio-temporal precipitation variability in arid-zone Australia and implications for groundwater recharge

Quantifying dryland groundwater recharge as a function of climate variability is becoming increasingly important in the face of a globally depleted resource, yet this remains a major challenge due to lack of adequate monitoring and the complexity of processes involved. A previously unpublished and unique dataset of high density and frequency rainfall measurements is presented, from the Fowlers Gap Arid Zone Research Station in western New South Wales (Australia). The dataset confirms extreme spatial and temporal variability in rainfall distribution which has been observed in other dryland areas and is generally explained by the dominance of individual storm cells. Contrary to previous observations, however, this dataset contains only a few localised storm cells despite the variability. The implications of spatiotemporal rainfall variability on the estimation of groundwater recharge is assessed and show that the most likely recharge mechanism is through indirect and localised, rather than direct, recharge. Examples of rainfall and stream gauge height illustrate runoff generation when a spatially averaged threshold of 15–25 mm (depending on the antecedent moisture conditions) is exceeded. Preliminary assessment of groundwater levels illustrates that the regional water table is much deeper than anticipated, especially considering the expected magnitude of indirect and localised recharge. A possible explanation is that pathways for indirect and localised recharge are inhibited by the large quantities of Aeolian dust observed at the site. Runoff readily occurs with water collecting in surface lakes which slowly dry and disappear. Assuming direct groundwater recharge under these conditions will significantly overestimate actual recharge.     

Delta13C and delta15N shifts in benthic invertebrates exposed to sewage from McMurdo Station, Antarctica

In an effort to identify biomonitors for contamination of Antarctic marine benthos by sewage, this study determines whether the US Antarctic Program’s McMurdo Station produces a benthic sewage footprint and whether resident megafauna are assimilating sewage-derived material. We identified strong C and N isotopic gradients in benthic sediment as a function of downstream distance from McMurdo Station’s point-source sewage addition. Sediment C and N isotope ratios approached marine background levels at the sampling end-point 612 m downcurrent. Based on isotope abundances in their tissues, at least some sewage C and N were assimilated by the sedentary, suspension feeding soft coral Alcyonium antarcticum, ascidian Cnemidocarpa verrucosa and bivalve Laternula elliptica. However, as inferred by tissue-sediment differences in downstream isotope trends, such assimilation was not in proportion to sewage exposure and input, therefore implying non-generalist feeding behavior by these species. In contrast, the motile, generalist feeding sea urchin Sterechinus neumayeri, sea star Odontaster validus and ribbon worm Parborlasia corrugatus showed isotopic evidence of sewage C and N assimilation roughly in proportion to sewage input. We recommend these generalist feeders for further use as biomonitors at this site now that sewage treatment has been implemented. As these species are circumpolar in distribution, they may also prove useful elsewhere in the Antarctic.     

Analytical methods that use natural heat as a tracer to quantify surface water–groundwater exchange, evaluated using field temperature records

Two methods applying natural heat as a tracer to quantify surface water–groundwater exchange were evaluated using field data. Arrays capable of monitoring and recording the streambed response to diurnal temperature variations in the surface water were deployed for a 2-month period in three locations in perennial pools at Maules Creek, New South Wales, Australia. Multi-level array design, field deployment and parameter estimation are discussed. The applicability of analytical solutions derived from the heat transport equation to the streambed environments was analysed using the recorded temperature time series. The stream was found to lose water to the aquifer, which was supported by simultaneously recorded hydraulic gradients. However, the one-dimensional (1D) analytical solutions did not adequately describe the observed streambed thermal response at two locations. The resulting artefacts in the estimated flow velocities are discussed. It was hypothesised that the artefacts originate from model limitation due to streambed heterogeneity and application of 1D solutions to multi-dimensional and dynamic streambed flow. This consequently imposes limitations on the field applicability of the methods. Nevertheless, in combination with time series of surface water and streambed water levels, the use of heat as a tracer provided a powerful tool for better understanding the shallow hydrogeological system.     

Snow zonation on Hielo Patagónico Sur, Southern Patagonia, derived from Landsat 5 TM data

Hielo Patagónico Sur (HPS), an icefield in Southern Patagonia, is the largest temperate ice mass in the southern hemisphere. Despite continued research efforts during the last decade many glaciological variables, especially mass balance, are still poorly known. This is partly because access to the icefield is difficult due to remoteness and persistent harsh weather conditions. Therefore, remote sensing appears to be a more suitable tool for the acquisition of data. In this work we present a remote sensing study of snow zonation on HPS using Landsat 5 Thematic Mapper data acquired on 12 March 2001. By using image processing and classification techniques, proved to be useful in other glaciated regions, we map for the first time the extent and occurrence of major snow zones on the whole HPS. We separate between two classes of ice and three classes of snow. Ice facies are classified as bare or debris (i.e. dirt) covered ice, covering 2454 km² and 777 km² respectively, or 18.4% and 5.8% of the icefield on the day of image acquisition. Snow types are classified according to spectral differences in the images, following the glacier facies concept. Two of the three snow cover types are interpreted to represent differences in snow grain size within a fairly homogeneous snow pack whereas the third one is interpreted to represent the slush zone. A first order altitudinal control on the distribution of these snow facies is evident. In addition, our results show that snow accumulation on HPS is markedly controlled by the interaction of strong west–northwest snow-bearing winds and the rough mountainous terrain. In order of decreasing altitude we find that the two snow facies and the slush facies occupy 3819 km², 3292 km² and 2295 km² respectively, or 28.6%, 24.6% and 17.2% of the icefield, on the day of image acquisition. Estimates of equilibrium line altitude using our results yield values of 800–900 m above sea level for the western side and 1500–1600 m above sea level for the eastern side, with an accumulation area ratio of 0.74.     

A 450-kyr record of hydrological conditions on the western Agulhas Bank Slope, south of Africa

The Agulhas Bank region, south of Africa, is an oceanographically important and complex area. The leakage of warm saline Indian Ocean water into the South Atlantic around the southern tip of Africa is a crucial factor in the global thermohaline circulation. Foraminiferal assemblage, stable isotope and sedimentological data from the top 10 m of core MD962080, recovered from the western Agulhas Bank Slope, are used to indicate changes in water mass circulation in the southeastern South Atlantic for the last 450 kyr. Sedimentological and planktonic foraminiferal data give clear signals of cold water intrusions. The benthic stable isotope record provides the stratigraphic framework and indicates that the last four climatic cycles are represented (i.e. down to marine isotope stage (MIS) 12). The planktonic foraminiferal assemblages bear a clear transitional to subantarctic character with Globorotalia inflata and Neogloboquadrina pachyderma (dextral) being the dominant taxa. Input of cold, subantarctic waters into the region by means of leakage through the Subtropical Convergence, as part of Agulhas ring shedding, and a general cooling of surface waters is suggested by increased occurrence of the subantarctic assemblage during glacial periods. Variable input of Indian Ocean waters via the Agulhas Current is indicated by the presence of tropical/subtropical planktonic foraminiferal species Globoquadrina dutertrei, Globigerinoides ruber (alba) and Globorotalia menardii with maximum leakage occurring at glacial terminations. The continuous presence of G. menardii throughout the core suggests that the exchange of water from the South Indian Ocean to the South Atlantic Ocean was never entirely obstructed in the last 450 kyr. The benthic carbon isotope record and sediment textural data reflect a change in bottom water masses over the core location from North Atlantic Deep Water to Upper Southern Component Water. Planktonic foraminiferal assemblages and sediment composition indicate a profound change in surface water conditions over the core site approximately 200–250 kyr BP, during MIS 7, from mixed subantarctic and transitional water masses to overall warmer surface water conditions.     

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.