Hydrograph structure informed calibration in the frequency domain with time localization

Abstract:

Complex models with large number of parameters are commonly used to estimate sediment yields and predict changes in sediment loads as a result of changes in management or conservation practice at large watershed (>2000 km2) scales. As sediment yield is a strongly non-linear function that responds to channel (peak or mean) velocity or flow depth, it is critical to accurately represent flows. The process of calibration in such models (e.g., SWAT) generally involves the adjustment of several parameters to obtain better estimates of goodness of fit metrics such as Nash Sutcliff Efficiency (NSE). However, such indicators only provide a global view of model performance, potentially obscuring accuracy of the timing or magnitude of specific flows of interest. We describe an approach for streamflow calibration that will greatly reduce the black-box nature of calibration, when response from a parameter adjustment is not clearly known. Fourier Transform or the Short Term Fourier Transform could be used to characterize model performance in the frequency domain as well, however, the ambiguity of a Fourier transform with regards to time localization renders its implementation in a model calibration setting rather useless. Brief and sudden changes (e.g. stream flow peaks) in signals carry the most interesting information from parameter adjustments, which are completely lost in the transform without time localization. Wavelet transform captures the frequency component in the signal without compromising time and is applied to contrast changes in signal response to parameter adjustments. Here we employ the mother wavelet called the Mexican hat wavelet and apply a Continuous Wavelet Transform to understand the signal in the frequency domain. Further, with the use of the cross-wavelet spectrum we examine the relationship between the two signals (prior or post parameter adjustment) in the time-scale plane (e.g., lower scales correspond to higher frequencies). The non-stationarity of the streamflow signal does not hinder this assessment and regions of change called boundaries of influence (seasons or time when such change occurs in the hydrograph) for each parameter are delineated. In addition, we can discover the structural component of the signal (e.g., shifts or amplitude change) that has changed.