Integrating lysimeter drainage and eddy covariance flux measurements in a groundwater recharge model

Abstract

Field-scale water balance is difficult to characterize because controls exerted by soils and vegetation are mostly inferred from local-scale measurements with relatively small support volumes. Eddy covariance flux and lysimeters have been used to infer and evaluate field-scale water balances because they have larger footprint areas than local soil moisture measurements. This study quantifies heterogeneity of soil deep drainage (D) in four 12.5-m² repacked lysimeters, compares evapotranspiration from eddy covariance (ET_EC) and mass balance residuals of lysimeters (ETwb_Lys), and models D to estimate groundwater recharge. Variation in measured D was attributed to redirection of snowmelt infiltration and differences in lysimeter hydraulic properties caused by surface soil treatment. During the growing seasons of 2010, 2011 and 2012, ETwb_Lys (278, 289 and 269 mm, respectively) was in good agreement with ET_EC (298, 301 and 335 mm). Annual recharge estimated from modelled D was 486, 624 and 613 mm for three calendar years 2010, 2011 and 2012, respectively. In summary, lysimeter D and ET_EC can be integrated to estimate and model groundwater recharge.

Editor D. Koutsoyiannis     

Geological and SHRIMP II U-Pb constraints on the age and origin of the Breves Cu-Au-(W-Bi-Sn) deposit,Carajás,Brazil

The Breves deposit in the Carajás Copper-Gold Belt, Brazil, a member of the Cu-Au-(W-Bi-Sn) group of deposits, contains about 50 Mt of 1.22% Cu, 0.75 g/t Au, 2.4 g/t Ag, 1,200 g/t W, 70 g/t Sn, 175 g/t Mo and 75 g/t Bi. It is hosted by sandstones and siltstones of the Águas Claras Formation (minimum age of 2,681±5 Ma) in the roof zone of a complex, highly altered granite intrusion. The mineralisation is disseminated in a greisenized zone, resulting from alteration of probable monzogranites and syenogranites. The ore-bearing greisen contains abundant xenomorphic quartz in association with Fe-chlorite and muscovite. The gangue assemblage also includes fluorite, tourmaline, and minor amounts of monazite, xenotime, chlorapatite, thorite, zircon, calcite, siderite and bastnäesite. Copper mineralisation is dominated by chalcopyrite associated with pyrite, arsenopyrite, pyrrhotite and molybdenite. Gold particles, in equilibrium with native bismuth, are common as inclusions in chalcopyrite. The greisen contains sub-economic concentrations of tungsten and niobium that are related to the presence of ferberite, qitianlingite and Nb-rutile. SHRIMP II zircon dating of the host granites gives ²⁰⁷Pb/²⁰⁶Pb ages of 1,878±8 and 1,880±9 Ma for two phases, and a combined age of 1,879±6 Ma. SHRIMP II dating of monazite and xenotime grains in late- to post-mineralisation veins gives a combined ²⁰⁷Pb/²⁰⁶Pb age of 1,872±7 Ma, indistinguishable from the ages of the granites. This provides a genetic connection between the Breves deposit and the ca. 1.88 Ga A-type granite magmatism that typifies the Carajás Belt as part of a much larger, intracratonic magmatic province that extends over much of the Amazonian Craton. The recognition of this association has exploration implications, not only for the geophysical signature of the granite roof zones, but also for likely geochemical dispersion around the deposits of this type.Editorial handling: G. Beaudoin     

Biomass expansion factor and root-to-shoot ratio for Pinus in Brazil

The Biomass Expansion Factor (BEF) and the Root-to-Shoot Ratio (R) are variables used to quantify carbon stock in forests. They are often considered as constant or species/area specific values in most studies. This study aimed at showing tree size and age dependence upon BEF and R and proposed equations to improve forest biomass and carbon stock. Data from 70 sample Pinus spp. grown in southern Brazil trees in different diameter classes and ages were used to demonstrate the correlation between BEF and R, and forest inventory data, such as DBH, tree height and age. Total dry biomass, carbon stock and CO₂ equivalent were simulated using the IPCC default values of BEF and R, corresponding average calculated from data used in this study, as well as the values estimated by regression equations. The mean values of BEF and R calculated in this study were 1.47 and 0.17, respectively. The relationship between BEF and R and the tree measurement variables were inversely related with negative exponential behavior. Simulations indicated that use of fixed values of BEF and R, either IPCC default or current average data, may lead to unreliable estimates of carbon stock inventories and CDM projects. It was concluded that accounting for the variations in BEF and R and using regression equations to relate them to DBH, tree height and age, is fundamental in obtaining reliable estimates of forest tree biomass, carbon sink and CO₂ equivalent.     

U---Pb and Rb---Sr systematics in a polyorogenic segment of the Precambrian shield, central southern Norway

U---Pb analyses of zircon, monazite and sphene as well as Rb---Sr analyses of whole rocks and minerals have been carried out in an attempt to elucidate the evolution of the Precambrian shield underlying the Caledonian Jotun-nappe in central southern Norway. The earliest event recognized in the area is a high grade metamorphism at 1518 ± 17 m.y. which followed intense magma formation, igneous activity, erosion and sedimentation over a period of maximum 300 m.y. The Sveconorwegian cycle is characterized by the intrusion of two anatectic magma generations at 1014 ± 35 m.y. and 930±10 m.y. that disturbed and partly reset the isotopic systems of the country rocks. A post-magmatic fracturing stage at 875 m.y. concludes this cycle. Rb---Sr mineral systems were not completely equilibrated in the undeformed shield during the Caledonian event in contrast to new grown minerals in strongly deformed overlying Lower Paleozoic sediments which record the main deformation at 384±18 m.y. Zircon lower intercept ages ranging between 330–370 m.y. show that lead loss of zircon in rocks subjected to low grade metamorphic conditions may be the result of such processes as annealing, alteration, dilation or a combination of them.     

Characteristics of chaotic pulse trains generated by lightning flashes

Pulse trains of irregular pulse characteristics and erratic amplitudes are observed in lightning flashes. These are termed “chaotic pulse trains” in this study. The average width of the individual pulses of these pulse trains is in the range of a few microseconds, with the lower limit of individual measurements extending into the sub-microsecond region. The typical pulse separation lies in the range of 2–. The pulse duration shows a consistent distribution with the most probable value at 400–. Chaotic pulse trains occur most often immediately prior to subsequent strokes. However, they occur also a few milliseconds to several tens of milliseconds prior to return strokes and sometimes do not precede return strokes at all. Twenty-six percent of the total subsequent strokes recorded in this study were preceded by chaotic pulse trains. A few of these pulse trains have occurred in a sequence of several tens of millisecond intervals in between return strokes. We have also found them to occur within very intense cloud flashes. These chaotic pulse trains have features that are different from those of other known pulse trains reported in the lightning literature. Hence, they may be caused by a physical process that has not yet been identified.     

Real-time cycle slip detection in triple-frequency GNSS

The modernization of the global positioning system and the advent of the European project Galileo will lead to a multifrequency global navigation satellite system (GNSS). The presence of new frequencies introduces more degrees of freedom in the GNSS data combination. We define linear combinations of GNSS observations with the aim to detect and correct cycle slips in real time. In particular, the detection is based on five geometry-free linear combinations used in three cascading steps. Most of the jumps are detected in the first step using three minimum-noise combinations of phase and code observations. The remaining jumps with very small amplitude are detected in the other two steps by means of two-tailored linear combinations of phase observations. Once the epoch of the slip has been detected, its amplitude is estimated using other linear combinations of phase observations. These combinations are defined with the aim of discriminating between the possible combinations of jump amplitudes in the three carriers. The method has been tested on simulated data and 1-second triple-frequency undifferenced GPS data coming from a friendly multipath environment. Results show that the proposed method is able to detect and repair all combinations of cycle slips in the three carriers.     

Segmentation and kinematics of the North America‐Caribbean plate boundary offshore Hispaniola

We explored the submarine portions of the Enriquillo–Plantain Garden Fault zone (EPGFZ) and the Septentrional–Oriente Fault zone (SOFZ) along the Northern Caribbean plate boundary using high‐resolution multibeam echo‐sounding and shallow seismic reflection. The bathymetric data shed light on poorly documented or previously unknown submarine fault zones running over 200 km between Haiti and Jamaica (EPGFZ) and 300 km between the Dominican Republic and Cuba (SOFZ). The primary plate‐boundary structures are a series of strike‐slip fault segments associated with pressure ridges, restraining bends, step overs and dogleg offsets indicating very active tectonics. Several distinct segments 50–100 km long cut across pre‐existing structures inherited from former tectonic regimes or bypass recent morphologies formed under the current strike‐slip regime. Along the most recent trace of the SOFZ, we measured a strike‐slip offset of 16.5 km, which indicates steady activity for the past ~1.8 Ma if its current GPS‐derived motion of 9.8 ± 2 mm a^?1 has remained stable during the entire Quaternary.     

Spatial and temporal variability of particulate matter in the benthic boundary layer at the N.W. European Continental Margin (Goban Spur)

Near bottom water samples and sediments were taken during five cruises to 6 stations forming a transect across the N.W. European Continental Margin at Goban Spur. Flow velocity spot measurements in the benthic boundary layer (BBL) always increased from the shelf to the upper slope (1470 m) from 5 to 9 cm s⁻¹ in spring/summer and from 15 to 37 cm s⁻¹ in autumn/winter. Decreasing values were detected at the lower slope (2000 m) and the lowest values of ca. 2 cm s⁻¹ at the continental rise at 4500 m water depth. Long term measurements with a benthic lander at 1470 m show that currents have a tidal component and reach maximum velocities up to 20 cm s⁻¹, sufficiently high periodically to resuspend and transport phytodetritus. During these long-term observations, currents were always weaker in spring/summer than in autumn/winter. Critical shear velocities of shelf/slope sediments increased with depth from 0.5 to 1.7 cm s⁻¹ and major resuspension events and Intermediate Nepheloid Layers (INLs) should occur around 1000 m. Chloroplastic Pigment Equivalents (CPE) ranged from 0.0 to 0.21 μg dm⁻³, Particulate Organic Carbon (POC) from 12 to 141 μg dm⁻³ and Total Particulate Matter (TPM) from 0.2 to 10.0 mg dm⁻³. Aggregates in the BBL occurred with a median diameter of 152 to 468 μm. Data on suspended particulate matter in the near-bottom waters showed that hydrodynamic sorting within the particulate organic fraction occurred. Phytodetritus was packaged in relatively large aggregates and contributed little to the total organic carbon pool in nearbottom waters (CPE/POC ca.0.2%). The main organic fraction has low settling velocities and high residence times within the benthic boundary layer. As POC was not concentrated in the near bed region the degree to which carbon is accessible to the benthic community depends on aggregate formation, subsequent settling and/or biodeposition of the POC. Close to the sea bed downslope transport may dominate. Under flow conditions high enough to resuspend fresh phythodetritus from sediments at the productive shelf edge, this could be transported to 1500 m (Goban Spur) or abyssal depth (Canyon site between Meriadzek and Goban Spur) within 21 days.     

Shale basins,sulfur-deficient ore brines and the formation of exhalative base metal deposits

The massive sulfide deposits of the Iberian Pyrite Belt are interbedded with felsic volcanic rocks and shale, and underlain by several thousand meters of siliciclastic sedimentary rocks known as the PQ Group. Isotope geochemistry and regional geology are both consistent with equilibration of the ore-forming fluids with the PQ Group, prior to ore deposition near the former seafloor. The average Cu:Zn:Pb ratio of the PQ Group rocks (ca. 26:55:19) is similar to the weighted average of all the massive sulfide orebodies combined (ca. 25:52:23).The genetic relationship between massive sulfide deposits and a siliciclastic sedimentary metal source is explained here by a thermodynamic model, proposing that mildly reducing redox conditions imposed by equilibration with the sedimentary rocks are most critical for the formation of an effective ore-forming fluid. Relatively metal-rich but organic-poor pyrite-bearing shale undergoing dewatering of saline pore fluids is an effective source for the generation of sulfur-deficient but relatively iron and base metal-rich brines. Thus, we propose that the giant deposits of the Iberian Pyrite Belt owe their existence not to exceptionally metal-enriched (e.g., magmatic) fluids, but to the existence of a fairly ordinary but large metal source in reactive siliciclastic sediments, combined with an underlying igneous heat source and a particularly efficient mechanism of sulfide precipitation by mixing with H₂S-rich fluids at or near the seafloor.Essentially similar mineral equilibria are imposed when saline fluids are buffered by typical continental basement rocks. Leaching of retrograde minerals and possibly residual salts from their magmatic or metamorphic prehistory is expected to generate similar, variably metal-rich but relatively sulfide-deficient fluids. Thus, the existence of mildly reducing rocks can be the dominant chemical control in the source of fluids generating many volcanogenic, Irish-type or sedex deposits, many of which are known to precipitate their metal load in response to biogenic sulfide addition at the ore deposition site.     

Rutile R632 – A New Natural Reference Material for U‐Pb and Zr Determination

A new natural rutile reference material is presented, suitable for U‐Pb dating and Zr‐in‐rutile thermometry by microbeam methods. U‐Pb dating of rutile R632 using laser ablation ICP‐MS with both magnetic sector field and quadrupole instruments as well as isotope dilution‐thermal ionisation mass spectrometry yielded a concordia age of 496 ± 2 Ma. The high U content (> 300 μg g^?1) enabled measurement of high‐precision U‐Pb ages despite its young age. The sample was found to have a Zr content of 4294 ± 196 μg g^?1, which makes it an excellent complementary reference material for Zr‐in‐rutile thermometry. Individual rutile grains have homogeneous compositions of a number of other trace elements including V, Cr, Fe, Nb, Mo, Sn, Sb, Hf, Ta and W. This newly characterised material significantly expands the range of available rutile reference materials relevant for age and temperature determinations.