Pyroclastic Flow Hazard at Volcán Citlaltépetl

Volcán Citlaltépetl (Pico de Orizaba) with an elevation of 5,675 m is the highest volcano in North America. Its most recent catastrophic events involved the production of pyroclastic flows that erupted approximately 4,000, 8,500, and 13,000 years ago. The distribution of mapped deposits from these eruptions gives an approximate guide to the extent of products from potential future eruptions. Because the topography of this volcano is constantly changing computer simulations were made on the present topography using three computer algorithms: energy cone, FLOW2D, and FLOW3D. The Heim Coefficient (), used as a code parameter for frictional sliding in all our algorithms, is the ratio of the assumed drop in elevation (H) divided by the lateral extent of the mapped deposits (L). The viscosity parameter for the FLOW2D and FLOW3D codes was adjusted so that the paths of the flows mimicked those inferred from the mapped deposits. We modeled two categories of pyroclastic flows modeled for the level I and level II events. Level I pyroclastic flows correspond to small but more frequent block-and-ash flows that remain on the main cone. Level II flows correspond to more widespread flows from catastrophic eruptions with an approximate 4,000-year repose period. We developed hazard maps from simulations based on a National Imagery and Mapping Agency (NIMA) DTED-1 DEM with a 90 m grid and a vertical accuracy of ±30 m. Because realistic visualization is an important aid to understanding the risks related to volcanic hazards we present the DEM as modeled by FLOW3D. The model shows that the pyroclastic flows extend for much greater distances to the east of the volcano summit where the topographic relief is nearly 4,300 m. This study was used to plot hazard zones for pyroclastic flows in the official hazard map that was published recently.     

Thickness estimates of the upper-mantle transition zone from bootstrapped velocity spectrum stacks of receiver functions

We modify the receiver-functions stacking technique known as velocity spectrum stacking (VSS) so as to estimate combinations of velocity model ( V_P and V_S ) and depth that stack the Ps conversion from upper-mantle discontinuities most coherently. We find that by estimating the differences in the depths to the 660 and 410 km discontinuities using velocities that maximize the stacked amplitudes of P410s and P660s phases we can estimate the thickness of the transition zone more accurately than the depths to either of these discontinuities. We present two examples indicating that the transition zone beneath Obninsk, Russia, is 252±6 km thick and that beneath Pasadena, California, is only 220±6 km thick.     

Sediment yield and storage variations in the Négron River catchment (south western Parisian Basin,France) during the Holocene period

In the Négron River catchment area (162 km²), surface‐sediment stores are composed of periglacial calcareous ‘grèze’ (5 × 10⁶ t) and loess (21 × 10⁶ t), and Holocene alluvium (12·6 × 10⁶ t), peat (0·6 × 10⁶ t) and colluvium (18·5 × 10⁶ t). Seventy‐five per cent of the Holocene sediments is stored along the thalwegs. Present net sediment yield, calculated from solid discharge at the Négron outlet, is low (0·6 t km^?2 a^?1) due to the dominance of carbonate rocks in the catchment. Mean sediment yield during the Holocene period is 7·0 t km^?2 a^?1 from alluvium stores and 7·6 t km^?2 a^?1 from colluvium stores. Thus, the gross sediment yield during the Holocene period is about 18·7 t km^?2 a^?1 and the sediment delivery ratio 3 per cent. The yield considerably varies from one sub‐basin to another (3·9 to 24·5 t km^?2 a^?1) according to lithology: about 25 per cent and 50 per cent of initial stores of periglacial grèze and loess respectively were reworked during the Holocene period. Sediment yield has increased by a factor of 6 in the last 1000 years, due to the development of agriculture. The very high rate of sediment storage on the slope during that period (88 per cent of the yield) can be accounted for by the formation of cultivation steps (‘rideaux’). It is predicted that the current destruction of these steps will result in a sediment wave reaching the valley floors in the coming decades. Subboreal and Subatlantic sediments and pollen assemblages in the Taligny marsh, where one‐third of the alluvium is stored, show the predominant influence of human activity during these periods in the Négron catchment. Copyright © 2002 John Wiley & Sons, Ltd.     

Absolute velocity of North America during the Mesozoic from paleomagnetic data

We use paleomagnetic data to map Mesozoic absolute motion of North America, using paleomagnetic Euler poles (PEP). First, we address two important questions: (1) How much clockwise rotation has been experienced by crustal blocks within and adjacent to the Colorado Plateau? (2) Why is there disagreement between the apparent polar wander (APW) path constructed using poles from southwestern North America and the alternative path based on poles from eastern North America? Regarding (1), a 10.5° clockwise rotation of the Colorado Plateau about a pole located near 35°N, 102°W seems to fit the evidence best. Regarding (2), it appears that some rock units from the Appalachian region retain a hard overprint acquired during the mid-Cretaceous, when the geomagnetic field had constant normal polarity and APW was negligible.We found three well-defined small-circle APW tracks: 245–200 Ma (PEP at 39.2°N, 245.2°E, R=81.1°, root mean square error (RMS)=1.82°), 200–160 Ma (38.5°N, 270.1°E, R=80.4°, RMS=1.06°), 160 to 125 Ma (45.1°N, 48.5°E, R=60.7°, RMS=1.84°). Intersections of these tracks (the “cusps” of Gordon et al. [Tectonics 3 (1984) 499]) are located at 59.6°N, 69.5°E (the 200 Ma or “J1” cusp) and 48.9°N, 144.0°E (the 160 Ma or “J2” cusp). At these times, the absolute velocity of North America appears to have changed abruptly.North America absolute motion also changed abruptly at the beginning and end of the Cretaceous APW stillstand, currently dated at about 125 and 88 Ma (J. Geophys. Res. 97 (1992b) 19651). During this interval, the APW path degenerates into a single point, implying rotation about an Euler pole coincident with the spin axis.Using our PEP and cusp locations, we calculate the absolute motion of seven points on the North American continent. Our intention is to provide a chronological framework for the analysis of Mesozoic tectonics. Clearly, if APW is caused by plate motion, abrupt changes in absolute motion should correlate with major tectonic events. This follows because large accelerations reflect important changes in the balance of forces acting on the plate, the most important of which are edge effects (subduction, terrane accretion, etc.). Some tectonic interpretations: (1) The J1 cusp may be associated with the inception of rifting of North America away from land masses to the east; the J2 cusp seems to mark the beginning of rapid spreading in the North Atlantic. (2) The J2 cusp signals the beginning of a period of rapid northwestward absolute motion of western North America; motion of tectonostratigraphic terranes in the westernmost Cordillera seems likely to have been directed toward the south during this interval. (3) The interval 88 to 80 Ma saw a rapid decrease in the paleolatitude of North America; unless this represents a period of true polar wander, terrane motion during this time should have been relatively northward.     

Évolution de l'onde semi-diurne M2 de la marée à Brest de 1846 à 2005

The work of searching, recovering and quality control of ancient sea-level measurements at Brest is presented. This work enables us to complete a study carried out by Cartwright in 1972, which showed a decrease in the tidal M2 semi-diurnal amplitude of 1% per century. After including these ancient data, as well as the last four decades of observations in the analysis, our results show an increase of the amplitude of M2 after 1960 and a decrease before 1880, suggesting a long-period oscillation rather than a steady secular trend. To cite this article: N. Pouvreau et al., C. R. Geoscience 338 (2006).     

Quantum-thermodynamic treatment of intrinsic anharmonicity; Wallace’s theorem revisited

Wallace (in Thermodynamics of crystals, 1972) developed a theorem, rooted in rigid lattice dynamics, which incorporates intrinsic anharmonic effects in solids. The practical application of this theorem in mineral physics is computationally involved and this is the main reason for the theorem not getting the attention it deserves. Because intrinsic anharmonicity is an important issue at the extreme conditions in planetary mantles, we derived a method which removes the computational obstacles in applying this theorem. We extended the theorem to incorporate details of the phonon spectrum and tested our algorithm on forsterite (Mg₂SiO₄). Using a least squares inversion technique applied to all available experimental data, we show that it results in an accurate representation of thermodynamic properties and sound wave velocities of Mg₂SiO₄ in its complete pressure–temperature stability range. We also show that the accuracy of our results is not significantly affected by the use of a different equation of state.     

Fitting the Linear Model of Coregionalization by Generalized Least Squares

In geostatistical studies, the fitting of the linear model of coregionalization (LMC) to direct and cross experimental semivariograms is usually performed with a weighted least-squares (WLS) procedure based on the number of pairs of observations at each lag. So far, no study has investigated the efficiency of other least-squares procedures, such as ordinary least squares (OLS), generalized least squares (GLS), and WLS with other weighing functions, in the context of the LMC. In this article, we compare the statistical properties of the sill estimators obtained with eight least-squares procedures for fitting the LMC: OLS, four WLS, and three GLS. The WLS procedures are based on approximations of the variance of semivariogram estimates at each distance lag. The GLS procedures use a variance–covariance matrix of semivariogram estimates that is (i) estimated using the fourth-order moments with sill estimates (GLS₁), (ii) calculated using the fourth-order moments with the theoretical sills (GLS₂), and (iii) based on an approximation using the correlation between semivariogram estimates in the case of spatial independence of the observations (GLS₃). The current algorithm for fitting the LMC by WLS while ensuring the positive semidefiniteness of sill matrix estimates is modified to include any least-squares procedure. A Monte Carlo study is performed for 16 scenarios corresponding to different combinations of the number of variables, number of spatial structures, values of ranges, and scale dependence of the correlations among variables. Simulation results show that the mean square error is accounted for mostly by the variance of the sill estimators instead of their squared bias. Overall, the estimated GLS₁ and theoretical GLS₂ are the most efficient, followed by the WLS procedure that is based on the number of pairs of observations and the average distance at each lag. On that basis, GLS₁ can be recommended for future studies using the LMC.