Time scales for pollution assessment
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 also often compared with the pollutant source time scales or biogeochemical processes to evaluate the relative importance of physical and water quality processes (Monsen et al., 2002; Salomon and Pommepuy, 1990). They are also used to understand population dynamics (Monsen et al., 2002) and to serve as indicators to classify and compare estuarine systems (Jay, 1994).
The first correlation between retention time scales and water quality was proposed in the classical empirical model of lake eutrophication of Vollenweider (1976), which describes algal biomass as a function of the phosphorus loading rate scaled hy 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, for a detailed list of applications).
In this article, three transport time scales commonly used to quantify the retention of water or pollutants, - flushing time, residence time and age - and their scope of applicability are described. The methodologies and tools to quantify each of these indicators are also presented.
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): 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 whole-system 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, hard-struture 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" (once-through 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 re-entrant 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 avalability 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
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
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 steady-state with a sinusoidal tidal signal, and
- the system is fully flushed in a single tidal cycle.
The last assumption can be alleviatted by a modified formulation (using (1-b)*P in the denominator) to account for an incomplete mixing, where b is the return flow factor (0-1).This method tends to underestimate the flushing time due to the steady-state hypothesis.
Continuosly stirred tank reactor (CSTR) method
where C is the concentration of a pollutant at time t, due to an instantaneous load at time t0 that leads to the initial concentration C0. 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 box-models, particle models and concentration models. Applications based on box-models are the simplest to set-up, requiring few information on the system and small computational resources. The residence times 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 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 is generally very complex in coastal systems. Particle models are by far the most popular approach to compute residence times. 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 Runge-Kutta 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), which may require very large computational resources.
Models based on concentration (transport models) have also been used for residence times calculations, based on two and three dimensional applications. For these models, residence times can be defined as the time necessary to reduce the initial pollutant concentration by 1/e (Sandery and Kampf, 2005) or to reduce total mass by 1/e (Wang et al., 2004). This approach can readily account for water quality processes, but the accuracy of the residence times 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 times. The readers are referred to Shen and Haas, 2004 and Kennedy et al., 2006 for details on the application of concentration and particle models for the calculation of age.
- 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:607-624.
- Bolin, B., A. Rodhe, 1973. A note on the concepts of age distribution and transit time in natural reservoirs. Tellus, 2558-62.
- 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. 107-117.
- 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. 177-206, 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), 328-340.
- Kennedy, M.G., D.P. Ahlfeld, D.P. Schmidt, J.E. Tobiason, 2006. Three-dimensional modeling for estimation of hydraulic retention time in a reservoir, Journal of Environmental Engineering, vol. 132/9, 976-984.
- Jay D., 1994. Residence times, box models and shear fluxes in tidal channel flows, Changes in Fluxes in Estuaries, Dyer and Orth (eds.), 3-12, 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), 1545-1553.
- Oliveira, A., A.M. Baptista, 1997. Diagnostic modeling of residence times in estuaries. Water Resources Research, 33(8), 1935-1946.
- Salomon, J.C., M. Pomepuy, 1990. Mathematical modeling of bacterial contaminationof the Morlaix estuary (France) Water Research, 24(8), 983-994.
- Sandery, P.A., J. Kampf, 2005. Winter-Spring flushing of Bass Strait, South-Eastern Australia: a numerical modeling study, Estuarine, Coastal and Shelf Science,63, 23-31.
- Shen, J., L. Haas, 2004. Calculating age and residence time in the tidal York River using three-dimensional model experiments, Estuarine, Coastal and Shelf Science, 61, 449-461.
- Vollenweider, R.A., 1976. Advances in defining critical loading levels of phosphorus in lake eutrophication. Mem. Ist. Ital. Idrobiol. 33, 53-83.
- Wang, C-F, M-H Hsu, A.Y. Kuo , 2004. Residence time of the Danshuei River estuary, Taiwan, Estuarine, Coastal and Shelf Science, 60, 381-393.
- 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), 149-191.