Quantifying and modelling urban stream temperature: a central US watershed study

Hydrologic models that rely on site specific linear and non‐linear regression water temperature (T_w) subroutines forced solely with observed air temperature (T_a) may not accurately estimate T_w in mixed‐use urbanizing watersheds where hydrogeological and land use complexity may confound common T_w regime assumptions. A nested‐scale experimental watershed study design was used to test T_w model predictions in a representative mixed‐use urbanizing watershed of the central USA. The linear regression T_w model used in the Soil and Water Assessment Tool (SWAT), a non‐linear regression T_w model, and a process‐based T_w model that accounts for watershed hydrology were evaluated. The non‐linear regression T_w model tested at a daily time step performed significantly (P < 0.01) better than the linear T_w model currently used in SWAT. Both regression T_w models overestimated T_w in lower temperature ranges (T_w < 10.0 °C) with percent bias (PBIAS) values ranging from ?28.2% (non‐linear T_w model) to ?66.1% (linear regression T_w model) and underestimated T_w in the higher temperature range (T_w > 25.0 °C) by 3.2%, and 7.2%, respectively. Conversely, the process‐based T_w model closely estimated T_w in lower temperature ranges (PBIAS = 4.5%) and only slightly underestimated T_w in the higher temperature range (PBIAS = 1.7%). Findings illustrate the benefit of integrating process‐based T_w models with hydrologic models to improve model transferability and T_w predictive confidence in urban mixed‐land use watersheds. The findings in this work are distinct geographically and in terms of mixed‐land use complexity and are therefore of immediate value to land‐use managers in similarly urbanizing watersheds globally. Copyright © 2015 John Wiley & Sons, Ltd.     

Storage,mixing, and fluxes of water in the critical zone across northern environments inferred by stable isotopes of soil water

Quantifying soil water storage, mixing, and release via recharge, transpiration, and evaporation is essential for a better understanding of critical zone processes. Here, we integrate stable isotope (²H and ¹⁸O of soil water, precipitation, and groundwater) and hydrometric (soil moisture) data from 5 long‐term experimental catchments along a hydroclimatic gradient across northern latitudes: Dry Creek (USA), Bruntland Burn (Scotland), Dorset (Canada), Krycklan (Sweden), and Wolf Creek (Canada). Within each catchment, 6 to 11 isotope sampling campaigns occurred at 2 to 4 sampling locations over at least 1 year. Analysis for ²H and ¹⁸O in the bulk pore water was done for >2,500 soil samples either by cryogenic extraction (Dry Creek) or by direct equilibration (other sites). The results showed a similar general pattern that soil water isotope variability reflected the seasonality of the precipitation input signal. However, pronounced differences among sampling locations occurred regarding the isotopic fractionation due to evaporation. We found that antecedent precipitation volumes mainly governed the fractionation signal, temperature and evaporation rates were of secondary importance, and soil moisture played only a minor role in the variability of soil water evaporation fractionation across the hydroclimatic gradient. We further observed that soil waters beneath conifer trees were more fractionated than beneath heather shrubs or red oak trees, indicating higher soil evaporation rates in coniferous forests. Sampling locations closer to streams were more damped and depleted in their stable isotopic composition than hillslope sites, revealing increased subsurface mixing towards the saturated zone and a preferential recharge of winter precipitation. Bulk soil waters generally comprised a high share of waters older than 14 days, which indicates that the water in soil pores are usually not fully replaced by recent infiltration events. The presented stable isotope data of soil water were, thus, a useful tool to track the spatial variability of water fluxes within and from the critical zone. Such data provide invaluable information to improve the representation of critical zone processes in spatially distributed hydrological models.     

Monitoring the water balance of Lake Victoria, East Africa, from space

Using satellite gravimetric and altimetric data, we examine trends in water storage and lake levels of multiple lakes in the Great Rift Valley region of East Africa for the years 2003–2008. GRACE total water storage estimates reveal that water storage declined in much of East Africa, by as much as , while altimetric data show that lake levels in some large lakes dropped by as much as 1–2 m. The largest declines occurred in Lake Victoria, the Earth’s second largest freshwater body. Because the discharge from the outlet of Lake Victoria is used to generate hydroelectric power, the role of human management in the lake’s decline has been questioned. By comparing catchment water storage trends to lake level trends, we confirm that climatic forcing explains only about 50decline. This analysis provides an independent means of assessing the relative impacts of climate and human management on the water balance of Lake Victoria that does not depend on observations of dam discharge, which may not be publically available. In the second part of the study, the individual components of the lake water balance are estimated. Satellite estimates of changes in lake level, precipitation, and evaporation are used with observed lake discharge to develop a parameterization for estimating subsurface inflows due to changes in groundwater storage estimated from satellite gravimetry. At seasonal timescales, this approach provides closure to Lake Victoria’s water balance to within . The third part of this study uses the water balance of a downstream water body, Lake Kyoga, to estimate the outflow from Lake Victoria remotely. Because Lake Kyoga is roughly 20 times smaller in area than Lake Victoria, its water balance is strongly influenced by inflow from Lake Victoria. Lake Kyoga has been shown to act as a linear reservoir, where its outflow is proportional to the height of the lake. This model can be used with satellite altimetric lake levels to estimate a time series of Lake Victoria discharge with an rms error of about .     

Synoptic sampling reveals spatial stability in flow and chemistry patterns despite large seasonal variability in a subarctic mountain catchment

Seasonality plays a critical role in cold mountain regions as variation in air temperature, ground thermal status, and precipitation phase alter the rate, timing and magnitude of hydrological and chemical transport. Additionally, cold mountain catchments can have highly variable topography, geology, permafrost, and landcover, which intrinsically add to this irregularity. Understanding how external and internal variability act to control mass fluxes requires sampling at a high spatial resolution over time, which rarely occurs in cold remote regions. In this work, we conduct five snapshot sampling surveys across 34 subcatchments during the ice-free period in Wolf Creek Research Basin (a mesoscale montane subarctic catchment) and two additional winter surveys across a subset of sites to assess the drivers of variability in stream chemistry and discharge. We sampled for specific conductance (SpC), major ions, and dissolved organic carbon (DOC) and used statistical metrics and Bayesian mixing analysis to quantify patterns of flow and chemistry across space and time. Our results indicate patterns in both flow and chemistry remain largely consistent across seasons for all solutes. However, there was weaker correlation of chemistry between sites, suggesting asynchronous behaviour within the catchment. There was evidence of increasing production of ions and DOC along the stream network during high spring flows but not during low flows. Although concentrations and flows exhibit high seasonality in subarctic mountains, this seasonal variability does not alter spatial patterns that arise from highly variable catchment characteristics.     

Microvariability in the early stages of cobble weathering by microenvironment on a glacier foreland,Storbreen, Jotunheimen,Norway

The initial stages of cobble weathering, measured as increasing percentage porosity, were calculated for sets of cobbles taken from the foreland of Storbreen, a cirque glacier in the Jotunheimen, Norway. Cobbles were taken from in front of the 1998 glacier snout, from the 1928, 1870, 1810 and 1750 moraine crests and from the ～10 000‐year‐old land surface beyond the Neoglacial foreland limit. Three microenvironments were examined within each site: (1) lichen‐free surfaces from exposed cobbles, (2) lichen‐covered surfaces from exposed cobbles and (3) buried cobbles taken from a soil depth of ～40 cm. Percentage porosity within plagioclase minerals was calculated using backscatter electron images of prepared thin sections. Porosity was calculated from five depth profiles into each cobble. Depth profiles were subdivided into a sequence of 50 µm × 50 µm rectangles running to at least a nominal depth of 500 µm within each cobble. Three cobbles from each dated land surface were sampled, except for the 1750 surface, which had five cobbles. Statistical analysis was by analysis of variance of rank‐order transformed data. Findings indicate that cobbles close to the glacier snout are largely unweathered; also, weathering is generally weak in the 1928, 1870 and 1810 positions, but statistically significantly higher in the 1750‐ and 10 000‐year‐old positions. Weathering of buried cobbles always exceeded weathering of exposed cobbles and may possibly reach a value beyond which it cannot progress while retaining surface cohesion. The degree of weathering on lichen‐free and lichen‐covered cobble surfaces is not initially distinguishable, but diverges sharply after ～250 years, when lichen‐covered surfaces experience significantly higher totals. Overall, the weathering trend in cobbles matches that found in soils at the same sites. Copyright © 2007 John Wiley & Sons, Ltd.     

Mercury crater statistics from MESSENGER flybys: Implications for stratigraphy and resurfacing history

The primary crater population on Mercury has been modified by volcanism and secondary craters. Two phases of volcanism are recognized. One volcanic episode that produced widespread intercrater plains occurred during the period of the Late Heavy Bombardment and markedly altered the surface in many areas. The second episode is typified by the smooth plains interior and exterior to the Caloris basin, both of which have a different crater size-frequency distribution than the intercrater plains, consistent with a cratering record dominated by a younger population of impactors. These two phases may have overlapped as parts of a continuous period of volcanism during which the volcanic flux tended to decrease with time. The youngest age of smooth plains volcanism cannot yet be determined, but at least small expanses of plains are substantially younger than the plains associated with the Caloris basin. The spatial and temporal variations of volcanic resurfacing events can be used to reconstruct Mercury's geologic history from images and compositional and topographic data to be acquired during the orbital phase of the MESSENGER mission.     

Radius and limb topography of Mercury obtained from images acquired during the MESSENGER flybys

Analysis of images obtained by the MESSENGER spacecraft during its three flybys of Mercury yields a new estimate for the planet's mean radius of 2439.25±0.69 km, in agreement with results from Mariner 10 and Earth-based observations, as well as with MESSENGER altimeter and occultation data. The mean equatorial radius and polar radius are identical to within error, suggesting that rotational oblateness is negligible when compared with other sources of topography. This result is consistent with the small gravitational oblateness of the planet. Minor differences in radius obtained at different locations reflect regional variations in topography. Residual topography along three limb profiles has a dynamic range of 7.4 km and a root-mean-square roughness of 0.8 km over hemispherical scales. Following MESSENGER's entry into orbit about Mercury in March 2011, we expect considerable additional improvements to our knowledge of Mercury's size and shape.     

Making other earths: dynamical simulations of terrestrial planet formation and water delivery

We present results from 44 simulations of late stage planetary accretion, focusing on the delivery of volatiles (primarily water) to the terrestrial planets. Our simulations include both planetary “embryos” (defined as Moon to Mars sized protoplanets) and planetesimals, assuming that the embryos formed via oligarchic growth. We investigate volatile delivery as a function of Jupiter's mass, position and eccentricity, the position of the snow line, and the density (in solids) of the solar nebula. In all simulations, we form 1-4 terrestrial planets inside 2 AU, which vary in mass and volatile content. In 44 simulations we have formed 43 planets between 0.8 and 1.5 AU, including 11 “habitable” planets between 0.9 and 1.1 AU. These planets range from dry worlds to “water worlds” with 100+oceans of water (1 ocean=1.5×10²⁴ g), and vary in mass between 0.23M_⊕ and 3.85M_⊕. There is a good deal of stochastic noise in these simulations, but the most important parameter is the planetesimal mass we choose, which reflects the surface density in solids past the snow line. A high density in this region results in the formation of a smaller number of terrestrial planets with larger masses and higher water content, as compared with planets which form in systems with lower densities. We find that an eccentric Jupiter produces drier terrestrial planets with higher eccentricities than a circular one. In cases with Jupiter at 7 AU, we form what we call “super embryos,” 1-2M_⊕ protoplanets which can serve as the accretion seeds for 2+M_⊕ planets with large water contents.     

Synthetic lead bromapatite: X-ray structure at ambient pressure and compressibility up to about 20 GPa

Lead bromapatite [Pb₁₀(PO₄)₆Br₂] has been synthesized via solid-state reaction at pressures up to 1.0 GPa, and its structure determined by single-crystal X-ray diffraction at ambient temperature and pressure. The large bromide anion is accommodated in the c-axis channel by lateral displacements of structural elements, particularly of Pb2 cations and PO₄ tetrahedra. The compressibility of bromapatite was also investigated up to about 20.7 GPa at ambient temperature, using a diamond-anvil cell and synchrotron X-ray radiation. The compressibility of lead bromapatite is significantly different from that of lead fluorapatite. The pressure–volume data of lead bromapatite (P < 10 GPa) fitted to the third-order Birch-Murnaghan equation yield an isothermal bulk modulus (K _T ) of 49.8(16) GPa and first pressure derivative (K_T^￠ K_{T}^{\prime } ) of 10.1(10). If K_T^￠ K_{T}^{\prime } is fixed at 4, the derived K _T is 60.8(11) GPa. The relative difference of the bulk moduli of these two lead apatites is thus about 12%, which is about two times the relative difference of the bulk moduli (~5%) of the calcium apatites fluorapatite [Ca₁₀(PO₄)₆F₂]_, chlorapatite [Ca₁₀(PO₄)₆Cl₂] and hydroxylapatite [Ca₁₀(PO₄)₆(OH)₂]. Another interesting feature apparently related to the replacement of F by Br in lead apatite is the switch in the principle axes of the strain ellipsoid: the c-axis is less compressible than the a-axis in lead bromapatite but more compressible in lead fluorapatite.