Geology Reports⌕ Search

Geology topics

Research about Virginia

Source-linked reports with geographic coverage including Virginia.

At least 109 records · Page 6Linked to original sources

Metabarcoding assays for the detection of freshwater mussels (Unionida) with environmental DNA

Freshwater mussels of the order Unionida are a widely distributed taxon that are important in maintaining freshwater ecosystems and are also highly imperiled throughout the world. Monitoring of mussel populations with environmental DNA (eDNA) is an attractive alternative to traditional methods because it is noninvasive and requires less labor and taxonomic knowledge from field personnel. We developed eDNA metabarcoding assays specific to freshwater mussels and tested them at six sites in the Clinch River, located in the southeastern United States. Our objective was to determine the utility of eDNA metabarcoding for future monitoring of mussel populations and restoration efforts in this watershed. Two metabarcoding assays that target the mitochondrial DNA regions of the cytochrome c oxidase subunit I (COI) and NADH dehydrogenase subunit (ND1) genes were developed and tested. Our assays appear to be order specific, amplifying members from the two families found in North America, Unionidae and Margaritiferidae, while not amplifying nontarget fish or other bivalve species. From the field collected samples, our assays together detected 19 species, eight of which are listed as federally endangered. The assays also detected 42%, 58%, and 54% of the species identified by recent quantitative visual mussel surveys at three sampling sites. Increased sampling effort by processing a greater water volume or number of samples will likely increase species detections. These eDNA metabarcoding assays may enable enhanced monitoring of freshwater mussel assemblages and subsequently inform conservation efforts.

Virginia↗

Dendritic reidite from the Chesapeake Bay impact horizon, Ocean Drilling Program Site 1073 (offshore northeastern USA): A fingerprint of distal ejecta?

High-pressure minerals provide records of processes not normally preserved in Earth’s crust. Reidite, a quenchable polymorph of zircon, forms at pressures >20 GPa during shock compression. However, there is no broad consensus among empirical, experimental, and theoretical studies on the nature of the polymorphic transformation. Here we decipher a multistage history of reidite growth recorded in a zircon grain in distal impact ejecta (offshore northeastern United States) from the ca. 35 Ma Chesapeake Bay impact event which, remarkably, experienced near-complete conversion (89%) to reidite. The grain displays two distinctive reidite habits: (1) intersecting sets of planar lamellae that are dark in cathodoluminescence (CL); and (2) dendritic epitaxial overgrowths on the lamellae that are luminescent in CL. While the former is similar to that described in literature, the latter has not been previously reported. A two-stage growth model is proposed for reidite formation at >40 GPa in Chesapeake Bay impact ejecta: formation of lamellar reidite by shearing during shock compression, followed by dendrite growth, also at high pressure, via recrystallization. The dendritic reidite is interpreted to nucleate on lamellae and replace damaged zircon adjacent to lamellae, which may be amorphous ZrSiO 4 or possibly an intermediate phase, all before quenching. These results provide new insights on the microstructural evolution of the high-pressure polymorphic transformation over the microseconds-long interval of reidite stability during meteorite impact. Given the formation conditions, dendritic reidite may be a unique indicator of distal ejecta.

Virginia↗

Application of a new species-richness based flow ecology framework for assessing flow reduction effects on aquatic communities

Water‐resources managers are challenged with maintaining a balance among beneficial uses throughout river networks and need robust means of assessing potential risks to aquatic life resulting from flow alterations. This study generated ecological limit functions from species‐streamflow relations to quantify potential fish richness response to flow alteration and compared results to currently accepted streamflow management guidelines. Modeled responses of absolute richness change were watershed specific and varied among sample sets derived from hydrologic unit classifications of different sizes (large HUC 6 basins to regional scale HUC 8). With a 20% flow reduction, 10% of HUC 8 predicted a richness decrease in one or more taxa. While absolute richness change was consistent across streams within a HUC, percent richness change was stream size dependent. Comparisons with Instream Flow Incremental Methodology habitat models predicted habitat loss greater than percent richness change; however, predictions for habitat and richness decreased similarly as stream size decreased. Watershed‐specific responses from flow reductions could allow water‐resources management decisions to be made locally based on the predicted richness change for certain sized streams. Quantitative results highlight the utility of a richness‐based framework for generating watershed‐specific risk assessments that validate and inform currently employed water‐resources management practices.

Virginia↗

Spatial and temporal patterns in streamflow, water chemistry, and aquatic macroinvertebrates of selected streams in Fairfax County, Virginia, 2007–18

Urbanization substantially alters the landscape in ways that can impact stream hydrology, water chemistry, and the health of aquatic communities. Stormwater best management practices (BMPs) are the primary tools used to mitigate the effects of urban stressors such as increased runoff, decreased baseflow, and increased nutrient and sediment transport. To date, Fairfax County Virginia’s stormwater management program has made substantial investments into the implementation of both structural and nonstructural BMPs aimed at restoring and protecting watersheds. The U.S. Geological Survey (USGS), in cooperation with Fairfax County, Virginia, established a long-term water-resources monitoring program to evaluate the watershed-scale effects of these investments. Monitoring began at 14 stations in 2007 and was expanded to 20 stations in 2013. This report utilized the first 10 years of data collection to (1) assess water quantity and quality, as well as ecological condition; (2) compute annual nutrient and sediment loads; and (3) evaluate trends in streamflow, water quality, and ecological condition. Efforts are underway to link the biotic and abiotic patterns described herein to watershed management practices as well as factors such as land use change, public works infrastructure, and climate. Hydrologic, chemical, and benthic macroinvertebrate community conditions in the streams monitored were similar to those observed in other studies of urban streams. Multidecadal trends in baseflow indices and runoff ratios at long-term Chesapeake Bay Non-tidal Network streamgages (CB-NTN) indicate a decrease in groundwater recharge and increase in storm runoff as a result of urbanization. Streamflow yields varied spatially with land cover, geology, and soil characteristics, whereas flashiness was positively related to impervious area. Dissolved oxygen typically was lowest in the Coastal Plain and across all Triassic Lowlands streams, and highest in the Piedmont. Dissolved oxygen concentrations generally were above Virginia’s minimum criterion of 4.0 milligrams per liter (mg/L), most violations occurred at Paul Spring Branch in the Coastal Plain during the warmest months of the year owing to increased chemical and biological oxygen demand. Typical pH values of the monitored streams centered on neutrality (pH = 7); however, diurnal fluctuations were most prevalent in the continuous pH data at Flatlick Branch (FLAT; a Triassic Lowlands station), as a result of increased photosynthesis catalyzed by phosphorus-rich geology. Specific conductance (SC) varied spatially owing to geology (highest at Triassic Lowlands stations) and anthropogenic disturbance (watersheds with high impervious land cover). Specific conductance typically was inversely related to streamflow except in winter months following deicing road salt applications, when values increased by several orders of magnitude. A significant increase in SC of about 2 percent per year was observed from the combined trend result of all monitoring stations over the 10-year period. Significant SC increases occurred at nearly all monitoring stations. Increasing trends were observed during winter and nonwinter months, which suggests that salts applied to deice roadways and other impervious surfaces are stored in the environment and released year-round. Suspended-sediment (SS) concentrations in monthly samples did not vary significantly between most stations, but typically were highest in the spring and lowest in the fall as a result of seasonal differences in streamflow and climate. Suspended-sediment yields ranged from 62 to 1,428 tons per square mile (ton/mi 2 ), with a median of 302 ton/mi 2 . Annual loads were greatest during the wettest water years (October 1-September 30; 2008, 2011, and 2014), with the greatest interannual variability occurring at Difficult Run above Fox Lake (DIFF) and South Fork Little Difficult Run (SFLIL). Suspended sediment was primarily composed of silts and clays; however, the proportion of sand in suspended sediment was related positively to streamflow. Cross-correlation analyses suggested the dominant sources of SS were streambank erosion and resuspension of in-channel material at DIFF and FLAT; whereas, upland sources and erosion of upper streambanks were more common at Dead Run (DEAD), Long Branch (LONG), and SFLIL. Median total phosphorus (TP) concentrations ranged from 0.016 to 0.077 mg/L, with a networkwide median of 0.022 mg/L, were highest in the warm season (April-September), and were composed primarily of dissolved phosphorous. Although TP concentrations were relatively low across the network, the highest concentrations were consistently at stations located in the Triassic Lowlands, owing to phosphorous-rich geology, and in the Coastal Plain, owing to the low-phosphorous sorptive capacity of those soils. A significant increase in TP concentration occurred in a few stations, but the combined trend results from all stations demonstrated a significant increase of about 4 percent per year. Networkwide increases were also observed in total dissolved phosphorus, orthophosphate, and total particulate phosphorus. The composition of TP shifted from dissolved to particulate as streamflow increased and for this reason loads primarily were composed of particulate phosphorous. Median annual TP loads were highest at FLAT and DEAD and ranged from 247 to 642 pounds per square mile (lbs/mi 2 ) networkwide. Interannual variability in phosphorous yields was apparent at most stations; the highest loading years were also the wettest years during the study period and coincident with the highest peak annual flows. Total nitrogen (TN) concentrations typically were low throughout the network with exceptions occurring at stations located in watersheds with a high density of septic infrastructure. Elevated TN concentrations also were observed in some watersheds without a high density of septic systems and may be attributable to geologic and soil properties that limit denitrification as well as other unknown anthropogenic inputs. Total nitrogen typically was dominated by nitrate during baseflows; however, the proportion of particulate nitrogen increased during stormflows. Total nitrogen yields were similar across stations, with medians ranging from about 3,600 to 6,300 lbs/mi 2 and were related to annual streamflow volume. Total nitrogen concentrations and flow-normalized concentrations decreased over the 10-year period at 7 stations, with median reductions of about 2.5 percent. Increasing trends were observed at the two stations with the highest median TN concentration (Captain Hickory Run and SFLIL, 3–5 mg/L), both watersheds contain a high density of septic infrastructure. The combined trend results from all stations revealed no trend in TN and a declining trend in nitrate of about 2 percent per year. Overall, benthic community metrics indicated that streams throughout Fairfax County were initially of poor health; however, many metrics show an improving trend (from poor to fair based on the Fairfax County Index of Biological Integrity [IBI]). Significant increasing trends in IBI occurred at the network-scale and at 4 individual stations; additionally, scores improved by at least 1 qualitative category (for example, poor to fair, fair to good) at 11 of the 14 stations between 2009 (the first year all 14 stations were sampled) and 2017. Changes in all metrics suggest that the biodiversity, function, and condition of streams in Fairfax County are improving, but some of these improvements are driven by increased diversity and percent composition of organisms that are tolerant of the urban environment.

Virginia↗

elfgen: A new instream flow framework for rapid generation and optimization of flow-ecology relations

Effective water resource management requires practical, data‐driven determination of instream flow needs. Newly developed, high‐resolution flow models and aquatic species databases provide enormous opportunity, but the volume of data can prove challenging to manage without automated tools. The objective of this study was to develop a framework of analytical methods and best practices to reduce costs of entry into flow–ecology analysis by integrating widely available hydrologic and ecological datasets. Ecological limit functions (ELFs) describing the relation between maximum species richness and stream size characteristics (streamflow or drainage area) were developed. Species richness is expected to increase with streamflow through a watershed up to a point where it either plateaus or transitions to a decreasing trend in larger streams. Our results show that identifying the location of this "breakpoint" is critical for producing optimal ELF model fit. We found that richness breakpoints can be estimated using automated low‐supervision methods, with high‐supervision providing negligible improvement in detection accuracy. Model fit (and predictive capability) was found to be superior in smaller hydrologic units. The ELF model ("elfgen" R package available on GitHub: https://github.com/HARPgroup/elfgen ) can be used to generate ELFs using built‐in datasets for the conterminous United States, or applied anywhere else streamflow and biodiversity data inputs are available.

Virginia↗

Science to support water-resource management in the upper Roanoke River watershed

Flooding, excessive sedimentation, and high bacteria counts are among the most challenging water resource issues affecting the Upper Roanoke River watershed. These issues threaten public safety, impair the watershed’s living resources, and threaten drinking water supplies, though mitigation is costly and difficult to manage. Urban development, land disturbance, and changing climatic patterns continue to challenge watershed managers who are tasked with protecting and improving the water quality of the Upper Roanoke River watershed. The U.S. Geological Survey helps watershed managers meet these demands by providing high-quality data and analyses designed to inform watershed restoration activities.

Virginia↗

Processes influencing marsh elevation change in low- and high-elevation zones of a temperate salt marsh

The movement of salt marshes into uplands and marsh submergence as sea level rises is well documented; however, predicting how coastal marshes will respond to rising sea levels is constrained by a lack of process-based understanding of how various marsh zones adjust to changes in sea level. To assess the way in which salt marsh zones differ in their elevation response to sea-level change, and to evaluate how potential hydrologic drivers influence the response, surface elevation tables, marker horizons, and shallow rod surface elevation tables were installed in a Virginia salt marsh in three zones that differed in elevation and vegetation type. Decadal rates of elevation change, surface accretion, and shallow subsidence or expansion were examined in the context of hydrologic drivers that included local sea-level rise, flooding frequency, hurricane storm surge, and precipitation. Surface elevation increases were fastest in the low-elevation zone, intermediate in the middle-elevation zone, and slowest in the high-elevation zone. These rates are similar to (low and middle marsh) or less than (high marsh) local rates of sea-level rise. Root zone expansion, presumably due to root growth and organic matter accumulation, varied among the three salt marsh zones and accounted for 37%, but probably more, of the increase in marsh surface elevation. We infer that, during marsh transgression, soil-forming processes shift from biogenic (high marsh) to minerogenic (low marsh) in response, either directly or indirectly, to changing hydrologic drivers.

Virginia↗

Chemical constituent concentrations in stream water, streambed sediment, and soils of Fort Belvoir, Virginia—A characterization of ambient conditions in 2019

Introduction The U.S. Army Fort Belvoir (FTBL) installation is on the banks of the Potomac River in Fairfax County, northeastern Virginia. The installation was founded by the U.S. Army during World War I. It has been home to a variety of military organizations over the course of its more than 100-year history and currently houses more than 145 mission partners. The installation consists of two noncontiguous units, the Main Post, and a smaller area to the northwest, Fort Belvoir North Area (FTNA). FTBL encompasses 8.91 square miles. There is concern that activities on FTBL, including a long history of training, operations, and maintenance, may have resulted in contamination of stream water, streambed sediment, and (or) soils. Of particular concern is the U.S. Environmental Protection Agency (EPA) Target Analyte List (TAL). TAL refers to “the list of inorganic compounds/elements designated for analysis as contained in the version of the EPA Contract Laboratory Program Statement of Work for Inorganics Analysis, Multi-Media, Multi-Concentration in effect as of the date on which the laboratory is performing the analysis” ( https://www.nj.gov/dep/srp/guidance/tcl_tal/ ). Because of the potential for TAL contamination at FTBL, the U.S. Geological Survey (USGS), in cooperation with U.S. Army Fort Belvoir, conducted a survey of FTBL’s stream water, streambed sediment, and soils during calendar year 2019. The terminology “ambient concentrations” is used in this report to represent the concentrations of the TAL and other constituents at the time of sampling. This is in contrast to “background concentrations,” a term that “refers to areas in which the concentrations of chemicals have not been elevated by site activities”. Although some of the samples collected for this project may represent “background concentrations,” there is no assurance that they do, so all data collected are described as having “ambient concentrations.” The purpose of the study was to obtain environmental data to characterize ambient concentrations of EPA TAL constituents in stream water, streambed sediment, and soils in FTBL, Virginia. This report describes methods and results of sampling stream water, streambed sediment, and soils during 2019. The purpose of this report is four-fold: (1) to describe the field sampling methods used to collect stream water, streambed sediment, and soils; (2) to describe the laboratory methods used to analyze the samples; (3) to report summaries of the field and laboratory results; and (4) to report the quality assurance and quality control results.

Virginia↗

Evaluating elevation change thresholds between structure-from-motion DEMs derived from historical aerial photos and 3DEP LiDAR data

This study created digital terrain models (DTMs) from historical aerial images using Structure from Motion (SfM) for a variety of image dates, resolutions, and photo scales. Accuracy assessments were performed on the SfM DTMs, and they were compared to the United States Geological Survey’s three-dimensional digital elevation program (3DEP) light detection and ranging (LiDAR) DTMs to evaluate geomorphic change thresholds based on vertical accuracy assessments and elevation change methodologies. The results of this study document a relationship between historical aerial photo scales and predicted vertical accuracy of the resultant DTMs. The results may be used to assess geomorphic change thresholds over multi-decadal timescales depending on spatial scale, resolution, and accuracy requirements. This study shows that if elevation changes of approximately ±1 m are to be mapped, historical aerial photography collected at 1:20,000 scale or larger would be required for comparison to contemporary LiDAR derived DTMs.

Virginia↗

Geomorphological mapping and anthropogenic landform change in an urbanizing watershed using structure-from-motion photogrammetry and geospatial modeling techniques

Increasing urbanization and suburban growth in cities globally has highlighted the importance of land planning using detailed geomorphologic maps that depict anthropogenic landform changes. Such mapping provides information crucial for land management, hazard identification, and the management of the challenges arising from urbanization. The development and use of quantitative and repeatable methods to map anthropogenic and natural processes are required to advance the science of urban geomorphological mapping. This study investigated the application of geospatial modeling, structure-from-motion (SfM) photogrammetric methods and DEM differencing as means of quantifying anthropogenic landform changes from archival aerial imagery. Anthropogenic landforms were incorporated into a detailed geomorphologic map in an urbanizing watershed located in the Washington, D.C. metropolitan suburb of Vienna, Virginia.

Virginia↗

Geochronologic age constraints on tectonostratigraphic units of the central Virginia Piedmont, USA

New geologic mapping coupled with uranium-lead (U-Pb) zircon geochronology (sensitive high-resolution ion microprobe-reverse geometry [SHRIMP-RG] and laser ablation-inductively coupled plasma-mass spectrometry [LA-ICP-MS]) analyses of 10 samples, provides new constraints on the tectonostratigraphic framework of the central Virginia Piedmont. Detrital zircon analysis confirms that the Silurian-Devonian Quantico Formation is a postorogenic successor basin, with zircons derived primarily from Ordovician Chopawamsic Formation volcanic rocks. Detrital zircons from strata of the Long Island syncline, previously mapped as a separate successor basin, have a peri-Gondwanan component distinct from Laurentian-sourced rocks of the Potomac terrane to the west. Volcanism of the Chopawamsic Formation spanned at least 14 million years during the Ordovician. The Chopawamsic Formation contains sheet-like Late Ordovician-Silurian granodioritic and tonalitic intrusions that were once mapped as Carboniferous. Biotite-muscovite migmatitic paragneiss, which borders the Chopawamsic Formation on its southeast side and also occurs east of the Lakeside fault, preserves evidence of Silurian deformation and metamorphism, with a Carboniferous (Alleghanian) overprint. Limited SHRIMP-RG analysis of detrital zircons from this paragneiss yields a Laurentian (Mesoproterozoic) signature, which suggests that the structurally concordant contact between volcanic rocks of the Chopawamsic Formation and paragneiss is either a pre-Alleghanian fault or an unconformity.

Virginia↗

Impacts of seagrass dynamics on the coupled long-term evolution of barrier-marsh-bay systems

Seagrass provides a wide range of economically and ecologically valuable ecosystem services, with shoreline erosion control often listed as a key service, but can also alter the sediment dynamics and waves within back-barrier bays. Here we incorporate seagrass dynamics into an existing barrier-marsh exploratory model, GEOMBEST++, to examine the coupled interactions of the back-barrier bay with both adjacent (marsh) and non-adjacent (barrier island) subsystems. While seagrass reduces marsh edge erosion rates and increases progradation rates in many of our 288 model simulations, seagrass surprisingly increases marsh edge erosion rates when sediment export from the back-barrier basin is negligible because the ability of seagrass to reduce the volume of marsh sediment eroded matters little for back-barrier basins in which all sediment is conserved. Our model simulations also suggest that adding seagrass to the bay subsystem leads to increased deposition in the bay, reduced sediment available to the marsh, and enhanced marsh edge erosion until the bay reaches a new, shallower equilibrium depth. In contrast, removing seagrass liberates previously-sequestered sediment that is then delivered to the marsh, leading to enhanced marsh progradation. Lastly, we find that seagrass reduces barrier island migration rates in the absence of back-barrier marsh by filling accommodation space in the bay. These model observations suggest that seagrass meadows operate as dynamic sources and sinks of sediment that can influence the evolution of coupled marsh and barrier island landforms in unanticipated ways.

Virginia↗

Preliminary bedrock geologic map of the Lahore 7.5-minute quadrangle, Orange, Spotsylvania, and Louisa Counties, Virginia

Introduction Bedrock geologic mapping of the Lahore, Va., 7.5-minute quadrangle was completed as part of a broader project, undertaken jointly between the U.S. Geological Survey, the Virginia Division of Geology and Mineral Resources, and other Federal and State agencies to better understand the causative mechanisms of the magnitude-5.8 (M5.8) earthquake that occurred near Mineral, Va., on August 23, 2011. This project involved detailed mapping of eight quadrangles in the epicentral region of the Mineral, Va., earthquake in order to improve our understanding of the geologic framework of the central Virginia seismic zone, which has a long record of historical and prehistoric seismicity. The Lahore 7.5-minute quadrangle contains the contact between Ordovician to Silurian, dioritic and granodioritic rocks of the Lahore and Ellisville plutons and older metasedimentary and metavolcanic rocks. The Lahore quadrangle is northeast of the Ferncliff and Louisa, Va., quadrangles, where the Shores complex is intruded by the Ellisville pluton along the pluton’s southwestern margin. The new mapping in the Lahore quadrangle shows that the Shores complex continues northeast of the Ellisville pluton. A northeast-trending mafic- and ultramafic-bearing belt within the Shores complex is a fault-bounded accretionary zone (accretionary wedge) between rocks of the Chopawamsic Formation and Laurentian slope-and-rise deposits. In the Lahore quadrangle, this belt contains several mappable, northeast- to southwest-trending mafic bodies and also includes small exposures of gabbro and talc schist. The Lahore quadrangle contains structures of both early- and late-Paleozoic age that correspond to the Taconic and Alleghanian orogenies. Taconic (Late Ordovician) S 1 schistosity in layered rocks is typically fine-grained and parallel to compositional layering, when present. Alleghanian (Pennsylvanian) S 2 schistosity is coarser and more micaceous than S 1 and is locally accompanied by a lineation that is represented by mineral lineations, micro-crenulations, or mullion fabric, and represents the hinges of F 2 folds. A foliation in the plutonic rocks is represented by an equilibrium assemblage of aligned mafic minerals and is early Paleozoic in age. Metamorphic grade in the non-plutonic rocks of the Lahore quadrangle ranges from lower-greenschist to the northwest to upper-greenschist to the southeast, as represented by mineral assemblages in non-plutonic rocks. The biotite isograd may be, in part, lithologically controlled by the contact between the informally-named Hardware and Byrd Mill formations, and locally affected by contact metamorphism by the Lahore pluton (western part of map). The Taconic garnet isograd is defined by the sparse presence of small (less than 1 millimeter), euhedral garnet crystals. Both the biotite and garnet isograds continue along strike to the southwest into the Ferncliff and Louisa quadrangles, where the isograds have been identified as Ordovician age (Taconic orogeny) based on muscovite, biotite, and amphibole 40 Ar/ 39 Ar cooling ages. Regionally, the most common trend and plunge of joints is northwest and subvertical, respectively, and orthogonal to the regional strike of foliation. Early Mesozoic extension may have reactivated the Harris Creek fault (southeast corner of map), a late Paleozoic (Alleghanian orogeny) transpressional fault that marks the contact between granodiorite of the Ellisville pluton and the Chopawamsic Formation.

Virginia↗

Patterns of denitrification potential in tidal freshwater forested wetlands

Limited evidence for spatial patterns of denitrification in tidal freshwater forested wetlands (TFFWs), seemingly due to high spatial variability in the process, is surprising considering the various spatial gradients of its biogeochemical and hydrogeomorphic controls in these ecosystems. Because certain physical environmental gradients may be useful for the prediction of denitrification in TFFWs, we measured denitrification and ecosystem attributes in hummock-hollow microtopography of TFFWs along longitudinal riverine positions (upper, middle, and lower tidal river sites, and nearby upstream nontidal forested floodplains) of the adjoining Pamunkey and Mattaponi Rivers, Virginia. We tested differences by river, site, and plot in denitrification enzyme activity (DEA) and substrate limitations of denitrification potential (DP). The Pamunkey River carries greater river nitrate concentrations, and we found less nitrate limitation of DP and greater soil nitrate in hollows of this river. DEA in tidal hummocks was positively correlated with soil organic matter, nitrogen, and carbon, with the highest rates in lower tidal sites. Hummocks also promoted greater oxygen-controlled substrate limitation of DP, whereby experimental aeration stimulated DP under subsequent inundation more in hummocks than hollows. Additionally, tidal sites had greater DEA than nontidal sites, inferred to be caused by a combination of higher moisture, organic, and nutrient content. Our results indicate that the increasing nitrogen concentrations in these rivers will increase denitrification more on the Mattaponi River by alleviating its greater nitrogen limitation compared to the Pamunkey River, and modification to sedimentation, inundation, or microtopography from sea level rise may alter denitrification gradients in TFFWs and upstream low-elevation nontidal floodplains.

Virginia↗

Geologic evidence for an icehouse Earth before the Sturtian global glaciation

Snowball Earth episodes, times when the planet was covered in ice, represent the most extreme climate events in Earth’s history. Yet, the mechanisms that drive their initiation remain poorly constrained. Current climate models require a cool Earth to enter a Snowball state. However, existing geologic evidence suggests that Earth had a stable, warm, and ice-free climate before the Neoproterozoic Sturtian global glaciation [ca. 717 million years (Ma) ago]. Here, we present eruption ages for three felsic volcanic units interbedded with glaciolacustrine sedimentary rocks from southwest Virginia, USA, that demonstrate that glacially influenced sedimentation occurred at tropical latitudes ca. 751 Ma ago. Our findings are the first geologic evidence of a cool climate teetering on the edge of global glaciation several million years before the Sturtian Snowball Earth.

Virginia↗

Bot fly parasitism of Allegheny woodrats (Neotoma magister) in Virginia

The Allegheny woodrat ( Neotoma magister ) is a species of high conservation concern and relatively well-studied with respect to habitat use/associations, food habits, conservation genetics, and population trends. However, with the exception of raccoon roundworm ( Baylisascaris procyonis ) occurrence and etiology in woodrats, most disease and parasite ecology aspects for the woodrat are unknown. Herein, we examined the prevalence of bot flies ( Cuterebra ) over nearly three decades of woodrat surveys (1990–2018) in the central Appalachian Mountains of western Virginia. We use genetic analyses to identify recent bot fly specimen collections from a woodrat captured in 2017. Though highly variable from year to year, the overall prevalence of parasitism was low (typically < 4% of captures). As such, bot flies do not appear to be a widespread parasitic burden to Allegheny woodrats in Virginia. Genetic analysis of four collected bot fly larvae was inconclusive, as the genetic signature of these woodrat bots did not match any of the six bot species known to parasitize rodents and lagomorphs in the eastern United States. Further collections and genetic analyses will be needed to determine if the genetic database is incomplete or incorrect, or if our find is a new species of bot fly not yet taxonomically recognized.

Virginia↗

Evaluating artificial shelter arrays as a minimally invasive monitoring tool for the hellbender Cryptobranchus alleganiensis

Hellbenders Cryptobranchus alleganiensis are critically imperiled amphibians throughout the eastern USA. Rock-lifting is widely used to monitor hellbenders but can severely disturb habitat. We asked whether artificial shelter occupancy (the proportion of occupied shelters in an array) would function as a proxy for hellbender abundance and thereby serve as a viable alternative to rock-lifting. We hypothesized that shelter occupancy would vary spatially in response to hellbender density, natural shelter density, or both, and would vary temporally with hellbender seasonal activity patterns and time since shelter deployment. We established shelter arrays (n = 30 shelters each) in 6 stream reaches and monitored them monthly for up to 2 yr. We used Bayesian mixed logistic regression and model ranking criteria to assess support for hypotheses concerning drivers of shelter occupancy. In all reaches, shelter occupancy was highest from June-August each year and was higher in Year 2 relative to Year 1. Our best-supported model indicated that the extent of boulder and bedrock (hereafter, natural shelter) in a reach mediated the relationship between hellbender abundance and shelter occupancy. More explicitly, shelter occupancy was positively correlated with abundance when natural shelter covered <20% of a reach, but uncorrelated with abundance when natural shelter was more abundant. While shelter occupancy should not be used to infer variation in hellbender relative abundance when substrate composition varies among reaches, we showed that artificial shelters can function as valuable monitoring tools when reaches meet certain criteria, though regular shelter maintenance is critical.

Virginia↗

Hydrogeologic framework of the Virginia Eastern Shore

The Yorktown-Eastover aquifer system of the Virginia Eastern Shore consists of upper, middle, and lower confined aquifers overlain by correspondingly named confining units and underlain by the Saint Marys confining unit. Miocene- to Pliocene-age marine-shelf sediments observed in 205 boreholes include medium- to coarse-grained sand and shells that compose the aquifers and fine-grained sand, silt, and clay that compose the confining units. The upper confining unit also includes fine-grained and organic-rich back-barrier and estuarine sediments of Pleistocene age. An overlying surficial aquifer is composed mostly of Pleistocene-age nearshore sand and gravel with smaller amounts of cobbles and boulders. In addition, Pleistocene-age sediments that fill three buried paleochannels are for the first time explicitly delineated here as distinct hydrogeologic units. Two aquifers are composed of medium- to coarse-grained fluvial sand and gravel, and an intervening confining unit is composed of fine-grained estuarine sand, silt, clay, and organic material. Aquifer and confining-unit sediments are also mixed with reworked marine-shelf sediments eroded from the sides of the paleochannels. Hydrogeologic units of the Yorktown-Eastover aquifer system generally dip eastward, are as much as several tens of feet thick, and have an undulating configuration possibly resulting from the underlying Chesapeake Bay impact crater. Aquifers and confining units are incised by the three paleochannels along an upward-widening and eastward-lengthening series of structural “windows.” Hydrogeologic units within mainstems and branching tributaries of the paleochannels dip southeastward parallel to slopes of the paleochannels, are as much as several tens of feet thick, and laterally abut the Yorktown-Eastover aquifer system along paleochannel sidewalls. The Yorktown-Eastover aquifer system is thereby hydraulically breached by the paleochannels to alternately create barriers to or conduits for groundwater flow. Results of previously documented aquifer tests at 58 wells indicate that transmissivity is generally greatest in young, shallow, and coarse-grained nearshore and fluvial sediments of the surficial aquifer and paleochannels. Transmissivity progressively decreases with depth in older, deeper, and finer grained marine-shelf sediments of the Yorktown-Eastover aquifer system, probably because they have undergone compaction as a result of greater overburden pressure over longer periods of time. Compiled chloride concentrations in samples from 330 wells generally increase downward, with most of the samples collected at altitudes above −300 feet and with most concentrations less than 250 milligrams per liter. The saltwater-transition zone has a broad trough-like shape aligned with the peninsula, being relatively shallow along the coastline and deeper along the central “spine.” Because movement of the saltwater is slow, the configuration largely reflects groundwater flow prior to widespread groundwater withdrawals. Fresh groundwater has leaked downward along deep parts of the saltwater-transition zone and leaked upward along shallower parts to discharge at the coast. The saltwater-transition zone also exhibits an anomalous ridge across the center of the peninsula. Groundwater levels indicate that the saltwater ridge formed primarily by the Exmore paleochannel acting as a large lateral collector drain. Groundwater levels were lowered, and the position of saltwater-transition zone was elevated, by a flow conduit that intercepted groundwater that otherwise would have flowed toward and discharged along the coastline. Nearly all freshwater on the Virginia Eastern Shore is supplied by groundwater withdrawals, which have lowered water levels, altered hydraulic gradients, and created a concern for saltwater intrusion. Previous characterizations of groundwater conditions that are relied on to manage groundwater development have been limited by a lack of hydrogeologic information, particularly data on buried paleochannels that are critical to safeguarding the groundwater supply. Using recently available expanded information, the U.S. Geological Survey undertook a study in cooperation with the Virginia Department of Environmental Quality during 2016–19 to develop an improved description of the groundwater system called a “hydrogeologic framework.” The hydrogeologic framework can aid water-supply planning and development by providing information on broad trends in aquifer configurations, hydraulic properties, and proximity to saltwater to avoid chloride contamination. Digital models to evaluate effects of groundwater withdrawals can also be improved with expanded data and capabilities to evaluate paleochannel hydraulic connections and the potential for saltwater movement. The hydrogeologic framework is limited by the nonuniform distribution of boreholes and the subjective delineation of aquifers and confining units, including those within paleochannels that are regarded as preliminary. The configuration of the saltwater-transition zone is also regarded as preliminary because of the nonuniform distribution of groundwater samples. Low well-sampling frequency precludes characterizing movement of the saltwater-transition zone. A monitoring strategy of sampling and possibly electromagnetic-induction well logging could be used to detect saltwater movement.

Virginia↗