Predictive GIS-Based Model of Rockfall Activity in Mountain Cliffs

Rockfall susceptibility has been analysed in mountain cliffs of the Cantabrian Range, North Spain. The main aim of this analysis has been to build a predictive model of rockfall activity from a low number of environmental and geological variables. The rockfall activity has been quantified in a GIS. The cartographic information used shows the spatial distribution of all the recent talus screes as well as their associated source areas in the rock-slopes. The area relation At/Ar (recent talus scree polygon/source basins) in the rock slopes has been used as the rockfall activity indicator. This relation has been validated in 50 pilot rock-slopes and compared with the relation number of recent rock fragments/source basin, obtained from field work. The environmental factors causing rockfall depend on the rock slope situation, and these are: altitude and sun radiation on the rock cliff. The geological factors considered are: lithology, relative position of the main discontinuities with respect to the topographic surface and two morphologic parameters: the roughness and slope gradient. A logistic regression analysis has been applied to a population of 442 limestone and quartzite rock cliffs. The dependent variable is the rockfall activity indicator, which allows the definition of two classes of rock cliff units: low and high activity. The independent variables are altitude, sun radiation (equinox radiation, summer solstice radiation, winter solstice radiation), slope roughness, slope gradient,anisotropy and lithology. Results suggest that it is possible tobuild a valid cartographic predictive model for rockfall activity in mountain rock cliffs from a limited number of easily obtainable variables. The method is especially applicable in massive rock slopes or in regions with uniform rock mass characteristics.     

Computational study of tetrahedral Al–Si and octahedral Al–Mg ordering in phengite

 As part of a wider study of the nature and origins of cation order–disorder in micas, a variety of computational techniques have been used to investigate the nature of tetrahedral and octahedral ordering in phengite, K₂ ^[6](Al₃Mg)^[4](Si₇Al)O₂₀(OH)₄. Values of the atomic exchange interaction parameters J _n used to model the energies of order–disorder were calculated. Both tetrahedral Al–Si and octahedral Al–Mg ordering were studied and hence three types of interaction parameter were necessary: for T–T, O–O and T–O interactions (where T denotes tetrahedral sites and O denotes octahedral sites). Values for the T–T and O–O interactions were taken from results on other systems, whilst we calculated new values for the T–O interactions. We have demonstrated that modelling the octahedral and tetrahedral sheets alone and independently produces different results from modelling a whole T–O–T layer, hence justifying the inclusion of the T–O interactions. Simulations of a whole T–O–T layer of phengite indicated the presence of short-range order, but no long-range order was observed. Received: 8 August 2002 / Accepted: 14 February 2003 Acknowledgements The authors are grateful to EPSRC (EJP) and the Royal Society (CIS) for financial support. Monte Carlo simulations were performed on the Mineral Physics Group's Beowulf cluster and the University of Cambridge's High Performance Computing Facility.     

Validation of Landslide Susceptibility Maps; Examples and Applications from a Case Study in Northern Spain

A procedure for validating landslide susceptibility maps wasapplied in a study area in northern Spain and the results obtained compared. Validationwas used to carry out sensitivity analysis for individual variables and combinationsof variables. The validity of different map-making methods was tested, as well as theutility of different types of Favourability Functions. The results obtained show thatvalidation is essential to determine the predictive value of susceptibility maps. Italso helps to better select the most suitable function and significant variables, thus improving the efficiency of the mapping process. Validation based on a temporal strategy makes it possible to derive hazard maps from susceptibility maps.     

The Influence of Geomorphologic Heritage on Present Nival Erosion: Peñalara, Spain

This article examines how snow plays a role in current erosive processes in a high mountain area (1800—2400 m a.s.l.) known as Peñalara, located in Spain's Central Range (40°50' N; 3°58' W). The hypothesis maintains that snow becomes an important erosive factor when it accumulates over sedimentary or weathered materials, therefore geomorphological heritage is a key factor in nival erosion. To test this hypothesis, the authors identified the landforms in the study area and determined their relative ages by weathering and lichenometry ( Rizocarpon geographicum ag. ), differentiating between preglacial, glacial (Recent Pleistocene) and postglacial (Holocene) forms. The information was used to plot a reticulate pattern of observation sites for the study area. Snow depth and the movement of selected blocks at each site were recorded from October 1991 to June 1995. The relationship between late-lying snowpatches, geomorphological heritage and current erosive processes was determined. Between 1800 and 2000 m a.s.l., there is an indirect relationship between snowpatches and predominant processes (stream incision and gelifluction) on terminal moraines. Between 2000 and 2200 m, direct action is present where there are late-laying snowpatches on lateral moraines and some glacial steps. Between 2200 and 2400 m, gelifraction and gravity processes are also in direct relation to snowpatches.     

Influence of the Load Wavelength on the Permeability of a Viscosity Interface in the Mantle

Assuming a radially stratified Newtonian mantle in a steady-state approximation, we demonstrate that the permeability of a viscosity interface at 660-km depth strongly depends on the wavelength of buoyancy forces driving the flow. The flow induced by long-wavelength loads penetrates through the boundary freely even if the viscosity increases by two orders. In contrast, the boundary is practically impermeable for short-wavelength loads located in the upper mantle. Thus, a stepwise increase of viscosity is a significant obstacle for small descending features in the upper mantle, but huge upper mantle downwellings, or upwellings formed in the-lower mantle can overcome it easily. This indicates that certain care is necessary in interpreting the seismic structure of the mantle by means of flow models. The global tomographic image includes only the first few degrees of the harmonic series and, consequently, its interpretation in terms of a present-day flow field results in a predominantly whole-mantle circulation even for extreme viscosity contrasts.     

Converting Radial Distortions Due to Geopotential Models to Geopotential Stokes Coefficients

A methodology for improving geopotential models has been developed. Theoretical relations have been derived converting coefficients in harmonic expansions for radial distortions due to geopotential models into geopotential Stokes coefficients. Terms of the order of 10 ^–10 in magnitude have been retained.     

Greenhouse Gas Emissions from Hydroelectric Reservoirs in Tropical Regions

This paper discusses emissions by power-dams in the tropics. Greenhouse gas emissions from tropical power-dams are produced underwater through biomass decomposition by bacteria. The gases produced in these dams are mainly nitrogen, carbon dioxide and methane. A methodology was established for measuring greenhouse gases emitted by various power-dams in Brazil. Experimental measurements of gas emissions by dams were made to determine accurately their emissions of methane (CH₄) and carbon dioxide (CO₂) gases through bubbles formed on the lake bottom by decomposing organic matter, as well as rising up the lake gradient by molecular diffusion.The main source of gas in power-dams reservoirs is the bacterial decomposition (aerobic and anaerobic) of autochthonous and allochthonous organic matter that basically produces CO₂ and CH₄. The types and modes of gas production and release in the tropics are reviewed.     

Remote sensing detection of atmospheric pollutants by differential absorption LIDAR 510M/SODAR PA2 mobile system

Summary Air pollution problem in the Czech Republic is very complex due to large number of factors as turbulent dispersion and chemical reactions concurring in the status of the PBL. For this 3D distribution of pollutants is a key information that cannot be given by conventional monitoring stations. Combined use of DIAL lidar and sodar can give access to such information. Case studies of air pollutants transport are presented using 3D concentration fields of NO₂ and ozone measured by lidar and meteorological conditions monitored by Doppler Sodar.     

GRASS GIS-embedded Decision Support Framework for Flood Simulation and Forecasting

This study presents a spatiotemporal analysis tool, called Shyska. This tool allows the simulation and prediction of flash floods in semiarid basins. Shyska has been developed by Geographical Information System (GIS)‐embedded functions, allowing the integration of hydrometeorological information from modern technologies of data acquisition in real time. A Digital Elevation Model (DEM) is used in order to obtain the relevant parameters from the integrated rainfall‐runoff models. Some of its most relevant modules and methodology employed for its development are described. Case studies in basins of south‐east Spain illustrate the applicability of the proposed techniques.     

Seismic studies of the Brasília fold belt at the western border of the São Francisco Craton, Central Brazil, using receiver function, surface-wave dispersion and teleseismic tomography

The Tocantins Province in Central Brazil is composed of a series of SSW–NNE trending terranes of mainly Proterozoic ages, which stabilized in the Neoproterozoic in the final collision between the Amazon and São Francisco cratons. No previous information on crustal seismic properties was available for this region. Several broadband stations were used to study the regional patterns of crustal and upper mantle structure, extending the results of a recent E–W seismic refraction profile. Receiver functions and surface wave dispersion showed a thin crust (33–37 km) in the Neoproterozoic Magmatic Arc terrane. High average crustal Vp/Vs ratios (1.74–1.76) were consistently observed in this unit. The foreland domain of the Brasília foldbelt, on the other hand, is characterized by thicker crust (42–43 km). Low Vp/Vs ratios (1.70–1.72) were observed in the low-grade foreland fold and thrust zone of the Brasília belt adjacent to the São Francisco craton. Teleseismic P-wave tomography shows that the lithospheric upper mantle has lower velocities beneath the Magmatic Arc and Goiás Massif compared with the foreland zone of the belt and São Francisco craton. The variations in crustal thickness and upper mantle velocities observed with the broadband stations correlate well with the measurements along the seismic refraction profile. The integration of all seismic observations and gravity data indicates a strong lithospheric contrast between the Goiás Massif and the foreland domain of the Brasília belt, whereas little variation was found across the foldbelt/craton surface boundary. These results support the hypothesis that the Brasília foreland domain and the São Francisco craton were part of a larger São Francisco-Congo continental plate in the final collision with the Amazon plate.