Scaling aspects of river flow routing

Scaling aspects of river flow routing are studied by comparing two flow routing schemes, one designed for use in coupled general circulation models (GCMs) and operated at large spatial scales (～350 km), and the other designed for use in typical hydrological applications at small spatial scales (～25 km). The same runoff data are used as input into the two routing schemes, and comparisons are made between mean annual, mean monthly and daily streamflow simulated at four locations within the Mackenzie River Basin. The results suggest that for the purpose of realistically modelling monthly streamflow at the mouth of the rivers in GCMs, flow routing at large spatial scales gives similar results. However, the amplitude of the annual streamflow cycle is slightly but characteristically larger, when routing is performed at large spatial scales. Flow routing at large spatial scales also results in overestimation of high flows, while low flows are underestimated. Copyright © 2001 John Wiley & Sons, Ltd.     

A physically based distributed subsurface–surface flow dynamics model for forested mountainous catchments

This study was designed to develop a physically based hydrological model to describe the hydrological processes within forested mountainous river basins. The model describes the relationships between hydrological fluxes and catchment characteristics that are influenced by topography and land cover. Hydrological processes representative of temperate basins in steep terrain that are incorporated in the model include intercepted rainfall, evaporation, transpiration, infiltration into macropores, partitioning between preferential flow and soil matrix flow, percolation, capillary rise, surface flow (saturation‐excess and return flow), subsurface flow (preferential subsurface flow and baseflow) and spatial water‐table dynamics. The soil–vegetation–atmosphere transfer scheme used was the single‐layer Penman–Monteith model, although a two‐layer model was also provided. The catchment characteristics include topography (elevation, topographic indices), slope and contributing area, where a digital elevation model provided flow direction on the steepest gradient flow path. The hydrological fluxes and catchment characteristics are modelled based on the variable source‐area concept, which defines the dynamics of the watershed response. Flow generated on land for each sub‐basin is routed to the river channel by a kinematic wave model. In the river channel, the combined flows from sub‐basins are routed by the Muskingum–Cunge model to the river outlet; these comprise inputs to the river downstream. The model was applied to the Hikimi river basin in Japan. Spatial decadal values of the normalized difference vegetation index and leaf area index were used for the yearly simulations. Results were satisfactory, as indicated by model efficiency criteria, and analysis showed that the rainfall input is not representative of the orographic lifting induced rainfall in the mountainous Hikimi river basin. Also, a simple representation of the effects of preferential flow within the soil matrix flow has a slight significance for soil moisture status, but is insignificant for river flow estimations. Copyright © 2005 John Wiley & Sons, Ltd.     

Simulating the three-dimensional circulation and hydrography of Halifax Harbour using a multi-nested coastal ocean circulation model

Halifax Harbour is located on the Atlantic coast of Nova Scotia, Canada. It is one of the world’s largest, ice-free natural harbours and of great economic importance to the region. A good understanding of the physical processes controlling tides, flooding, transport and dispersion, and hydrographic variability is required for pollution control and sustainable development of the Harbour. For the first time, a multi-nested, finite difference coastal ocean circulation model is used to reconstruct the three-dimensional circulation and hydrography of the Harbour and its variability on timescales of hours to months for 2006. The model is driven by tides, wind and sea level pressure, air-sea fluxes of heat, and terrestrial buoyancy fluxes associated with river and sewage discharge. The predictive skill of the model is assessed by comparing the model simulations with independent observations of sea level from coastal tide gauges and currents from moored instruments. The simulated hydrography is also compared against a new monthly climatology created from all available temperature and salinity observations made in the Harbour over the last century. It is shown that the model can reproduce accurately the main features of the observed tides and storm surge, seasonal mean circulation and hydrography, and wind driven variations. The model is next used to examine the main physical processes controlling the circulation and hydrography of the Harbour. It is shown that non-linear interaction between tidal currents and complex topography occurs over the Narrows. The overall circulation can be characterized as a two-layer estuarine circulation with seaward flow in the thin upper layer and landward flow in the broad lower layer. An important component of this estuarine circulation is a relatively strong, vertically sheared jet situated over a narrow sill connecting the inner Harbour to the deep and relatively quiescent Bedford Basin. Local wind driven variability is strongest in winter as expected but it is also shown that a significant part of the temperature and salinity variability is driven by physical processes occurring on the adjacent inner continental shelf, especially during storm and coastal upwelling events.     

Development of a hydrodynamic model of the Bay of Biscay. Validation of hydrology

The hydrology of the Bay of Biscay was investigated using the regional ocean model MARS3D (Model for Application at Regional Scale). The simulated hydrology is compared to a set of various data encompassing monthly climatology, remote sensing SST, CTD casts, and coastal salinity measurements. Special focus was put on the validation over the continental shelf. This paper reports that despite some misfits, the climatological hydrology and its seasonal variability are correctly simulated. Various statistics computed over the period from 1999–2004 highlight different aspects of the hydrology. The biases and root mean square errors (RMSE) remain very weak at all depths when comparing salinity (<0.1 and <0.6 psu respectively). The predicted temperature shows a global overestimation of temperature (bias of around 0.8 °C) and the maximum errors are located near the thermocline (rmse of 1 °C at 20–40 m). The model is shown to properly reproduce the annual dynamics of sea surface temperature, as well as the dynamics of large river plumes observed by high frequency time series from coastal salinity gauges. The misfits highlighted by these various comparisons between model and observations are attributed to heat fluxes and mixing parameterisation.     

Investigating soil moisture sensitivity to precipitation and evapotranspiration errors using SiB2 model and ensemble Kalman filter

Accurate soil moisture information is useful in agricultural practice, weather forecasting, and various hydrological applications. Although land surface modeling provides a viable approach to simulating soil moisture, many factors such as errors in the precipitation can affect the accuracy of soil moisture simulations. This paper examined how precipitation rate and evapotranspiration rate affect the accuracy of soil moisture simulation using simple biosphere model with and without data assimilation through ensemble Kalman filter (EnKF). For each of the two variables, seven levels of relative errors (?20, ?10, ?5, 0, 5, 10 and 20 %) were introduced independently, thus a total of 49 combined cases were investigated. Observations from Wudaogou Hydrology Experimental site in the Huaihe River basin, China, were used to drive and verify the simulations. Results indicate that when the error of precipitation rate is within 10 % of the observations, the resulting error in soil moisture simulations is less significant and manageable, thus the simulated precipitation can be used to drive hydrological models in poorly gauged catchments when observations are not available. When the error of evapotranspiration rate is within 20 % of the observations, which is partly caused by model structural and parameterization errors, its impact on soil moisture simulation is less significant and can be acceptable. This study also demonstrated that the EnKF can perform consistently well to improve soil moisture simulation with less sensitivity to precipitation errors.     

Impact of ocean–wave coupling on typhoon-induced waves and surge levels around the Korean Peninsula: a case study of Typhoon Bolaven

Typhoon-induced waves and surges are important when predicting potential hazards near coastal regions. In this paper, we applied a coupled modeling system for ocean–wave interaction to examine prediction capabilities for typhoon-induced waves and surges around the Korean Peninsula. To identify how ocean–wave coupling impacts wave and surge simulations during typhoon conditions, a set of comparative experiments was performed during Typhoon Bolaven (2012): (1) a fully coupled ocean–wave model, (2) a one-way coupled ocean–wave model without surface current feedback and ocean-to-wave water levels, and (3) a stand-alone ocean model without considering wave-based sea surface roughness (SSR). When coupled with the ocean model, the surface current reduced significantly the wave height on the right-hand side of the advancing typhoon track and improved prediction accuracy along the southern coast of Korea. Compared with the observed surge levels, the simulated surge height yielded improved results for peak height magnitude and timing compared with the uncoupled model. For wave-to-surge feedback, we found that wave-induced SSR plays an important role by modulating wind stress in the surface layer. The modulated wind stress directly affected the surge height, which improved surge peak prediction during the typhoon.     

Correcting bias in rainfall inputs to a semi-distributed hydrological model using downstream flow simulation errors

The focus of this study is on bias correcting semi-distributed rainfall inputs into a hydrological model applied in the Okavango River basin in southern Africa, where there are very few local observations and heavy reliance is placed on global rainfall datasets. While the hydrological model, before rainfall bias correction, is able to represent the broad characteristics of the sub-basin streamflow responses, as demonstrated by good agreement between observed and simulated flow duration curves, there are many years where the annual volumes are over- or underestimated. The long records of observed flow at downstream stations are successfully used to bias correct the rainfall inputs to the upstream sub-basins using an analysis of their individual contributions to downstream flow and their annual rainfall–runoff response ratios. The results show improved simulations for the relatively shorter observation periods at the upstream gauging stations.     

Numerical simulation of ice-ocean variability in the Barents Sea region

A dynamic–thermodynamic sea ice model has been coupled to a three-dimensional ocean general circulation model for the purpose of conducting ocean climate dynamical downscaling experiments for the Barents Sea region. To assess model performance and suitability for such an application, the coupled model has been used to conduct a hindcast for the period 1990–2002. A comparison with available observations shows that the model successfully tracks seasonal and inter-annual variability in the ocean temperature field and that the simulated horizontal and vertical distribution of temperature are in good agreement with observations. The model results follow the seasonal and inter-annual variability in sea ice cover in the region, with the exception that the model results show too much ice melting in the northern Barents Sea during summer. The spatial distribution of the winter simulated sea ice cover is in close agreement with observations. Modelled temperatures and ice concentrations in the central Barents Sea are biased too high and too low, respectively. The probable cause is too high inflow of Atlantic Water into the Barents. The seasonal and inter-annual fluctuations in temperature and sea ice cover in the central Barents are, however, in excellent agreement with observations. Salt release during the freezing process in the numerical simulation exhibits considerable inter-annual variability and tends to vary in an opposite manner to the net inflow volume flux at the western entrance of the Barents Sea. Overall, the model produces realistic ice-ocean seasonal and inter-annual variability and should prove to be a useful tool for dynamical downscaling applications.     

Linking 1D coastal ocean modelling to environmental management: an ensemble approach

The use of a one-dimensional interdisciplinary numerical model of the coastal ocean as a tool contributing to the formulation of ecosystem-based management (EBM) is explored. The focus is on the definition of an experimental design based on ensemble simulations, integrating variability linked to scenarios (characterised by changes in the system forcing) and to the concurrent variation of selected, and poorly constrained, model parameters. The modelling system used was previously specifically designed for the use in “data-rich” areas, so that horizontal dynamics can be resolved by a diagnostic approach and external inputs can be parameterised by nudging schemes properly calibrated. Ensembles determined by changes in the simulated environmental (physical and biogeochemical) dynamics, under joint forcing and parameterisation variations, highlight the uncertainties associated to the application of specific scenarios that are relevant to EBM, providing an assessment of the reliability of the predicted changes. The work has been carried out by implementing the coupled modelling system BFM-POM1D in an area of Gulf of Trieste (northern Adriatic Sea), considered homogeneous from the point of view of hydrological properties, and forcing it by changing climatic (warming) and anthropogenic (reduction of the land-based nutrient input) pressure. Model parameters affected by considerable uncertainties (due to the lack of relevant observations) were varied jointly with the scenarios of change. The resulting large set of ensemble simulations provided a general estimation of the model uncertainties related to the joint variation of pressures and model parameters. The information of the model result variability aimed at conveying efficiently and comprehensibly the information on the uncertainties/reliability of the model results to non-technical EBM planners and stakeholders, in order to have the model-based information effectively contributing to EBM.     

Identifying storm flow pathways in a rainforest catchment using hydrological and geochemical modelling

The hydrological model TOPMODEL is used to assess the water balance and describe flow paths for the 9·73 ha Lutz Creek Catchment in Central Panama. Monte Carlo results are evaluated based on their fit to the observed hydrograph, catchment‐averaged soil moisture and stream chemistry. TOPMODEL, with a direct‐flow mechanism that is intended to route water through rapid shallow‐soil flow, matched observed chemistry and discharge better than the basic version of TOPMODEL and provided a reasonable fit to observed soil moisture and wet‐season discharge at both 15‐min and daily time‐steps. The improvement of simulations with the implementation of a direct‐flow component indicates that a storm flow path not represented in the original version of TOPMODEL plays a primary role in the response of Lutz Creek Catchment. This flow path may be consistent with the active and abundant pipeflow that is observed or delayed saturation overland flow. The ‘best‐accepted’ simulations from 1991 to 1997 indicate that around 41% of precipitation becomes direct flow and around 10% is saturation overland flow. Other field observations are needed to constrain evaporative and groundwater losses in the model and to characterize chemical end‐members posited in this paper. Published in 2004 by John Wiley & Sons, Ltd.     

Intercomparison of oceanic and atmospheric forced and coupled mesoscale simulations. Part I: Surface fluxes

A mesoscale non-hydrostatic atmospheric model has been coupled with a mesoscale oceanic model. The case study is a four-day simulation of a strong storm event observed during the SEMAPHORE experiment over a 500 × 500 km² domain. This domain encompasses a thermohaline front associated with the Azores current. In order to analyze the effect of mesoscale coupling, three simulations are compared: the first one with the atmospheric model forced by realistic sea surface temperature analyses; the second one with the ocean model forced by atmospheric fields, derived from weather forecast re-analyses; the third one with the models being coupled. For these three simulations the surface fluxes were computed with the same bulk parametrization. All three simulations succeed well in representing the main oceanic or atmospheric features observed during the storm. Comparison of surface fields with in situ observations reveals that the winds of the fine mesh atmospheric model are more realistic than those of the weather forecast re-analyses. The low-level winds simulated with the atmospheric model in the forced and coupled simulations are appreciably stronger than the re-analyzed winds. They also generate stronger fluxes. The coupled simulation has the strongest surface heat fluxes: the difference in the net heat budget with the oceanic forced simulation reaches on average 50 Wm⁻² over the simulation period. Sea surface-temperature cooling is too weak in both simulations, but is improved in the coupled run and matches better the cooling observed with drifters. The spatial distributions of sea surface-temperature cooling and surface fluxes are strongly inhomogeneous over the simulation domain. The amplitude of the flux variation is maximum in the coupled run. Moreover the weak correlation between the cooling and heat flux patterns indicates that the surface fluxes are not responsible for the whole cooling and suggests that the response of the ocean mixed layer to the atmosphere is highly non-local and enhanced in the coupled simulation.     

Improvement of surface run‐off in the hydrological model ParFlow by a scale‐consistent river parameterization

We propose an improvement of the overland‐flow parameterization in a distributed hydrological model, which uses a constant horizontal grid resolution and employs the kinematic wave approximation for both hillslope and river channel flow. The standard parameterization lacks any channel flow characteristics for rivers, which results in reduced river flow velocities for streams narrower than the horizontal grid resolution. Moreover, the surface areas, through which these wider model rivers may exchange water with the subsurface, are larger than the real river channels potentially leading to unrealistic vertical flows. We propose an approximation of the subscale channel flow by scaling Manning's roughness in the kinematic wave formulation via a relationship between river width and grid cell size, following a simplified version of the Barré de Saint‐Venant equations (Manning–Strickler equations). The too large exchange areas between model rivers and the subsurface are compensated by a grid resolution‐dependent scaling of the infiltration/exfiltration rate across river beds. We test both scaling approaches in the integrated hydrological model ParFlow. An empirical relation is used for estimating the true river width from the mean annual discharge. Our simulations show that the scaling of the roughness coefficient and the hydraulic conductivity effectively corrects overland flow velocities calculated on the coarse grid leading to a better representation of flood waves in the river channels.     

On the Coupling Between Channel Level and Surface Ground-Water Flows

This paper is devoted to a mathematical analysis of some general models of mass transport and other coupled physical processes developed in simultaneous flows of surface, soil and ground waters. Such models are widely used for forecasting (numerical simulation) of a hydrological cycle for concrete territories. The mathematical models that proved a more realistic approach are obtained by combining several mathematical models for local processes. The water-exchange models take into account the following factors: Water flows in confined and unconfined aquifers, vertical moisture migration allowing earth surface evaporation, open-channel flow simulated by one-dimensional hydraulic equations, transport of contamination, etc. These models may have different levels of sophistication. We illustrate the type of mathematical singularities which may appear by considering a simple model on the coupling of a surface flow of surface and ground waters with the flow of a line channel or river.     

Hydrological response of a medium-sized mountainous catchment to climate changes

Abstract

The long term hydrological response of a medium-sized mountainous catchment to climate changes has been examined, The climate changes were represented by a set of hypothetical scenarios of temperature increases coupled with precipitation and potential evapotranspiration changes. Snow accumulation and ablation, plus runoff from the study catchment (the Mesochora catchment in central Greece) were simulated under present (historical) and altered climate conditions using the US National Weather Service snowmelt and soil moisture accounting models. The results of this research obtained through alternative scenarios suggest strongly that all the hypothetical climate change scenarios would cause major decreases in winter snow accumulation and hence increases in winter runoff, as well as decreases in spring and summer runoff. The simulated changes in annual runoff were minor compared with the changes in the monthly distribution of runoff. Attendant changes in the monthly distribution of soil moisture and actual evapotranspiration would also occur. Such hydrological results would have significant implications on future water resources design and management.     

Water level observations from unmanned aerial vehicles for improving estimates of surface water–groundwater interaction

Integrated hydrological models are usually calibrated against observations of river discharge and piezometric head in groundwater aquifers. Calibration of such models against spatially distributed observations of river water level can potentially improve their reliability and predictive skill. However, traditional river gauging stations are normally spaced too far apart to capture spatial patterns in the water surface, whereas spaceborne observations have limited spatial and temporal resolution. Unmanned aerial vehicles can retrieve river water level measurements, providing (a) high spatial resolution; (b) spatially continuous profiles along or across the water body, and (c) flexible timing of sampling. A semisynthetic study was conducted to analyse the value of the new unmanned aerial vehicle‐borne datatype for improving hydrological models, in particular estimates of groundwater–surface water (GW–SW) interaction. Mølleåen River (Denmark) and its catchment were simulated using an integrated hydrological model (MIKE 11–MIKE SHE). Calibration against distributed surface water levels using the Differential Evolution Adaptive Metropolis algorithm demonstrated a significant improvement in estimating spatial patterns and time series of GW–SW interaction. After water level calibration, the sharpness of the estimates of GW–SW time series improves by ~50% and root mean square error decreases by ~75% compared with those of a model calibrated against discharge only.     

Sequential assimilation of soil moisture and streamflow data in a conceptual rainfall–runoff model

Soil moisture is a key hydrological variable in flood forecasting: it largely influences the partition of rain between runoff and infiltration and thus controls the flow at the outlet of a catchment. The methodology developed in this paper aims at improving the commonly used hydrological tools in an operational forecasting context by introducing soil moisture data into streamflow modelling. A sequential assimilation procedure, based on an extended Kalman filter, is developed and coupled with a lumped conceptual rainfall–runoff model. It updates the internal states of the model (soil and routing reservoirs) by assimilating daily soil moisture and streamflow data in order to better fit these external observations. We present in this paper the results obtained on the Serein, a Seine sub-catchment (France), during a period of about 2 years and using Time Domain Reflectivity probe soil moisture measurements from 0–10 to 0–100 cm and stream gauged data. Streamflow prediction is improved by assimilation of both soil moisture and streamflow individually and by coupled assimilation. Assimilation of soil moisture data is particularly effective during flood events while assimilation of streamflow data is more effective for low flows. Combined assimilation is therefore more adequate on the entire forecasting period. Finally, we discuss the adequacy of this methodology coupled with Remote Sensing data.     

Numerical study on wave dynamics and wave-induced bed erosion characteristics in Potter Cove,Antarctica

Wave generation, propagation, and transformation from deep ocean over complex bathymetric terrains to coastal waters around Potter Cove (King George Island, South Shetland Islands, Antarctica) have been simulated for an austral summer month using the Simulating Waves Nearshore (SWAN) wave model. This study aims to examine and understand the wave patterns, energy fluxes, and dissipations in Potter Cove. Bed shear stress due to waves is also calculated to provide a general insight on the bed sediment erosion characteristics in Potter Cove.A nesting approach has been implemented from an oceanic scale to a high-resolution coastal scale around Potter Cove. The results of the simulations were compared with buoy observations obtained from the National Data Buoy Center, the WAVEWATCH III model results, and GlobWave altimeter data. The quality of the modelling results has been assessed using two statistical parameters, namely the Willmott’s index of agreement D and the bias index. Under various wave conditions, the significant wave heights at the inner cove were found to be about 40–50 % smaller than the ones near the mouth of Potter Cove. The wave power in Potter Cove is generally low. The spatial distributions of the wave-induced bed shear stress and active energy dissipation were found to be following the pattern of the bathymetry, and waves were identified as a potential major driving force for bed sediment erosion in Potter Cove, especially in shallow water regions. This study also gives some results on global ocean applications of SWAN.

     

Simulating wetland impacts on stream flow in southern Africa using a monthly hydrological model

The processes that occur in wetlands and natural lakes are often overlooked and not fully incorporated in the conceptual development of many hydrological models of basin runoff. These processes can exert a considerable influence on downstream flow regimes and are critical in understanding the general patterns of runoff generation at the basin scale. This is certainly the case for many river basins of southern Africa which contain large wetlands and natural lakes and for which downstream flow regimes are altered through attenuation, storage and slow release processes that occur within the water bodies. Initial hydrological modelling studies conducted in some of these areas identified the need to explicitly account for wetland storage processes in the conceptual development of models. This study presents an attempt to incorporate wetland processes into an existing hydrological model, with the aim of reducing model structural uncertainties and improving model simulations where the impacts of wetlands or natural lakes on stream flow are evident. The approach is based on relatively flexible functions that account for the input–storage–output relationships between the river channel and the wetland. The simulation results suggest that incorporating lake and wetland storage processes into modelling can provide improved representation (the right results for the right reason) of the hydrological behaviour of some large river basins, as well as reducing some of the uncertainties in the quantification of the original model parameters used for generating the basin runoff. Copyright © 2013 John Wiley & Sons, Ltd.     

On the Use of Hydrological Models and Satellite Data to Study the Water Budget of River Basins Affected by Human Activities: Examples from the Garonne Basin of France

Natural and anthropogenic forcing factors and their changes significantly impact water resources in many river basins around the world. Information on such changes can be derived from fine scale in situ and satellite observations, used in combination with hydrological models. The latter need to account for hydrological changes caused by human activities to correctly estimate the actual water resource. In this study, we consider the catchment area of the Garonne river (in France) to investigate the capabilities of space-based observations and up-to-date hydrological modeling in estimating water resources of a river basin modified by human activities and a changing climate. Using the ISBA–MODCOU and SWAT hydrological models, we find that the water resources of the Garonne basin display a negative climate trend since 1960. The snow component of the two models is validated using the moderate-resolution imaging spectroradiometer snow cover extent climatology. Crop sowing dates based on remote sensing studies are also considered in the validation procedure. Use of this dataset improves the simulated evapotranspiration and river discharge amounts when compared to conventional data. Finally, we investigate the benefit of using the MAELIA multi-agent model that accounts for a realistic agricultural and management scenario. Among other results, we find that changes in crop systems have significant impacts on water uptake for agriculture. This work constitutes a basis for the construction of a future modeling framework of the sociological and hydrological system of the Garonne river region.