The El Salvador earthquakes of January and February 2001: context, characteristics and implications for seismic risk

The small Central American republic of El Salvador has experienced, on average, one destructive earthquake per decade during the last hundred years. The latest events occurred on 13 January and 13 February 2001, with magnitudes M_w 7.7 and 6.6, respectively. The two events, which were of different tectonic origin, follow the patterns of the seismicity of the region although neither event has a known precedent in the earthquake catalogue in terms of size and location. The earthquakes caused damage to thousands of traditionally built houses and triggered hundreds of landslides, which were the main causes of fatalities. The earthquakes have clearly demonstrated trends of increasing seismic risk in El Salvador due to rapid population expansion in areas of high shaking and landslide hazard, exacerbated by deforestation and uncontrolled urbanisation. The institutional mechanisms required for the control of land use and building practice are very weak and present a major obstacle to risk mitigation.     

Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data

This paper presents the extension of the self-calibrating method to the coupled inverse modelling of groundwater flow and mass transport. The method generates equally likely solutions to the inverse problem that display the variability as observed in the field and are not affected by a linearisation of the state equations. Conditioning to the state variables is measured by an objective function including, among others, the mismatch between the simulated and measured concentrations. Conditioning is achieved by minimising the objective function by gradient-based methods. The gradient contains the partial derivatives of the objective function with respect to: log conductivities, log storativities, prescribed heads at boundaries, retardation coefficients and mass sources. The derivatives of the objective function with respect to log conductivity are the most cumbersome and need the most CPU-time to be evaluated. For this reason, to compute this derivative only advective transport is considered. The gradient is calculated by the adjoint-state method. The method is demonstrated in a controlled, synthetic study, in which the worth of concentration data is analysed. It is shown that concentration data are essential to improve transport predictions and also help to improve aquifer characterisation and flow predictions, especially in the upstream part of the aquifer, even in the case that a considerable amount of other experimental data like conductivities and heads are available. Besides, conditioning to concentration data reduces the ensemble variances of estimated transmissivity, hydraulic head and concentration.     

Digital model of corrected accumulated flow for peak discharge data acquisition and drainage system design

The design of a drainage system for a roofing slate quarry was implemented by the enhancement of discharge peak estimation, and the uncertainty inevitably associated with the engineering model was reduced.

The development of a topographical, geological, and vegetation cover database developed from a Geographical Information System (GIS) allowed for the definition of the drainage network for a hydraulic system, along with the calculation of the runoff coefficient. This is applied to the digital model of accumulated flow (DMF) as a weight correction coefficient, using a matrix-based model at 5×5 m resolution. The new digital model of corrected accumulated flow (DMCF) is the result of combining the thematic maps with the map of slope <3%, which was previously created from the slope model. It is demonstrated that this new model allows to apply the “Rational Method” on cartographic units defined by the GIS.

The DMCF is compared with other traditional applications of the Rational Method based on the calculation of the discharge peak considering: (1) the drainage basin as a single watershed or (2) defining an average runoff coefficient in each sub-watershed. Both approaches have bigger discharge peaks than those obtained by the DMCF since the slope, lithology, and vegetation cover have average values, and the runoff coefficient is poorly defined, increasing the uncertainty in the discharge peak.     

Preliminary Rock-Slope-Susceptibility Assessment Using GIS and the SMR Classification

This paper aims to examine the impact of large-scale structuraladjustments (like the Greater Dhaka Flood Protection Project, GDFPP) on local living environment.It focuses the importance of environmental factors in flood hazard mitigation, and examines theenvironmental attitudes of the floodplain residents arising from the large-scale structural adjustments.Based on `perceived natural hazard research perspectives', this paper examines: (i) the reasons for persistentfloodplain occupation, and (ii) the importance of environmental factors in the choice, motivations and decision-makingof floodplain residents.This research used data collected from 300 households situated inthe eastern part of Dhaka. The face-to-face household survey data provided individuals' responses to a structuredquestionnaire on hazards and environment. Survey concerned urban floodplains, and looked fordata on housing, household characteristics, and residents' attitudes. Results of interview surveys wereused to: (i) explore the reasons of floodplain occupation, and (ii) residents' attitudes to tolerable levelof flood risk and willingness to accept environmental change resulting from the proposed structural embankments inthe eastern perimeter of Dhaka City, Bangladesh.Findings revealed that floodplain occupation (by theindividuals' decision-making) was a result of overall reaction to the Government's structuraladjustment policies that resulted from institutional, locational and socio-economic factors. The attitude survey results provided residents' perception to hazards and environment to be dependenton the socio-economic factors – but in a complex manner, many factors are interrelated.In addition to support for structural embankments, the study sample displayed a common concernand widespread environmental awareness. In terms of any `trade-off' between thebenefits (resources) from the embankments and costs (hazards) due to the detrimental impact on environment, the residents of Dhaka, despite some concern forsacrificing embankments for environment, tended to show a generalconsensus for embankments.     

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.     

Paleomagnetism of the Cretaceous Morelos and Mezcala Formations, southern Mexico

A paleomagnetic study of platform-facies carbonate rocks of the mid-Cretaceous Morelos Formation and deep-water carbonate rocks of the overlying Upper Cretaceous Mezcala Formation, sampled at Zopilote canyon, in Guerrero State, southern Mexico, indicates that their characteristic magnetization was acquired contemporaneously with folding of these rocks during the Late Cretaceous Laramide orogeny. The remanence carrier is interpreted to be magnetite, although other mineral phases of high coercivity carry recent secondary overprints. The overall mean is of Dec=323.1° and Inc=36.5° (k=162.7; α₉₅=2.7°; N=18 sites; 64% unfolding). Comparison with the North America reference direction indicates that this area has experienced a small, yet statistically significant, counterclockwise direction of 19.2±4.0°. Similar rotations are documented in other localities from southern Mexico; rotations are linked to mid-Tertiary deformation associated with the left-lateral strike-slip fault system that accommodated motion of the Chortis and Xolapa blocks.     

The anisotropy of magnetic susceptibility in biotite, muscovite and chlorite single crystals

The anisotropy of magnetic susceptibility (AMS) of single crystals of biotite, muscovite and chlorite has been measured in order to provide accurate values of the magnetic anisotropy properties for these common rock-forming minerals. The low-field AMS and the high-field paramagnetic susceptibility are defined. For the high-field values, it is necessary to combine the paramagnetic deviatoric tensor obtained from the high-field torque magnetometer with the paramagnetic bulk susceptibility measured from magnetization curves of the crystals. This leads to the full paramagnetic susceptibility ellipsoid due to the anisotropic distribution of iron cations in the silicate lattice. The ellipsoid of paramagnetic susceptibility, which was obtained for the three phyllosilicates, is highly oblate in shape and the minimum susceptibility direction is subparallel to the crystallographic c-axes. The anisotropy of the susceptibility within the basal plane of the biotite has been evaluated and found to be isotropic within the accuracy of the instrumental measurements. The degree of anisotropy of biotite and chlorite is compatible with previously reported values while for muscovite the smaller than previously published values. The shape of the chlorite AMS ellipsoid for all the samples is near-perfect oblate in contrast with a wide distribution of oblate and prolate values reported in earlier studies. Reliable values are important for deriving models of the magnetic anisotropy where it reflects mineral fabrics and deformation of rocks.     

Simulation of a long-lived meso-β scale convective system over the mediterranean coast of Spain. Part II: Sensitivity to external forcings

Summary The numerical simulation of a long-lived, stationary mesoscale convective system (MCS) already described in a previous paper (Fernández et al., 1995) is analyzed in greater detail. The influence of various external forcings, such as sea surface temperature, local orography or terrain roughness, upon the characteristics of the system is studied. This analysis makes it possible not only to identify the most important factors, but also to deduce the importance of some other internal forcings and to propose explanations for some dynamic features of the system that were difficult to understand. Hence, the sensitivity test methodology applied seems to be a useful tool to clarify the complex dynamics of some moist convective events. In the modelled MCS, sea surface temperature and orography are identified as key factors. The results also indicate that the upstream triggering of convection provoked by an orographic blocking effect is the main cause of the development of the system, while upslope triggering plays a secondary role.With 18 Figures     

Assessing land use/cover changes: a nationwide multidate spatial database for Mexico

A nationwide multidate GIS database was generated in order to carry out the quantification and spatial characterization of land use/cover changes (LUCC) in Mexico. Existing cartography on land use/cover at a 1:250,000 scale was revised to select compatible inputs regarding the scale, the classification scheme and the mapping method. Digital maps from three different dates (the late 1970s, 1993 and 2000) were revised, evaluated, corrected and integrated into a GIS database. In order to improve the reliability of the database, an attempt was made to assess the accuracy of the digitalisation procedure and to detect and correct unlikely changes due to thematic errors in the maps. Digital maps were overlaid in order to generate LUCC maps, transition matrices and to calculate rates of conversion. Based upon this database, rates of deforestation between 1976 and 2000 were evaluated as 0.25 and 0.76% per year for temperate and tropical forests, respectively.     

A constrained 2D gravity model of the Sebastián Vizcaíno Basin, Baja California Sur, Mexico

The subsurface geometry of the Sebastián Vizcaíno Basin is obtained from the 2D inversion of gravity data, constrained by a density-versus-depth relationship derived from an oil exploration deep hole. The basin accumulated a thick pile of marine sediments that evolved in the fore-arc region of the compressive margin prevalent along western North America during Mesozoic and Tertiary times. Our interpretation indicates that the sedimentary infill in the Sebastián Vizcaíno Basin reaches a maximum thickness of about 4 km at the centre of a relatively symmetric basin. At the location of the Suaro-1 hole, the depth to the basement derived from this work agrees with the drilled interface between calcareous and volcaniclastic members of the Alisitos Formation. A sensitivity analysis strongly suggests that the assumed density function leads to a nearly unique solution of the inverse problem.