Daniel H. Doctor

Sources, transformations, and hydrological processes that control stream nitrate and dissolved organic matter concentrations during snowmelt in an upland forest

We explored catchment processes that control stream nutrient concentrations at an upland forest in northeastern Vermont, USA, where inputs of nitrogen via atmospheric deposition are among the highest in the nation and affect ecosystem functioning. We traced sources of water, nitrate, and dissolved organic matter (DOM) using stream water samples collected at high frequency during spring snowmelt. Hydrochemistry, isotopic tracers, and end-member mixing analyses suggested the timing, sources, and source areas from which water and nutrients entered the stream. Although stream-dissolved organic carbon (DOC) and dissolved organic nitrogen (DON) both originated from leaching of soluble organic matter, flushing responses between these two DOM components varied because of dynamic shifts of hydrological flow paths and sources that supply the highest concentrations of DOC and DON.

Coupled hydrological and biogeochemical processes controlling variability of nitrogen species in streamflow during autumn in an upland forest

Autumn is a season of dynamic change in forest streams of the northeastern United States due to effects of leaf fall on both hydrology and biogeochemistry. Few studies have explored how interactions of biogeochemical transformations, various nitrogen sources, and catchment flow paths affect stream nitrogen variation during autumn. To provide more information on this critical period, we studied (1) the timing, duration, and magnitude of changes to stream nitrate, dissolved organic nitrogen (DON), and ammonium concentrations; (2) changes in nitrate sources and cycling; and (3) source areas of the landscape that most influence stream nitrogen. We collected samples at higher temporal resolution for a longer duration than typical studies of stream nitrogen during autumn. This sampling scheme encompassed the patterns and extremes that occurred during base flow and stormflow events of autumn. Base flow nitrate concentrations decreased by an order of magnitude from 5.4 to 0.7 µmol L-1 during the week when most leaves fell from deciduous trees. Changes to rates of biogeochemical transformations during autumn base flow explained the low nitrate concentrations; in-stream transformations retained up to 72% of the nitrate that entered a stream reach. A decrease of in-stream nitrification coupled with heterotrophic nitrate cycling were primary factors in the seasonal nitrate decline. The period of low nitrate concentrations ended with a storm event in which stream nitrate concentrations increased by 25-fold. In the ensuing weeks, peak stormflow nitrate concentrations progressively decreased over closely spaced, yet similarly sized events. Most stormflow nitrate originated from nitrification in near-stream areas with occasional, large inputs of unprocessed atmospheric nitrate, which has rarely been reported for nonsnowmelt events. A maximum input of 33% unprocessed atmospheric nitrate to the stream occurred during one event. Large inputs of unprocessed atmospheric nitrate show direct and rapid effects on forest streams that may be widespread, although undocumented, throughout nitrogen-polluted temperate forests. In contrast to a week-long nitrate decline during peak autumn litterfall, base flow DON concentrations increased after leaf fall and remained high for 2 months. Dissolved organic nitrogen was hydrologically flushed to the stream from riparian soils during stormflow. In contrast to distinct seasonal changes in base flow nitrate and DON concentrations, ammonium concentrations were typically at or below the detection limit, similar to the rest of the year. Our findings reveal couplings among catchment flow paths, nutrient sources, and transformations that control seasonal extremes of stream nitrogen in forested landscapes.

Predicting the Occurrence of Cave-Inhabiting Fauna Based on Features of the Earth Surface Environment

One of the most challenging fauna to study in situ is the obligate cave fauna because of the difficulty of sampling. Cave-limited species display patchy and restricted distributions, but it is often unclear whether the observed distribution is a sampling artifact or a true restriction in range. Further, the drivers of the distribution could be local environmental conditions, such as cave humidity, or they could be associated with surface features that are surrogates for cave conditions. If surface features can be used to predict the distribution of important cave taxa, then conservation management is more easily obtained. We examined the hypothesis that the presence of major faunal groups of cave obligate species could be predicted based on features of the earth surface. Georeferenced records of cave obligate amphipods, crayfish, fish, isopods, beetles, millipedes, pseudoscorpions, spiders, and springtails within the area of Appalachian Landscape Conservation Cooperative in the eastern United States (Illinois to Virginia and New York to Alabama) were assigned to 20 x 20 km grid cells. Habitat suitability for these faunal groups was modeled using logistic regression with twenty predictor variables within each grid cell, such as percent karst, soil features, temperature, precipitation, and elevation. Models successfully predicted the presence of a group greater than 65% of the time (mean = 88%) for the presence of single grid cell endemics, and for all faunal groups except pseudoscorpions. The most common predictor variables were latitude, percent karst, and the standard deviation of the Topographic Position Index (TPI), a measure of landscape rugosity within each grid cell. The overall success of these models points to a number of important connections between the surface and cave environments, and some of these, especially soil features and topographic variability, suggest new research directions. These models should prove to be useful tools in predicting the presence of species in understudied areas.

Spatial analysis of geologic and hydrologic features relating to sinkhole occurrence in Jefferson County, West Virginia

In this study the influence of geologic features related to sinkhole susceptibility was analyzed and the results were mapped for the region of Jefferson County, West Virginia. A model of sinkhole density was constructed using Geographically Weighted Regression (GWR) that estimated the relations among discrete geologic or hydrologic features and sinkhole density at each sinkhole location. Nine conditioning factors on sinkhole occurrence were considered as independent variables: distance to faults, fold axes, fracture traces oriented along bedrock strike, fracture traces oriented across bedrock strike, ponds, streams, springs, quarries, and interpolated depth to groundwater. GWR model parameter estimates for each variable were evaluated for significance, and the results were mapped. The results provide visual insight into the influence of these variables on localized sinkhole density, and can be used to provide an objective means of weighting conditioning factors in models of sinkhole susceptibility or hazard risk.

Using Geophysics to Map Bedrock Faults, Dikes, and Surficial Geology in Relation to Karst Features in the Briery Branch Quadrangle, Rockingham County, Virginia

at the surface through outcrop observation. Sinkholes and other karst features were also identified in the ERI profiles, generally expressed as low-resistivity anomalies. A pseudo-3D ERI survey was completed over a closed depression in highly weathered alluvium overlying limestone bedrock that clearly images the outline of the depression in the subsurface, as well as a possible water or mud-filled solutional conduit. The surface expression of the depression was approximately 165 ft (50 m) in diameter, while the low-resistivity anomaly was 15 to 50 ft (5 to 15 m) wide and appeared at depth of 15 to 80 ft (5 to 25 m).

Hypogene caves of the central Appalachian Shenandoah Valley in Virginia

Several caves in the Shenandoah Valley in Virginia show evidence for early hypogenic conduit development with later-enhanced solution under partly confined phreatic conditions guided by geologic structures. Many (but not all) of these caves have been subsequently invaded by surface waters as a result of erosion and exhumation. Those not so affected are relict phreatic caves, bearing no relation to modern drainage patterns. Field and petrographic evidence shows that carbonate rocks hosting certain relict phreatic caves were dolomitized and/or silicified by early hydrothermal fluid migration in zones that served to locally enhance rock porosity, thus providing preferential pathways for later solution by groundwater flow, and making the surrounding bedrock more resistant to surficial weathering to result in caves that reside within isolated hills on the land surface. Features suggesting that deep phreatic processes dominated the development of these relict caves include (1) cave passage morphologies indicative of ascending fluids, (2) cave plans of irregular pattern, reflecting early maze or anastomosing development, (3) a general lack of cave breakdown and cave streams or cave stream deposits, and (4) calcite wall and pool coatings within isolated caves intersecting the local water table, and within unroofed caves at topographic locations elevated well above the local base level. Episodes of deep karstification were likely separated by long periods of geologic time, encompassing multiple phases of sedimentary fill and excavation within caves, and reflect a complex history of deep fluid migration that set the stage for later shallow speleogenesis that continues today.

New insight into the origin of manganese oxide ore deposits in the Appalachian Valley and Ridge of northeastern Tennessee and northern Virginia, USA

Manganese oxide deposits have long been observed in association with carbonates within the Appalachian Mountains, but their origin has remained enigmatic for well over a century. Ore deposits of Mn oxides from several productive sites located in eastern Tennessee and northern Virginia display morphologies that include botryoidal and branching forms, massive nodules, breccia matrix cements, and fracture fills. The primary ore minerals include hollandite, cryptomelane, and romanechite. Samples of Mn oxides from multiple localities in these regions were analyzed using electron microscopy, X-ray analysis, Fourier transform infrared spectroscopy, and trace/rare earth element geochemistry. The samples from eastern Tennessee have biological morphologies, contain residual biopolymers, and exhibit REE signatures that suggest the ore formation was due to supergene enrichment (likely coupled with microbial activity). In contrast, several northern Virginia ores hosted within quartz-sandstone breccias exhibit petrographic relations, mineral morphologies, and REE signatures indicating inorganic precipitation, and a likely hydrothermal origin with supergene overprinting. Nodular accumulations of Mn oxides within weathered alluvial deposits that occur near to breccia-hosted Mn deposits in Virginia show geochemical signatures that are distinct from the breccia matrices, and appear to reflect remobilization of earlier-emplaced Mn and concentration within supergene traps. Based on the proximity of all of the productive ore deposits to mapped faults or other zones of deformation, we suggest that the primary source of all of the Mn may have been deep-seated, and that Mn oxides with supergene and/or biological characteristics result from the local remobilization and concentration of this primary Mn.

Instrumenting Caves to Collect Hydrologic and Geochemical Data: Case Study from James Cave, Virginia

Karst aquifers are productive groundwater systems, supplying approximately 25 % of the world’s drinking water. Sustainable use of this critical water supply requires information about rates of recharge to karst aquifers. The overall goal of this project is to collect long-term, high-resolution hydrologic and geochemical datasets at James Cave, Virginia, to evaluate the quantity and quality of recharge to the karst system. To achieve this goal, the cave has been instrumented for continuous (10-min interval) measurement of the (1) temperature and rate of precipitation; (2) temperature, specific conductance, and rate of epikarst dripwater; (3) temperature of the cave air; and (4) temperature, conductivity, and discharge of the cave stream. Instrumentation has also been installed to collect both composite and grab samples of precipitation, soil water, the cave stream, and dripwater for geochemical analysis. This chapter provides detailed information about the instrumentation, data processing, and data management; shows examples of collected datasets; and discusses recommendations for other researchers interested in hydrologic and geochemical monitoring of cave systems. Results from the research, briefly described here and discussed in more detail in other publications, document a strong seasonality of the start of the recharge season, the extent of the recharge season, and the geochemistry of recharge.