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.     

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.     

Hydrochemical properties of transition zone between fresh groundwater and seawater in karst environment of the Adriatic islands,Croatia

Eastern coast of the Adriatic Sea consists of karstified carbonates. It belongs to the well-known Dinaric karst region. The coast is extremely indented and there are 718 islands with numerous rock crags and reefs. Some of the inhabited islands use their own water resources for public water supply, or plan to do it in the future. Since karst rocks are extremely permeable, the seawater intrudes into underground water resources, thereby forming the wedge. A wide transition zone occurs between this seawater wedge and fresh water aquifers. Consequently, island groundwater reserves turn brackish to a certain extent. In this study, 77 water samples were collected from a wide variety of water resources. Comprehensive statistical and mathematical multivariate analysis of these data was performed. Simple statistical approach showed several useful correlations among some parameters, and more complex multivariate techniques extracted three factors in connection with three natural processes: (1) mixing with the seawater, (2) carbonate dissolution and (3) human influence (pollution) and nitrogen transformation processes. The results of this study demonstrate that in situ measurement of electrical conductivity is adequate for the very rough field estimation of numerous parameters.     

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.     

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 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).     

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.     

Water flow and solute transport model of potentially toxic elements through unsaturated zone at regional wellfield Kosnica

The objective of this work was to build a prognostic water flow model and potentially toxic elements (lead, cadmium, zinc) transport model in the unsaturated zone. Research was conducted in the catchment area of Kosnica regional wellfield, where the unsaturated zone is characterised by Fluvisol. Lower sorption capacities were determined in the first horizons for all three potentially toxic elements. Correlation coefficient of the measured and simulated values of tracer concentration is 0.58 for the AC horizon and 0.84 for the 2C/C1 horizon. Based on calibrated water flow and transport parameters, a prognostic water flow model and potentially toxic elements (lead, cadmium, zinc) transport model in the unsaturated zone was built. In case of an accidental spill of potentially toxic elements with concentrations of 1000 mg/l, the risk of contamination of the aquifer is present. Copyright © 2016 John Wiley & Sons, Ltd.     

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.