· · · Institutional Repository

An analysis of field measurements recorded over a tidal cycle in the Rotterdam Waterway is presented. These measurements are the first to elucidate the processes influencing the along-channel current structure and the excursion of the salt wedge in this estuary. The salt wedge structure remained stable throughout the measuring period. The velocity measurements indicate decoupling effects between the layers and that bed-generated turbulence is confined below the pycnocline. The barotropic M4 overtide structure is imposed at the mouth of the estuary, and the generation of M4 overtides within the estuary is found to be relatively small. Internal tidal asymmetry does not make a significant contribution to the M4 velocity frequency band. Instead, the combination of barotropic and baroclinic forcing, in conjunction with the suppression of turbulence at the interface, provides the main explanation for the time dependence and mean structure of the flow in the Rotterdam Waterway. This gives rise to the observed differences in the length of the flood and ebb, in the magnitudes of the flood and ebb velocities, in the length of the slack water periods, and in the timing of the onset of slack water at the surface and near the bed. It results in the formation of distinct exchange flow profiles at the head of the salt wedge around slack water and the creation of maximal velocities at the pycnocline during flood. Advection governs the displacement and structure of the salt wedge since turbulent mixing is suppressed. The tidal displacement of the salt wedge controls the height of the pycnocline above the bed at a particular site. Hence, it controls the height to which bed-generated turbulence can protrude into the water column. Consequently, the authors find asymmetries in the structure of the internal flow, turbulent mixing, and bed stresses that are not related to classical internal tidal asymmetry.

The generation of residual circulation in a tidally energetic estuary with constant longitudinal salinity gradient and parabolic cross section is examined by means of a two-dimensional cross-sectional numerical model, neglecting river runoff and Stokes drift. It is shown how the longitudinal and lateral residual circulation can be decomposed into contributions from various processes such as tidal straining circulation, gravitational circulation, advectively driven circulation, and horizontal mixing circulation. The sensitivity of the residual circulation and its components from various processes to changes in forcing is investigated by varying the Simpson number (nondimensional longitudinal buoyancy gradient) and the unsteadiness parameter (nondimensional tidal frequency), as well as the bed roughness and the width of the estuary. For relatively weak salinity gradient forcing, the tidal straining circulation dominates the residual exchange circulation in support of classical estuarine circulation (up-estuary flow near the bed and down-estuary flow near the surface). The strength of the longitudinal estuarine circulation clearly increases with increased salinity gradient forcing. However, when the Simpson number exceeds 0.15, the relative contributions of both gravitational circulation and advectively driven circulation to estuarine circulation increase substantially. Lateral residual circulation is relatively weak for small Simpson numbers and becomes flood oriented (divergent flow near the bed and convergent flow near the surface) for larger Simpson numbers because of increasing contributions from gravitational and advectively driven circulation. Increasing the unsteadiness number leads to decreased longitudinal and lateral residual circulation. Although changes in bed roughness result in relatively small changes in residual circulation, results are sensitive to the width of the estuary, mainly because of changes in residual exchange circulation driven by tidal straining.

Tidal straining, which can mathematically be described as the covariance between eddy viscosity and vertical shear of the along-channel velocity component, has been acknowledged as one of the major drivers for estuarine circulation in channelized tidally energetic estuaries. In this paper, the authors investigate the role of lateral circulation for generating this covariance. Five numerical experiments are carried out, starting with a reference scenario including the full physics and four scenarios in which specific key physical processes are neglected. These processes are longitudinal internal pressure gradient forcing, lateral internal pressure gradient forcing, lateral advection, and the neglect of temporal variation of eddy viscosity. The results for the viscosity–shear covariance are correlated across different experiments to quantify the change due to neglect of these key processes. It is found that the lateral advection of vertical shear of the along-channel velocity component and its interaction with the tidally asymmetric eddy viscosity (which is also modified by the lateral circulation) is the major driving force for estuarine circulation in well-mixed tidal estuaries.

In this study the response of dry convective boundary layers to nonstationary surface heat fluxes is systematically investigated. This is relevant not only during sunset and sunrise but also, for example, when clouds modulate incoming solar radiation. Because the time scale of the associated change in surface heat fluxes may differ from case to case, the authors consider the generic situation of oscillatory surface heat fluxes with different frequencies and amplitudes and study the response of the boundary layer in terms of transfer functions. To this end both a mixed layer model (MLM) and a large-eddy simulation (LES) model are used; the latter is used to evaluate the predictive quality of the mixed layer model. The mixed layer model performs generally quite well for slow changes in the surface heat flux and provides analytical understanding of the transfer characteristics of the boundary layer such as amplitude and phase lag. For rapidly changing surface fluxes (i.e., changes within a time frame comparable to the large eddy turnover time), it proves important to account for the time it takes for the information to travel from the surface to higher levels of the boundary layer such as the inversion zone. As a follow-up to a 1997 study by Sorbjan, who showed that the conventional convective velocity scale is inadequate as a scaling quantity during the decay phase, this paper addresses the issue of defining, in (generic) transitional situations, a velocity scale that is solely based on the surface heat flux and its history.

Modeling of the systematic errors in the second-order moments of wind speeds measured by continuouswave (ZephIR) and pulsed (WindCube) lidars is presented. These lidars use the conical scanning technique to measure the velocity field. The model captures the effect of volume illumination and conical scanning. The predictions are compared with the measurements from the ZephIR, WindCube, and sonic anemometers at a flat terrain test site under different atmospheric stability conditions. The sonic measurements are used at several heights on a meteorological mast in combination with lidars that are placed on the ground. Results show that the systematic errors are up to 90% for the vertical velocity variance, whereas they are up to 70% for the horizontal velocity variance. For the ZephIR, the systematic errors increase with height, whereas for the WindCube, they decrease with height. The systematic errors also vary with atmospheric stability and are low for unstable conditions. In general, for both lidars, the model agrees well with the measurements at all heights and under different atmospheric stability conditions. For the ZephIR, the model results are improved when an additional low-pass filter for the 3-s scan is also modeled. It is concluded that with the current measurement configuration, these lidars cannot be used to measure turbulence precisely.

Field measurements are presented, which are the first to quantify the processes influencing the entrapment of suspended particulate matter (SPM) at the limit of saltwater intrusion in the Rotterdam Waterway. The estuarine turbidity maximum (ETM) is shown to be maintained by the trapping of fluvial SPM at the head of the salt wedge. The trapping process is associated with the raining out of fluvial SPM from the upper, fresher part of the water column, into the layer below the pycnocline. The dominant mechanisms responsible are baroclinic shear flows and the abrupt change in turbulent mixing characteristics due to damping of turbulence at the pycnocline. This view contrasts with the assumption of landward transport of marine SPM by asymmetries in bed stress. The SPM transport capacity of the tidal flow is not fully utilized in the ETM, and the ETM is independent of a bed-based supply of mud. This is explained by regular exchange of part of the ETM with harbor basins, which act as efficient sinks, and that the Rotterdam Waterway is not a complete fluvial SPM trap. The supply of SPM by the freshwater discharge ensures that the ETM is maintained over time. Hence, theETMis an advective phenomenon. Relative motion between SPM and saltwater occurs because of lags introduced by resuspension. Moreover,SPM that lags behind the salt wedge after high water slack (HWS) is eventually recollected at the head. Hence, SPM follows complex transport pathways and the mechanisms involved in trapping and transport of SPM are inherently three-dimensional.

Amethod has been developed that largely automates the labor-intensive extraction work for large amounts of rainfall strip charts and paper rolls. The method consists of the following five basic steps: 1) scanning the charts and rolls to high-resolution digital images, 2) manually and visually registering relevant meta information from charts and rolls and preprocessing rolls to locate day transitions, 3) applying automatic curve extraction software in a batch process to determine the coordinates of cumulative rainfall lines on the images, 4) postprocessing the curves that were not correctly determined in step 3, and 5) aggregating the cumulative rainfall in pixel coordinates to the desired time resolution. The core of the method is in step 3. Here a color detection procedure is introduced that automatically separates the background of the charts and rolls from the grid and subsequently the rainfall curve. The rainfall curve is detected by minimization of a cost function. In total, 321 station years of locations in the Netherlands have successfully been digitized and transformed to long-term rainfall time series with 5-min resolution. In about 30% of the cases, semiautomatic postprocessing of the results was needed using a purpose-built graphical interface application. This percentage, however, strongly depends on the quality of the recorded curves and the charts and rolls. Although developed for rainfall, the method can be applied to other elements as well.

A new method for retrieving air velocity fluctuations in the cloud-capped boundary layer (BL) using radar reflectivity and the Doppler velocity fields is proposed. The method was developed on the basis of data obtained by the Transportable Atmospheric Radar (TARA) located in Cabauw, Netherlands, at 0500–0812 UTC 8 May 2004, and tested using a detailed trajectory ensemble model of the cloud-capped BL. During the observations, the BL depth was 1200 m, and the cloud base (measured by a lidar) was at 500–550 m. No preliminary assumptions concerning the shapes of drop size distributions were made. On the basis of the TARAradar data, vertical profiles of the vertical air velocity standard deviation, of turbulent dissipation rate, etc. were estimated. The correlation functions indicate the existence of large eddies in the BL with a characteristic horizontal scale of about 600 m. Analysis of the slope (the scaling parameter) of the structure functions indicates that turbulence above 400 m can be considered to be isotropic. Below this level, the turbulence becomes anisotropic. The rate of anisotropy increases with the decrease of the height above the surface. The averaged values of the dissipation rate were evaluated as 1–2 cm2 s23. The importance of using the cloud-capped BL model as a link between different types of observed data (radar, lidar, aircraft, etc.) is discussed. More data should be analyzed to understand the changes in the turbulent structure of the BL during its growth, as well as during cloud and drizzle formation.

The mechanism behind the collapse of turbulence in the evening as a precursor to the onset of the very stable boundary layer is investigated. To this end a cooled, pressure-driven flow is investigated by means of a local similarity model. Simulations reveal a temporary collapse of turbulence whenever the surface heat extraction, expressed in its nondimensional form h/L, exceeds a critical value. As any temporary reduction of turbulent friction is followed by flow acceleration, the long-term state is unconditionally turbulent. In contrast, the temporary cessation of turbulence, which may actually last for several hours in the nocturnal boundary layer, can be understood from the fact that the time scale for boundary layer diffusion is much smaller than the time scale for flow acceleration. This limits the available momentum that can be used for downward heat transport. In case the surface heat extraction exceeds the so-called maximum sustainable heat flux (MSHF), the near-surface inversion rapidly increases. Finally, turbulent activity is largely suppressed by the intense density stratification that supports the emergence of a different, calmer boundary layer regime.

The collapse of turbulence in the nocturnal boundary layer is studied by means of a simple bulk model that describes the basic physical interactions in the surface energy balance. It is shown that for a given mechanical forcing, the amount of turbulent heat that can be transported downward is limited to a certain maximum. In the case of weak winds and clear skies, this maximum can be significantly smaller than the net radiative loss minus soil heat transport. In the case when the surface has low heat capacity, this imbalance generates rapid surface cooling that further suppresses the turbulent heat transport, so that eventually turbulence largely ceases (positive feedback mechanism). The model predicts the minimum wind speed for sustainable turbulence for the so-called crossing level. At this level, some decameters above the surface, the wind is relatively stationary compared to lower and higher levels. The critical speed is predicted in the range of about 5–7 m s21, depending on radiative forcing and surface properties, and is in agreement with observations at Cabauw. The critical value appears not very sensitive to model details or to the exact values of the input parameters. Finally, results are interpreted in terms of external forcings, such as geostrophic wind. As it is generally larger than the speed at crossing height, a 5 m s21 geostrophic wind may be considered as the typical limit below which sustainable, continuous turbulence under clear-sky conditions is unlikely to exist. Below this threshold emergence of the very stable nocturnal boundary layer is anticipated.