Time scales for pollution assessment
Time scales for pollution assessment can be considered using three commonly used methods: flushing time, residence time and age. In this article the methodologies and tools to quantify each of these methods is explained as well as the applicability of each method.
Inhoud
Introduction
Transport time scales are useful tools to quantify the importance of hydrodynamic processes in the transport and fate of pollutants in coastal and estuarine water systems. Indeed, the water quality of a system depends crucially on the retention of pollutants within the system and its ability to flush them out. Transport time scales are often compared with the pollutant source time scales or biogeochemical processes to evaluate the relative importance of physical and water quality processes. ^{[1]} ^{[2]} They are also used to understand population dynamics ^{[1]} and to serve as indicators to classify and compare estuarine systems.^{[3]}
The first correlation between retention time scales and water quality was proposed in the classical empirical model of lake eutrophication of Vollenweider (1976),^{[4]} which describes algal biomass as a function of the phosphorus loading rate scaled by the hydraulic residence time. Many other applications followed, from heavy metals problems, mineralization rates of organic matter and primary production (the reader is refered to Monsen et al. 2002 ^{[1]} for a detailed list of applications).
Terminology
Many retention time scales can be found in the literature, often with distinct definitions for the same concept (see for instance Oliveira and Baptista, 1997, Shen and Haas, 2004 for more details). Here, three fundamentally different concepts are presented, following the classification presented in Monsen et al. (2002) ^{[1]} :flushing time, residence time and age.
Flushing time has been defined as the "time to replace the freshwater volume of the estuary by the total freshwater input flux (river, discharges, rainfall,...)" (Officer and Kester, 1991) or, in a more general way, as "the ratio of the mass of a scalar in a reservoir to the rate of renewal of the scalar" (Geyer et al., 2000). This transport time scale is a wholesystem indicator of the renewal capacity, but does not allow for the distinction between several forcing mechanisms (e.g., the influence of tidal motion to flush out the system) or the spatial and time variability of the renewal capacity.
Residence time has been defined as "the time it takes for any waterparcel of the sample to leave the lagoon through its outlet to the sea" (Dronkers and Zimmerman, 1982). They can be used to analyse the spatial variation of flushing properties and the variation of renewal for different environmental conditions (effect of tidal amplitude and phase, relative importance of different forcings – waves, currents, …). These detailed properties makes residence times very useful for comparative analyses of the effects of several engineering interventions (dredging, hardstruture building), and also to characterize changes in the system’s contaminant inputs (changes in the nature, location and frequency of the sources of contaminants). The specific way to define "the time to leave the system" can also lead to different concepts of residence times which can be very important for water quality analyses. If residence times are defined as the "time for a water parcel to leave the system once" (oncethrough residence times, Oliveira and Baptista, 1997), the concept is very useful to characterize the flushing of pollutants that are significantly altered once outside the system (for instance due to strong variations in salinity and/or temperature). An opposite definition is "the time for a water parcel to leave the system without returning at a later tide" (denoted reentrant residence times in Oliveira and Baptista, 1997), which is useful for the analysis of the retention of conservative tracers in a system. Finally, residence times can also be defined as the "time spent in the domain of interest" (denoted as exposure time in Deleersnjider and Delhez, 2005). This definition is particularly important in coastal pollution since it quantifies the time of exposure of a system to a specific contaminant.
Using a reservoir concept, the age of a contaminant can be defined as "the time elapsed since it entered the system" (Bolin and Rodhe, 1973). Zimmerman, (1976) proposed a definition that explicitly accounts for the spatial variability of this time scale: “age of a water parcel is the time elapsed since the particle departed the region where its age is zero”. A general theory for age, based on tracer concentrations, can be found in Deleersnijder et al. (2001). Age is the complement of residence time and can be used to understand the pathways of contaminants and organisms within the system, in particular in the influence of external sources/sinks of contaminants to coastal systems.
Methodologies and tools
Each transport time scale can be estimated by several methodologies based on models or field data. The choice of the methodology depends on the specific goals and on the avaiability of data to characterize the system.
The Flushing Time (FT) is the least demanding parameter. It can be calculated in several ways, using either data or numerical model results:
Freshwater faction method
$ FT = \frac{s_0s}{s_0}*\frac{V}{Q} $
where S0 – ocean salinity, s – average salinity in domain, Q – freshwater input This method assumes steady state conditions and can only be applied to systems with a significant freshwater input. For complex systems, it requires a considerable amount of salinity data or the use of a numerical model.
Tidal prism method
$ FT = V*\frac{T}{P} $
where V  water volume at high tide, T  tidal period, P  tidal prism of a representative flood tide. This method only requires basin geometry and tidal range information. It assumes that
 the system is well mixed,
 tidal flow is the dominant flushing mechanism,
 the system is at steadystate with a sinusoidal tidal signal, and
 the system is fully flushed in a single tidal cycle.
The last assumption can be alleviated by a modified formulation (using (1b)*P in the denominator) to account for an incomplete mixing, where b is the return flow factor (01).This method tends to underestimate the flushing time due to the steadystate hypothesis.
Continuously stirred tank reactor (CSTR) method
$ C(t) = C_0*e^{\frac{t}{FT}} $
where C is the concentration of a pollutant at time t, due to an instantaneous load at time t_{0} that leads to the initial concentration C_{0}. This method assumes that
 no further mass is introduced in the system,
 flow and volume in the CSTR are constant in time and
 instantaneous and complete mixing of the pollutant in the system.
Residence times are typically evaluated with numerical model results. Applications have been done with models of increasing complexity, including boxmodels, particle models and concentration models.
 Boxmodels are the simplest to setup, requiring little information on the system and small computational resources. The residence time of the pollutants in each box can be computed based on the knowledge of the geometry of the system, the hydrodynamics and the pollution loads (the reader is referred to Hagy et al. (2000) ^{[5]} for a detailed application of this approach). Calculations of residence times based on this method are limited by the capability of box models to solve the underlying hydrodynamics, which are generally very complex in coastal systems.
 Particle models are by far the most popular approach to compute residence time. Residence times are computed based on the release of large numbers of particles, scattered throughout the domain of interest, at several release times within the tidal cycle and for different tidal amplitudes. This approach generally accounts for both advection and diffusion processes, through the use of Euler or RungeKutta methods for the advective term and random walk methods for the diffusive term. The advantages of this approach are its ability to handle the spatial and temporal variability of residence times, and its applicability to several types of tracers and multiple sources of contamination. Limitations include the limited hability to solve complex water quality processes using particles and the need to use very large numbers of particles for comprehensive studies (e.g. O(500,000) in the Gulf of Maine, Bilgili et al. 2005) ^{[6]}, which may require very large computational resources.
 Concentration models (transport models) have also been used for residence time calculations, based on two and three dimensional applications. For these models, residence time can be defined as the time necessary to reduce the initial pollutant concentration by 1/e ^{[7]} or to reduce total mass by 1/e ^{[8]}. This approach can readily account for water quality processes, but the accuracy of the residence time can be limited by numerical errors (mass imbalances, numerical diffusion,..) and cannot be applied for several simultaneous sources of the same contaminant.
The methodologies to compute age are the same as those used for residence time. The readers are referred to Shen and Haas, (2004)^{[9]} and Kennedy et al. (2006)^{[10]} for details on the application of concentration and particle models for the calculation of age.
References
 ↑ ^{1,0} ^{1,1} ^{1,2} ^{1,3} Monsen N.E., Cloern J.E. and Lucas L.V (2002). A comment on the use of flushing time, residence time, and age as transport time scales. Limnology and Oceanography, 47(5) 15451553.
 ↑ Salomon J.C. and Pomepuy M. (1990). Mathematical modeling of bacterial contaminationof the Morlaix estuary (France). Water Research 24(8) 983994.
 ↑ Jay D. (1994). Residence times, box models and shear fluxes in tidal channel flows, Changes in Fluxes in Estuaries, Dyer and Orth (eds.), 312, Olsen and Olsen, Fredensborg, Denmark.
 ↑ Vollenweider R.A. (1976). Advances in defining critical loading levels of phosphorus in lake eutrophication. Mem. Ist. Ital. Idrobiol. 33 5383.
 ↑ Hagy J.D., Boynton W.R. and Sanford L.P. (2000). Estimation of net physical transport and hydraulic residencetimes for a coastal plain estuary using box models, Estuaries, 23(3) 328340.
 ↑ Bilgili A., Proehl J., Lynch D., Smith K. and Swift R. (2005). Estuary/Ocean Exchange and Tidal Mixing in a Gulf of Maine Estuary: A Lagrangian Modeling Study, Estuarine Coastal and Shelf Science, 65 607624.
 ↑ Sandery P.A. and Kampf J. (2005). WinterSpring flushing of Bass Strait, SouthEastern Australia: a numerical modeling study, Estuarine, Coastal and Shelf Science, 63 2331.
 ↑ Wang CF., Hsu MH. and Kuo A.Y. (2004). Residence time of the Danshuei River estuary, Taiwan, Estuarine, Coastal and Shelf Science, 60 381393.
 ↑ Shen J. and Haas L. (2004). Calculating age and residence time in the tidal York River using threedimensional model experiments, Estuarine, Coastal and Shelf Science, 61 449461.
 ↑ Kennedy M.G., Ahlfeld D.P., Schmidt D.P. and Tobiason J.E. (2006). Threedimensional modeling for estimation of hydraulic retention time in a reservoir, Journal of Environmental Engineering, 132(9) 976984.
 Bilgili, A., J. Proehl, D. Lynch, K. Smith, R. Swift. 2005. Estuary/Ocean Exchange and Tidal Mixing in a Gulf of Maine Estuary: A Lagrangian Modeling Study, Estuarine Coastal and Shelf Science, 65:607624.
 Bolin, B., A. Rodhe, 1973. A note on the concepts of age distribution and transit time in natural reservoirs. Tellus, 255862.
 Dronkers, J. and J.T.F. Zimmerman, 1982. Some principles of mixing in tidal lagoons. Oceanologica Acta. Proceedings of the International Symposium on Coastal Lagoons, p. 107117.
 Geyer, W.R., J.T. Morris, F.G. Pahl, and D.A. Jay, 2000. Interaction between physical processes and ecosystem structure. A comparative approach. p. 177206, In EStuarine Science: a synthetic approach to research and practice, Hobbie, J.E. (ed.)., Island Press.
 Hagy, J.D., W.R. Boynton, L.P. Sanford, 2000. Estimation of net physical transport and hydraulic residencetimes for a coastal plain estuary using box models, Estuaries, 23(3), 328340.
 Kennedy, M.G., D.P. Ahlfeld, D.P. Schmidt, J.E. Tobiason, 2006. Threedimensional modeling for estimation of hydraulic retention time in a reservoir, Journal of Environmental Engineering, vol. 132/9, 976984.
 Jay D., 1994. Residence times, box models and shear fluxes in tidal channel flows, Changes in Fluxes in Estuaries, Dyer and Orth (eds.), 312, Olsen and Olsen, Fredensborg, Denmark.
 Monsen, N.E., J.E.Cloern, L.V. Lucas, 2002. A comment on the use of flushing time, residence time, and age as transport time scales. Limnology and Oceanography, 47(5), 15451553.
 Oliveira, A., A.M. Baptista, 1997. Diagnostic modeling of residence times in estuaries. Water Resources Research, 33(8), 19351946.
 Salomon, J.C., M. Pomepuy, 1990. Mathematical modeling of bacterial contaminationof the Morlaix estuary (France) Water Research, 24(8), 983994.
 Sandery, P.A., J. Kampf, 2005. WinterSpring flushing of Bass Strait, SouthEastern Australia: a numerical modeling study, Estuarine, Coastal and Shelf Science,63, 2331.
 Shen, J., L. Haas, 2004. Calculating age and residence time in the tidal York River using threedimensional model experiments, Estuarine, Coastal and Shelf Science, 61, 449461.
 Vollenweider, R.A., 1976. Advances in defining critical loading levels of phosphorus in lake eutrophication. Mem. Ist. Ital. Idrobiol. 33, 5383.
 Wang, CF, MH Hsu, A.Y. Kuo , 2004. Residence time of the Danshuei River estuary, Taiwan, Estuarine, Coastal and Shelf Science, 60, 381393.
 Zimmerman, J.T.F., 1976.Mixing and flushing of tidal embayments in the westren Dutch Wadden Sea. Part I: Distribution of salinity and calculation of mixing time scales. Netherlands Journal of Sea Research, 10(2), 149191.
Please note that others may also have edited the contents of this article.
