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

Abstract

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.