Groundwater recharge and evapotranspiration for two natural ecosystems covered with oak and heather

The evapotranspiration and groundwater recharge from two natural areas with high (oak) and low (heath) vegetation were estimated by calibrating a semi-physical numerical soil water and heat model to fit 8 and 7 years of TDR-measurements of water content, respectively. The measurements were made between the surface and 7 m depth. For the oak stand, the estimated annual recharge for the years 1992–1999 is 390 mm, the evaporation from soil and interception is 205 mm, and the transpiration is 285 mm. For the heath area estimation was carried out for the years 1993–1999. However, the heath was struck by a heavy beetle attack in 1994, which strongly affected the vegetation and thus the water balance for the following 3 years. For years not affected, the estimated recharge is 733 mm (about 50% larger than for the oak stand for the same years), the evaporation is 316 mm, and the transpiration is 128 mm. The estimated recharge values compare fairly well to estimates obtained from bromide tracer experiments. However, the recharge estimates obtained from the tracer experiments are very uncertain. The uncertainty is mainly due to spatial heterogeneity making the three replicate samples taken here for each time and depth insufficient.

The analyses of TDR-measurements and tracer data showed that water front movement depends on the antecedent soil water content. Some layers are bypassed, especially at low water contents, and at high soil water contents preferential flow was observed at the heath site.     

Subduction of the Pacific plate under Kamchatka: Seismic velocity of underthrust

The results of studying the process of underthrusting in the Pacific plate under Kamchatka are presented, and the phenomena associated with this process are described. The seismic component of the velocity of the plate underthrust is estimated on the basis of (1) data from the CMT (Centroid Moment Tensor) catalog and (2) the sequence of the strongest Kamchatka earthquakes. A flat site with a strike azimuth of 217° and a dip angle of 25° located at depths of 30–70 km is assumed to be the interface between the plunging plate and Kamchatka. From CMT focal mechanisms, the underthrust velocity is estimated at V = 0.7 cm/yr for the southern zone (located south of Cape Shipunski) and at V = 1.1 cm/yr for the central zone (from Cape Shipunski to Cape Kronotski). From pairs of the strongest earthquakes that have occurred successively since 1737, the underthrust velocity for the southern zone is estimated at V = 6.6–7.1 cm/yr (from two pairs) and for the central zone, at V = 6.6 cm/yr. The creep portion of the underthrust amounts to 5–15% of the total velocity (the velocity of motion of the Pacific plate is 8 cm/yr).     

The Effect of a Mining and Smelting Plant on the Dynamics of Heavy Metals in Small River Basins in the Zone of Kursk–Belgorod Magnetic Anomaly

Heavy metal distribution under the impact of integrated mining and smelting works is studied. Factors affecting small river pollution are analyzed. The capability of small streams to create and preserve a barrier function in respect of water contamination with toxic heavy metals is estimated.     

Are extrusive rhyolites produced from permeable foam eruptions?

The permeable foam hypothesis is suggested by Eichelberger et al. (1986) to explain a major loss of water from rhyolitic magmas in the volcanic conduit. Evidence for the high-water content of the major portion of the magma is herein examined and rejected. Eichelberger's hypothesis does not take into account the large (2 orders of magnitude) viscosity change that would occur in the conduit as a result of water loss. It also requires that the permeable foam collapse and weld to form an obsidian that in thin section displays no evidence of the foam. An alternate hypothesis to explain the existence of small amounts of high water content rhyolite glasses in acid volcanoes is that rhyolite magmas are relatively dry (0.1–0.3% H₂O) and that water enters the magma from the environment to produce a water-rich selvage which then is kneaded into the body of the magma.     

Geokinematics of Central Europe: New insights from the CERGOP-2/Environment Project

The Central European Geodynamics Project CERGOP-2, funded by the European Union from 2003 to 2006 under the 5th Framework Programme, benefited from repeated measurements of the coordinates of epoch and permanent GPS stations of the Central European GPS Reference Network (CEGRN), starting in 1994. Here we report on the results of the systematic processing of available data up to 2005. The analysis has yielded velocities for some 60 sites, covering a variety of Central European tectonic provinces, from the Adria Indenter to the Tauern Window, the Dinarides, the Pannonian Basin, the Vrancea Seismic Zone and the Carpathian Mountains. The estimated velocities define kinematical patterns which outline, with varying spatial resolution depending on the station density and history, the present-day surface kinematics in Central Europe. Horizontal velocities are analyzed after removal from the ITRF2000 estimated velocities of a rigid rotation accounting for the mean motion of Europe: a 2.3 mm/year north–south oriented convergence rate between Adria and the Southern Alps that can be considered to be the present-day velocity of the Adria Indenter relative to the European Foreland. An eastward extrusion zone initiates at the Tauern Window. The lateral eastward flow towards the Pannonian Basin exhibits a gentle gradient from 1 to 1.5 mm/year immediately east of the Tauern Window to zero in the Pannonian Basin. This kinematic continuity implies that the Pannonian plate fragment recently suggested by seismic data does not require a specific Eulerian pole. On the southeastern boundary of the Adria microplate, we report a velocity drop from 4 to 4.5 mm/year motion near Matera to 1 mm/year north of the Dinarides, in the southwestern part of the Pannonian Basin. A positive velocity gradient as one moves south from West Ukraine across Rumania and Bulgaria is estimated to be 2 mm/year on a scale of 600–800 km, as if the crust were dragged by the counterclockwise rotation along the North Anatolian Fault Zone. This regime apparently does not interfere with the Vrancea Seismic Zone: earthquakes there are sufficiently deep (>100 km) that the brittle deformation at depth can be considered as decoupled from the creep at the surface. We conclude that models of the Quaternary tectonics of Central and Eastern Europe should not neglect the long wavelength, nearly aseismic deformation affecting the upper crust in the Romanian and Bulgarian regions.     

Spatiotemporal variations of actual evapotranspiration over the Lake Selin Co and surrounding small lakes (Tibetan Plateau) during 2003–2012

Actual evapotranspiration(ET_a) over the Tibetan Plateau(TP) is an important component of the water cycle,and greatly influences the water budgets of the TP lake basins.Quantitative estimation of ET_a within lake basins is fundamental to physically understanding ET_a changes,and thus will improve the understanding of the hydro logical processes and energy balance throughout the lake basins.In this study,the spatiotemporal dynamic changes of ET_α within the Lake Selin Co(the TP's largest lake) and its surrounding small lakes and land area during 2003-2012 are examined at the basin scale.This was carried out using the well-established Water and Energy Budget-based Distributed Hydrological Model(WEB-DHM) for the land area,the Penman method for the water area when unfrozen,and a simple sublimation estimation approach for the water area when frozen.The relationships between ET_a changes and controlling factors are also discussed.Results indicate that the simulated land ET_a from the WEB-DHM reasonably agrees with the estimated ET_a values from the nonlinear complementary relationship model using appropriately calibrated parameter values at a point scale.Land ET_a displayed a non-significant increase of 7.03 mm year~(-1),and largely depends on precipitation.For the water area,the combined effects of reduced wind speed and net radiation offset the effect of rising temperature and vapor pressure deficit,and contributed to a non-significant decrease in evaporation of 4.17 mm year~(-1).Sensitivity analysis shows that vapor pressure deficit and wind speed are the most sensitive variables to the changes of evaporation from the water area.     

Estimation of site effects using strong motion data of BYTNet array in Turkey

Simultaneous estimation of effects of source, propagation path, and local site amplification was carried out using observed strong motion records in a frequency range from 0.8 to 20 Hz for the purpose of empirical evaluation of the local site effects in different geological conditions in the northwestern part of Turkey. The analyzed data are S-wave portions of 162 accelerograms from 39 shallow events observed at 14 sites of BYTNet array. A spectral separation method was applied to the observed S-wave spectra. The solutions for source spectra, inelasticity factor of propagation path for S-waves (Q _s-value), and factor of site amplification at each site were obtained in a least squares sense. In the analysis, we assumed that the factor of the site amplification at a reference site is the same as that of theoretical amplification of S-waves to the soil model whose bottom layer has an S-wave velocity around 2.15 km/s. The estimated Q _s-value of the propagation path is modeled as Q _s(f)?=?87.4f^0.78. The estimated site amplifications are characterized into three groups. The sites in the first group belong to rock site with no dominant peaks at a frequency range of 2 to 10 Hz. The second group of hard soil sites is characterized with moderately dominant peaks at a frequency of 5 Hz. The last group for soft soil sites has common peaks at a frequency of 4 Hz with larger amplitudes than those in the hard soil group. We, then, compare the amplifications with average S-wave velocity in top 30 m of the shallow S-wave profiles and proposed linear empirical formula between them at each frequency. We, furthermore, inverted the observed amplification factors into S-wave velocity and Q _s-value profiles of the deep soil over the basement.     

Estimation of tree transpiration and response of tree conductance to meteorological variables in desert-oasis system of Northwest China

Transpirations of three dominated tree species, namely Mongol Scotch Pine (Pinus sylvestris var. mongolica Litvin), White elm (Ulmus pumila) and Gansu Poplar (Populus gansuensis Wang et Yang) in oasis shelter forest (Linze site) and of two dominated tree species, namely Euphrates Poplar (Populus euphratica Oliv.) and Russia olive (Elaeagnus angustifolia Linn.) in lowland desert (Erjinaqi site) have been estimated using measured sapflow in summer, autumn and winter, 2002 and in spring, 2003. An ENVIS System was used for each site to measure microclimate variables, soil moisture and sapflow every half an hour, and the study time scale is one day. In the 104 days of observation during the growing season at the Linze site, the average daily sapflow of Gansu Poplar is 9.93 L· d^-1, and the average transpiration per unit leaf area is 1.99 mm · d^-1. For White elm tree, the daily average sapflow is 4.08 L· d^-1, while the daily average transpiration per unit leaf area is 0.49 mm · d^-1. The values for Mongol Scotch Pine are 3.91 L· d^-1 and 0.25 mm · d^-1, respectively. In the total 73 days of observation during the growing season at the Erjinaqi site, the daily average sapflows of Russia olive and Euphrates Poplar are 12.1 and 20.97 L· d^-1, respectively, and the average transpirations per unit leaf area are 0.22 and 0.31 mm · d^-1, respectively. In the observation period of the growing season, tree conductances of Mongol Scotch Pine, White elm, Gansu Poplar or Russia olive show an exponential relationship with the daily average air temperature or vapour pressure deficit, but the relationship is not so obvious between tree conductance and global radiation. The transpiration process of each tree species is affected by all the observed four environmental variables. The response of tree conductance to different climatic factors changes with tree species. The effect of the same factor to the same tree species is also variable in different growing stages. The sapflow of every tree species is relatively large in later spring to early summer, and low in summer, and then reaches its largest value in later September. In the mid-November, the sapflow is relatively large, especially the deciduous tree species. This may be characteristic of the tree species in Arid Regions of Northwest China.     

Gravity and elevation changes at Askja,Iceland

Ground tilt measurements demonstrate that Askja is in a state of unrest, and that in the period 1988–1991 a maximum 48±3 rad tilt occurred down towards the centre of the caldera. This is consistent with 126 mm of deflation at the centre of the caldera with a 2.5–3.0 km depth to the source of deformation. The volume of the subsidence bowl is 6.2x10⁶ m³. When combined with high precision microgravity measurements, the overall change in sub-surface mass may be quantified. After correction for the observed elevation change using the free air gradient of gravity measured for each station, the total decrease in mass is estimated to be less than 10⁹ kg. A small residual ground inflation and net gravity increase in the southeastern part of the caldera may be caused by dyke intrusion in this region. The minimum dimensions of such an intrusion or complex of intrusions are 1 m width, up to 100 m deep and up to several hundred metres thick.     

Near-bed mass flux profiles in aeolian sand transport: high-resolution measurements in a wind tunnel

Vertical profiles of the streamwise mass flux of blown sand in the near-bed (< 17 mm) region are analysed from high-resolution measurements made using an optical sensor in a wind tunnel. This analysis is complemented by detailed measurements of mass flux and mean velocity profiles throughout the boundary layer depth (0·17 m) using passive, chambered sand traps of small dimensions and armoured thermal anemometers, respectively. The data permit a preliminary analysis of the relations between the observed forms of the profiles of near-bed fluid stress and horizontal mass flux within a carefully conditioned boundary layer. Profiles of mass flux density are found to be characterized by three regions of differing gradient with transitions at about 2 mm and 19 mm above the bed. The exponential decay of mass flux with height is confirmed for elevations above 19 mm, and when plotted as a function of u_*²/g (a parameter of mean vertical trajectory height in saltation), the gradient of mass flux in this region scales with the wake-corrected friction velocity (u), where u > 0·30 m s⁻¹. A separate near-bed region of more intense transport below 19 mm is identified which carries 80 per cent of the total mass flux. This region is evident in some previous field and wind tunnel data but not in profiles simulated by numerical models. Ventilated passive sand traps underestimate mass flux in this region by 37 per cent. At slow or moderate wind speeds a third significant region below 2 mm is observed. These regions are likely to be related to grain populations in successive saltation, low-energy ejections and intermittent bed contact, respectively. Optical measurements reveal locally high grain concentrations at some elevations below 5 mm; these heights scale with transport rate, mass flux gradient and wind speed. Copyright © 1999 John Wiley & Sons, Ltd.     

Recognition of environmentally caused variations in radon time series

Four techniques (pattern recognition, Fourier transformation, cross-correlation, and multiple linear regression) have been used to recognize nontectonic environmental factors affecting groundwater radon data collected by the Caltech automated geochemical network. Several factors, including air temperature, rainfall, water level, carbon dioxide concentration in the water, and barometric pressure, were found to correlate with radon level at some stations in the network. For example, approximately 60% of the variance in radon signal at one site studied (Pacoima) could be accounted for by nontectonic influences. However, it was not possible to correlate all the observed changes in radon concentration with nontectonic environmental variables. Some of the observed radon anomalies did appear to be related to either regional changes in tectonic strain or some individual earthquakes. A model is developed to account for spike-like increases in radon caused by carbon dioxide emission.     

The magnetic observatory of Pendeli (Athens,Greece)

Summary A new magnetic observatory, named the Magnetic Observatory of Pendeli, was established and put in operation in 1958 near Athens (Greece). This Observatory was organized by and belongs to the Greek «Institute for Geology and Subsurface Research». The geographical position of the Pendeli Observatory is given by =38° 02.8, =23°51.8 andh=495 m (above sea level). The gemagnetic coordinates of the same are =36°.2, =102°.0. The Observatory is situated near Pendeli Mt. (18 km NEE of Athens). The site of the Observatory consists of marmor underlain by mica schists, both magnetically inactive.The building of the Observatory is constructed of stone and its roof made of tiles. The magnetograph room is in the underground of the building. The magnetic and thermal conditions in the variometer room are satisfactory enough.The variometers of the Observatory forH andZ are of the magnetic balance type. ForD a fibre suspension declinometer is used. The scale values of the variometers are _H =7.2 /mm, _D =1.0/mm (7.6 /mm), _Z =11.5 /mm. The speed of recording amounts to 12 mm/h and the width of the record is 9 cm.The Pendeli Observatory has been operating since April 1958. The record is changed every day. The scale and base-line values are determined every 10 days. The room for the absolute measurements is found in the ground floor of the building. A field magnetic theodolite is being used in the absolute measurements.D is measured with two magnets in a fibre suspension declinometer.H is measured by means of the deflection oscillation method and with a QHM as well. For the measurement ofI an earthinductor is available. The values of the magnetic elements are properly corrected in order to represent the external normal field.     

Background aspects of lake restoration: water balance,heavy metal content,phosphorus homeostasis

Perturbations of the lake water balance, inputs of heavy metals to lakes, and intensifying fertilization of lakes through input and accumulation of phosphorus—these three classes of phenomena are among the more important background processes in lake restoration. Lake restoration consists of a series of measures animed at producing a homeostatic response of a lake system to external perturbations. The success of its implementation is affected by the morphometric and edaphic parameters of different types of lakes. The relationship between the volume (V) and mean-depth of fresh-water lakes indicates a trend of . Glacial lakes occuring on or near crystalline shields have relatively shallow depths, whereas volcanic lakes, rift valley and deep valley lakes have relatively greater depths for the same volume. For saline lakes (21 lakes, V>1 km³) that undergo cycles of expansion and shrinkage, the V to relationship is closer to power 1. Water residence times (τ) of small and big fresh-water lakes show a trend of τ approximately linear in or τ∝V^0.3. Volcanic lakes and Maare have longer residence times in comparison to other lakes of similar volumes. For the major inorganic chemical species and heavy metals, the regulatory upper-limit concentrations in drinking water in the USA and EEC are from several times to more than 100 times higher than their concentrations in a global mean river water. Only three elements (Fe, P, and Al) occur in river water at concentrations approaching such upper-limit recommendations. Rates of accumulation of phosphorus in lake water and sediments, computed as the difference between input and ouflow removal rates for 23 fresh-water lakes, are generally lower for lakes of longer water residence time. The rate of accumulation is a measure of homeostatic response of the lake system to input load: it is equivalent to the rate of all the removal processes needed to maintain phosphorus concentration in lake water at a steady state.     

Mining-induced microseismic event location errors: Accuracy and precision of two location systems

The accuracy and precision of microseismic event locations were measured, analyzed, and compared for two types of location systems: anolog and digital. In the first system, relative times of first arrival were estimated from analog signals using automated hardware circuitry; station positions were estimated from mine map coordinates; and event locations were determined using the BLD (Blake, Leighton, and Duvall) direct solution method. In the second system, arrival times were manually measured during interactive displays of digital waveforms; station coordinates were surveyed; and the SW-GBM (Salamon and Wiebols; Godson, Bridges, and McKavanagh) direct basis function was used to solve for locations. Both systems assume constant isotropic seismic velocity of slightly different signals data sets, calibration blast signals with known source site and origin time, and microseismic event signals, were recorded by each location system employing the same array of high-frequency (5 kHz) accelerometers with 150 m maximum dimension. The calibration blast tests indicated a location precision of ±2 m and accuracy of ±10 m for the analog system. Location precision and accuracy for the digital system measured ±1 m and ±8 m, respectively. Numerical experiments were used to assess the contributions of errors in velocity, arrival times, and station positions on the location accuracy and precision for each system. Measured and estimated errors appropriate to each system for microseismic events were simulated in computing source locations for comparison with exact synthetic event locations. Discrepancy vectors between exact locations and locations calculated with known data errors averaged 7.7 and 1.4 m for the analog and digital systems, respectively. These averages are probably more representative of the location precision of microseismic events, since the calibration blast tests produce impulsive seismic arrivals resulting in smaller arrival-time pick errors in the analog system. For both systems, location accuracy is limited by inadequate modeling of the velocity structure. Consequently, when isotropic velocity models are used in the travel-time inversions, the increased effort expended with the digital location system does not, for the particular systems studied, result in increased accuracy.     

A daily rainfall occurrence process

A model for the periodic (annual cycle), discrete rainfall occurrence process is presented. Using this model the probabilistic properties of the process in -day intervals can be investigated. In such an interval the rainfall occurrence process is approximated by some stationary processa _t ,tIN. The processa _t ,tIN is described by the distributions of the lengths of wet and dry sequences. It is assumed that the lengths of successive wet and dry sequences are independent. For this process the distribution of the number of wet days in -day intervals is calculated. The model is fitted to 50-year rainfall data from Wroclaw, Poland. Rainfall amounts of 0.1, 1.0 and 2.0 mm are considered as thresholds defining a wet day. To estimate the distribution of the length of wet and dry sequences the family of Pascal distribution is chosen.     

The transport of water in subduction zones

The transport of water from subducting crust into the mantle is mainly dictated by the stability of hydrous minerals in subduction zones. The thermal structure of subduction zones is a key to dehydration of the subducting crust at different depths. Oceanic subduction zones show a large variation in the geotherm, but seismicity and arc volcanism are only prominent in cold subduction zones where geothermal gradients are low. In contrast, continental subduction zones have low geothermal gradients, resulting in metamorphism in cold subduction zones and the absence of arc volcanism during subduction. In very cold subduction zone where the geothermal gradient is very low(?5?C/km), lawsonite may carry water into great depths of ?300 km. In the hot subduction zone where the geothermal gradient is high(25?C/km), the subducting crust dehydrates significantly at shallow depths and may partially melt at depths of 80 km to form felsic melts, into which water is highly dissolved. In this case, only a minor amount of water can be transported into great depths. A number of intermediate modes are present between these two end-member dehydration modes, making subduction-zone dehydration various. Low-T/low-P hydrous minerals are not stable in warm subduction zones with increasing subduction depths and thus break down at forearc depths of ?60–80 km to release large amounts of water. In contrast, the low-T/low-P hydrous minerals are replaced by low-T/high-P hydrous minerals in cold subduction zones with increasing subduction depths, allowing the water to be transported to subarc depths of 80–160 km. In either case, dehydration reactions not only trigger seismicity in the subducting crust but also cause hydration of the mantle wedge. Nevertheless, there are still minor amounts of water to be transported by ultrahigh-pressure hydrous minerals and nominally anhydrous minerals into the deeper mantle. The mantle wedge overlying the subducting slab does not partially melt upon water influx for volcanic arc magmatism, but it is hydrated at first with the lowest temperature at the slab-mantle interface, several hundreds of degree lower than the wet solidus of hydrated peridotites. The hydrated peridotites may undergo partial melting upon heating at a later time. Therefore, the water flux from the subducting crust into the overlying mantle wedge does not trigger the volcanic arc magmatism immediately.     

Preliminary15N studies on atmospheric nitrogenous trace gases

Preliminary nitrogen isotope data for ammonia from animal urine, fuel combustion, fertilizer use and fertilizer factories have been measured or estimated. It turns out that direct nitrogen isotope measurements of atmospheric ammonia at Jülich are in the expected range calculated from the ranges of different sources. For deposition of atmospheric ammonium in Jülich-rain a depletion in¹⁵N with respect to atmospheric ammonia has been found which is explained by isotope fractionations during rainout and washout. In correspondence with this fractionation model are nitrogen isotope data of rain-ammonium from coastal areas, which are enriched in¹⁵N due to the fact that sea water acts as a sink for atmospheric ammonia.For Jülich rain-nitrate a pronounced seasonal trend has been detected with lower¹⁵N data in spring and summer than in autumn and winter. This trend is interpreted by different nitrogen isotope data of anthropogenic and natural nitric oxides which have been measured or estimated from isotope fractionation effects during nitrification and denitrification reactions in soils. It should be possible to get better global estimations for anthropogenic and natural nitric oxides from nitrogen isotope measurements.     

Phosphorus Forms and Their Changes in Water of Three Subalpine Lakes in Italy

Total phosphorus and its main forms: dissolved mineral, dissolved organic, particulate organic and particulate mineral in the vertical water column of three subalpine lakes of various types in Italy, has been estimated during the winter-vernal season. The range of variation in the phosphorus content in these waters was as follows: total phosphorus 16 ± 2860 μg/1 PO₄, dissolved mineral phosphorus 4 ± 1040 μg/l PO₄, dissolved organic phosphorus 1 ± 160μg/l PO₄, particulate organic phosphorus 0 ± 290 μg/l PO₄ and particulate mineral phosphorus 1 ± 100 μg/l PO₄, Generally the content of total phosphorus and dissolved mineral phosphorus (phosphates) increased with the degree of eutrophy with the depth and with the progress of the vernal season towards the summer stagnation time. The amount of phosphates increased in water with the depletion of oxygen, both in the verical water column and with the progress of stagnation time. The amounts of dissolved organic phosphorus decressed with the depth of the vertical water column whereas the dissolved mineral phosphorus increased. The development of the particulate organic phosphorus stratification in the vertical water column was clearly visible in the eutrophic lake. The quantities of total phosphorus and its main component, dissolved mineral phosphorus, decreased evidently from January to May in all three lakes, mostly in the eutrophic lake. The reason of this decrease is sorption by lake sediments and to a certain degree sedimentation of phosphorus sorbed by ferric hydroxide. The increase of dissolved mineral phosphorus and that of total phosphorus in the vertical water column and with the progress of summer stagnation had as a reason the liberation of phosphorus from sediments, and not so much decomposition of sedimentating plankton or dissolved organic phosphorus. The share of single (mean) values of phosphorus forms in the total phosphorus was as follows: In the oligotrophic lake the share of particulate mineral phosphorus was extremely high in March (21% of the total), probably because of the inflow of the melting waters from the drainage area. The development of vertical stratification in waters of three subalpine Italian lakes at the end of the vernal season (May) indicates the quantitative prevailing of dissolved mineral phosphorus with its increase with the depth and domination of dissolved organic phosphorus in the trophogenic zone.