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     

Hyperspectral remote sensing analysis of short rotation woody crops grown with controlled nutrient and irrigation treatments

Hyperspectral remote sensing research was conducted to document the biophysical and biochemical characteristics of controlled forest plots subjected to various nutrient and irrigation treatments. The experimental plots were located on the Savannah River Site near Aiken, SC. AISA hyperspectral imagery were analysed using three approaches, including: (1) normalized difference vegetation index based simple linear regression (NSLR), (2) partial least squares regression (PLSR) and (3) machine-learning regression trees (MLRT) to predict the biophysical and biochemical characteristics of the crops (leaf area index, stem biomass and five leaf nutrients concentrations). The calibration and cross-validation results were compared between the three techniques. The PLSR approach generally resulted in good predictive performance. The MLRT approach appeared to be a useful method to predict characteristics in a complex environment (i.e. many tree species and numerous fertilization and/or irrigation treatments) due to its powerful adaptability.     

Effective sea ice area based on a thickness threshold

As Arctic sea ice declines in response to climate change, a shift from thick multiyear ice to a thinner ice cover is occurring. With this transition, ice thicknesses approach a threshold below which ice no longer insulates the atmosphere from oceanic surface fluxes. While this is well known, there are no estimates of the magnitude of this threshold, nor of the proportion of sea ice area that is below this threshold as ice thins. We determine this threshold by simulating the atmospheric response to varying thicknesses, ranging from 0.0 to 2.0 m and determine that threshold to be 0.40–0.50 m. The resulting “effective” ice area is 4–14% lower than reported total ice area, as 0.39–0.97 × 10⁶ km² of the total ice area falls below the threshold throughout the twentieth century, including during notable ice minima. The atmosphere above large non-insulating ice-covered regions is susceptible to more than 2 °C of warming despite ice presence. Observed mean Arctic Ocean ice thickness is projected to fall below this threshold as early as the mid-2020s. Studies on ocean–atmosphere interactions in relation to sea ice area should focus on this insulating sea ice area, where ice is at least 0.40–0.50 m thick, and treat ice regions below 0.40–0.50 m thickness with caution.

     

Bifurcation of the Pacific North Equatorial Current in a wind-driven model: response to climatological winds

The sensitivity of the bifurcation of the North Equatorial Current in the Pacific to different wind products is investigated. Variations of the bifurcation latitude with season is simulated in a purely wind-driven model and is found to be in agreement with recent observations. The seasonal cycle is nearly independent of the wind climatology, but the annual average latitude depends on the wind stress curl. It is also shown that in the upper ocean, the poleward shift in bifurcation latitude with depth is realistic in our simple model. This implies that given a stratification close to the observed, it is primarily the wind forcing that determines the location of the bifurcation and its seasonal variation.     

Petroleum-related hydrocarbons in deep and subsurface sediments from South-Western Barents Sea

Subsurface sediments from a pockmark area in South-Western Barents Sea have been earlier found to contain elevated levels of petroleum-related polycyclic aromatic hydrocarbons. This work describes a comprehensive analysis of various biomarkers, including the highly source-specific hopanes, in a 4.5?m long gravity core from the same area, together with subsurface sediment samples from other areas in the region without pockmarks present ("background samples"). A clear difference between the pockmark gravity core and the background sediment cores was found, both with regard to genesis and the level of transformation of organic matter. A number of indicator parameters, such as methylphenanthrene index (MPI-1), point towards a significantly higher maturity of hydrocarbons in the pockmark core throughout its length as compared to the other sampled locations. Higher contents of microbial hopanoids (hopenes) may indicate the former presence of petroleum. These findings confirm the hypothesis of a natural hydrocarbon source in the deeper strata present in the studied location with pockmarks.     

The first record of a neonatal ornithopod dinosaur from Gondwana

Discrete post-embryonic teeth and bone fragments have been recovered from the matrix with the holotype skeleton (MPM‐10001) of the ornithopod dinosaur, Talenkauen santacrucensis Novas et al., 2004 (Upper Cretaceous, Argentina). The minute tooth crowns are 1 mm apicobasally tall and 1.7 mm mesodistally wide. The crowns are symmetrical and have a centrally located primary ridge on the lingual surface. Secondary ridges lead to five marginal denticles on both teeth. The tooth morphology is consistent with dentary teeth in euiguanodontids. There is no evidence of transport, suggesting that the material is autochthonous with respect to the adult body block of T. santacrucensis (MPM‐10001). Steeply inclined wear facets on the lingual surface and associated microstriae support the conclusion that the minute teeth were from a post-embryonic euiguanodontid dinosaur rather than early stage replacement teeth. The morphology, size, and wear of the teeth and small bone fragments found in the body block of MPM‐10001 suggest that this material belongs to a neonatal T. santacrucensis. This is the first record of neonatal ornithopod remains from Gondwana.     

Reionization and the Cosmic Dawn with the Square Kilometre Array

The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21 cm line from the earliest phases of star and galaxy formation in the Universe. This 21 cm signal provides a new and unique window both on the time of the formation of the first stars and accreting black holes and the subsequent period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.     

Berm formation and dynamics on a gently sloping beach; the effect of water level and swash overtopping

Berm formation and morphological development of the beach face have been observed during a neap–neap tidal cycle on the gently sloping and accreting beach at Vejers, Denmark. During the field campaign, an intertidal bar migrated onshore and stabilized as a berm on the foreshore. A new intertidal bar occurred on the lower beach face, migrated onshore on the rising tide and finally merged with the pre‐existing berm. As the tide continued to rise, the new berm translated further onshore as an intertidal bar to the uppermost part of the foreshore. The sediment transport during the berm transition was onshore directed in the upper swash and offshore directed in the lower swash. This berm development can be described through both the neap‐berm, ridge‐and‐runnel and berm‐ridge development concepts proposed by Hine (Sedimentology 1979; 26: 333–351), and all three stages were observed during only three tidal cycles. The main factors controlling this fast transformation were the gentle slope of the cross‐shore profile, rapid water level translation rates, substantial swash overtopping of the berm, and low infiltration rates. Despite the onshore migration of intertidal bars and berm formation, no net foreshore accretion took place during the field campaign. This was largely due to the formation of rip channels with strong rip currents cutting through the intertidal bars and the berm, which acted as a sediment drain in the profile. Copyright © 2009 John Wiley & Sons, Ltd.     

The spatial ecology of deer–vehicle collisions

Deer–vehicle collisions (DVCs) are a constant problem in North America that cost millions of dollars of damage and loss of human life. This study investigated the spatial relationship of DVCs for a county in western Indiana using nearest neighbor analysis, chi-square test of independence, and a landscape metric: the percentage of like adjacencies. The study was able to show that DVCs are not randomly located, and that habitat type and structure may play an important role in the location of DVCs.     

Object and Pixel-Based Reservoir Modeling of a Braided Fluvial Reservoir

To assess differences between object and pixel-based reservoir modeling techniques, ten realizations of a UK Continental Shelf braided fluvial reservoir were produced using Boolean Simulation (BS) and Sequential Indicator Simulation (SIS). Various sensitivities associated with geological input data as well as with technique-specific modeling parameters were analyzed for both techniques. The resulting realizations from the object-based and pixel-based modeling efforts were assessed by visual inspection and by evaluation of the values and ranges of the single-phase effective permeability tensors, obtained through upscaling. The BS method performed well for the modeling of two types of fluvial channels, yielding well-confined channels, but failed to represent the complex interaction of these with sheetflood and other deposits present in the reservoir. SIS gave less confined channels and had great difficulty in representing the large-scale geometries of one type of channel while maintaining its appropriate proportions. Adding an SIS background to the Boolean channels, as opposed to a Boolean background, resulted in an improved distribution of sheetflood bodies. The permeability results indicated that the SIS method yielded models with much higher horizontal permeability values (20–100%) and lower horizontal anisotropy than the BS versions. By widening the channel distribution and increasing the range of azimuths, however, the BS-produced models gave results approaching the SIS behavior. For this reservoir, we chose to combine the two methods by using object-based channels and a pixel-based heterogeneous background, resulting in moderate permeability and anisotropy levels.