Application of causal discovery of factors driving dissolved oxygen in estuarine environments.
Dissolved oxygen (DO) concentrations in estuarine bottom waters are a manifestation of multiple, interacting physical and biogeochemical processes, yet identifying their independent contributions remains challenging. Here, we analyze monthly water quality monitoring data from eight stations across Long Island Sound from 1994 to 2022 using a causal discovery framework (PCMCI+) and transformation of forcing variables. Our goal is to identify and isolate variables that causally influence bottom DO and improve predictive models by minimizing overfitting and multicollinearity. PCMCI+ reveals surface-layer temperature as the most important and consistent negative driver of bottom DO, followed by stratification. Wind events exhibit only brief relief by advection and mixing, while river discharge shows no direct causal link to DO, making it less influential than previously thought. Biogeochemical variables, including chlorophyll-a (Chl-a), nitrate and nitrite, and particulate carbon, influence DO through both contemporaneous and time-lagged pathways, often with signs that shift depending on the process. The derived models were evaluated by comparing skill scores, mean squared error, and Akaike Information Criterion. Both model types perform well, with coefficient of determination values exceeding 0.90 at multiple stations using only 3-5 predictors. Our analysis reveals that the best causal predictors are surface-layer temperature, stratification, Chl-a, and particle carbon. This approach provides a scalable framework for improving prediction models and understanding the mechanistic links that control the seasonal variability of DO in estuarine systems.