Forest operations,tree species composition and decline in rainfall explain runoff changes in the Nacimiento experimental catchments,south central Chile

Few long-term studies have explored how intensively managed short rotation forest plantations interact with climate variability. We examine how prolonged severe drought and forest operations affect runoff in 11 experimental catchments on private corporate forest land near Nacimiento in south central Chile over the period 2008–2019. The catchments (7.7–414 ha) contain forest plantations of exotic fast-growing species (Pinus radiata, Eucalyptus spp.) at various stages of growth in a Mediterranean climate (mean long-term annual rainfall = 1381 mm). Since 2010, a drought, unprecedented in recent history, has reduced rainfall at Nacimiento by 20%, relative to the long-term mean. Pre-drought runoff ratios were <0.2 under 8-year-old Eucalyptus; >0.4 under 21-year-old Radiata pine and >0.8 where herbicide treatments had controlled vegetation for 2 years in 38% of the catchment area. Early in the study period, clearcutting of Radiata pine (85%–95% of catchment area) increased streamflow by 150 mm as compared with the year before harvest, while clearcutting and partial cuts of Eucalyptus did not increase streamflow. During 2008–2019, the combination of emerging drought and forestry treatments (replanting with Eucalyptus after clearcutting of Radiata pine and Eucalyptus) reduced streamflow by 400–500 mm, and regeneration of previously herbicide-treated vegetation combined with growth of Eucalyptus plantations reduced streamflow by 1125 mm (87% of mean annual precipitation 2010–2019). These results from one of the most comprehensive forest catchment studies in the world on private industrial forest land indicate that multiple decades of forest management have reduced deep soil moisture reservoirs. This effect has been exacerbated by drought and conversion from Radiata pine to Eucalyptus, apparently largely eliminating subsurface supply to streamflow. The findings reveal tradeoffs between wood production and water supply, provide lessons for adapting forest management to the projected future drier climate in Chile, and underscore the need for continued experimental work in managed forest plantations.     

Basin filling evolution of the central basins of Mallorca since the Pliocene

A new compilation of data from 436 drill cores using decompaction and backstripping techniques was used to reconstruct the basin filling history from the Pliocene until the present day in the Palma, Inca and Sa Pobla Basins on the island of Mallorca (Spain). Calcareous rocks dominate the source area and provide a limited amount of clastic input to the basins that has resulted in an average accumulation rate of between 5 and 20 m/Ma during the last 5.3 Ma. Carbonate sediment production dominated the basin filling history during early‐mid Pliocene, but during the Quaternary, the sedimentation processes in the Palma Basin were probably enhanced by an evolution in the drainage network that increased the sediment supply and the accumulated thickness caused by stream capture. However, the maximum sedimentation rate filling the depocentres of the three basins has been decreasing since the Pliocene, showing that not only the catchment transport efficiency but also the relative sea level have been controlling the sediment accumulation in these carbonate basins. The isopach cross‐sections support the idea that a palaeorelief was generated during the Messinian sea level drop and that heterogeneities were filled in from the Pliocene to the Quaternary. We conclude that the central basins of Mallorca were filled heterogeneously due to tectonic and geomorphic processes that controlled sediment transport and production, resulting in different average sedimentation thicknesses that decreased since the Pliocene as the accommodation space became filled and the relative sea level dropped.     

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Geothermal palaeogradients and metamorphic zonation from the conodont colour alteration index (CAI)

The conodont colour alteration index (CAI) is potentially valuable in thermal history investigations, but there are problems in its use. An approach is proposed herein to overcome the problem of establishing reliable geothermal palaeogradients in areas with burial metamorphism, based on knowledge of the stratigraphic succession thickness, the boundaries of chronostratigraphic units, and the age and duration of deformation episodes. Although also problematical that a correlation between CAI and Kübler index (KI) values does not exist, there are areas where CAI and KI values are notably consistent for the anchizone boundaries; however, anomalous correlations between CAI and KI values have also been found in some areas. A CAI anchizone (or ancaizone) is defined in order to enable CAI values to be used independently to establish a metamorphic zonation.     

Fundamental parameters and new variables of the galactic open cluster NGC 7128

CCD photometry in Johnson UBV and Strömgren uvby systems and medium-resolution spectroscopy of the galactic open cluster NGC 7128 are presented. Spectral types of the brightest 12 stars in the cluster field were determined based on equivalent widths of the H α and the He  i 6678-Å line. The spectroscopic observations also revealed two obvious and one probable Be-type stars showing H α emission. The analysis of the photometric diagrams gave a colour excess of E ( B − V )=1.03±0.06 mag, a distance modulus DM =13.0±0.2 mag and an age above 10 Myr. Time-resolved photometric observations obtained on one night resulted in the detection of short time-scale light variations of seven new and three already known variable stars in the cluster field.     

Lead isotope systematics of Late Cretaceous – Tertiary Andean arc magmas and associated ores between 8°N and 40°S: evidence for latitudinal mantle heterogeneity beneath the Andes

ABSTRACT Lead isotope variability of magmatic arc rocks and associated mineralization of the Central Andes is usually considered to be the result of mixing between a homogeneous mantle and heterogeneous continental crust. About 230 new lead isotope data on the Northern and Central Andes allow us to compare for the first time lead isotope systematics of the Late Cretaceous – Tertiary arc magmatism and associated mineralization along the Andean chain between 8°N and 40°S. Lead isotope compositions indicate mixing between mantle and upper crustal rocks along the whole Andean chain. Additionally, we have found that mantle end-members of the Late Cretaceous – Tertiary magmatism are heterogeneous and systematically shifted towards less radiogenic ²⁰⁶Pb/²⁰⁴Pb compositions from north to south along the Andes. This heterogeneity most likely results from mixing between a low radiogenic mantle, possibly carrying a DMM or EM I component, and a more radiogenic mantle, possibly carrying an HIMU component. Thus, our results imply that lead isotope variability of Andean magmas at the continental scale is caused not only by crustal but also by mantle heterogeneity.     

The effects of Population III stars on the chemical and photometrical evolution of ellipticals

We have studied the effects of a hypothetical initial generation containing very massive stars [   M > 100 M_⊙  , pair-creation supernovae] on the chemical and photometric evolution of elliptical galaxies. To this purpose, we have computed the evolution of a typical elliptical galaxy with luminous mass of the order of  10¹¹ M_⊙  and adopted chemical evolution models already tested to reproduce the main features of ellipticals. We have tested several sets of yields for very massive zero-metallicity stars: these stars should produce quite different amounts of heavy elements than lower-mass stars. We found that the effects of Population III stars on the chemical enrichment is negligible if only one or two generations of such stars occurred, whereas they produce quite different results from the standard models if they continuously formed for a period not shorter than 0.1 Gyr. In this case, the results are at variance with the main observational constraints of ellipticals such as the average  [〈α/ Fe〉_*]  ratio in stars and the integrated colours. Therefore, we conclude that if Population III stars ever existed they must have been present for a very short period of time and their effects on the following evolution of the parent galaxy must have been negligible. This effect is minimum if a more realistic model with initial infall of gas rather than the classic monolithic model is adopted. Ultimately, we conclude that there is no need to invoke a generation of very massive stars in ellipticals to explain their chemical and photometric properties.     

Chemical evolution of Elliptical Galaxies and the ICM

We present a new model for the chemical evolution of elliptical galaxies taking into account SN feedback, detailed nucleosynthesis and galactic winds. We discuss the effect of galactic winds on the chemical enrichment of the ICM and compute the energy per particle injected by the galaxies into the ICM. This revised version was published online in August 2006 with corrections to the Cover Date.