Impact of artificial subsurface drainage on groundwater travel times and baseflow discharge in an agricultural watershed,Iowa (USA)

Pollutant delivery through artificial subsurface drainage networks to streams is an important transport mechanism, yet the impact of drainage tiles on groundwater hydrology at the watershed scale has not been well documented. In this study, we developed a two‐dimensional, steady‐state groundwater flow model for a representative Iowa agricultural watershed to simulate the impact of tile drainage density and incision depth on groundwater travel times and proportion of baseflow contributed by tile drains. Varying tile drainage density from 0 to 0.0038 m^?1, while maintaining a constant tile incision depth at 1.2 m, resulted in the mean groundwater travel time to decrease exponentially from 40 years to 19 years and increased the tile contribution to baseflow from 0% to an upper bound of 37%. In contrast, varying tile depths from 0.3 to 2.7 m, while maintaining a constant tile drainage density of 0.0038 m^?1, caused mean travel times to decrease linearly from 22 to 18 years and increased the tile contribution to baseflow from 30% to 54% in a near‐linear manner. The decrease in the mean travel time was attributed to decrease in the saturated thickness of the aquifer with increasing drainage density and incision depth. Study results indicate that tile drainage affects fundamental watershed characteristics and should be taken into consideration when evaluating water and nitrate export from agricultural regions. Copyright © 2011 John Wiley & Sons, Ltd.     

Late Pleistocene glaciation of the Mt Giluwe volcano,Papua New Guinea

The Mt Giluwe shield volcano was the largest area glaciated in Papua New Guinea during the Pleistocene. Despite minimal cooling of the sea surface during the last glacial maximum, glaciers reached elevations as low as 3200 m. To investigate changes in the extent of ice through time we have re-mapped evidence for glaciation on the southwest flank of Mt Giluwe. We find that an ice cap has formed on the flanks of the mountain on at least three, and probably four, separate occasions. To constrain the ages of these glaciations we present 39 new cosmogenic ³⁶Cl exposure ages complemented by new radiocarbon dates. Direct dating of the moraines identifies that the maximum extent of glaciation on the mountain was not during the last glacial maximum as previously thought. In conjunction with existing potassium/argon and radiocarbon dating, we recognise four distinct glacial periods between 293–306 ka (Gogon Glaciation), 136–158 ka (Mengane Glaciation), centred at 62 ka (Komia Glaciation) and from >20.3–11.5 ka (Tongo Glaciation). The temperature difference relative to the present during the Tongo Glaciation is likely to be of the order of at least 5 °C which is a minimum difference for the previous glaciations. During the Tongo Glaciation, ice was briefly at its maximum for less than 1000 years, but stayed near maximum levels for nearly 4000 years, until about 15.4 ka. Over the next 4000 years there was more rapid retreat with ice free conditions by the early Holocene.     

On the modelling of sand bedforms using the semivariogram

This study shows the usefulness of the semivariogram for modelling sand ripples created by water flows of varied flow intensity. A combination of two mathematical functions is fitted to each sample semivariogram, that is an exponential (or stochastic) component and a periodic component. The parameters of each of these components have direct physical meaning. A non-dimensional ratio combining the two parameters of the exponential model is interpreted as a regularity index (which increases with the degree of regularity of bedform arrangement). This regularity index is inversely related to the Froude number of the flow. The non-dimensional wavelength, estimated from the dominant periodic function, is also inversely and closely related to the Froude number. The wave height, accurately estimated from properties of the two fitted components, is a direct function of flow velocity and is also proportional to the standard deviation of bed elevations. The bedform shape introduces a considerable discrepancy between the generally assumed normal frequency distribution and the empirical distributions of bed height. The series of bed elevations are generally characterized by a mixture of normal distributions having the same variance but different means. The calculation of a covariance assuming a constant and single mean (as in spectral analysis) can therefore be misleading and the problem may be avoided by using the semivariogram.     

Relationship between waterfowl and American coot abundance with submersed macrophytic vegetation in Currituck Sound,North Carolina

Waterfowl and American coot data collected from aerial censusing on Currituck Sound were compared during three periods: before Eurasian watermilfoil (Myriophyllum spicatum) colonization (1961–1965); during Eurasian watermilfoil infestation (1968–1977); and following the near absence of submersed macrophytic vegetation after Eurasian watermilfoil decline (1981–1990). Numbers of dabbling ducks (Anas rubripes, Anas platyrhynchos, Anas acuta, Anas americana, Anas crecca), ring-necked ducks (Aythya collaris), and American coot (Fulica americana) were significantly different in Currituck Sound (p<0.05) in at least two of these periods and were highest during Eurasian watermilfoil infestation. These species also constituted higher percentages of the Atlantic Flyway counts during Eurasian watermilfoil infestation. than in periods before or after. Numbers of snow geese (Chen caerulescens), Canada geese (Branta canadensis), and canvasbacks (Aythya valisneria) were significantly different between periods; however, these birds exhibited a consistent decline from the early to the late period. This study suggests that submersed macrophytes, specifically Eurasian watermilfoil, support increased numbers of dabbling ducks, ring-necked ducks, and American coot.     

Shape preferred orientation of object populations: automatic analysis of digitized images

An automatic technique is presented for systematic shape preferred orientation (SPO) analysis of crystal fabrics in rocks. This technique is based on digitization of an image with a CCD camera and separation of an object population from its matrix. The image is analysed globally using a rotating grid that is superimposed on the square pixel grid of the computer screen. The parametric method is based on counting intercepts in all directions. The interceps rose provides information on both the orientation and strength of the SPO. Problems of hardware-related counting noise anisotropy are avoided using an original method of filtering the intercept function with a linear, sliding mask on the counting grid. The technique was tested both on single objects of different shapes and on different populations of identical objects. A very high degree of accuracy was obtained for orientation measurements. Fabric shape parameters are introduced which are specific to this intercept technique; they compare favourably with shape parameters used in axial fabric analysis.