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.     

Determination of subcatchment and watershed boundaries in a complex and highly urbanized landscape

Urban development significantly alters the landscape by introducing widespread impervious surfaces, which quickly convey surface run‐off to streams via stormwater sewer networks, resulting in “flashy” hydrological responses. Here, we present the inadequacies of using raster‐based digital elevation models and flow‐direction algorithms to delineate large and highly urbanized watersheds and propose an alternative approach that accounts for the influence of anthropogenically modified land cover. We use a semi‐automated approach that incorporates conventional drainage networks into overland flow paths and define the maximal run‐off contributing area. In this approach, stormwater pipes are clustered according to their slope attributes, which define flow direction. Land areas drained by each cluster and contributing (or exporting) flow to a topographically delineated catchment were determined. These land masses were subsequently added or removed from the catchment, modifying both the shape and the size. Our results in a highly urbanized Toronto, Canada, area watershed indicate a moderate net increase in the directly connected watershed area by 3% relative to a topographically forced method; however, differences across three smaller scale subcatchments are greater. Compared to topographic delineation, the directly connected watershed areas of both the upper and middle subcatchments decrease by 5% and 8%, respectively, whereas the lower subcatchment area increases by 15%. This is directly related to subsurface storm sewer pipes that cross topographic boundaries. When directly connected subcatchment area is plotted against total streamflow and flashiness indices using this method, the coefficients of variation are greater (0.93 to 0.97) compared to the use of digital elevation model‐derived subcatchment areas (0.78 to 0.85). The accurate identification of watershed and subcatchment boundaries should incorporate ancillary data such as stormwater sewer networks and retention basin drainage areas to reduce water budget errors in urban systems.     

Assessment of sandstone deterioration at Ta Keo temple (Angkor): first results and future prospects

A first application of geomorphological methods to the assessment of sandstone deterioration at Angkor is presented. Damage diagnosis was carried out on the first eastern tier of the central pyramid of the 1,000 year-old Ta Keo temple. Methods combine field observations and measurements at 230 sampling points, high-resolution lasergrammetry and stereophotogrammetry on a 2-m² test zone, and SEM observations. The first results indicate that decay operates through a synergistic combination of weathering phenomena dominated by scaling and solution, and exhibits a high spatial variability. Percentages of deteriorated surfaces vary from 17.6 to 93.8%, and average stone recession values from 0.00 to 2.71 cm (minimum) and 0.34 to 5.49 cm (maximum). On the test-zone, stereophotogrammetry and 3D-mapping of the present and reconstructed initial states using lasergrammetry indicate that erosion scars up to 6 cm deep have formed since 1963. On the whole, the amount of deteriorated surfaces more than tripled between 1963 and 2008. The degree of implication of salts in stone decay remains unclear for most efflorescences are composed of calcite (CaCO3), with secondary importance of barite (BaSO4) and gypsum (CaSO₄·2H₂O). Future prospects aim to evaluate the impact on stone decay of the clearing out of the temple from the forest in the 1920s.     

Otolith shape as a stock discrimination tool for ling (Genypterus blacodes) in the fjords of Chilean Patagonia

ABSTRACT

Genypterus blacodes, in terms of its fishing history and local economic importance, is an emblematic species harvested in Chilean Patagonia (41°00’–57°00’S). Most of the current fisheries and biological knowledge of this species come from the open ocean, whereas information about the species in fjords and inner channels is fragmentary. In 2018, two research surveys targeting G. blacodes were conducted in the fjords and inner channels of Chilean Patagonia. A total of 253 pairs of sagittal otoliths were sampled at three different localities, and their contours were modelled using wavelet analysis as a tool for stock discrimination. Contours were compared using canonical analysis, and classification was performed using linear discriminant and Random Forest analyses. The results indicated that the wavelet method is efficient in modelling otolith contours, and the discriminant analyses showed differences among fishing grounds across the latitudinal gradient, thus confirming the hypothesis that G. blacodes conform to at least two separate stock units in Chilean Patagonia. Fishing grounds that were closer in space showed higher levels of misclassification. The discussion focuses on how environmental variables and the geography of fjords shape stock differences and how this information can be used for the sustainable management of G. blacodes.     

Oxidative stress in Laeonereis acuta (Polychaeta, Nereididae): environmental and seasonal effects

Laeonereis acuta was seasonally collected in an industrially polluted site (P) and in an unpolluted site (UP) at the Patos Lagoon estuary (southern Brazil). Glutathione S-transferase (GST), superoxide dismutase (SOD) and catalase (CAT) activity (U/mg protein) was determined in five groups of worms from each site. Metallothionein (MT - mol GSH/g ww) and lipid peroxides content (LPO - nmoles of cumene hydroperoxide/g ww) were also measured. Annual mean values for CAT (UP=3.7+/-0.3; P=5.7+/-0.6), GST (UP=0.034+/-0.003; P=0.045+/-0.004) and MT (UP=0.15+/-0.02; P=0.23+/-0.03) were higher (p<0.05) in worms from the P site. In autumn, CAT activity was higher (p<0.05) in worms from the P site (7.6 +/- 1.3) than in those from the UP site (3.6 +/- 0.4). In summer, MT concentration was higher in worms from the P site (0.37 +/- 0.03) than in those from the UP site (0.19 +/- 0.01). No significant difference (p>0.05) in the LPO content was observed in worms from the different sites or collected in different seasons. These results indicate that worms from the polluted site showed higher antioxidant responses than those from the unpolluted site, sufficient to prevent oxidative damage in terms of LPO.     

Bayesian estimation of seismic hazard for two sites in Switzerland

Seismic hazard analysis is based on data and models, which both are imprecise and uncertain. Especially the interpretation of historical information into earthquake parameters, e.g. earthquake size and location, yields ambiguous and imprecise data. Models based on probability distributions have been developed in order to quantify and represent these uncertainties. Nevertheless, the majority of the procedures applied in seismic hazard assessment do not take into account these uncertainties, nor do they show the variance of the results. Therefore, a procedure based on Bayesian statistics was developed to estimate return periods for different ground motion intensities (MSK scale).Bayesian techniques provide a mathematical model to estimate the distribution of random variables in presence of uncertainties. The developed method estimates the probability distribution of the number of occurrences in a Poisson process described by the parameter . The input data are the historical occurrences of intensities for a particular site, represented by a discrete probability distribution for each earthquake. The calculation of these historical occurrences requires a careful preparation of all input parameters, i.e. a modelling of their uncertainties. The obtained results show that the variance of the recurrence rate is smaller in regions with higher seismic activity than in less active regions. It can also be demonstrated that long return periods cannot be estimated with confidence, because the time period of observation is too short. This indicates that the long return periods obtained by seismic source methods only reflects the delineated seismic sources and the chosen earthquake size distribution law.     

Red dust rain within the Spanish Mediterranean area

Dust rain belongs to the climatological conditions of the Iberian Peninsula's Mediterranean seaboard. Traditionally known as muddy' or bloody rains, red dust rain has become a topical issue as a result of more frequent rainfall of this kind in recent years. In spite of the difficulties that studying this phenomenon involves, owing to the lack of systematically kept records at meteorological observatories, details are provided of the high frequency of such phenomena within the Spanish Mediterranean area, using records from other phenological sources. The study analyzes the chemical composition of red dust rain and the atmospheric factors giving rise to such phenomena. Whether there is a relationship betwen the higher frequency of dust rain and possible changes in global atmospheric circulation is as yet unclear.     

About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

Molecular indicators for palaeoenvironmental change in a Messinian evaporitic sequence (Vena del Gesso, Italy). II: High-resolution variations in abundances and 13C contents of free and sulphur-bound carbon skeletons in a single marl bed

The extractable organic matter of 10 immature samples from a marl bed of one evaporitic cycle of the Vena del Gesso sediments (Gessoso-solfifera Fm., Messinian, Italy) was analyzed quantitatively for free hydrocarbons and organic sulphur compounds. Nickel boride was used as a desulphurizing agent to recover sulphur-bound lipids from the polar and asphaltene fractions. Carbon isotopic compositions (delta vs PDB) of free hydrocarbons and of S-bound hydrocarbons were also measured. Relationships between these carbon skeletons, precursor biolipids, and the organisms producing them could then be examined. Concentrations of S-bound lipids and free hydrocarbons and their delta values were plotted vs depth in the marl bed and the profiles were interpreted in terms of variations in source organisms, 13 C contents of the carbon source, and environmentally induced changes in isotopic fractionation. The overall range of delta values measured was 24.7%, from -11.6% for a component derived from green sulphur bacteria (Chlorobiaceae) to -36.3% for a lipid derived from purple sulphur bacteria (Chromatiaceae). Deconvolution of mixtures of components deriving from multiple sources (green and purple sulphur bacteria, coccolithophorids, microalgae and higher plants) was sometimes possible because both quantitative and isotopic data were available and because either the free or S-bound pool sometimes appeared to contain material from a single source. Several free n-alkanes and S-bound lipids appeared to be specific products of upper-water-column primary producers (i.e. algae and cyanobacteria). Others derived from anaerobic photoautotrophs and from heterotrophic protozoa (ciliates), which apparently fed partly on Chlorobiaceae. Four groups of n-alkanes produced by algae or cyanobacteria were also recognized based on systematic variations of abundance and isotopic composition with depth. For hydrocarbons probably derived from microalgae, isotopic variations are well correlated with those of total organic carbon. A resistant aliphatic biomacromolecule produced by microalgae is, therefore, probably an important component of the kerogen. These variations reflect changes in the depositional environment and early diagenetic transformations. Changes in the concentrations of S-bound lipids induced by variations in conditions favourable for sulphurization were discriminated from those related to variations in primary producer assemblages. The water column of the lagoonal basin was stratified and photic zone anoxia occurred during the early and middle stages of marl deposition. During the last stage of the marl deposition the stratification collapsed due to a significant shallowing of the water column. Contributions from anaerobic photoautotrophs were apparently associated with variations in depth of the chemocline.     

Dynamical evolution of NEAs: Close encounters, secular perturbations and resonances

We discuss the main mechanisms affecting the dynamical evolution of Near-Earth Asteroids (NEAs) by analyzing the results of three numerical integrations over 1 Myr of the NEA (4179) Toutatis. In the first integration the only perturbing planet is the Earth. So the evolution is dominated by close encounters and looks like a random walk in semimajor axis and a correlated random walk in eccentricity, keeping almost constant the perihelion distance and the Tisserand invariant. In the second integration Jupiter and Saturn are present instead of the Earth, and the 3/1 (mean motion) and v ₆ (secular) resonances substantially change the eccentricity but not the semimajor axis. The third, most realistic, integration including all the three planets together shows a complex interplay of effects, with close encounters switching the orbit between different resonant states and no approximate conservation of the Tisserand invariant. This shows that simplified 3-body or 4-body models cannot be used to predict the typical evolution patterns and time scales of NEAs, and in particular that resonances provide some fast-track dynamical routes from low-eccentricity to very eccentric, planet-crossing orbits.On leave from the Department of Mathematics, University of Pisa, Via Buonarroti 2, 56127 Pisa, Italy, thanks to the G. Colombo fellowships of the European Space Agency.