Changes of erosive rainfall for El Niño and La Niña years in the northern Andean highlands of Peru

Information related to rainfall erosivity in the Andes is scarce. This study was carried out to determine the characteristics of rainfall events at the La Encañada watershed, northern Peru, using daily rainfall data from the 1995 to 2000 period that included all the El Niño and Southern Oscillation (ENSO) phases. Three weather stations were installed within the study area, at the top, middle and bottom of the watershed. We analysed the total amount, duration, intensity, kinetic energy and probability of return of rainfall events. In general, 80% of the rainfall events at watershed level had an average rainfall intensity lower than 2.5 mm h^?1 and only 4% had an average intensity larger than 7.5 mm h^?1. Rainfall erosivity registered at the bottom of the watershed was slightly higher than in the rest of the area. The highest intensities were observed during an El Niño year whereas a La Niña year was characterized by the highest amount of total rainfall compared to the other ENSO phases and by the low intensity rain events. Simulations using the WEPP model estimated higher sediment yield and runoff for the bottom of the watershed during a La Niña year versus El Niño or Neutral years. Even when the analysed rainfall data was too limited to conclude erosion and runoff during any ENSO phase, the simulated results showed us the trend of the behaviour of rainfall erosivity under the ENSO phases at different locations.     

The impact of recent climate change on flooding and sediment supply within a Mediterranean mountain catchment,southwestern Crete,Greece

This paper presents work from a geomorphological investigation carried out in the Aradena Gorge, southwestern Crete, Greece. The gorge is typical of many steepland fluvial systems in the Mediterranean, with steep relief, coarse‐gravel sediments and high rates of sedimentation generated during intense winter storm events. Hillslope deposits and coarse‐gravel flood units within a 5 km section of the gorge have been mapped, dated (using lichenometry and dendrochronology), and their sedimentological characteristics recorded to establish a c. 200‐year record of flood frequency/magnitude and hillslope/channel sediment supply variability. This record has been compared with instrumented and previously published records of climate change from Crete and the Mediterranean region and used to establish the major controls on flooding and sediment dynamics within the Aradena Gorge. Rates of colluviation and sediment delivery to the channel appear to have been greater than the present sometime before c. AD 1800 and may be related to cooler climates with a more seasonal precipitation regime during the Little Ice Age (c. AD 1450 to 1850). In gorge sections where the present rate of sediment supply from hillslope colluvium is very low, the channel has incised into older alluvial and colluvial deposits. Conversely, in the few sections where sediment supply is currently very high, the channel is aggrading with a braided pattern. Major rock‐fall deposits at certain locations in the gorge have restricted any major downstream sediment transfer. Twelve periods of increased flooding during the last 150 years have been identified and these correlate quite well with negative or declining phases of the North Atlantic Oscillation (NAO). Analysis of daily precipitation data from Crete suggests negative phases of the winter NAO are characterized by an increase in the number of long‐duration, high‐intensity storms. These storms, particularly those with five‐day and greater duration, appear to be significant in triggering major floods in the Aradena Gorge. During the last 40 years the NAO index has been increasing and become locked into a positive phase. As a consequence of this, major flooding appears to have declined during the same period. Copyright © 2002 John Wiley & Sons, Ltd.     

Divariant and trivariant reaction line slopes in FMAS and CFMAS: theory and applications

The theory of divariant reactions is investigated in some detail and approaches involving simultaneously operating reactions and X _Mg isopleths are combined. Using the Gibbs' method (Spear et al. 1982), a general equation is derived for X _Mg isopleth slopes in FMAS and CFMAS, which may be applied to any divariant assemblage in these systems. Isopleth slopes are calculated for a number of divariant FMAS and CFMAS reactions inferred to have operated in Sri Lankan granulites, giving more rigorous constraints on the local P-T path. Isopleth slopes are also calculated for the trivariant CFMAS reaction orthopyroxene + plagioclase = garnet + quartz. The resulting dP/dT range from +8 to +18 bar/K. Application to Sri Lankan granulites confirms an earlier hypothesis (Faulhaber and Raith 1991) that a retrograde path with a dP/dT of ～10–15 bar/K may cause garnet growth in charnockites and garnet breakdown in metabasites. It is shown here that dP/dT depends primarily on the Ca contents of garnet and plagioclase. Large variations in isopleth and reaction line slopes are possible, due to the thermodynamic effect of additional components in solid solutions, e.g. albite in plagioclase. This has a good potential for the assessment of the exact form of P-T paths. Equations to calculate isopleth slopes could be included in software packages aimed at defining P-T estimates, enabling geologists to derive realistic P-T paths.     

Ambiguity resolution strategies using the results of the International GPS Geodynamics Service (IGS)

Resolving the initial phase ambiguities of GPS carrier phase observations was always considered an important aspect of GPS processing techniques. Resolution of the so-called wide-lane ambiguities using a special linear combination of theL ₁ andL ₂ carrier and code observations has become standard. New aspects have to be considered today: (1) Soon AS, the so-called Anti-Spoofing, will be turned on for all Block II spacecrafts. This means that precise code observations will be no longer available, which in turn means that the mentioned approach to resolve the wide-lane ambiguities will fail. (2) Most encouraging is the establishment of the new International GPS Geodynamics Service (IGS), from where high quality orbits, earth rotation parameters, and eventually also ionospheric models will be available. We are reviewing the ambiguity resolution problem under these new aspects: We look for methods to resolve the initial phase ambiguities without using code observations but using high quality orbits and ionospheric models from IGS, and we study the resolution of the narrow-lane ambiguities (after wide-lane ambiguity resolution) using IGS orbits.     

The Hohonu Batholith of North Westland, New Zealand: granitoid compositions controlled by source H2O contents and generated during tectonic transition

Geochemical studies on the Hohonu Batholith, of the West Coast, South Island, New Zealand, have recognised two distinct but chemically related suites of mid-Cretaceous granitoids. The suites are characterised by restricted radiogenic isotopic compositions (Sr_(i) = 0.7062 to 0.7085; ɛNd_(i) = −4.4 to −6.1), and represent melting of a mafic lithosphere source followed by interaction with Ordovician metasediments. The two suites (Te Kinga Suite and Deutgam Suite) are distinguished by contrasting contents of Al₂O₃, Na₂O, Sr, Ba, Eu and HREE, attributable to different residual asssemblages controlled by differing H₂O contents during melting of a metabasaltic source. The relatively mafic, metaluminous, I-type Deutgam Suite represents magmas derived by dehydration melting in equilibrium with an amphibolitic (plagioclase + amphibole) residue. In contrast, the peraluminous, high silica compositions of the Te Kinga Suite were produced by melting at higher H₂O contents, reducing the stability of plagioclase and resulting in a melt in equilibrium with a plagioclase-free eclogitic (garnet + amphibole) residue. Residual plagioclase during generation of the Deutgam Suite resulted in lower Al₂O₃, Na₂O, Sr, Ba and Eu contents, whereas residual garnet during generation of the Te Kinga suite resulted in depleted HREE contents. The mid-Cretaceous granitoids of the Hohonu Batholith were generated during a period of rapid tectonic transition from crustal thickening during collision to crustal thinning and core complex formation during extension. Received: 23 July 1996 / Accepted: 21 August 1997     

Analytic modeling of groundwater dynamics with an approximate impulse response function for areal recharge

An analytic approach is presented for the simulation of variations in the groundwater level due to temporal variations of recharge in surficial aquifers. Such variations, called groundwater dynamics, are computed through convolution of the response function due to an impulse of recharge with a measured time series of recharge. It is proposed to approximate the impulse response function with an exponential function of time which has two parameters that are functions of space only. These parameters are computed by setting the zeroth and first temporal moments of the approximate impulse response function equal to the corresponding moments of the true impulse response function. The zeroth and first moments are modeled with the analytic element method. The zeroth moment may be modeled with existing analytic elements, while new analytic elements are derived for the modeling of the first moment. Moment matching may be applied in the same fashion with other approximate impulse response functions. It is shown that the proposed approach gives accurate results for a circular island through comparison with an exact solution; both a step recharge function and a measured series of 10 years of recharge were used. The presented approach is specifically useful for modeling groundwater dynamics in aquifers with shallow groundwater tables as is demonstrated in a practical application. The analytic element method is a gridless method that allows for the precise placement of ditches and streams that regulate groundwater levels in such aquifers; heads may be computed analytically at any point and at any time. The presented approach may be extended to simulate the effect of other transient stresses (such as fluctuating surface water levels or pumping rates), and to simulate transient effects in multi-aquifer systems.