Effects of Short Time Variation in the River Discharge on the Salt Wedge Intrusion in the Yura Estuary, a Micro Tidal Estuary, Japan ()
1. Introduction
Estuaries are formed around the river mouths, where lighter fresh water meets denser sea water. There is a great deal of variety in mixing processes and salinity distributions, affected by the balance between tides and river discharges [1]. For instance, during periods of strong tidal forces and low runoff, the Columbia River Estuary becomes well mixed, where salinity structure is vertically homogeneous [2]. However, estuaries along the Mediterranean coast are known as highly stratified estuaries, where sharp vertical salinity gradients are observed because of their small tidal ranges [3]. Reference [3] categorized highly stratified estuaries into two types based on the river discharge, estuarine cross section, and tidal range: 1) the salt wedge regime is established during low discharge periods, whereas it is washed away during high river discharge; 2) The salt wedge regime is established during high river discharge periods, but the water is partially mixed and the obvious salt wedge is not observed during low discharge periods.
Fresh water is utilized for irrigation and drinking purposes in many rivers. Therefore, salt wedge intrusions have been studied for management of estuarine water quality [4]. Reference [5] used a two-layer hydrodynamic model to evaluate effects of the change in the river discharge on the salt wedge dynamics of the Ebro River estuary, which belongs to the first type of highly stratified estuary according to Ibañez’s classification [3]. Model results suggest that the decrease in the river discharge enlarges and prolongs the salt wedge. However, these studies dealt with long term (more than few months) trends, and there are few studies on the short term variation in the slat wedge intrusion in micro tidal estuaries.
Compared to American and Southeastern Asian continental countries, Japan is characterized by a relatively narrow strip of mountainous land. This implies that the root of its rivers bear a strong altitude gradient from the origin to the mouth, leading to short residence time of rainwater, as well as strong and rapid variations in river discharge [6]. These sudden changes in turn impact estuarine dynamics. This study was conducted in the Yura Estuary, the first type of micro tidal estuary according to Ibañez’s classification, which has the salt wedge in summer. A previous study revealed that the salt wedge in the Yura River changes seasonally and is affected by variations in river discharge and sea level [7]. It is also shown that a chlorophyll maximum was generated around a halocline during summer salt wedge regime, indicating that the salt wedge plays an important role in the ecology of the Yura Estuary. During summer, discharge of Yura River also varies drastically over short time scales due to heavy rains. In this study, therefore, we evaluate the effect of short temporal variations in the river discharge on the salt wedge intrusion in the Yura Estuary using a three dimensional numerical model. First, the model was developed to reproduce the real salt wedge intrusion, which was validated with field observations [7]. The model was then applied to an experiment that examines the response of the salt wedge to the rapid change in the river discharge in the Yura Estuary as representative of micro tidal estuaries.
2. Materials and Methods
2.1. Study Site
The Yura River is a 146 km long and flows into Wakasa bay, Sea of Japan (Figure 1). The catchment area is about 1880 km2. The mean annual discharge measured at Fukuchiyama is ~51.5 m3·s–1 in the last 54 years. During winter, northwesterly winds off the Asian continent bring heavy snow on Sea of Japan side. The average river discharge is large (~ 55 m3·s–1) in winter and early spring because of snow melting. However, the average discharge decreases to 21 m3·s–1, because summer to autumn is relatively dry season. The catchment area is still subjected to sudden heavy rain in summer so that the discharges sometimes reach over 100 m3·s–1. Therefore the salt wedge intrusion varies seasonally [7], and it is expected to change suddenly in summer.
The typical tidal range of the Yura Estuary is less than 0.5 m, and it is thus classified as a micro tidal estuary. Therefore the effects of tidal currents on the physical conditions are negligible. The Yura Estuary enters the salt wedge regime during low discharge periods, with a maximum salt wedge length of 18 km from the river mouth [7]. The estuarine width is 100 m at the upper part and increases to 500 m at the mouth. Average water depth is 3.7 m and there is a deep part (~10 m) at 5 - 6 km from the mouth.
2.2. Numerical Model
2.2.1. Model Design
The salt wedge dynamics were calculated by Delft3DFlow [8]. This model is a three dimensional hydrodynamic model developed by Delft Hydraulics. The model domain extends 25 km from the upper part of the Yura Estuary to 3 km offshore from the river mouth (Figures 1 and 2). Curvilinear grids were applied for the width of river part domain so that the model domain corresponds to the real topography. The original bathymetry data of the upper Yura estuarine part were obtained from the Ministry of Land, Infrastructure, Transport and Tourism,
Figure 1. Location and bathymetry of the Yura Estuary.
Figure 2. Model domain and curvilinear grids used in the simulations.
and the coastal area from the Japan Coast Guard. The model system has 14 σ-levels in the vertical; five layers of 3%, three of 5%, four of 15%, two of 5% from the surface to bottom. The horizontal grid scale ranges from 15 m by 200 m in the river to 150 m by 200 m in the sea. The vertical eddy viscosity and diffusivity are calculated by a k-ε model. The horizontal eddy viscosity and diffusivity are calculated by a 3D-turbulence closure model [8]. Both the river end (y = 0 km) and the sea side (y = 28 km) are open boundaries. Monthly observations of salinity and temperature [7] are used to force the open boundary. Daily averaged sea level at Maizuru and daily average river discharge at Fukuchiyama are also used for the sea side open boundary and river side open boundary, respectively (Figures 2 and 3).
The salt wedge intrusion was simulated for two years; from April 2006 to March 2008. The model was spun up and reached steady state after one month, and then real calculation started with observed boundary conditions. Calculated salinity distribution and salt wedge length by the model were compared with those obtained by field observations [7]. In this study, salt wedge length is defined as the distance from the river mouth to the tip of salt wedge, which salinity is 10 at the bottom.
2.2.2. Responses of the Salt Wedge to the Flood
In order to study the response of the salt wedge to variations in river discharge, the river discharge condition was changed while the boundary conditions of temperature, salinity and sea level were kept constant. With the intention of reproduce the summer salinity distribution, the
Figure 3. Time series of (a) Sea level in Maizurum, and (b) River discharge in Fukuchiyama.
Figure 4. Time series of river discharge with each scenarios of Q0.
salinity and temperature of the river open boundary were set to 0 and 15.4˚C respectively, and the sea side boundary conditions were set to 33.4, 17.6˚C and 0.2 m for salinity, temperature and sea level condition respectively. The river discharge (Q) of a flood event was empirically calculated as follows;
, (1)
where Q0 represents maximum discharge (m3·s–1) and T is date from the flood event (day). Equation (1) is derived from typical flood event conditions and river flow decreases gradually without any increasing of the river discharge after the flood.
First, the model was run for one month with a constant river discharge (15 m3·s–1) in order to obtain a steady state as representative for dry season. The second stage of the experiment was forced with the steady state as the initial condition and then the river discharge was changed according to Equation (1). Eight flood scenarios of Q0 equals 50, 100, 150, 200, 250, 500, 1000 and 1500 were conducted. The minimum river discharge was set to 15 m3·s–1 Examples of variation of discharge are 50, 150, 250, 500 and 1000 m3·s–1 are shown in Figure 4. This shows that the difference of Q among the scenarios was large in the early stages, but reduced to 15 m3·s–1 20 days after the flood regardless of Q0. This is consistent with the observed river discharge at Fukuchiyama. The response of the salt wedge was evaluated by comparing the minimum salt wedge length and replacing time, which represents the time that takes salt wedge to re-intrude from the minimum position to 12.5 km (T12.5) and 15 km (T15). In this study, the minimum salt wedge length is defined as the minimum position during the experiment.
3. Results
3.1. Reproduction of the Observed Salinity Distribution
Figure 5 shows the salinity distributions obtained by field observations and model simulations. The river discharge was low on 11 December 2006 and 22 August 2007, while high on 29 January 2008 and 25 February 2008. The salt wedge intruded into the river, which was strongly stratified under low river discharge conditions (Figures 5(a) and (c)). The upper fresh water layer was thicker in December 2006 than that in August 2007, depending on the river discharge and sea level (Figure 3), although the salt wedge intruded into the river in both dates. The model reproduced the salinity distribution well, showing the salt wedge intrusion through the bottom layer with the strong halocline (Figures 5(b) and (d)). In contrast, the river was occupied by the fresh water and the salt wedge was washed away on 29 January and 25 February 2008 (Figures 5(f) and (h)). Slightly