Characterizing spatial and temporal variation in 18O and 2H content of New Zealand river water for better understanding of hydrologic processes

Time series of hydrogen and oxygen stable isotope ratios (δ²H and δ¹⁸O) in rivers can be used to quantify groundwater contributions to streamflow, and timescales of catchment storage. However, these isotope hydrology techniques rely on distinct spatial or temporal patterns of δ²H and δ¹⁸O within the hydrologic cycle. In New Zealand, lack of understanding of spatial and temporal patterns of δ²H and δ¹⁸O of river water hinders development of regional and national-scale hydrological models. We measured δ²H and δ¹⁸O monthly, together with river flow rates at 58 locations across New Zealand over a two-year period. Results show: (a) general patterns of decreasing δ²H and δ¹⁸O with increasing latitude were altered by New Zealand's major mountain ranges; δ²H and δ¹⁸O were distinctly lower in rivers fed from higher elevation catchments, and in eastern rain-shadow areas of both islands; (b) river water δ²H and δ¹⁸O values were partly controlled by local catchment characteristics (catchment slope, PET, catchment elevation, and upstream lake area) that influence evaporation processes; (c) regional differences in evaporation caused the slope of the river water line (i.e., the relationship between δ²H and δ¹⁸O in river water) for the (warmer) North Island to be lower than that of the (cooler, mountain-dominated) South Island; (d) δ²H seasonal offsets (i.e., the difference between seasonal peak and mean values) for individual sites ranged from 0.50‰ to 5.07‰. Peak values of δ¹⁸O and δ²H were in late summer, but values peaked 1 month later at the South Island sites, likely due to greater snow-melt contributions to streamflow. Strong spatial differences in river water δ²H and δ¹⁸O caused by orographic rainfall effects and evaporation may inform studies of water mixing across landscapes. Generally distinct seasonal isotope cycles, despite the large catchment sizes of rivers studied, are encouraging for transit time analysis applications.     

Do southern Appalachian Mountain summer stream temperatures respond to removal of understory rhododendron thickets?

Dense understory thickets of the native evergreen shrub Rhododendron maximum expanded initially following elimination of American chestnut by the chestnut blight, and later in response to loss of the eastern hemlock due to hemlock woolly adelgid invasion. Rhododendron thickets often blanket streams and their riparian zones, creating cool, low-light microclimates. To determine the effect of such understory thickets on summer stream temperatures, we removed riparian rhododendron understory on 300 m reaches of two southern Appalachian Mountain headwater streams, while leaving two 300 m reference reaches undisturbed. Overhead canopy was left intact in all four streams, but all streams were selected to have a significant component of dead or dying eastern hemlock in the overstory, creating time-varying canopy gaps throughout the reach. We continuously monitored temperatures upstream, within and downstream of treatment and reference reaches. Temperatures were monitored in all four streams in the summer before treatments were imposed (2014), and for two summers following treatment (2015, 2016). Temperatures varied significantly across and within streams prior to treatment and across years for the reference streams. After rhododendron removal, increases in summer stream temperatures were observed at some locations within the treatment reaches, but these increases did not persist downstream and varied by watershed, sensor, and year. Significant increases in daily maxima in treatment reaches ranged from 0.9 to 2.6°C. Overhead canopy provided enough shade to prevent rhododendron removal from increasing summer temperatures to levels deleterious to native cold-water fauna (average summer temperatures remained below 16°C), and local temperature effects were not persistent.     

The relationship between abundances of meiofauna and their suspected microbial food (diatoms and bacteria)

Samples were taken bi-weekly for one year at a sand site and a mud site in the North Inlet Estuary, Georgetown, South Carolina, for meiofauna, their suspected microbial food (bacteria and diatoms), and associated physical factors. Linear regression techniques were used to correlate food abundance and physical factors with the density of meiofaunal taxa. At both sites diatoms positively correlated with meiofauna taxa, but bacteria did not. Physical factors were not correlated with meiofaunal or microbial abundances at the sand site. Whereas, at the mud site meiofauna and diatom abundances were positively correlated with the depth of the redox layer and inversely correlated with temperature. Peaks of meiofaunal abundance did not follow peaks of food abundance. Analysis of copepods at the species level indicated that taxa response was due to the response of the dominant species. Even though some correlations existed, this study suggests that copepod species and meiofauna at the gross taxonomic level do not respond to changes in potential food abundance. Physical factors apparently influence both meiofauna and diatoms in the same fashion. However, bacterial abundance was not positively correlated with any of the factors studied.     

Emplacement age and thermal footprint of the diamondiferous Ellendale E9 lamproite pipe, Western Australia

The diamondiferous Ellendale 9 (E9) pipe is a funnel-shaped maar-diatreme volcano consisting of inward-dipping tuff sequences intruded by lamproite plugs and dykes. The host rocks for the E9 pipe are Permian sandstones. The multiple lithological contacts exposed within the mined maar volcano provide a natural laboratory in which to study the effect of volcanic processes on U–Th–Pb–He systematics. Zircon from the regional sandstone and E9 lamproite display a bimodal distribution of ages on (U–Th)/He–U/Pb plots. The zircon U/Pb ages for the E9 pipe (n?=?52) range from 440 to 2,725 Ma, while the cluster of (U–Th)/He ages for the lamproite dyke zircon indicate that dyke emplacement occurred at 20.6?±?2.8 Ma, concordant with a maximum emplacement age of about ≤22 Ma from phlogopite ⁴⁰Ar/³⁹Ar. These ages indicate a xenocrystic origin for the zircon entrained in the E9 dyke. The U/Pb ages of detrital zircon from the regional sandstone host (373–3,248 Ma; n?=?41) are indistinguishable from those of the lamproite zircon xenocrysts, whereas the detrital zircon in the host sandstone yield (U–Th)/He ages from 260 to 1,500 Ma. A thermochronology traverse across the E9 lamproite dyke reveals that the zircon (U–Th)/He ages in the host sandstone have not been significantly thermally reset during dyke emplacement, even at the contact. The capability of the zircon (U–Th)/He method to distinguish deep, mantle source lithologies from upper crustal source lithologies could be used in geochemical exploration for diamonds. Pre-screening of detrital samples using etching and helium assay methods will improve the efficiency and decrease the cost of greenfields exploration.     

Historical Trend in River Ice Thickness and Coherence in Hydroclimatological Trends in Maine

We analyzed long-term records of ice thickness on the Piscataquis River in central Maine and air temperature in Maine to determine whether there were temporal trends that were associated with climate warming. The trend in ice thickness was compared and correlated with regional time series of winter air temperature, heating degree days (HDD), date of river ice-out, seasonal center-of-volume date (SCVD) (date on which half of the stream runoff volume during the period 1 Jan. to 31 May has occurred), water temperature, and lake ice-out date. All of these variables except lake ice-out date showed significant temporal trends during the 20th century. Average ice thickness around 28 February decreased by about 23 cm from 1912 to 2001. Over the period 1900 to 1999, winter air temperature increased by 1.7 °C and HDD decreased by about 7.5%.Final ice-out date on the Piscataquis River occurred earlier (advanced), by 0.21 days yr^–1 over the period 1931 to 2002, and the SCVD advancedby 0.11 days yr^–1 over the period 1903 to 2001. Ice thickness was significantly correlated (P-value <0.01) with winter air temperature, HDD, river ice-out, and SCVD. These systematic temporal trends in multiple hydrologic indicator variables indicate a coherent response to climate forcing.     

Global observations of nonlinear mesoscale eddies

Sixteen years of sea-surface height (SSH) fields constructed by merging the measurements from two simultaneously operating altimeters are analyzed to investigate mesoscale variability in the global ocean. The prevalence of coherent mesoscale features (referred to here as “eddies”) with radius scales of O(100 km) is readily apparent in these high-resolution SSH fields. An automated procedure for identifying and tracking mesoscale features based on their SSH signatures yields 35,891 eddies with lifetimes ?16 weeks. These long-lived eddies, comprising approximately 1.15 million individual eddy observations, have an average lifetime of 32 weeks and an average propagation distance of 550 km. Their mean amplitude and a speed-based radius scale as defined by the automated procedure are 8 cm and 90 km, respectively.The tracked eddies are found to originate nearly everywhere in the World Ocean, consistent with previous conclusions that virtually all of the World Ocean is baroclinically unstable. Overall, there is a slight preference for cyclonic eddies. However, there is a preference for the eddies with long lifetimes and large propagation distances to be anticyclonic. In the southern hemisphere, the distributions of the amplitudes and rotational speeds of eddies are more skewed toward large values for cyclonic eddies than for anticyclonic eddies. As a result, eddies with amplitudes >10 cm and rotational speeds >20 cm s⁻¹ are preferentially cyclonic in the southern hemisphere. By contrast, there is a slight preference for anticyclonic eddies for nearly all amplitudes and rotational speeds in the northern hemisphere.On average, there is no evidence of anisotropy of these eddies. Their average shape is well represented as Gaussian within the central 2/3 of the eddy, but the implied radius of maximum rotational speed is 64% smaller than the observed radius of maximum speed. In part because of this mismatch between the radii of maximum axial speed in the observations and the Gaussian approximation, a case is made that a quadratic function that is a very close approximation of the mode profile of the eddy (i.e., the most frequently occurring value at each radius) is a better representation of the composite shape of the eddies. This would imply that the relative vorticity is nearly constant within the interiors of most eddies, i.e., the fluid motion consists approximately of solid-body rotation.Perhaps the most significant conclusion of this study is that essentially all of the observed mesoscale features outside of the tropical band 20°S-20°N are nonlinear by the metric U/c, where U is the maximum circum-average geostrophic speed within the eddy interior and c is the translation speed of the eddy. A value of U/c > 1 implies that there is trapped fluid within the eddy interior. Many of the extratropical eddies are highly nonlinear, with 48% having U/c > 5 and 21% having U/c > 10. Even in the tropics, approximately 90% of the observed mesoscale features are nonlinear by this measure.Two other nondimensional parameters also indicate strong degrees of nonlinearity in the tracked eddies. The distributions of all three measures of nonlinearity are more skewed toward large values for cyclonic eddies than for anticyclonic eddies in the southern hemisphere extratropics but the opposite is found in the northern hemisphere extratropics. There is thus a preference for highly nonlinear extratropical eddies to be cyclonic in the southern hemisphere but anticyclonic in the northern hemisphere.Further evidence in support of the interpretation of the observed features as nonlinear eddies is the fact that they propagate nearly due west with small opposing meridional deflections of cyclones and anticyclones (poleward and equatorward, respectively) and with propagation speeds that are nearly equal to the long baroclinic Rossby wave phase speed. These characteristics are consistent with theoretical expectations for large, nonlinear eddies. While there is no apparent dependence of propagation speed on eddy polarity, the eddy speeds relative to the local long Rossby wave phase speeds are found to be about 20% faster in the southern hemisphere than in the northern hemisphere. The distributions of the propagation directions of cyclones and anticyclones are essentially the same, except mirrored about a central azimuth angle of about 1.5° equatorward. This small, but we believe statistically significant, equatorward rotation of the central azimuth may be evidence of the effects of ambient currents (meridional advection or the effects of vertical shear on the potential vorticity gradient vector) on the propagation directions of the eddies.While the results presented here are persuasive evidence that most of the observed westward-propagating SSH variability consists of isolated nonlinear mesoscale eddies, it is shown that the eddy propagation speeds are about 25% slower than the westward propagation speeds of features in the SSH field that have scales larger than those of the tracked eddies. This scale dependence of the propagation speed may be evidence for the existence of dispersion and the presence of features that obey linear Rossby wave dynamics and have larger scales and faster propagation speeds than the nonlinear eddies. The amplitudes of these larger-scale signals are evidently smaller than those of the mesoscale eddy field since they are not easily isolated from the energetic nonlinear eddies.