Boundary impulse response functions in a century-long eddying global ocean simulation

Results are presented from a century-long 1/10A degrees global ocean simulation that included a suite of age-related passive tracers. In particular, an ensemble of five global Boundary Impulse Response functions (BIRs, which are statistically related to the more fundamental Transit Time Distributions, TTDs) was included to quantify the character of the TTD when mesoscale eddies are explicitly simulated rather than parameterized. We also seek to characterize the level of variability in water mass ventilation timescales arising from eddy motions. The statistics of the BIR timeseries are described, and it is shown that the greatest variability occurs at early times, followed by a remarkable conformity between ensemble members at longer timescales. The statistics of the first moment of the BIRs are presented, and the upper-ocean spatial distribution of the standard deviation of the first moment of the BIRs discussed. It is shown that variations in the BIR first moment with respect to the ensemble average are typically only a few percent, and that the variability slightly decreases with increasing ensemble size, implying that only a few ensemble members may be necessary for a reasonable estimate of the TTD. The completeness of the estimated TTD, i.e., the degree to which the century long BIRs capture the range of global ocean ventilation timescales is discussed, and the potential for extrapolation of the BIR to longer times is briefly explored. Several regional BIRs were also simulated in order to quantify the relative abundance of fluid parcels that originate in specific geographical locations.