Water losses risk assessment: an example from Carpathian karst

The multipurpose dam and reservoir “Bogovina” in the Crni Timok valley is one of several projects that have been initiated to improve water availability in Serbia. The site and height of the dam were initially proposed without consideration of risk factors such as water leakage or stability of the banks. Later, complex hydrogeological and speleological surveys and tests identified the potential of significant water losses from the reservoir area connected mainly to the karstic aquifer, and the nearby Bogovina cave system. The analyses show that once the reservoir is filled up, groundwater flow currently oriented towards the future reservoir would saturate the upper part of the karstified rocks, reactivate currently unsaturated pathways and form a reverse discharge outside of the reservoir area. In response to these findings, the dam design and technical details have been adapted accordingly: the dam height has been reduced by 9 m, and it is proposed that grouting and consolidation work be conducted both at the foundation of the dam and extensively on the embankments. It was concluded that although the proposed remedial measures cannot guarantee reservoir tightness, they can reduce the risk of large-scale leaking.     

Areal versus linear evaluation of relationship between drainage basin lithology and geochemistry of stream and overbank sediments in low-order mountainous drainage basins

The comparison of two multiple regression models is based on the assumption that geochemical composition of the drainage basin alluvial sediments is derived primarily from the underlying bedrock lithology. The parent material is integrated with both stream sediment and overbank sediment geochemistry via the two essentially different approaches as regards the drainage basin geomorphological data: (1) as the relative area of influence representing a portion of the catchment basin occupied by a specific rock type and (2) as the relative “line” of influence representing a narrow zone of the underlying bedrock traversed by the perennial streams which form the active stream network. The model comparison is established on the goodness-of-fit test for both experimental designs and for the same set of data. Both experiments converge on the linear approach as the more appropriate model in evaluating the lithologic influence on the analysed sample media in small mountainous watersheds.     

The Požarevac loess–paleosol sequence: a record of increased aridity in the south-eastern margin of the Carpathian Basin during the last 350 ka

The loess sequence preserved in the Požarevac brickyard in north-eastern Serbia comprises eight loess units separated by seven paleosols. Geochronological investigation using amino acid racemization and luminescence dating support stratigraphic correlations of loess units L3, S2LL1 and L1 at the Požarevac section with loess of glacial cycles E [Marine Isotope Stage (MIS) 10], D (MIS 9–8), C (MIS 7–6) and B (MIS 5–2) across central Europe. Correlation with the marine oxygen-isotope stratigraphy and associated paleoclimatic inferences are further supported by magnetic susceptibility, particle size and carbonate content measured in Požarevac sediments. Malacological investigations at the Požarevac section reveal the continuous presence of the Chondrula tridens and Helicopsis striata faunal assemblages throughout the last 350 ka. The loess malacological fauna, which is characterized by the complete absence of cold-resistant and cold-preferring species, suggests a stable, dry and relatively warm glacial and interglacial climate, compared with other central European loess localities. Together these data suggest that the south-eastern part of the Carpathian (Pannonian, Middle Danube) Basin was a refugium for warm-preferring and xerophilous land-snails during the generally unfavorable glacial climates of the late Middle and Late Pleistocene.     

Evapotranspiration and canopy characteristics of two lodgepole pine stands following mountain pine beetle attack

Over the past decade, British Columbia (BC), has experienced the largest mountain pine beetle (MPB) outbreak on record. This study used the eddy‐covariance (EC) technique to examine the impact of the MPB attack on evapotranspiration (E) and associated canopy characteristics of two lodgepole pine stands with secondary structure (trees, saplings and seedlings surviving the attack) located in central BC. MPB‐06, an 85‐year‐old almost pure stand of pine trees, was first attacked in 2006, and by 2010, ~80% of the trees had been killed. MPB‐03, a 110‐year‐old stand with an overstory consisting of over 90% pine and a developed sub‐canopy, was first attacked in 2003 and by 2007 had > 95% pine canopy mortality. EC measurements began in August 2006 at MPB‐06 and in March 2007 at MPB‐03, and continued for four years. Annual total E ranged from 226 mm to 237 mm at MPB‐06, and from 280 to 297 mm at MPB‐03, showing relatively little year‐to‐year change at both sites over the four years. Increased E from the accelerated growth of the surviving vegetation (secondary structure, shrubs and herbs) compensated for reduction in E due to the death of the overstory. Monthly average daytime canopy conductance, the Priestley–Taylor (α), and the canopy–atmosphere decoupling coefficient (Ω) steadily increased during the growing season reaching approximate maximum values of 5 mm s^?1, 0.75 and 0.12, respectively. Potential evapotranspiration was approximated using a vapour pressure deficit‐dependent α obtained at high soil water content. Calculated water deficits indicated some water‐supply limitation to the surviving trees and understory at both sites. Rates of root zone drainage during the growing season were low relative to precipitation. Copyright © 2013 John Wiley & Sons, Ltd.     

Energy Balance Closure Using Eddy Covariance Above Two Different Land Surfaces and Implications for CO<Subscript>2</Subscript> Flux Measurements

Components of the surface energy balance of a mature boreal jack pine forest and a jack pine clearcut were analysed to determine the causes of the imbalance that is commonly observed in micrometeorological measurements. At the clearcut site (HJP02), a significant portion of the imbalance was caused by: (i) the overestimation of net radiation (R _n ) due to the inclusion of the tower in the field of view of the downward facing radiometers, and (ii) the underestimation of the latent heat flux (λE) due to the damping of high frequency fluctuations in the water vapour mixing ratio by the sample tube of the closed-path infrared gas analyzer. Loss of low-frequency covariance induced by insufficient averaging time as well as systematic advection of fluxes away from the eddy-covariance (EC) tower were discounted as significant issues. Spatial and temporal distributions of the total surface-layer heat flux (T), i.e. the sum of sensible heat flux (H) and λE, were well behaved and differences between the relative magnitudes of the turbulent fluxes for several investigated energy balance closure (C) classes were observed. Therefore, it can be assumed that micrometeorological processes that affected all turbulent fluxes similarly did not cause the variation in C. Turbulent fluxes measured at the clearcut site should not be forced to close the energy balance. However, at the mature forest site (OJP), loss of low-frequency covariance contributed significantly to the systematic imbalance when a 30-min averaging time was used, but the application of averaging times that were long enough to capture all of the low-frequency covariance was inadequate to resolve all of the high-frequency covariance. Although we found qualitative similarity between T and the net ecosystem exchange (NEE) of carbon dioxide (CO₂), forcing T to closure while retaining the Bowen ratio and applying the same factor to CO₂ fluxes (F _C ) cannot be generally recommended since it remains uncertain to what extent long wavelength contributions affect the relationship between T, F _C and C.     

Dynamics of the Hungaria asteroids

To try to understand the dynamical and collisional evolution of the Hungaria asteroids we have built a large catalog of accurate synthetic proper elements. Using the distribution of the Hungaria, in the spaces of proper elements and of proper frequencies, we can study the dynamical boundaries and the internal structure of the Hungaria region, both within a purely gravitational model and also showing the signature of the non-gravitational effects. We find a complex interaction between secular resonances, mean motion resonances, chaotic behavior and Yarkovsky-driven drift in semimajor axis. We also find a rare occurrence of large scale instabilities, leading to escape from the region. This allows to explain the complex shape of a grouping which we suggest is a collisional family, including most Hungaria but by no means all; we provide an explicit list of non-members of the family. There are finer structures, of which the most significant is a set of very close asteroid couples, with extremely similar proper elements. Some of these could have had, in a comparatively recent past, very close approaches with low relative velocity. We argue that the Hungaria, because of the favorable observing conditions, may soon become the best known sub-group of the asteroid population.     

Pyrite Formation in Organic-rich Clay, Calcitic and Coal-Forming Environments

The early diagenetic characteristics of pyrite formation processes in a Miocene freshwater sequence of mixed sediments （coal fragments in clays, sandstones or shales） alternating with continuous brown coal layers was investigated. Based on abundant minerals, the following main sedimentary environments were distinguished： the illite-montmorillonitic （I-M）, calcitic （Ct） and coal-forming environment （CL）. For these hydrogeochemically differing environments the effects of limiting factors on the pyrite formation process （availability of sulphate and Fe, amount of organic matter and participation of organic sulphur） were assessed by correlation analysis. Significant differences in the effects of these limiting factors in the particular environments were observed. These differences were explained taking in account the different oxidative activity, Fe-complex and surface complex forming properties of humic substances in dependence of pH of environment and the abundance of sorptionally active clay minerals. In environments having a relatively low pH and containing clay minerals （I-M- and CL-environments） the oxidative activity of humic substances （Hs） on pyrite precursors was greatly prevented however pyrite formation depended on reactive Fe availability as the consequence of complex formation. On the contrary, in environments with a relatively high pH, as it was the calcitic, the oxidative activity of Hs was greatly enhanced, thus oxidizing the sulfur precursors of pyrite. The oxidation degree of organic matter was probably also a consequence of the differing activity of the humic electron-acceptors.     

Asteroid Mean Elements: Higher Order and Iterative Theories

Mean orbital elements are obtained from osculating ones by removing the short periodic perturbations. Large catalogues of asteroid mean elements need to be computed, as a first step in the computation of proper elements, used to study asteroid families. The algorithms for this purpose available so far are only accurate to first order in the masses of the perturbing planet; the mean elements have satisfactory accuracy for most of the asteroid belt, but degraded accuracy in the neighbourhoods of the main mean motion resonances, especially the 2:1. We investigate a number of algorithms capable of improving this approximation; they belong to the two classes of Breiter-type methods and iterative methods. The former are obtained by applying some higher order numerical integration scheme, such as Runge–Kutta, to the differential equation whose solution is a transformation removing the fast angular variables from the equations; they can be used to compute a full second order theory, however, only if the full second order determining function is explicitly computed, and this is computationally too cumbersome for a complicated problem such as the N-body. The latter are fixed point iterative schemes, with the first order theory as an iteration step, used to compute the inverse map from mean to osculating elements; formally the method is first order, but because they implement a fixed frequency perturbation theory, they are more accurate than conventional single iteration methods; a similar method is already in use in our computation of proper from mean elements. Many of these methods are tested on a sample of asteroid orbits taken from the Themis family, up to the edge of the 2:1 resonance, and the dispersion of the values of the computed mean semimajor axis over 100 000 years is used as quality control. The results of these tests indicate that the iterative methods are superior, in this specific application, to the Breiter methods, in accuracy and reliability. This is understood as the result of the cancellations occurring between second order perturbation terms: the incomplete second order theory, resulting from the use of a Breiter method with the first order determining function only, can be less accurate than complete, fixed frequency theories of the first order. We have therefore computed new catalogues of asteroid mean and proper elements, incorporating an iterative algorithm in both steps (osculating to mean and mean to proper elements). This new data set, significantly more reliable even in the previously degraded regions of Themis and Cybele, is in the public domain. This revised version was published online in July 2006 with corrections to the Cover Date.     

Asteroid long-periodic perturbations: The second order Hamiltonian

The Hamiltonian of the second order with respect to the disturbing mass, as defined in the higher order-higher degree theory of asteroid secular perturbations by Yuasa (1973), is expressed in the heliocentric, ecliptic coordinate system. Errors found in the original paper with terms coming from the principal part of the disturbing function are removed, and corrected values of the coefficients are computed. The importance of second-order perturbations and the improvement in the accuracy of proper element determination, achieved by using the newly-obtained coefficients, are demonstrated. Finally, a table of the secular frequencies as functions of the semimajor axis is given, and compared with the analogous one by Kozai (1979).     

Families among high-inclination asteroids

We present a new classification of families identified among the population of high-inclination asteroids. We computed synthetic proper elements for a sample of 18,560 numbered and multi-opposition objects having sine of proper inclination greater than 0.295. We considered three zones at different heliocentric distances (inner, intermediate and outer region) and used the standard approach based on the Hierarchical Clustering Method (HCM) to identify families in each zone. In doing so, we used slightly different approach with respect to previously published methodologies, to achieve a more reliable and robust classification. We also used available SDSS color data to improve membership and identify likely family interlopers. We found a total of 38 families, as well as a significant number of clumps and clusters deserving further investigation.