Dynamic faulting under rate-dependent friction

We discuss the effects of rate-dependent friction on the propagation of seismic rupture on active faults. Several physicists using Burridge and Knopoff's box and spring model of faulting have proposed that fault complexity may arise from the spontaneous development of a self-similar stress distribution on the fault plane. If this model proves to be correct, it has important consequences for the origin of the complexity of seismic sources. In order to test these ideas on a more realistic earthquake model, we developed a new boundary integral equation method for studying rupture propagation along an antiplane fault in the presence of nonlinear rate-dependent friction. We study rupture dynamics of models with single and twin asperities. In our models, asperities are places on the fault with a higher value of prestress. Othewise all fault parameters are homogeneous. We show that for models with such asperities, a slip velocity weakening friction leads to the propagation of supersonic healing phases and to the spontaneous arrest of fracture if the prestress outside the asperities is low enough. For models with asperities, we can also observe narrow slip velocity pulses, qualitatively similar to the so-called Heaton pulses observed in some earthquake accelerograms. We also observe a complex distribution of stress after the rupture that depends on details of the initial distribution of asperities and on the details of the friction law.     

Seismic activity related to the 1991 eruption of Colima Volcano,Mexico

Ten years after the last effusive eruption and at least 15 years of seismic quiescence, volcanic seismic activity started at Colima volcano on 14 February 1991, with a seismic crisis which reached counts of more than 100 per day and showed a diversity of earthquake types. Four other distinct seismic crises followed, before a mild effusive eruption in April 1991. The second crisis preceded the extrusion of an andesitic scoriaceous lava lobe, first reported on 1 March; during this crisis an interesting temporary concentration of seismic foci below the crater was observed shortly before the extrusion was detected. The third crisis was constituted by shallow seismicity, featuring possible mild degassing explosion-induced activity in the form of hiccups (episodes of simple wavelets that repeat with diminishing amplitude), and accompanied by increased fumarolic activity. The growth of the new lava dome was accompanied by changing seismicity. On 16 April during the fifth crisis which consisted of some relatively large, shallow, volcanic earthquakes and numerous avalanches of older dome material, part of the newly extruded dome, which had grown towards the edge of the old dome, collapsed, producing the largest avalanches and ash flows. Afterwards, block lava began to flow slowly along the SW flank of the volcano, generating frequent small incandescent avalanches. The seismicity associated with the stages of this eruptive activity shows some interesting features: most earthquake foci were located north of the summit, some of them relatively deep (7–11 km below the summit level), underneath the saddle between the Colima and the older Nevado volcanoes. An apparently seismic quiet region appears between 4 and 7 km below the summit level. In June, harmonic tremors were detected for the first time, but no changes in the eruptive activity could be correlated with them. After June, the seismicity decreasing trend was established, and the effusive activity stopped on September 1991.     

Space in geography

The concept of space is one of the most fundamental of geographical concepts. There is no work in geography that does not certain it. Nevertheless, geography has not as yet formulated an explicit and unambiguous definition of geographical space. This fact has had negative consequences for geographical theory, methodology and application. In the effort of contributing to its elimination we will try to outline the basic connotations of the concept.The concept of geographical space is a relational one. It acquires meaning and sense only when related to other concepts. The concept of space may be conceived as a supplement to things, i.e. substantively conceived objects. Space conceived in this way is the synonym of emptiness. The concept of space may be also conceived in relation ti individual landscape elements as their environments. Space conceived in this way has the character of a field of force. And, finally, space may be conceived also with respect to the totality of landscape elements, i.e. the system expressed by the term synergic. It is only this third variant of space which should be understood as the geographical space in the full menaing of the term. It is only this conception of geographical space as the space filled with qualities in relations and proportions that is considered as one of the basic prerequisites for the formulation of the theory of geography as science capable of prediction and thus also of practical utilization.     

The Relegation Algorithm

The relegation algorithm extends the method of normalization by Lie transformations. Given a Hamiltonian that is a power series = ₀+ ₁+ ... of a small parameter , normalization constructs a map which converts the principal part ₀into an integral of the transformed system — relegation does the same for an arbitrary function [G]. If the Lie derivative induced by [G] is semi-simple, a double recursion produces the generator of the relegating transformation. The relegation algorithm is illustrated with an elementary example borrowed from galactic dynamics; the exercise serves as a standard against which to test software implementations. Relegation is also applied to the more substantial example of a Keplerian system perturbed by radiation pressure emanating from a rotating source.This revised version was published online in October 2005 with corrections to the Cover Date.     

Entropy‐based equation to assess hillslope sediment production

A new equation to assess hillslope sediment production, based on physical and probabilistic approaches, is presented. The equation, which allows the computation of the delivery ratio for every event, considers the physical variables of travel distance, stream power, settling velocity and gross erosion. The probability density function that arises from the new formulation is solved using the principle of maximum entropy. Based on data from five watersheds in both tropical and temperate zones, the new delivery parameter K_v is calibrated and associated with vegetation cover and conservation practice. The proposed equation is rationally based in relation to parameter K_v. The entropy‐based equation was applied to assess sediment yield in two other experimental watersheds, showing good predictability for the set (mean absolute error of 20·8%). No systematic error was found in the analysed data. The entropy‐based equation showed good predictability for long‐term sets of data and for high‐erosivity events, but did not perform well for the low‐erosivity ones. Copyright © 2007 John Wiley & Sons, Ltd.     

Finding regions of influence on SO2 and SO= 4 daily concentration measurements at four sites in Spain

Contamination by the pollutants SO₂ and SO⁼₄ was analyzed for the 1989–1992 period at four regional stations in Spain under the auspices of the EMEP-BAPMON program. The evolution of the time series of the daily pollution has also been assessed, and high mean concentrations at La Cartuja and Logroño observed, with values of 3.8 and 4.5 g m⁻³ for SO₂, respectively. Maximum annual concentrations were recorded in 1989, when SO₂ reached values of 6.24, 5.39, 5.71, and 9.30 g m⁻³ for the stations of La Cartuja, San Pablo de los Montes, Roquetas, and Logroño, respectively. This work attempts to establish a relationship between the concentrations of the pollutants - both SO₂ gas and SO⁼₄ aerosol - and the zones of emission or persistence of these long-range transported pollutants. In this way, those regions showing a greater impact on the air quality in each season have been determined. To achieve this, the trajectories of the air masses carrying away the pollution to each of the receiving stations were considered and followed by a sectorial analysis. Nonparametric statistical methods were implemented to contrast the chemical homogeneity among the different sectors. The criterion that several homogeneous sectors form a chemically homogeneous region was used. To improve this sectorial analysis, we have proposed a new technique based on the Potential-Source-Contribution Function (PSCF). Starting out from a set of specified regions, considered to be chemically homogeneous domains, it is possible to determine the likelihood that an air mass with particular characteristics (e.g., that a value of the daily concentration higher than the mean recorded at the station has been obtained) will arrive at a given station after having crossed one of the previously defined regions. Using this technique, it is possible to determine the source regions through which the air masses circulate and bring high pollution concentrations to the studied stations. Thanks to the PSCF, these statistical methods offer, through a sectorial analysis, the possibility to pass from a qualitative to a more quantitative view.     

A statistical estimator of the spatial distribution of the water-table altitude

An algorithm was designed to statistically estimate the areal distribution of water-table altitude. The altitude of the water table was bounded below by the minimum water-table surface and above by the land surface. Using lake elevations and stream stages, and interpolating between lakes and streams, the minimum water-table surface was generated. A multiple linear regression among the minimum water-table altitude, the differerence between land-surface and minimum water-table altitudes, and the water-level measurements from surficial aquifier system wells resulted in a consistently high correlation for all groups of physiographic regions in Florida. A simple linear regression between land-surface and water-level measurements resulted in a root-mean-square residual of 4.23 m, with residuals ranging from -8.78 to 41.54 m. A simple linear regression between the minimum water table and the water-level measurements resulted in a root-mean-square residual of 1.45 m, with residuals ranging from -7.39 to 4.10 m. The application of the multiple linear regression presented herein resulted in a root-mean-square residual of 1.05 m, with residuals ranging from -5.24 to 5.63 m. Results from complete and partial F tests rejected the hypothesis of eliminating any of the regressors in the multiple linear regression presented in this study.     

Vulnerability assessment in a volcanic risk evaluation in Central Mexico through a multi-criteria-GIS approach

The Valley of Toluca is a major industrial and agricultural area in Central Mexico, especially the City of Toluca, the capital of The State of Mexico. The Nevado de Toluca volcano is located to the southwest of The Toluca Basin. Results obtained from the vulnerability assessment phase of the study area (5,040 km² and 42 municipalities) are presented here as a part of a comprehensive volcanic risk assessment of The Toluca Basin. Information has been gathered and processed at a municipal level including thematic maps at 1:250,000 scale. A database has been built, classified and analyzed within a GIS environment; additionally, a Multi-Criteria Evaluation (MCE) approach was applied as an aid for the decision-making process. Cartographic results were five vulnerability maps: (1) Total Population, (2) Land Use/Cover, (3) Infrastructure, (4) Economic Units and (5) Total Vulnerability. Our main results suggest that the Toluca and Tianguistenco urban and industrial areas, to the north and northeast of The Valley of Toluca, are the most vulnerable areas, for their high concentration of population, infrastructure, economic activity, and exposure to volcanic events.     

Thermochronology and exhumation rates of granitic intrusions at Mesa Central,Mexico

ABSTRACT

The Mesa Central of Mexico (MC) is an elevated plateau located 2000 m above sea level in central Mexico, where intrusions outcrop that register the history of exhumation-erosion occurring during the Late Cretaceous-Paleogene. The tectonic history of the region records formation of the Late Cretaceous-Paleogene ‘Mexican orogen’; this was followed by extension of the entire region and several plutons were then exhumed. The age and magnitude of the crustal uplift and erosion occurring during exhumation has not been addressed to date. Therefore, this study reports the crystallization and cooling ages of two plutons, the Tesorera Granodiorite and the Comanja Granite, and estimates their emplacement depths. Based on these data, the exhumation age of the Tesorera Granodiorite is estimated to be between ~73 Ma and ~63 Ma at an exhumation rate of ~528 m/m. y. and that of the Comanja Granite is 52 Ma and 48 Ma at an exhumation rate of ~2500 m/m. y. Exhumation-erosion event of the Tesorera Granodiorite was located on the trace of the San Luis-Tepehuanes Fault System and that of the Comanja Granite on the a trace of the El Bajío Fault System. Furthermore, the high exhumation rate in the Comanja Granite suggests that gravitational collapse played an important role during exhumation.     

Analytical study of the energy rate balance equation for the magnetospheric storm-ring current

We present some results of the analytical integration of the energy rate balance equation, assuming that the input energy rate is proportional to the azimuthal interplanetary electric field, E_y, and can be described by simple rectangular or triangular functions, as approximations to the frequently observed shapes of E_y, especially during the passage of magnetic clouds. The input function is also parametrized by a reconnection-transfer efficiency factor (which is assumed to vary between 0.1 and 1). Our aim is to solve the balance equation and derive values for the decay parameter compatible with the observed Dst peak values. To facilitate the analytical integration we assume a constant value for through the main phase of the storm. The model is tested for two isolated and well-monitored intense storms. For these storms the analytical results are compared to those obtained by the numerical integration of the balance equation, based on the interplanetary data collected by the ISEE-3 satellite, with the values parametrized close to those obtained by the analytical study. From the best fit between this numerical integration and the observed Dst the most appropriate values of are then determined. Although we specifically focus on the main phase of the storms, this numerical integration has been also extended to the recovery phase by an independent adjust. The results of the best fit for the recovery phase show that the values of may differ drastically from those corresponding to the main phase. The values of the decay parameter for the main phase of each event, _m, are found to be very sensitive to the adopted efficiency factor, , decreasing as this factor increases. For the recovery phase, which is characterized by very low values of the power input, the response function becomes almost independent of the value of and the resulting values for the decay time parameter, _r, do not vary greatly as varies. As a consequence, the relative values of between the main and the recovery phase, _m/_r, can be greater or smaller than one as varies from 0.1 to 1.