Temporal coupling of subsurface and surface soil CO2 fluxes: Insights from a nonsteady state model and cross-wavelet coherence analysis

Inferences about subsurface CO2 fluxes often rely on surface soil respiration (Rsoil) estimates because directly measuring subsurface microbial and root respiration (collectively, CO2 production, STotal) is difficult. To evaluate how well Rsoil serves as a proxy for STotal, we applied the nonsteady state DEconvolution of Temporally varying Ecosystem Carbon componenTs model (0.01-m vertical resolution), using 6-hourly data from a Wyoming grassland, in six simulations that cross three soil types (clay, sandy loam, and sandy) with two depth distributions of subsurface biota. We used cross-wavelet coherence analysis to examine temporal coherence (localized linear correlation) and offsets (lags) between STotal and Rsoil and fluxes and drivers (e.g., soil temperature and moisture). Cross-wavelet coherence revealed higher coherence between fluxes and drivers than linear regressions between concurrent variables. Soil texture and moisture exerted the strongest controls over coherence between CO2 fluxes. Coherence between CO2 fluxes in all soil types was strong at short (~1 day) and long periods (>8 days), but soil type controlled lags, and rainfall events decoupled the fluxes at periods of 1?8 days for several days in sandy soil, up to 1 week in sandy loam, and for a month or more in clay soil. Concentrating root and microbial biomass nearer the surface decreased lags in all soil types and increased coherence up to 10% in clay soil. The assumption of high temporal coherence between Rsoil and STotal is likely valid in dry, sandy soil, but may lead to underestimates of short-term STotal in semiarid grasslands with fine-grained and/or wet soil.