Quantified vegetation change over 42 years at Cape Hallett, East Antarctica

This paper reports on the remapping of a carefully documented vegetation plot at Cape Hallett (72°19′S 170°16′E) to provide an assessment of the rates of vegetation change over decadal time scales. E.D. Rudolph, in 1962, mapped in detail the vegetation of a site approximately 28 m by 120 m at Cape Hallett, Victoria Land, Antarctica. This site was relocated and remapped in January 2004 and changes were assessed using GIS techniques. This appears to be the longest available time period for assessing vegetation change in Antarctica. The analysis indicated that considerable change had occurred in moss and algae distribution patterns and this seems to have been caused by increased water supply, particularly in wetter areas. There was also evidence of some change in lichen distribution. The extent of the change indicates that vegetation cover can be used for monitoring change in areas as extreme as the Ross Sea region. For this analysis to be successful it was important that the mapping techniques used were totally explicit and could easily be replicated. Fortunately, Rudolph had defined his cover classes and the site was also clearly marked. The application of GIS mapping techniques allows the mapping to be more explicitly defined and easily replicated

Source-tracking cadmium in New Zealand agricultural soils: a stable isotope approach

Cadmium (Cd) is a toxic heavy metal, which is accumulated by plants and animals and therefore enters the human food chain. In New Zealand (NZ), where Cd mainly originates from the application of phosphate fertilisers, stable isotopes can be used to trace the fate of Cd in soils and potentially the wider environment due to the limited number of sources in this setting. Prior to 1997, extraneous Cd added to soils in P fertilisers was essentially limited to a single source, the small pacific island of Nauru. Analysis of Cd isotope ratios (ɛ114/110Cd) in Nauru rock phosphate, pre-1997 superphosphate fertilisers, and Canterbury (Lismore Stony Silt Loam) topsoils (Winchmore Research Farm) has demonstrated their close similarity with respect to ɛ114/110Cd. We report a consistent ɛ114/110Cd signature in fertiliser-derived Cd throughout the latter twentieth century. This finding is useful because it allows the application of mixing models to determine the proportions of fertiliser-derived Cd in the wider environment. We believe this approach has good potential because we also found the ɛ114/110Cd in fertilisers to be distinct from unfertilised Canterbury subsoils. In our analysis of the Winchmore topsoil series (1949-2015), the ɛ114/110Cd remained quite constant following the change from Nauru to other rock phosphate sources in 1997, despite a corresponding shift in fertiliser ɛ114/110Cd at this time. We can conclude that to the present day, the Cd in topsoil at Winchmore still mainly originates from historical phosphate fertilisers. One implication of this finding is that the current applications of P fertiliser are not resulting in further Cd accumulation. We aim to continue our research into Cd fate, mobility and transformations in the NZ environment by applying Cd isotopes in soils and aquatic environments across the country

Low-temperature thermochronology and thermokinematic modeling of deformation, exhumation, and development of topography in the central Southern Alps, New Zealand

Apatite and zircon (U-Th)/He and fission track ages were obtained from ridge transects across the central Southern Alps, New Zealand. Interpretation of local profiles is difficult because relationships between ages and topography or local faults are complex and the data contain large uncertainties, with poor reproducibility between sample duplicates. Data do form regional patterns, however, consistent with theoretical systematics and corroborating previous observations: young Neogene ages occur immediately southeast of the Alpine Fault (the main plate boundary structure on which rocks are exhumed); partially reset ages occur in the central Southern Alps; and older Mesozoic ages occur further toward the southeast. Zircon apparent ages are older than apatite apparent ages for the equivalent method. Three-dimensional thermokinematic modeling of plate convergence incorporates advection of the upper Pacific plate along a low-angle detachment then up an Alpine Fault ramp, adopting a generally accepted tectonic scenario for the Southern Alps. The modeling incorporates heat flow, evolving topography, and the detailed kinetics of different thermochronometric systems and explains both complex local variations and regional patterns. Inclusion of the effects of radiation damage on He diffusion in detrital apatite is shown to have dramatic effects on results. Geometric and velocity parameters are tuned to fit model ages to observed data. Best fit is achieved at 9 mm a−1 plate convergence, with Pacific plate delamination on a gentle 10°SE dipping detachment and more rapid uplift on a 45–60° dipping Alpine Fault ramp from 15 km depth. Thermokinematic modeling suggests dip-slip motion on reverse faults within the Southern Alps should be highest ∼22 km from the Alpine Fault and much lower toward the southeast

Investigation of the relationships between basin morphology, tectonic uplift, and denudation from the study of an active fold belt in the Siwalik Hills, central Nepal

The present study investigates correlations between an extensive range of geomorphic properties that can be estimated from a digital elevation model and the uplift rate on geological timescales. The analysis focuses on an area in the Siwalik Hills (central Nepal), where lithology and climate can be considered as uniform. This area undergoes rapid tectonic uplift at rates of up to 15 mm yr^(−1), which are derived from the geometric pattern of a fault-bend model of fold growth. The selected geomorphic properties can be divided in two categories, depending on whether or not the vertical dimension is taken into account. None of the planar properties are significantly correlated to uplift rate, unlike those that include the vertical dimension, such as the mean elevation of basins, hypsometric curve, and hypsometric integral, and relief defined by the amplitude factor of length scaling analysis. Correlation between relief and uplift rate is observed for all length scales of topography shorter than 600 m, which suggests that all orders of the streams are able to adjust to the tectonic signal. Simple mass balance considerations imply that the average elevation is only 10% of surface uplift, suggesting that a dynamic equilibrium has been reached quite rapidly. Using a simple two-process model for erosion, we find that fairly high diffusion coefficients (order of 10 m^2 yr^(−1)) and efficient transport of the material by rivers are required. This unusually high value for mass diffusivity at small length scales may be obtained by either a very efficient linear diffusion or by landsliding. Actually, both processes may be active, which appears likely given the nature of the unconsolidated substratum and the favorable climatic conditions. Local relief in the study area may therefore be used to predict either uplift or denudation, but the prediction is calibrated only for that specific climatic and lithologic conditions and cannot be systematically applied to other contexts

Soil erosion in the Alps : causes and risk assessment

The issue of soil erosion in the Alps has long been neglected due to the low economic value of the agricultural land. However, soil stability is a key parameter which affects ecosystem services like slope stability, water budgets (drinking water reservoirs as well as flood prevention), vegetation productivity, ecosystem biodiversity and nutrient production. In alpine regions, spatial estimates on soil erosion are difficult to derive because the highly heterogeneous biogeophysical structure impedes measurement of soil erosion and the applicability of soil erosion models. However, remote sensing and geographic information system (GIS) methods allow for spatial estimation of soil erosion by direct detection of erosion features and supply of input data for soil erosion models. Thus, the main objective of this work is to address the problem of soil erosion risk assessment in the Alps on catchment scale with remote sensing and GIS tools. Regarding soil erosion processes the focus is on soil erosion by water (here sheet erosion) and gravity (here landslides). For these two processes we address i) the monitoring and mapping of the erosion features and related causal factors ii) soil erosion risk assessment with special emphasis on iii) the validation of existing models for alpine areas. All investigations were accomplished in the Urseren Valley (Central Swiss Alps) where the valley slopes are dramatically affected by sheet erosion and landslides. For landslides, a natural susceptibility of the catchment has been indicated by bivariate and multivariate statistical analysis. Geology, slope and stream density are the most significant static landslide causal factors. Static factors are here defined as factors that do not change their attributes during the considered time span of the study (45 years), e.g. geology, stream network. The occurrence of landslides might be significantly increased by the combined effects of global climate and land use change. Thus, our hypothesis is that more recent changes in land use and climate affected the spatial and temporal occurrence of landslides. The increase of the landslide area of 92% within 45 years in the study site confirmed our hypothesis. In order to identify the cause for the trend in landslide occurrence time-series of landslide causal factors were analysed. The analysis revealed increasing trends in the frequency and intensity of extreme rainfall events and stocking of pasture animals. These developments presumably enhanced landslide hazard. Moreover, changes in land-cover and land use were shown to have affected landslide occurrence. For instance, abandoned areas and areas with recently emerging shrub vegetation show very low landslide densities. Detailed spatial analysis of the land use with GIS and interviews with farmers confirmed the strong influence of the land use management practises on slope stability. The definite identification and quantification of the impact of these non-stationary landslide causal factors (dynamic factors) on the landslide trend was not possible due to the simultaneous change of several factors. The consideration of dynamic factors in statistical landslide susceptibility assessments is still unsolved. The latter may lead to erroneous model predictions, especially in times of dramatic environmental change. Thus, we evaluated the effect of dynamic landslide causal factors on the validity of landslide susceptibility maps for spatial and temporal predictions. For this purpose, a logistic regression model based on data of the year 2000 was set up. The resulting landslide susceptibility map was valid for spatial predictions. However, the model failed to predict the landslides that occurred in a subsequent event. In order to handle this weakness of statistic landslide modelling a multitemporal approach was developed. It is based on establishing logistic regression models for two points in time (here 1959 and 2000). Both models could correctly classify >70% of the independent spatial validation dataset. By subtracting the 1959 susceptibility map from the 2000 susceptibility map a deviation susceptibility map was obtained. Our interpretation was that these susceptibility deviations indicate the effect of dynamic causal factors on the landslide probability. The deviation map explained 85% of new independent landslides occurring after 2000. Thus, we believe it to be a suitable tool to add a time element to a susceptibility map pointing to areas with changing susceptibility due to recently changing environmental conditions or human interactions. In contrast to landslides that are a direct threat to buildings and infrastructure, sheet erosion attracts less attention because it is often an unseen process. Nonetheless, sheet erosion may account for a major proportion of soil loss. Soil loss by sheet erosion is related to high spatial variability, however, in contrast to arable fields for alpine grasslands erosion damages are long lasting and visible over longer time periods. A crucial erosion triggering parameter that can be derived from satellite imagery is fractional vegetation cover (FVC). Measurements of the radiogenic isotope Cs-137, which is a common tracer for soil erosion, confirm the importance of FVC for soil erosion yield in alpine areas. Linear spectral unmixing (LSU), mixture tuned matched filtering (MTMF) and the spectral index NDVI are applied for estimating fractional abundance of vegetation and bare soil. To account for the small scale heterogeneity of the alpine landscape very high resolved multispectral QuickBird imagery is used. The performance of LSU and MTMF for estimating percent vegetation cover is good (r²=0.85, r²=0.71 respectively). A poorer performance is achieved for bare soil (r²=0.28, r²=0.39 respectively) because compared to vegetation, bare soil has a less characteristic spectral signature in the wavelength domain detected by the QuickBird sensor. Apart from monitoring erosion controlling factors, quantification of soil erosion by applying soil erosion risk models is done. The performance of the two established models Universal Soil Loss Equation (USLE) and Pan-European Soil Erosion Risk Assessment (PESERA) for their suitability to model erosion for mountain environments is tested. Cs-137 is used to verify the resulting erosion rates from USLE and PESERA. PESERA yields no correlation to measured Cs-137 long term erosion rates and shows lower sensitivity to FVC. Thus, USLE is used to model the entire study site. The LSU-derived FVC map is used to adapt the C factor of the USLE. Compared to the low erosion rates computed with the former available low resolution dataset (1:25000) the satellite supported USLE map shows “hotspots” of soil erosion of up to 16 t ha-1 a-1. In general, Cs-137 in combination with the USLE is a very suitable method to assess soil erosion for larger areas, as both give estimates on long-term soil erosion. Especially for inaccessible alpine areas, GIS and remote sensing proved to be powerful tools that can be used for repetitive measurements of erosion features and causal factors. In times of global change it is of crucial importance to account for temporal developments. However, the evaluation of the applied soil erosion risk models revealed that the implementation of temporal aspects, such as varying climate, land use and vegetation cover is still insufficient. Thus, the proposed validation strategies (spatial, temporal and via Cs-137) are essential. Further case studies in alpine regions are needed to test the methods elaborated for the Urseren Valley. However, the presented approaches are promising with respect to improve the monitoring and identification of soil erosion risk areas in alpine regions

Impacts of anthropogenic and environmental factors on the occurrence of shallow landslides in an alpine catchment (Urseren Valley, Switzerland)

Changes in climate and land use pose a risk to stability of alpine soils, but the direction and magnitude of the impact is still discussed controversially with respect to the various alpine regions. In this study, we explicitly consider the influence of dynamic human-induced changes on the occurrence of landslides in addition to natural factors. Our hypothesis was that if changes in land use and climate have a significant influence on the occurrence of landslides we would see a trend in the incidence of landslides over time. We chose the Urseren Valley in the Central Swiss Alps as investigation site because the valley is dramatically affected by landslides and the land use history is well documented. Maps of several environmental factors were used to analyse the spatial landslide pattern. In order to explain the causation of the temporal variation, time-series (45 years) of precipitation characteristics, cattle stocking and pasture maps were compared to a series of seven landslide investigation maps between 1959 and 2004. We found that the area affected by landslides increased by 92% from 1959 to 2004. Even though catchment characteristics like geology and slope largely explain the spatial variation in landslide susceptibility (68%), this cannot explain the temporal trend in landslide activity. The increase in stocking numbers and the increased intensity of torrential rain events had most likely an influence on landslide incidence. In addition, our data and interviews with farmers pointed to the importance of management practice

Seismic response to evolving injection at the Rotokawa geothermal field, New Zealand

Catalogs of microseismicity are routinely compiled at geothermal reservoirs and provide valuable insights into reservoir structure and fluid movement. Hypocentral locations are typically used to infer the orientations of structures and constrain the extent of the permeable reservoir. However, frequency-magnitude distributions may contain additional, and underused, information about the distribution of pressure. Here, we present a four-year catalog of seismicity for the Rotokawa geothermal field in the central Taupō Volcanic Zone, New Zealand starting two years after the commissioning of the 140 MWe Nga Awa Purua power station. Using waveform-correlation-based signal detection we double the size of the previous earthquake catalog, refine the location and orientation of two reservoir faults and identify a new structure. We find the rate of seismicity to be insensitive to major changes in injection strategy during the study period, including the injectivity decline and shift of injection away from the dominant injector, RK24. We also map the spatial distribution of the earthquake frequency-magnitude distribution, or b-value, and show that it increases from ∼1.0 to ∼1.5 with increasing depth below the reservoir. As has been proposed at other reservoirs, we infer that these spatial variations reflect the distribution of pressure in the reservoir, where areas of high b-value correspond to areas of high pore-fluid pressure and a broad distribution of activated fractures. This analysis is not routinely conducted by geothermal operators but shows promise for using earthquake b-value as an additional tool for reservoir monitoring and management

Rapid methods of landslide hazard mapping : Fiji case study

A landslide hazard probability map can help planners (1) prepare for, and/or mitigate against, the effects of landsliding on communities and infrastructure, and (2) avoid or minimise the risks associated with new developments. The aims of the project were to establish, by means of studies in a few test areas, a generic method by which remote sensing and data analysis using a geographic information system (GIS) could provide a provisional landslide hazard zonation map. The provision of basic hazard information is an underpinning theme of the UN’s International Decade for Natural Disaster Reduction (IDNDR). It is an essential requirement for disaster preparedness and mitigation planning. This report forms part of BGS project 92/7 (R5554) ‘Rapid assessment of landslip hazards’ Carried out under the ODA/BGS Technology Development and Research Programme as part of the British Government’s provision of aid to developing countries. It provides a detailed technical account of work undertaken in a test area in Viti Levu in collaboration with Fiji Mineral Resources Department. The study represents a demonstration of a methodology that is applicable to many developing countries. The underlying principle is that relationships between past landsliding events, interpreted from remote sensing, and factors such as the geology, relief, soils etc provide the basis for modelling where future landslides are most likely to occur. This is achieved using a GIS by ‘weighting’ each class of each variable (e.g. each lithology ‘class’ of the variable ‘geology’) according to the proportion of landslides occurring within it compared to the regional average. Combinations of variables, produced by summing the weights in individual classes, provide ‘models’ of landslide probability. The approach is empirical but has the advantage of potentially being able to provide regional scale hazard maps over large areas quickly and cheaply; this is unlikely to be achieved using conventional ground-based geotechnical methods. In Fiji, landslides are usually triggered by intense rain storms commonly associated with tropical cyclones. However, the regional distribution of landslides has not been mapped nor is it known how far geology and landscape influence the location and severity of landsliding events. The report discusses the remote sensing and GIS methodology, and describes the results of the pilot study over an area of 713 km2 in south east Viti Levu. The landslide model uses geology, elevation, slope angle, slope aspect, soil type, and forest cover as inputs. The resulting provisional landslide hazard zonation map, divided into high, medium and low zones of landslide hazard probability, suggests that whilst rainfall is the immediate cause, others controls do exert a significant influence. It is recommended that consideration be given in Fiji to implementing the techniques as part of a national strategic plan for landslide hazard zonation mapping

Patterns of Sponge Abundance Across a Gradient of Habitat Quality in the Wakatobi Marine National Park, Indonesia

Sponges are important components of reef communities worldwide, fulfilling a number of important functional roles. Habitat degradation caused by the loss of hard corals has the potential to cause increases in sponge abundance and percentage cover as they gain access to resources such as space and food. In this study we compared sponge densities and percentage cover at sites with varying hard coral cover in the Wakatobi Marine National Park, Indonesia (WMNP). We found significant differences in sponge densities at the study sites but no significant difference in sponge densities on different surface angles. Unexpectedly, we also found a weak positive correlation between coral cover and sponge density. This indicates that spatial competition is unlikely to be the most important factor determining sponge abundance in the WMNP. In contrast to sponge density data, we found that sponge percentage cover and hard coral cover were weakly negatively correlated, but found no significant difference in sponge percentage cover between the study sites. Finally, multivariate analysis of the benthic communities at the study sites indicated that while sites with higher coral cover were characterised by coral (proportionally), lower coral cover sites were characterised by algae and sponges. This suggests that although there was no significant difference in sponge percentage cover between the study sites conditions that led to the loss of hard coral at lower quality sites mean that these sites are characterised by sponges and algae rather than by any other groups of benthic organisms

A tribute to Professor William Dearman : new small-scale engineering geological maps of the United Kingdom

Professor William Robert (“Bill”) Dearman was the first British Professor of Engineering Geology and a world leader in the development of engineering geological mapping techniques and methodologies, recognised by being awarded the IAEG’s Hans Cloos Medal. The maps described here would not have been completed without his initial ideas and interpretative work. It is a great regret that the maps were not finished in time for him to see them before his death in January 2009. This paper and the maps described in it are, therefore, dedicated to his memory