Improving the spatial representation of soil properties and hydrology using topographically derived initialization processes in the SWAT model

Topography exerts critical controls on many hydrologic, geomorphologic and biophysical processes. However, many watershed modelling systems use topographic data only to define basin boundaries and stream channels, neglecting opportunities to account for topographic controls on processes such as soil genesis, soil moisture distributions and hydrological response. Here, we demonstrate a method that uses topographic data to adjust spatial soil morphologic and hydrologic attributes: texture, depth to the C‐horizon, saturated conductivity, bulk density, porosity and the water capacities at field (33 kpa) and wilting point (1500 kpa) tensions. As a proof of concept and initial performance test, the values of the topographically adjusted soil parameters and those from the Soil Survey Geographic Database (SSURGO; available at 1 : 20 000 scale) were compared with measured soil pedon pit data in the Grasslands Soil and Water Research Lab watershed in Riesel, TX. The topographically adjusted soils were better correlated with the pit measurements than were the SSURGO values. We then incorporated the topographically adjusted soils into an initialization of the Soil and Water Assessment Tool model for 15 Riesel research watersheds to investigate how changes in soil properties influence modelled hydrological responses at the field scale. The results showed that the topographically adjusted soils produced better runoff predictions in 50% of the fields, with the SSURGO soils performing better in the remainder. In addition, the a priori adjusted soils result in fewer calibrated model parameters. These results indicate that adjusting soil properties based on topography can result in more accurate soil characterization and, in some cases, improve model performance. Copyright © 2016 John Wiley & Sons, Ltd.     

Development and application of a physically based landscape water balance in the SWAT model

Watershed scale hydrological and biogeochemical models rely on the correct spatial‐temporal prediction of processes governing water and contaminant movement. The Soil and Water Assessment Tool (SWAT) model, one of the most commonly used watershed scale models, uses the popular curve number (CN) method to determine the respective amounts of infiltration and surface runoff. Although appropriate for flood forecasting in temperate climates, the CN method has been shown to be less than ideal in many situations (e.g. monsoonal climates and areas dominated by variable source area hydrology). The CN model is based on the assumption that there is a unique relationship between the average moisture content and the CN for all hydrologic response units (HRUs), and that the moisture content distribution is similar for each runoff event, which is not the case in many regions. Presented here is a physically based water balance that was coded in the SWAT model to replace the CN method of runoff generation. To compare this new water balance SWAT (SWAT‐WB) to the original CN‐based SWAT (SWAT‐CN), two watersheds were initialized; one in the headwaters of the Blue Nile in Ethiopia and one in the Catskill Mountains of New York. In the Ethiopian watershed, streamflow predictions were better using SWAT‐WB than SWAT‐CN [Nash–Sutcliffe efficiencies (NSE) of 0·79 and 0·67, respectively]. In the temperate Catskills, SWAT‐WB and SWAT‐CN predictions were approximately equivalent (NSE > 0·70). The spatial distribution of runoff‐generating areas differed greatly between the two models, with SWAT‐WB reflecting the topographical controls imposed on the model. Results show that a water balance provides results equal to or better than the CN, but with a more physically based approach. Copyright © 2010 John Wiley & Sons, Ltd.     

Capture orbits and melnikov integrals in the planar three-body problem

The separatrix between bounded and unbounded orbits in the three-body problem is formed by the manifolds of forward and backward parabolic orbits. In an ideal problem these manifolds coincide and form the boundary between the sets of bounded and unbounded orbits. As a mass parameter increases the movement of the parabolic manifolds is approximated numerically and by Melnikov's method. The evidence indicates that for positive values of the mass parameter these manifolds no longer coincide, and that capture and oscillatory orbits exist.     

Coupling thermodynamic modeling and high-resolution in situ LA-ICP-MS monazite geochronology: evidence for Barrovian metamorphism late in the Grenvillian history of southeastern Ontario

The Flinton Group is a greenschist to upper amphibolite facies package of metasediments in southeastern Ontario that was metamorphosed during the Ottawan Orogeny. Thermodynamic modeling of metapelitic mineral assemblages suggests an increase in peak conditions of metamorphism across the 40 km wide study area from 3.5 to 7.9 kbar and 540 to 715 °C. Garnet isopleth thermobarometry applied to the cores of compositionally zoned porphyroblasts reveals remarkably similar P-T conditions of initial crystallization at approximately 3.7–4.0 kbar and 512–520 °C, corresponding to a relatively high geothermal gradient of ca. 34–45 °C km^?1. It is inferred from modeling and reaction textures that metamorphism was along Barrovian P-T paths. Major and trace element zoning in garnet from one sample records a complex growth history as evidenced by major and trace element zoning and the distribution of xenotime, allanite and monazite inclusions. High-resolution (6 μm) LA-ICP-MS U-Pb geochronology performed on monazite in the rock matrix and included in the outer 150 μm of garnet rim-ward of a Y annulus revealed an age of 976?±?4 Ma. The age is interpreted to reflect monazite growth at the expense of allanite and apatite late in garnet’s growth history over the P-T interval 4.5–6.8 kbar and 540–640 °C. This new age estimate for near peak metamorphism fits well into the regional framework but is significantly younger than previously reported ages for Ottawan metamorphism. Based on microstructures this new age suggests that compressional tectonics were operating much later in the history of the Grenville of southeastern Ontario than previously thought.     

Predicting phosphorus dynamics in complex terrains using a variable source area hydrology model

Phosphorus (P) loss from agricultural watersheds has long been a critical water quality problem, the control of which has been the focus of considerable research and investment. Preventing P loss depends on accurately representing the hydrological and chemical processes governing P mobilization and transport. The Soil and Water Assessment Tool (SWAT) is a watershed model commonly used to predict run‐off and non‐point source pollution transport. SWAT simulates run‐off employing either the curve number (CN) or the Green and Ampt methods, both assume infiltration‐excess run‐off, although shallow soils underlain by a restricting layer commonly generate saturation‐excess run‐off from variable source areas (VSA). In this study, we compared traditional SWAT with a re‐conceptualized version, SWAT‐VSA, that represents VSA hydrology, in a complex agricultural watershed in east central Pennsylvania. The objectives of this research were to provide further evidence of SWAT‐VSA's integrated and distributed predictive capabilities against measured surface run‐off and stream P loads and to highlight the model's ability to drive sub‐field management of P. Thus, we relied on a detailed field management database to parameterize the models. SWAT and SWAT‐VSA predicted discharge similarly well (daily Nash–Sutcliffe efficiencies of 0.61 and 0.66, respectively), but SWAT‐VSA outperformed SWAT in predicting P export from the watershed. SWAT estimated lower P loss (0.0–0.25 kg ha^?1) from agricultural fields than SWAT‐VSA (0.0–1.0+ kg ha^?1), which also identified critical source areas – those areas generating large run‐off and P losses at the sub‐field level. These results support the use of SWAT‐VSA in predicting watershed‐scale P losses and identifying critical source areas of P loss in landscapes with VSA hydrology. Copyright © 2014 John Wiley & Sons, Ltd.     

Petrology of the hilina formation,Kilauea volcano,Hawaii

The Hilina Formation comprises the oldest sequence of lava flows and tuffs exposed on Kilauea Volcano. These rocks are only exposed in kipukas in younger Puna Formation lavas along cliffs on the south flank of Kilauea Volcano. Locally, tuffs and flows of the Pahala Formation separate the underlying Hilina Formation rocks rom the overlying Puna Formation rocks. Charcoal collected from the base of the Pahala Formation yielded a C¹⁴ age of 22.800±340 years B.P. which defines a minimum age for the Hilina Formation. Hilina Formation lavas crop out over a wide region and probably originated from the summit area and from both rift zones. The Hilina Formation contains both olivine-controlled and differentiated lavas (using the terminology ofWright, 1971). The olivine-controlled lavas of the Hilina Formation are distinguishable mineralogically and geochemically from younger olivine-controlled Kilauea lavas. The younger lavas generally contain discrete low-calcium pyroxene grains. greater glass contents, higher K₂O/P₂O₅ ratios and lower total iron contents. Similar geochemical trends prevail for Manuna Loa lavas, and may typify the early lavas of Hawaiian shield volcanoes. Despite these similarities, the Hilina Formation (and all Kilauea) lavas have higher TiO₂ and CaO, and lower SiO₂ and Al₂O₃ contents than Mauna Loa Lavas. These differences have existed for over 30,000 years. Therefore, it is unlikely that the older lavas of Kilauea are compositionally similar to recent Mauna Loa lavas as was previously suggested. K₂O, TiO₂, Na₂ and Zr contents of lavas from a stratigraphic sequence of Hilina Formation lavas are variable. These variations may be utilized to subdivide the sequence into geochemical groups. These groups are not magma batches. Rather, they represent lavas from batches whose compositions may have been modified by crystal fractionation and magma mixing.     

Holocene sea-level changes in Palau,West Caroline Islands

¹⁴C ages supplemented by ${}^{230}\text{Th}{}^{234}\text{U}$ determinations have been obtained for calcareous deposits on Koror, Babelthuap, Auluptagel, and Adorius Islands in the Palau Group, West Caroline Islands. Test borings for a bridge between Koror and Babelthuap reveal shoreface terraces consisting largely of bioclastic sand, but resembling fringing reefs. The base of the shoreface terraces dates to between 7000 and 8000 ¹⁴C yr B.P. and the upper surfaces are slightly younger than 4000 yr B.P. Spacing of the subsurface isochrons indicates that the rate of sedimentation increased up-section on Koror terrace and decreased up-section on Babelthuap terrace. The average rate of deposition in the terraces was 0.5 cm/yr. Isochrons in the Babelthuap shoreface terrace are 7 m higher than those in the Koror terrace. If deposition was sufficient to keep the surfaces of the terraces at low tide level, then the Babelthuap side was essentially stable and the Koror side was uplifted between about 8000 and 6300 yr ago, and then subsided until 4000 yr ago; since then there has been about 2 m of uplift. “Top hat” microatolls on Koror terrace indicate that relative sea level has dropped about 30 cm in the past 75 yr. In the absence of good evidence for changes of level in water ponded in a moat, it is likely that the microatolls indicate uplift of the terrace. On the other hand, if tectonic activity was minimal, then differences in the two terraces are due to differences in sedimentation with the Koror side of the channel being substantially subtidal between about 7000 and 5000 yr B.P. Taking the composite eustatic sea-level curve of Hawaii and elsewhere as a reference standard, it is deduced that Auluptagel Island has risen 0.8 m in the last 2900 yr, and Adorius Island has risen approximately 8 m in the last 5000–6000 yr.     

GRAIN-SIZE DISTRIBUTION AND MORPHOLOGY OF METAL IN E-CHONDRITES

The size distribution and morphology of metal grains have been examined in 11 sections of types I and II E-chondrites. The changes in the grain-size distribution and morphology of metal grains correspond with the petrologic types and define a series that reflects increase in thermal metamorphism in the following order: type I, Kota Kota-Indarch-South Oman-St. Mark's; and type II, Jajh deh Kot Lalu-Atlanta-Daniel's Kuil-Hvittis-Pillistfer-Khairpur-Blithfield. Concentrations of metal grains adjacent to the perimeters of chondrules are observable throughout the sequence and delineate relic chondritic structure in six of the seven type II E-chondrites; relic structures are absent from Blithfield.     

IRON METEORITES OF GROUP IA AND IIA: STUDY OF PHOSPHIDE MORPHOLOGY ENHANCED BY DRY CI2 ETCHING

A method is described for etching iron meteorites by heating at 350 °C in dry chlorine gas. Iron is volatilized from low-nickel iron phases, e.g. kamacite, and the residual non-volatile chlorides are removed by solution. The study of the morphological character of the exposed phosphides on 16 specimens from Group IA and IIA showed no distinctive habit restricted to either group.