Elucidating the dynamics and mixing agents of a shallow fjord through age tracer modelling

The mixing agents and their role in the dynamics of a shallow fjord are elucidated through an Eulerian implementation of artificial tracers in a three-dimensional hydrodynamic model. The time scales of vertical mixing in this shallow estuary are short, and the artificial tracers are utilized in order to reveal information not detectable in the temperature or salinity fields. The fjord's response to external forcing is investigated through a series of model experiments in which we quantify vertical mixing, transport time scales of fresh water runoff and estuarine circulation in relation to external forcing.Using age tracers released at surface and bottom, we quantify the time scales of downward mixing of surface water and upward mixing of bottom water. Wind is shown to be the major agent for vertical mixing at nearly all depth levels in the fjord, whereas the tide or external sea level forcing is a minor agent and only occasionally more important just close to the bottom. The time scale of vertical mixing of surface water to the bottom or ventilation time scale of bottom water is estimated to be in the range 0.7 h to 9.0 days, with an average age of 2.7 days for the year 2004.The fjord receives fresh water from two streams entering the innermost part of the fjord, and the distribution and age of this water are studied using both ageing and conservative tracers. The salinity variations outside this fjord are large, and in contrast to the salinity, the artificial tracers provide a straight forward analysis of river water content. The ageing tracer is used to estimate transport time scales of river water (i.e. the time elapsed since the water left the river mouth). In May 2004, the typical age of river water leaving the fjord mouth is 5 days. As the major vertical mixing agent is wind, it controls the estuarine circulation and export of river water. When the wind stress is set to zero, the vertical mixing is reduced and the vertical salinity stratification is increased, and the river water can be effectively exported out of the fjord.We also analyse the river tracer fields and salinity field in relation to along estuary winds in order to detect signs of wind-induced straining of the along estuary density gradient. We find that events of down estuary winds are primarily associated with a reduced along estuary salinity gradient due to increased surface salinity in the innermost part of the fjord, and with an overall decrease in vertical stratification and river water content at the surface. Thus, our results show no apparent signs of wind-induced straining in this shallow fjord but instead they indicate increased levels of vertical mixing or upwelling during down estuary wind events.     

On the importance of soil moisture in calibration of rainfall–runoff models: two case studies

Streamflow modelling results from the GR4H and PDM hydrological models were evaluated in two Australian sub-catchments, using (1) calibration to streamflow and (2) joint-calibration to streamflow and soil moisture. Soil moisture storage in the models was evaluated against soil moisture observations from field measurements. The PDM had the best performance in terms of both streamflow and soil moisture estimations during the calibration period, but was outperformed by GR4H during validation. It was also shown that the soil moisture estimation was improved significantly by joint-calibration for the case where streamflow and soil moisture estimations were poor. In other cases, addition of the soil moisture constraint did not degrade the results. Consequently, it is recommended that GR4H be used, in preference to the PDM, in the foothills of the Murrumbidgee catchment or other Australian catchments with semi-arid to sub-humid climate, and that soil moisture data be used in the calibration process.     

The sheltering effect of the Balearic Islands in the hindcast wave field

The sheltering effect of the Balearic Islands in the hindcast wave field was studied for typical Mediterranean wave situations of Llevant, Tramuntana and Mestral and for mild conditions such as the Garbí and Ponent winds. For this purpose, a third generation wave model was applied to the Mediterranean Sea and different patterns of the sheltered areas were found for the various representative situations depending on the wind variability and on the magnitude of the wind speed. From the analysis it was concluded that the sheltered zones created during storms generally persist for short periods of time of the order of 6 h, possibly reaching a maximum of 12 h. In contrast with earlier results obtained for swell dominated ocean areas, it was observed that in this area, due to the short fetches the sea states are mainly local wind seas and thus the wave field behind the islands depends on the local wind.     

Implications of subsurface thermal–hydraulic–mechanical processes associated with glaciation on Shield flow system evolution and performance assessment

A deep geologic repository (DGR) situated on the Canadian Shield will be subject to long-term climate change that will markedly alter surface conditions as a result of glaciation and permafrost penetration. Systematic, two-dimensional and three-dimensional coupled thermal–hydraulic–mechanical finite-element simulations with varying degrees of coupling, including depth-dependent salinity (represented as a change in groundwater density) and temperature-dependent density and viscosity, were undertaken to address the implications of glaciation on groundwater flow system dynamics as it could affect DGR performance. The modelling domain consisted of a 1.6-km deep sub-regional scale (≈100 km²) fractured Shield flow system. Initial and transient thermal, hydraulic and mechanical boundary conditions were developed from two realizations of the University of Toronto Glacial Systems Model of the last Laurentide glaciation. Results indicate that during the glacial loading/unloading cycle, for this particular conceptual model, there is limited penetration of glacial meltwaters to depth and small residual anomalous hydraulic head. During glacial coverage, the mechanical factor of safety increases in the moderately fractured and sparsely fractured rock mass, but principal effective stress reorientation also occurs. Given the assumed nonglacial in situ state of stress and mechanical properties, the fracture zones were predicted to be less stable under glacial conditions.

Natural flood management,land use and climate change trade-offs: the case of Tarland catchment,Scotland

A distributed hydrological model (WaSiM-ETH) was applied to a mesoscale catchment to investigate natural flood management as a nonstructural approach to tackle flood risks from climate change. Peak flows were modelled using climate projections (UKCP09) combined with afforestation-based land-use change options. A significant increase in peak flows was modelled from climate change. Afforestation could reduce some of the increased flow, with greatest benefit from coniferous afforestation, especially replacing lowland farmland. Nevertheless, large-scale woodland expansion was required to maintain peak flows similar to present and beneficial effects were significantly reduced for larger “winter-type” extreme floods. Afforestation was also modelled to increase low-flow risks. Land-use scenarios showed catchment-scale trade-offs across multiple objectives meant “optimal” flood risk solutions were unlikely, especially for afforestation replacing lowland farmland. Hence, combined structural/nonstructural measures may be required in such situations, with integrated catchment management to synergize multiple objectives.     

Numerical Study Of The Impact Of Coherent Structures On Vertical Transfers In The Atmospheric Boundary Layer

In two preceding papers, coherent structures of theatmospheric boundary layer (ABL), such as rollvortices or cells, were investigated through radar andaircraft observations collected during the TRAC-93(Turbulence Radar Aircraft Cells) experiment held inFrance in June 1993. The analysis of this experimentaldata set provided information on the spatialcharacteristics of these organisations (length scale,orientation, type ... ), their temporal and verticalevolution, and their relation with the dynamic andthermodynamic conditions of the ABL. For the thirdpaper in this series, a large eddy simulation model is used to examine the impact of thecoherent structures on the ABL vertical fluxes. Theanalysis of the simulated horizontal fields is madewith two-dimensional auto and cross-correlationsapplied on different pertinent ABL variables. Theresults emphasise a directional anisotropy of theseorganised fields throughout the ABL, much morepronounced in the heat flux fields, not only at thelength scale of organisations but also at theturbulence scales. This finding has an importantconsequence for traditional ABL flux measurementsbased on the hypothesis of isotropic and homogeneousturbulence. It can explain part of the underestimationof the surface fluxes often mentioned in theliterature. This approach makes it possible tomodify the concept of diffusion time (in chemicalmodelling) and could also lead to revised ABLparameterisations in Range Scale models.