A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

Study of geothermal water intrusion due to groundwater exploitation in the Puebla Valley aquifer system, Mexico

Significant intrusion of geothermal water into fresh groundwater takes place in the Puebla Valley aquifer system, Mexico. The decline in the potentiometric surface due to the overexploitation of the groundwater induces this intrusion. This hydrological system comprises three aquifers located in Plio-Quaternary volcanic sediments and Mesozoic calcareous rocks. The hydraulic balance of the aquifer shows that the annual output exceeds the natural inputs by 12 million m³. Between 1973 and 2002, a drop in the potentiometric surface, with an 80 m cone of depression, was identified in a 5-km-wide area located southwest of the city of Puebla. Chemical analyses performed on water samples since 1990 have shown an increase in total dissolved solids (TDS) of more than 500 mg/L, coinciding with the region showing a cone of depression in the potentiometric surface. A three-dimensional flow and transport model, based on the hydrogeological and geophysical studies, was computed by using the MODFLOW and MT3D software. This model reproduces the evolution of the aquifer system during the last 30 years and predicts for 2010 an additional drawdown in the potentiometric surface of 15 m, and an increase in the geothermal water intrusion.     

Catastrophic debris flows near Machu Picchu village (Aguas Calientes), Peru

Slope movements together with intensive river erosion and the following accumulation are the leading processes in the landscape evolution in the area of Machu Picchu village (former Aguas Calientes), which is located close to the Machu Picchu Sanctuary. Debris flows affect not only the bottoms of valleys or canyons, but also debris fans at the termini of the drainage basins, which are heavily inhabited at some places. The most recent event in the Machu Picchu village occurred in April 2004, but several others were documented in a broader area in the last 50 years. The field inspections at Machu Picchu (May and September 2004; June and September 2005) together with oral testimony revealed the nature and behavior of the debris flow. Machu Picchu village can be assessed as a zone with high landslide risk in relation to its urban development. Despite that, the village recorded a rapid growth (threefold population increase) without urban control within the past two decades. Precipitation, which is the main triggering factor of the debris flows, and natural hazard management of the Machu Picchu village are discussed in this paper.     

The Spačice eclogite: constraints on the P–T–t history of the Gföhl granulite terrane, Moldanubian Zone, Bohemian Massif

Summary In the Kutná Hora Complex, the Běstvina Formation, which is similar to Gf?hl granulite, contains eclogite that has escaped widespread retrograde recrystallization. The eclogite assemblage, garnet + omphacite + quartz + rutile ± plagioclase, yields an estimate for peak metamorphic conditions of 18–20 kbar and 835–935 °C, which is comparable to that determined from felsic granulite, 14–20 kbar and 900–1000 °C. Garnet in eclogite exhibits both prograde and retrograde compositional zoning, from which constraints on thermal history of the Gf?hl terrane can be derived by diffusion modelling. At 900 °C, a garnet grain of 800–1000 μm radius would homogenize in 7.5–11.7 million years, but the existence of compositional gradients on a length scale of 100–200 μm suggests that the duration of peak metamorphism may have been limited to ∼500,000 years. Diffusion modelling of retrograde zoning in garnet yields a cooling rate of 150–100 °C/m.y. for a radius of 800–1000 μm and initial temperature of 900 °C. The relatively brief duration of high-pressure/high-temperature metamorphism and rapid cooling and exhumation of the Gf?hl terrane may be a consequence of lithospheric delamination during Early Carboniferous collision of Bohemia (Teplá-Barrandia) and Moldanubia (Franke, 2000).     

Factors influencing the presence or absence of tributary-junction fans in the Iberian Range, Spain

This paper analyses the factors which influence the presence or absence of tributary-junction fans in the Iberian Range, northern Spain. Two valleys were selected, both characterised by wide variations in lithology, altitude, land use and plant cover. Two groups of factors were studied: those related to the internal characteristics of the drainage basins, which particularly control sediment generation; and those related to the characteristics of the depositional area which control accommodation space and main river power. Among the internal factors, the development of alluvial fans was related to: (i) the capacity of the basin to yield large volumes of sediment, (ii) the occurrence of intense human pressure until recent times, a good indicator of sediment yield, and (iii) the capacity of the basin to quickly increase discharge during rainstorms (discharge density and torrentiality). It is suggested that the areas that were intensively cultivated in the past, and have therefore been affected by intense erosion, have played a decisive role on the development of alluvial fans. This would imply that many of these alluvial fans have a relatively recent origin, perhaps related to the beginning of a widespread deforestation. The basins without alluvial fans are characterised by relatively steep hillslope gradients (that is, slopes that never were subjected to historical cultivation), low drainage densities and dense forest and shrub cover, mostly coinciding with high altitude basins composed of quartzite and shale bedrocks. Regarding the external factors, the shape, size and longitudinal gradient of the main river to which the fans are tributary are the most relevant conditioning factors determining the development of alluvial fans.     

Avalanche Risk During Backcountry Skiing – An Analysis of Risk Factors

Skier-triggered avalanches are the main cause of avalanche accidents in backcountry skiing. The risk of accidents during backcountry skiing was analysed statistically and related to factors such as elevation level, aspect, stability rating and the time of the year. The analysis is based on a database about terrain usage and avalanche accidents from a large heli-skiing operator in Canada, which makes it possible to study the conditional probability of accidents given the recorded pattern of terrain usage. This study shows that the historical risk of accidentally triggering an avalanche greater than size 1 depends highly on the stability rating, with the highest risk occurring during “poor” stability. The risk is greater at high elevations, and it is lower during the late season than earlier on. Skier risk does not depend as much on aspect as may be indicated from avalanche data alone. However, it is relatively high in the N–NE–E sector. These factors are not independent of each other and therefore analyses of combined factors were also performed. Questionnaires and interviews were used to gain knowledge about the terrain selection of professional mountain guides. These results indicate that when selecting terrain, guides first look at the overall shape and size of the terrain, but avalanche history of terrain and inclination are also important factors. Finally, remarks in avalanche reports were analysed, and common human factors identified.     

Effect of sampling and linkage on fault length and length–displacement relationship

The power-law exponent (n) in the equation: D=cL ⁿ , with D = maximum displacement and L = fault length, would be affected by deviations of fault trace length. (1) Assuming n=1, numerical simulations on the effect of sampling and linkage on fault length and length–displacement relationship are done in this paper. The results show that: (a) uniform relative deviations, which means all faults within a dataset have the same relative deviation, do not affect the value of n; (b) deviations of the fault length due to unresolved fault tip decrease the values of n and the deviations of n increase with the increasing length deviations; (c) fault linkage and observed dimensions either increase or decrease the value of n depending on the distribution of deviations within a dataset; (d) mixed deviations of the fault lengths are either negative or positive and cause the values of n to either decrease or increase; (e) a dataset combined from two or more datasets with different values of c and orders of magnitude also cause the values of n to deviate. (2) Data including 19 datasets and spanning more than eight orders of fault length magnitudes (10⁻²–10⁵ m) collected from the published literature indicate that the values of n range from 0.55 to 1.5, the average value being 1.0813, and the peak value of n _d (double regression) is 1.0–1.1. Based on above results from the simulations and published data, we propose that the relationship between the maximum displacement and fault length in a single tectonic environment with uniform mechanical properties is linear, and the value of n deviated from 1 is mainly caused by the sampling and linkage effects.     

Geometry and accuracy of reflecting points in bistatic satellite altimetry

The analysis of the time and space distribution of specular (reflecting) points in bistatic altimetry between GPS and CHAMP satellites or SAC-C (taken as examples) is extended from Wagner and Klokočník (2003 J. Geod 77: 128–138). We demonstrate a significantly higher number and density of reflecting points in bistatic altimetry in comparison with traditional monostatic altimetry. After an outline of our older accuracy assessment for the vertical position of the reflecting point, we add a new independent derivation and compare both approaches. We account for orbit errors of both the transmitters (GPS) and receiver (CHAMP) satellites, and the measurement (delay) error. We found that the accuracy of the vertical position of the reflecting point decreases only slowly with increasing off-nadir angle and that the orbit errors must be accounted for if decimeter and better accuracy is required. In this paper, we do not study errors such as state of the ocean, technical parameters of the receiving system, and atmospheric corrections.     

Daily Meteorological Observations in Cádiz – San Fernando. Analysis of the Documentary Sources and the Instrumental Data Content (1786–1996)

Meteorological observations in the city of Cádiz are acknowledged as having been made from the middle of the 18th century onwards although they were only recorded and preserved in documentary form from 1789 onwards. Data readings were taken at the new Naval Observatory in San Fernando, ten kilometres from Cádiz, from 1797 onwards. Continuous series for temperature and atmospheric pressure at a daily resolution have been compiled and constructed from 1817 until 1996. The data series is composed of thrice daily observations made at Cádiz (1821–1880) by local observers and hourly data from Naval Observatory at San Fernando (1870–1996).     

Simultaneous use of geochemical and geophysical methods to characterise abandoned landfills

. Groundwater sampling and geophysical methods determined a serious contamination problem associated with refilled exploitation sites at the Cal Dimoni area, Llobregat delta, Barcelona, Spain. To characterise this process, hydrogeochemical analyses were performed and showed the following modifications to groundwater chemical composition: increasing pH values, changing redox conditions, significant increases in total organic carbon (TOC) and certain trace elements, and high groundwater conductivity values. Major ion content accumulations were found under the refilled area. In contrast, elements involved in the oxidation–reduction processes, such as iron, manganese and nitrates, clearly diminished. Electromagnetic prospecting methods were also performed and delineated the contamination plume extent. These methods also showed separate sources of contamination, one clearly related to the groundwater–refilled zone leachate interaction, another as a consequence of the manure–accumulation surface site. Geochemical and geophysical methods have shown similar results for locating groundwater contamination sources, and for determining leachate generation mechanisms and flow paths.