Sources of Suspended-Sediment Flux in Streams of the Chesapeake Bay Watershed: A Regional Application of the SPARROW Model1

Brakebill, John W., Scott W. Ator, and Gregory E. Schwarz, 2010. Sources of Suspended-Sediment Flux in Streams of the Chesapeake Bay Watershed: A Regional Application of the SPARROW Model. Journal of the American Water Resources Association (JAWRA) 46(4): 757-776. DOI: 10.1111/j.1752-1688.2010.00450.x Abstract: We describe the sources and transport of fluvial suspended sediment in nontidal streams of the Chesapeake Bay watershed and vicinity. We applied SPAtially Referenced Regressions on Watershed attributes, which spatially correlates estimated mean annual flux of suspended sediment in nontidal streams with sources of suspended sediment and transport factors. According to our model, urban development generates on average the greatest amount of suspended sediment per unit area (3,928 Mg/km²/year), although agriculture is much more widespread and is the greatest overall source of suspended sediment (57 Mg/km²/year). Factors affecting sediment transport from uplands to streams include mean basin slope, reservoirs, physiography, and soil permeability. On average, 59% of upland suspended sediment generated is temporarily stored along large rivers draining the Coastal Plain or in reservoirs throughout the watershed. Applying erosion and sediment controls from agriculture and urban development in areas of the northern Piedmont close to the upper Bay, where the combined effects of watershed characteristics on sediment transport have the greatest influence may be most helpful in mitigating sedimentation in the bay and its tributaries. Stream restoration efforts addressing floodplain and bank stabilization and incision may be more effective in smaller, headwater streams outside of the Coastal Plain.     

Monitoring total suspended solids by using nephelometry

Correlation curves were developed relating nephelometric turbidity units (NTU) with total suspended solids (TSS) for diked upland dredged material placement site effluents of three US Army Corps of Engineers (COE) maintenance dredging projects in the Chesapeake Bay, Maryland. The procedure was developed in an effort to ensure compliance with Maryland's 400 milligrams per liter (mg/l) TSS standard for COE dredging projects. Samples of the sediments to be dredged were collected and analyzed, correlating turbidity readings with TSS determined by standard gravimetric techniques. The correlation curves were provided to the COE inspectors to measure the effluent with a turbidity meter and to extract a TSS concentration from the correlation curve. Samples collected and analyzed after initiation of the dredging indicated that the correlation curves were an overestimate of the actual TSS concentrations of the effluent discharges. The procedure, endorsed by the State of Maryland, provided immediate on-site TSS analysis eliminating the previously encountered delays in obtaining gravimetric analysis of effluent discharges and potential contract management problems.     

USING MARYLAND'S STREAM CORRIDOR ASSESSMENT SURVEY TO PRIORITIZE WATERSHED RESTORATION EFFORTS1

ABSTRACT: The Stream Corridor Assessment survey has been developed by the Maryland Department of Natural Resources as a watershed management tool to identify environmental problems and to help prioritize restoration opportunities on a watershed basis. Potential environmental problems commonly identified during the survey include: stream channel alterations, excessive bank erosion, exposed pipes, inadequate stream buffers, fish migration blockages, trash dumping sites, near stream construction, pipe outfalls, and unusual conditions. In addition, the survey records information on the location of potential wetlands creation sites and collects data on the general condition of instream and riparian habitats. Over the past several years, in collaboration with the Maryland Conservation Corps, watershed associations, and local governments, more than 3,293 km (2,046 miles) of Maryland streams have been surveyed. Overall, the survey has proven to be a cost effective starting point for many watershed restoration efforts.     

Impact of floodplain and Stage 0 stream restoration on flood attenuation and floodplain exchange during small frequent storms

River flooding impacts human life and infrastructure, yet provides habitat and ecosystem services. Traditional flood control (e.g., levees, dams) reduces habitat and ecosystem services, and exacerbates flooding elsewhere. Floodplain restoration (i.e., bankfull floodplain reconnection and Stage 0) can also provide flood management, but has not been sufficiently evaluated for small frequent storms. We used 1D unsteady Hydrologic Engineering Center's River Analysis System to simulate small storms in a 5 km-long, second-order generic stream from the Chesapeake Bay watershed, and varied % channel restored (starting at the upstream end), restoration location, restoration bank height (distinguishes bankfull from Stage 0 restoration), and floodplain width/Manning's n. Stream restoration decreased (attenuated) peak flow up to 37% and increased floodplain exchange by up to 46%. Floodplain width and % channel restored had the largest impact on flood attenuation. The incremental effects of new restoration projects on flood attenuation were greatest when little prior restoration had occurred. By contrast, incremental effects on floodplain exchange were greatest in the presence of substantial prior restoration, setting up a tradeoff. A similar tradeoff was revealed between attenuation and exchange for project location, but not bank height or floodplain width. In particular, attenuation and exchange were always greater for Stage 0 than for bankfull floodplain restoration. Stage 0 thus may counteract human impacts such as urbanization.     

Restoration of a Canadian Prairie Wetland with Agricultural and Municipal Wastewater

3 /day (800,000 US gallons) of municipal wastewater and beef processing wastewater. A large nongovernmental organization hastened restoration with a development process that outlined restoration goals and management objectives to satisfy a dual mandate of wastewater treatment and wildlife habitat creation. In 1995, after five years of wastewater additions, the basins had been refilled and the surrounding uplands had been acquired and restored. The Frank Lake Conservation Area currently provides high-quality habitat for a variety of wildlife in a region where many of the native plants and animals species have been lost due to habitat loss and fragmentation. The success of upland and water management strategies is reflected in the increase of target species' abundance and richness: 50 shorebird species, 44 waterfowl species, 15 raptor species, and 28 other new bird species have returned to the marsh since restoration. As well, significant N and P reduction occurs as waters flow through the first basin of the marsh. The management strategies of this project that satisfied a dual mandate serve as a model to guide managers of other large-scale wetland restoration projects.     

Deriving Chesapeake Bay Water Quality Standards

Achieving and maintaining the water quality conditions necessary to protect the aquatic living resources of the Chesapeake Bay and its tidal tributaries has required a foundation of quantifiable water quality criteria. Quantitative criteria serve as a critical basis for assessing the attainment of designated uses and measuring progress toward meeting water quality goals of the Chesapeake Bay Program partnership. In 1987, the Chesapeake Bay Program partnership committed to defining the water quality conditions necessary to protect aquatic living resources. Under section 303(c) of the Clean Water Act, States and authorized tribes have the primary responsibility for adopting water quality standards into law or regulation. The Chesapeake Bay Program partnership worked with U.S. Environmental Protection Agency to develop and publish a guidance framework of ambient water quality criteria with designated uses and assessment procedures for dissolved oxygen, water clarity, and chlorophyll a for Chesapeake Bay and its tidal tributaries in 2003. This article reviews the derivation of the water quality criteria, criteria assessment protocols, designated use boundaries, and their refinements published in six addendum documents since 2003 and successfully adopted into each jurisdiction's water quality standards used in developing the Chesapeake Bay Total Maximum Daily Load.     

Legacy Effects of Colonial Millponds on Floodplain Sedimentation,Bank Erosion,and Channel Morphology,Mid‐Atlantic,USA1

Abstract: Many rivers and streams of the Mid‐Atlantic Region, United States (U.S.) have been altered by postcolonial floodplain sedimentation (legacy sediment) associated with numerous milldams. Little Conestoga Creek, Pennsylvania, a tributary to the Susquehanna River and the Chesapeake Bay, is one of these streams. Floodplain sedimentation rates, bank erosion rates, and channel morphology were measured annually during 2004‐2007 at five sites along a 28‐km length of Little Conestoga Creek with nine colonial era milldams (one dam was still in place in 2007). This study was part of a larger cooperative effort to quantify floodplain sedimentation, bank erosion, and channel morphology in a high sediment yielding region of the Chesapeake Bay watershed. Data from the five sites were used to estimate the annual volume and mass of sediment stored on the floodplain and eroded from the banks for 14 segments along the 28‐km length of creek. A bank and floodplain reach based sediment budget (sediment budget) was constructed for the 28 km by summing the net volume of sediment deposited and eroded from each segment. Mean floodplain sedimentation rates for Little Conestoga Creek were variable, with erosion at one upstream site (?5 mm/year) to deposition at the other four sites (highest = 11 mm/year) despite over a meter of floodplain aggradation from postcolonial sedimentation. Mean bank erosion rates range between 29 and 163 mm/year among the five sites. Bank height increased 1 m for every 10.6 m of channel width, from upstream to downstream (R² = 0.79, p < 0.0001) resulting in progressively lowered hydraulic connectivity between the channel and the floodplain. Floodplain sedimentation and bank erosion rates also appear to be affected by the proximity of the segments to one existing milldam, which promotes deposition upstream and scouring downstream. The floodplain and bank along the 28‐km reach produced a net mean sediment loss of 5,634 Mg/year for 2004‐2007, indicating that bank erosion was exceeding floodplain sedimentation. In particular, the three segments between the existing dam and the confluence with the Conestoga River (32% of the studied reach) account for 97% of the measured net sediment budget. Future research directed at understanding channel equilibria should facilitate efforts to reduce the sediment impacts of dam removal and legacy sediment.     

ENVIRONMENTAL FEATURES FOR FLOOD CONTROL CHANNELS1

ABSTRACT: The environmental effects of flood control channel modifications such as clearing and snagging, straightening, enlargement, and/or paving can be quite severe in some cases. Information review reveals that several environmental features have been incorporated into the design, construction, operation, or maintenance of recent flood control channel projects to avoid adverse environmental impacts and enhance environmental quality. Typically, these features have been proposed by conservation agencies and designed with minimal quantitative analysis. Environmental features for channel projects include selective clearing and snagging techniques, channel designs with nonuniform geometry such as single bank modification and floodways, restoration and enhancement of aquatic habitat, improved techniques for placement of excavated material, and revegetation.     

Water Quality Functions of Riparian Forest Buffers in Chesapeake Bay Watersheds

/ Maryland, Virginia, and Pennsylvania, USA, have agreed to reduce nutrient loadings to Chesapeake Bay by 40% by the year 2000. This requires control of nonpoint sources of nutrients, much of which comes from agriculture. Riparian forest buffer systems (RFBS) provide effective control of nonpoint source (NPS) pollution in some types of agricultural watersheds. Control of NPS pollution is dependent on the type of pollutant and the hydrologic connection between pollution sources, the RFBS, and the stream. Water quality improvements are most likely in areas of where most of the excess precipitation moves across, in, or near the root zone of the RFBS. In areas such as the Inner Coastal Plain and Piedmont watersheds with thin soils, RFBS should retain 50%-90% of the total loading of nitrate in shallow groundwater, sediment in surface runoff, and total N in both surface runoff and groundwater. Retention of phosphorus is generally much less. In regions with deeper soils and/or greater regional groundwater recharge (such as parts of the Piedmont and the Valley and Ridge), RFBS water quality improvements are probably much less. The expected levels of pollutant control by RFBS are identified for each of nine physiographic provinces of the Chesapeake Bay Watershed. Issues related to of establishment, sustainability, and management are also discussed.KEY WORDS: Riparian forest buffers; Chesapeake Bay; Nonpoint source pollution; Nitrogen; Phosphorus; Sediment     

Monitoring Regional Riparian Forest Cover Change Using Stratified Sampling and Multiresolution Imagery1

Claggett, Peter R., Judy A. Okay, and Stephen V. Stehman, 2010. Monitoring Regional Riparian Forest Cover Change Using Stratified Sampling and Multiresolution Imagery. Journal of the American Water Resources Association (JAWRA) 46(2):334-343. DOI: 10.1111/j.1752-1688.2010.00424.x Abstract: The Chesapeake Bay watershed encompasses 165,760 km² of land area with 464,098 km of rivers and streams. As part of the Chesapeake Bay restoration effort, state and federal partners have committed to restoring 26,000 miles (41,843 km) of riparian forest buffers. Monitoring trends in riparian forest buffers over large areas is necessary to evaluate the efficacy of these restoration efforts. A sampling approach for estimating change in riparian forest cover from 1993/1994 to 2005 was developed and implemented in Anne Arundel County, Maryland, to exemplify a method that could be applied throughout the Bay watershed. All stream reaches in the county were stratified using forest cover change derived from Landsat imagery. A stratified random sample of 219 reaches was selected and forest cover change within the riparian buffer of each sampled reach was interpreted from high-resolution aerial photography. The estimated footprint of gross change in riparian forest cover (i.e., the sum of gross gain and gross loss) for the county was 1.83% (SE = 0.22%). Stratified sampling taking advantage of a priori knowledge of locations of change proved to be a practical and efficient protocol for estimating riparian forest buffer change at the county scale and the protocol would readily extend to much broader scale monitoring.     

Prediction of Annual Streambank Erosion for Sequoia National Forest,California

Many bank erosion models have limitations that restrict their use in wildland settings. Scientists and land managers at the Sequoia National Forest would like to understand the mechanisms and rates of streambank erosion to evaluate management issues and post‐wildfire effects. This study uses bank erosion hazard index (BEHI) and near‐bank stress (NBS) methods developed in Rosgen (2006 Watershed Assessment of River Stability and Sediment Supply [WARSSS]) for predicting streambank erosion in a geographic area that is dominated by colluvium and in which streambank erosion modeling has not been previously evaluated. BEHI evaluates bank susceptibility to erosion based on bank angle, bank and bankfull height, rooting depth and density, surface protection, and stratification of material within the banks. NBS assesses energy distribution against the bank measured as a ratio of bankfull near‐bank maximum depth to mean bankfull depth. We compared BEHI classes and NBS to actual bank erosion measured from 2008 to 2012. This index predicted streambank erosion with clear separation among BEHI ratings with R² values of 0.76 for extreme, 0.37 for high/very high, 0.49 for moderate, and 0.70 for low BEHI. The relationships between measured erosion and BEHI extend the application of BEHI/NBS to a new region where they can inform management priorities, afforestation, stream/riparian restoration projects, and potentially burned area rehabilitation.     

CHARGES FOR URBAN RUNOFF: ISSUES IN IMPLEMENTATION1

ABSTRACT: Maryland officials have identified stormwater utilities as a potential method of financing programs to control nutrients in urban stormwater runoff that are proposed in Maryland's Chesapeake Bay Nutrient Reduction Plan. This paper reviews a number of issues related to the equity, efficiency and acceptability of user charge schemes. Overall, charges are found to be preferable to property taxes from both equity and efficiency perspectives. In addition, evidence suggests that elected officials will support creation of utilities. Obstacles to the implementation of utilities are identified.     

Contrasting Landscape Influences on Sediment Supply and Stream Restoration Priorities in Northern Fennoscandia (Sweden and Finland) and Coastal British Columbia

Sediment size and supply exert a dominant control on channel structure. We review the role of sediment supply in channel structure, and how regional differences in sediment supply and landuse affect stream restoration priorities. We show how stream restoration goals are best understood within a common fluvial geomorphology framework defined by sediment supply, storage, and transport. Landuse impacts in geologically young landscapes with high sediment yields (e.g., coastal British Columbia) typically result in loss of instream wood and accelerated sediment inputs from bank erosion, logging roads, hillslopes and gullies. In contrast, northern Sweden and Finland are landscapes with naturally low sediment yields caused by low relief, resistant bedrock, and abundant mainstem lakes that act as sediment traps. Landuse impacts involved extensive channel narrowing, removal of obstructions, and bank armouring with boulders to facilitate timber floating, thereby reducing sediment supply from bank erosion while increasing export through higher channel velocities. These contrasting landuse impacts have pushed stream channels in opposite directions (aggradation versus degradation) within a phase-space defined by sediment transport and supply. Restoration in coastal British Columbia has focused on reducing sediment supply (through bank and hillslope stabilization) and restoring wood inputs. In contrast, restoration in northern Fennoscandia (Sweden and Finland) has focused on channel widening and removal of bank-armouring boulders to increase sediment supply and retention. These contrasting restoration priorities illustrate the consequences of divergent regional landuse impacts on sediment supply, and the utility of planning restoration activities within a mechanistic sediment supply-transport framework.     

Backfilling canals to mitigate Wetland dredging in Louisiana coastal marshes

Returning canal spoil banks into canals, or backfilling, is used in Louisiana marshes to mitigate damage caused by dredging for oil and gas extraction. We evaluated 33 canals backfilled through July 1984 to assess the success of habitat restoration. We determined restoration success by examining canal depth, vegetation recolonization, and regraded spoil bank soils after backfilling. Restoration success depended on: marsh type, canal location, canal age, marsh soil characteristics, the presence or absence of a plug at the canal mouth, whether mitigation was on- or off-site, and dredge operator performance.Backfilling reduced median canal depth from 2.4 to 1.1 m, restored marsh vegetation on the backfilled spoil bank, but did not restore emergent marsh vegetation in the canal because of the lack of sufficient spoil material to fill the canal and time. Median percentage of cover of marsh vegetation on the canal spoil banks was 51.6%. Median percentage of cover in the canal was 0.7%. The organic matter and water content of spoil bank soils were restored to values intermediate between spoil bank levels and predredging marsh conditions.The average percentage of cover of marsh vegetation on backfilled spoil banks was highest in intermediate marshes (68.6%) and lowest in fresh (34.7%) and salt marshes (33.9%). Average canal depth was greatest in intermediate marshes (1.50 m) and least in fresh marshes (0.85 m). Canals backfilled in the Chenier Plain of western Louisiana were shallower (average depth = 0.61 m) than in the eastern Deltaic Plain (mean depth range = 1.08 to 1.30 m), probably because of differences in sediment type, lower subsidence rate, and lower tidal exchange in the Chenier Plain. Canals backfilled in marshes with more organic soils were deeper, probably as a result of greater loss of spoil volume caused by oxidation of soil organic matter. Canals ten or more years old at the time of backfilling had shallower depths after backfilling. Depths varied widely among canals backfilled within ten years of dredging. Canal size showed no relationship to canal depth or amount of vegetation reestablished. Plugged canals contained more marsh reestablished in the canal and much greater chance of colonization by submerged aquatic vegetation compared with unplugged canals. Dredge operator skill was important in leveling spoil banks to allow vegetation reestablishment. Wide variation in dredge performance led to differing success of vegetation restoration.Complete reestablishment of the vegetation was not a necessary condition for successful restoration. In addition to providing vegetation reestablishment, backfilling canals resulted in shallow water areas with higher habitat value for benthos, fish, and waterfowl than unfilled canals. Spoil bank removal also may help restore water flow patterns over the marsh surface. Increased backfilling for wetland mitigation and restoration is recommended.     

ENVIRONMENTAL DESIGNS FOR STREAMBANK PROTECTION PROJECTS1

ABSTRACT: Streambank protection projects are intended to prevent streambank erosion, thereby preventing streambank failure and maintaining a desirable channel alignment. Streambank erosion is a natural process of unaltered, dynamic river systems, and protection projects seek to impose stability on this natural system. The environmental impacts of such projects are primarily changes to terrestrial and aquatic habitats and to aesthetics. Adverse environmental impacts have been minimized and enhancement of existing habitat and aesthetics have been achieved through the development of new, innovative designs or modifications to existing designs and through use of construction and maintenance practices that promote habitat and aesthetics. Designs based on channel flow characteristics, e.g., revetments using a variety of structural materials, can result in preservation of wildlife habitat by reducing the use of structural protection by matching the erosion potential of flow at the bank with the protection capability of the materials used. Designs based on streambed stabilization prevent bank failure caused by bank undermining, result in preservation or establishment of streamside vegetation, and enhance aesthetics. Protection schemes that manage and preserve floodplains, berms, and riparian areas preserve the natural condition of the floodplain area. Designs based on deflection of erosive flows, e.g., dikes, minimize disturbance to the bank vegetation and create low-velocity aquatic habitats. Use of vegetation for bank protection is most effective when used in combination with structural components. Construction and maintenance practices can be scheduled and modified to minimize impacts to floodplain areas and to enhance wildlife habitat while preserving the integrity of the protection structure.     

Integration of SWAT and HSPF for Simulation of Sediment Sources in Legacy Sediment‐Impacted Agricultural Watersheds

A total maximum daily load for the Chesapeake Bay requires reduction in pollutant load from sources within the Bay watersheds. The Conestoga River watershed has been identified as a major source of sediment load to the Bay. Upland loads of sediment from agriculture are a concern; however, a large proportion of the sediment load in the Conestoga River has been linked to scour of legacy sediment associated with historic millpond sites. Clarifying this distinction and identifying specific segments associated with upland vs. channel sources has important implications for future management. In order to address this important question, we combined the strengths of two widely accepted watershed management models — Soil and Water Assessment Tool (SWAT) for upland agricultural processes, and Hydrologic Simulation Program FORTRAN (HSPF) for instream fate and transport — to create a novel linked modeling system to predict sediment loading from critical sources in the watershed including upland and channel sources, and to aid in targeted implementation of management practices. The model indicates approximately 66% of the total sediment load is derived from instream sources, in agreement with other studies in the region and can be used to support identification of these channel source segments vs. upland source segments, further improving targeted management. The innovated linked SWAT‐HSPF model implemented in this study is useful for other watersheds where both upland agriculture and instream processes are important sources of sediment load.     

Critique of present wetlands mitigation policies in the united states based on an analysis of past restoration projects in San Francisco Bay

A detailed evaluation of past wetland restoration projects in San Francisco Bay was undertaken to determine their present status and degree of success. Many of the projects never reached the level of success purported and others have been plagued by serious problems. On the basis of these findings, it is debatable whether any sites in San Francisco Bay can be described as completed, active, or successful restoration projects at present. In spite of these limited accomplishments, wetland creation and restoration have been adopted in the coastal permit process as mitigation to offset environmental damage or loss of habitat. However, because the technology is still largely experimental, there is no guarantee that man-made wetlands will persist as permanent substitutes for sacrificed natural habitats. Existing permit policies should be reanalyzed to insure that they actually succeed in safeguarding diminishing wetlands resources rather than bartering them away for questionable habitat substitutes. Coastal managers must be more specific about project requirements and goals before approval is granted. Continued research on a regional basis is needed to advance marsh establishment techniques into a proven technology. In the meantime, policies encouraging or allowing quid pro quo exchanges of natural wetlands with man-made replacements should proceed with caution. The technology and management policies used at present are many steps ahead of the needed supporting documentation.     

Mitigating the Effects of Landscape Development on Streams in Urbanizing Watersheds

This collaborative study examined urbanization and impacts on area streams while using the best available sediment and erosion control (S&EC) practices in developing watersheds in Maryland, United States. During conversion of the agricultural and forested watersheds to urban land use, land surface topography was graded and vegetation was removed creating a high potential for sediment generation and release during storm events. The currently best available S&EC facilities were used during the development process to mitigate storm runoff water quality, quantity, and timing before entering area streams. Detailed Geographic Information System (GIS) maps were created to visualize changing land use and S&EC practices, five temporal collections of LiDAR (light detection and ranging) imagery were used to map the changing landscape topography, and streamflow, physical geomorphology, and habitat data were used to assess the ability of the S&EC facilities to protect receiving streams during development. Despite the use of the best available S&EC facilities, receiving streams experienced altered flow, geomorphology, and decreased biotic community health. These impacts on small streams during watershed development affect sediment and nutrient loads to larger downstream aquatic ecosystems such as the Chesapeake Bay.     

Application of SPARROW Modeling to Understanding Contaminant Fate and Transport from Uplands to Streams

Understanding spatial variability in contaminant fate and transport is critical to efficient regional water‐quality restoration. An approach to capitalize on previously calibrated spatially referenced regression (SPARROW) models to improve the understanding of contaminant fate and transport was developed and applied to the case of nitrogen in the 166,000 km² Chesapeake Bay watershed. A continuous function of four hydrogeologic, soil, and other landscape properties significant (α = 0.10) to nitrogen transport from uplands to streams was evaluated and compared among each of the more than 80,000 individual catchments (mean area, 2.1 km²) in the watershed. Budgets (including inputs, losses or net change in storage in uplands and stream corridors, and delivery to tidal waters) were also estimated for nitrogen applied to these catchments from selected upland sources. Most (81%) of such inputs are removed, retained, or otherwise processed in uplands rather than transported to surface waters. Combining SPARROW results with previous budget estimates suggests 55% of this processing is attributable to denitrification, 23% to crop or timber harvest, and 6% to volatilization. Remaining upland inputs represent a net annual increase in landscape storage in soils or biomass exceeding 10 kg per hectare in some areas. Such insights are important for planning watershed restoration and for improving future watershed models.     

Mapping Three‐Dimensional Water‐Quality Data in the Chesapeake Bay Using Geostatistics1

Abstract:  Data interpretation and visualization software tools with geostatistical capabilities were adapted, customized, and tested to assist the Chesapeake Bay Program in improving its water‐quality modeling protocols. Tools were required to interpolate, map, and visualize three‐dimensional (3D) water‐quality data, with the capability to determine estimation errors. Components of the software, originally developed for ground‐water modeling, were customized for application in estuaries. Additional software components were developed for retrieval, and for pre‐ and post‐ processing of data. The Chesapeake Bay Program uses the 3D mapped data for input to the Bay water‐quality model that projects the future health of the Bay and its tidal tributary system. In determining water‐quality attainment criteria, 3D kriging estimation errors are needed as a statistical measure of uncertainty. Furthermore, given the high cost of installing and operating new monitoring stations, geostatistical techniques can assist the Chesapeake Bay Program in the identification of suitable data collection locations. Following the evaluation, selection, and development of the software components phase, 3D ordinary kriging techniques with directional semi‐variograms to account for anisotropy were successfully demonstrated for mapping 3D fixed station water‐quality data, such as dissolved oxygen and salinity. Additionally, an improved delineation tool was implemented to simulate the upper and lower pycnocline boundary surfaces allowing the segregation of the interpolated 3D data into three separate zones for a better characterization of the pycnocline layer.