Modelling artificial channel and land‐use changes and their impact on floods and sediment yield to the Ishikari basin

The potential for flooding and sediment transport is greatly affected by river channel form and changes in land use. Therefore the modelling of channel morphology prior to canalization and of land‐use change is important with respect to the prediction of floods and sediment yield and their consequences. A combination of land‐use transformation maps and soil properties shows certain decision rules for the conversion of forest into arable or vice versa. The model proposed, from this study, was used to simulate possible past and/or future channel and land‐use patterns. Subsequently, the outcome of this simulation was used to assess the risk of flooding, sediment transport and soil‐erosion under different conditions. In this study, channel morphology prior to canalization and land‐use change in the Ishikari basin, Hokkaido, Japan, were analysed by comparing three scenarios using a physical based channel and slope model. The results indicate that pre‐canalization channel morphology has a significant impact on flood peak, but no significant effect on sediment yield. In contrast, land‐use change has a significant effect on soil eroded from hillslopes, but no significant effect on flooding for Ishikari basin. This study also illustrates the challenges that a simple model, such as a physical based channel and slope model, can simulate large‐scale river basin processes using fewer hydrological data resources. Copyright © 2004 John Wiley & Sons, Ltd.     

Preservation of meandering river channels in uniformly aggrading channel belts

Channel belt deposits from meandering river systems commonly display an internal architecture of stacked depositional features with scoured basal contacts due to channel and bedform migration across a range of scales. Recognition and correct interpretation of these bounding surfaces is essential to reconstruction of palaeochannel dimensions and to flow modelling for hydrocarbon exploration. It is therefore crucial to understand the suite of processes that form and transfer these surfaces into the fluvial sedimentary record. Here, the numerical model ‘NAYS2D’ is used to simulate a highly sinuous meandering river with synthetic stratigraphic architectures that can be compared directly to the sedimentary record. Model results highlight the importance of spatial and temporal variations in channel depth and migration rate to the generation of channel and bar deposits. Addition of net uniform bed aggradation (due to excess sediment input) allows quantification of the preservation of meander morphology for a wide range of depositional conditions. The authors find that the effect of vertical variation in scouring due to channel migration is generally orders of magnitude larger than the effect of bed aggradation, which explains the limited impact bed aggradation has on preservation of meander morphology. Moreover, lateral differences in stratigraphy within the meander belt are much larger than the stratigraphic imprint of bed aggradation. Repeatedly produced alternations of point bar growth followed by cut‐off result in a vertical trend in channel and scour feature stacking. Importantly, this vertical stacking trend differs laterally within the meander belt. In the centre of the meander belt, the high reworking intensity results in many bounding surfaces and disturbed deposits. Closer to the margins, reworking is infrequent and thick deposits with a limited number of bounding surfaces are preserved. These marginal areas therefore have the highest preservation potential for complete channel deposits and are thus best suited for palaeochannel reconstruction.     

An interpretation of the Nusselt-Rayleigh number relationship for convection in a spherical shell

Recent studies on the relationship between the Nusselt ( Nu ) and Rayleigh ( Ra ) numbers for base-heated convection in a spherical shell with a constant viscosity show that the power-law index is around 1/4, which is different from the value of 1/3 predicted by a simple boundary layer theory. We show that such a difference may be caused by the flow pattern due to the geometry. The flow pattern of the convection in a spherical shell at relatively low Ra , at least, less than 10⁶, is characterized by narrow upwelling and broad downwelling, which is similar to the opposite flow pattern of internally heated convection. Convection in the internally heated case predicts the power-law index of 1/4. We demonstrate this relationship based on the concept of 'local' Rayleigh ( Ra1 ) and Nusselt ( Nu1 ) numbers     

Numerical Investigations of Mean Winds Within Canopies of Regularly Arrayed Cubical Buildings Under Neutral Stability Conditions

Recently, several attempts have been made to model the wind velocity in an urban canopy in order to accurately predict the mixing and transport of momentum, heat, and pollutants within and above the canopy on an urban scale. For this purpose, unverified assumptions made by Macdonald (Boundary-Layer Meteorol 97:25–45, 2000) to develop a model for the profile of the mean wind velocity within an urban canopy have been used. In the present study, in order to provide foundations for improving the urban canopy models, the properties of the spatially-averaged mean quantities used to make these assumptions have been investigated by performing large-eddy simulations (LES) of the airflow around square and staggered arrays of cubical blocks with the following plan area densities: λ _p = 0.05, 0.11, 0.16, 0.20, 0.25, and 0.33. The LES results confirm that the discrepancy between the spatial average of wind velocity and Macdonald’s five-point average of wind velocity can be large in both types of arrays for large λ _p . It is also confirmed that Prandtl’s mixing length varies significantly with height within the canopy, contrary to Macdonald’s assumption for both types of arrays and for both small and large λ _p . On the other hand, in accordance with Macdonald’s assumption, the sectional drag coefficient is found to be almost constant with height except in the case of staggered arrays with high λ _p .     

Wind environment in the Lee of Kauai Island,Hawaii during trade wind conditions: weather setting for the Helios Mishap

On 26 June 2003 (approximately 1030 local time) the Helios ultralight aircraft broke apart off the west coast of Kauai Island, Hawaii as it was climbing out of the Kauai wind shadow. Following the aircraft mishap, a study was carried out to understand the conditions on the day of the crash and to better characterize the wind in the lee of Kauai. As part of this effort, both aircraft measurements and numerical modelling studies were carried out. Measurements and models showed the trade wind flow was enhanced around the island creating a region of wind shear surrounding the leeside calm zone. This wind shear region was found to be vertically oriented along the south side but tilted northward with height along the northern side of the calm zone. Several other factors on the day of the crash were investigated including water vapour gradients, diurnal Island heating, and gravity waves but their possible influences on the crash could not be confirmed. While the numerical model captured the general features of the Kauai leeside winds, the orientation of the calm zone was north of the observed one.     

Field measurement and modeling of near-bed sediment transport processes with fluid mud layer in Tokyo Bay

Tokyo Bay is one of the estuaries in Japan with a high population of almost 26 million people in the basin area. One of the major concerns for the environment in this water area is the decreasing ecosystem functions including the deterioration of water and sediment qualities caused by various anthropogenic activities. Since the bottom sediments around almost the entire area of the inner bay consist of fine materials with a high organic content, which cause the deterioration of water quality through processes such as hypoxia, an understanding of the fine sediment dynamics in the Bay is crucial for an environmental assessment of the water area. This paper proposes a model for the key processes of fine sediment dynamics, which reflects field data about muddy bed structures and their dynamics obtained during the monitoring campaign in 2007. One of the specific features of the sediment in the Bay at present is the persistent existence of fluid mud layers (water content over 300?%) with a thickness of around a few decimeters, which might be caused by deposition of abundant organic particles due to eutrophication. The present study shows that diffusion flux model delivers quite reliable results for estimating erosion flux from the top of fluid mud layers after calibrating the model parameter against the time series data of vertical flux measured by an acoustic Doppler velocimeter system. This study also derives analytical solutions, based on the Bingham fluid concept, of advection flux in the fluid mud layer on which external shear stress force is applied.     

The distribution of total and methylmercury concentrations in soils near the Idrija mercury mine, Slovenia, and the dependence of the mercury concentrations on the chemical composition and organic carbon levels of the soil

Although the mining activity of the Idrija mine in Slovenia ceased in 1995, a large amount of mining dregs containing high concentrations of mercury remains in the area. The mining dregs were transported with river flow and deposition along the Idrija River. To estimate the dispersion and change in the chemical form of mercury, a total of 28 soil core samples were taken around the river. The individual core samples were separated into layers for the analysis of their chemical composition, carbon contents, total mercury (T-Hg) and methylmercury (MeHg) concentrations. The chemical composition measured by X-ray fluorescence spectrometry was useful to estimate the dispersion of tailings: the fluvial terrace soil had a chemical composition similar to that of the tailings and could be distinguished clearly from the forest soil. The highest T-Hg concentration, 1,100 mg kg⁻¹, was observed in the fluvial terrace soil near the mine. Although the concentration decreased gradually along with distance from the mine, concentrations higher than 200 mg kg⁻¹ of T-Hg were still observed in the fluvial terrace soil approximately 20 km downstream from the mine. In the vertical distribution of T-Hg in the hillslope soil, a higher value was observed in the upper layers, which suggests the recent atmospheric deposition of mercury. The concentration of MeHg was the lowest at the riverside and higher in the hillslope soil, which was the opposite of the T-Hg distribution. The total organic carbon content tracked similarly with the distribution of MeHg and a linear relation with a significantly high correlation coefficient was obtained. The distinction may be related to the different dispersion process of mercury, and the organic carbon contents may be an important factor for MeHg formation.