Spatial scale effect on seasonal streamflows in permafrost catchments on the Qinghai–Tibet Plateau

The scale issue is of central concern in hydrological processes to understand the potential upscaling or downscaling methodologies, and to develop models for scaling the dominant processes at different scales and in different environments. In this study, a typical permafrost watershed in the Qinghai‐Tibet Plateau was selected. Its hydrological processes were monitored for 4 years from 2004 to 2008, measuring the effects of freezing and thawing depth of active soil layers on runoff processes. To identify the nature and cause of variation in the runoff response in different size catchments, catchments ranging from 1·07 to 112 km² were identified in the watershed. The results indicated that the variation of runoff coefficients showed a ‘V’ shape with increasing catchment size during the spring and autumn seasons, when the active soil was subjected to thawing or freezing processes. A two‐stage method was proposed to create runoff scaling models to indicate the effects of scale on runoff processes. In summer, the scaling transition model followed an exponential function for mean daily discharge, whereas the scaling model for flood flow exhibited a linear function. In autumn, the runoff process transition across multiple scales followed an exponential function with air temperature as the driving factor. These scaling models demonstrate relatively high simulation efficiency and precision, and provide a practical way for upscaling or downscaling runoff processes in a medium‐size permafrost watershed. For permafrost catchments of this scale, the results show that the synergistic effect of scale and vegetation cover is an important driving factor in the runoff response. Copyright © 2011 John Wiley & Sons, Ltd.     

Internal tides and waves near the continental shelf edge

Abstract

The subject is reviewed from the viewpoints of theory, internal tide and wave structure and their implications.

A wider theoretical context suggests scope for further investigation of natural or nearly-trapped forms above the inertial frequency.

Although internal tides in many locations are observed to have first-mode vertical structure, higher modes are seen offshore from shallow shelf-break forcing and for particular Froude numbers, and may be expected locally near generation. Bottom intensification is often observed where the sea floor matches the characteristic slope. Solitons form from internal tides of large amplitude or at large changes of depth.

Internal tides and solitons are observed also at many sills and in straits, and to intensify in canyons.

Non-linear effects of the waves, especially solitons, include the conveyance of water, nutrients, ‘‘mixing potential'’ etc. away from their source to other locations, and the generation of mean currents. The waves transfer energy and possibly heat between the ocean and shelf, may be a source of medium frequency waves on the shelf (periods of minutes) and can contribute to interior mixing and overturning, bottom stirring and sediment movement.     

Numerical simulation of fully saturated porous materials

Fully coupled, porous solid–fluid formulation, implementation and related modeling and simulation issues are presented in this work. To this end, coupled dynamic field equations with u?p?U formulation are used to simulate pore fluid and soil skeleton (elastic–plastic porous solid) responses. Present formulation allows, among other features, for water accelerations to be taken into account. This proves to be useful in modeling dynamic interaction of media of different stiffnesses (as in soil–foundation–structure interaction). Fluid compressibility is also explicitly taken into account, thus allowing excursions into modeling of limited cases of non‐saturated porous media. In addition to these features, present formulation and implementation models in a realistic way the physical damping, which dissipates energy. In particular, the velocity proportional damping is appropriately modeled and simulated by taking into account the interaction of pore fluid and solid skeleton. Similarly, the displacement proportional damping is physically modeled through elastic–plastic processes in soil skeleton. An advanced material model for sand is used in present work and is discussed at some length. Also explored in this paper are the verification and validation issues related to fully coupled modeling and simulations of porous media. Illustrative examples describing the dynamical behavior of porous media (saturated soils) are presented. The verified and validated methods and material models are used to predict the behavior of level and sloping grounds subjected to seismic shaking. Copyright © 2008 John Wiley & Sons, Ltd.     

Glacier fluctuations in the southern Peruvian Andes during the late‐glacial period,constrained with cosmogenic 3He

The occurrence of pronounced climate reversals during the last glacial termination has long been recognised in palaeoclimate records from both hemispheres and from high to low latitudes. Accurate constraint of both the timing and magnitude of events, such as the Younger Dryas and Antarctic Cold Reversal, is vital in order to test different hypotheses for the causes and propagation of abrupt climate change. However, in contrast to higher‐latitude regions, well‐dated records from the Tropics are rare and the structure of late‐glacial tropical climate remains uncertain. As a step toward addressing this problem, we present an in situ cosmogenic ³He surface exposure chronology from Nevado Coropuna, southern Peru, documenting a significant fluctuation of the ice margin during the late‐glacial period. Ten tightly clustered ages from a pair of moraines located halfway between the modern glacier and the Last Glacial Maximum terminus range from 11.9 to 13.9 ka and give an arithmetic mean age of 12.8 ± 0.7 ka (1σ). These data constitute direct evidence for a readvance, or prolonged stillstand, of glaciers in the arid Andes of southwestern Peru. Copyright © 2011 John Wiley & Sons, Ltd.     

Prediction of overland flow and seepage zones associated with the interaction of multiple infiltration devices (cascading infiltration devices)

Infiltration devices are traditionally evaluated as standalone entities that do not interact with each other. A model is outlined that will allow interactions between proposed infiltration devices to be predicted prior to a development commencing. The model allows prediction of seepage into downslope devices and the assessment of the locations where the combined ground‐water mound will reach the surface and result in overland flow. The volume of overland flow discharged by the seepage zone may exceed the overland and piped flow received by the infiltration devices. Copyright © 2007 John Wiley & Sons, Ltd.     

Internal testing of a numerical model of hillslope–channel coupling using laboratory flume experiments

The numerical model COUP 2D simulates the hydrological coupling between hillslopes and the river channel during a rainfall event. In order to test the numerical model, a 1:100 scaled laboratory flume which was modified to incorporate lateral hillslope elements, was used to run a series of experiments in which hillslope angle, channel angle, hillslope discharge and channel discharge were the varying parameters. Overall, there were 18 different experimental configurations with three replicates carried out for each condition, leading to a total of 54 experiments. These conditions were then used to parameterize and run COUP 2D. Internal model outputs of flow depth and flow velocity at four cross‐sections in the channel were compared to the measurements made in the physical model for the same parameter conditions. Statistical comparisons of the measured and modelled data were carried out for each experiment and across all experiments, using two goodness‐of‐fit measures—root mean square error and Nash–Sutcliffe coefficient of efficiency—in order to assess the performance of the model over an entire simulation as well as over all the simulations. The main effects on the goodness‐of‐fit measures for flow depth of each experimental variable, as well as the interactions between variables, were evaluated using statistical modelling. The results show that the model captures flow‐depth variations in response to changing channel and hillslope parameters. Statistical modelling suggests that the main effects on model error are cross‐section position, channel angle and channel discharge. Significant interactions also occur between all the channel variables and between the channel variables and hillslope discharge. The results of the testing procedure have significant implications for the consideration of different model components and for the interaction between data‐ and model evaluation. Copyright © 2007 John Wiley & Sons, Ltd.     

Climate variability,snow, and physiographic controls on storm hydrographs in small forested basins,western Cascades,Oregon

Large floods are often attributed to the melting of snow during a rain event. This study tested how climate variability, snowpack presence, and basin physiography were related to storm hydrograph shape in three small (<1 km²) basins with old‐growth forest in western Oregon. Relationships between hydrograph characteristics and precipitation were tested for approximately 800 storms over a nearly 30‐year period. Analyses controlled for (1) snowpack presence/absence, (2) antecedent soil moisture, and (3) hillslope length and gradient. For small storms (<150 mm precipitation), controlling for precipitation, the presence of a snowpack on near‐saturated soil increased the threshold of precipitation before hydrograph rise, extended the start lag, centroid lag, and duration of storm hydrographs, and increased the peak discharge. The presence of a snowpack on near‐saturated soil sped up and steepened storm hydrographs in a basin with short steep slopes, but delayed storm hydrographs in basins with longer or more gentle slopes. Hydrographs of the largest events, which were extreme regional rain and rain‐on‐snow floods, were not sensitive to landform characteristics or snowpack presence/absence. Although the presence of a snowpack did not increase peak discharge in small, forested basins during large storms, it had contrasting effects on storm timing in small basins, potentially synchronizing small basin contributions to the larger basin hydrograph during large rain‐on‐snow events. By altering the relative timing of hydrographs, snowpack melting could produce extreme floods from precipitation events whose size is not extreme. Further work is needed to examine effects of canopy openings, snowpack, and climate warming on extreme rain‐on‐snow floods at the large basin scale. Copyright © 2008 John Wiley & Sons, Ltd.     

The impact of land use change and check‐dams on catchment sediment yield

Extensive land use changes have occurred in many areas of SE Spain as a result of reforestation and the abandonment of agricultural activities. Parallel to this the Spanish Administration spends large funds on hydrological control works to reduce erosion and sediment transport. However, it remains untested how these large land use changes affect the erosion processes at the catchment scale and if the hydrological control works efficiently reduce sediment export. A combination of field work, mapping and modelling was used to test the influence of land use scenarios with and without sediment control structures (check‐dams) on sediment yield at the catchment scale. The study catchment is located in SE Spain and suffered important land use changes, increasing the forest cover 3‐fold and decreasing the agricultural land 2·5‐fold from 1956 to 1997. In addition 58 check‐dams were constructed in the catchment in the 1970s accompanying reforestation works. The erosion model WATEM‐SEDEM was applied using six land use scenarios: land use in 1956, 1981 and 1997, each with and without check‐dams. Calibration of the model provided a model efficiency of 0·84 for absolute sediment yield. Model application showed that in a scenario without check dams, the land use changes between 1956 and 1997 caused a progressive decrease in sediment yield of 54%. In a scenario without land use changes but with check‐dams, about 77% of the sediment yield was retained behind the dams. Check‐dams can be efficient sediment control measures, but with a short‐lived effect. They have important side‐effects, such as inducing channel erosion downstream. While also having side‐effects, land use changes can have important long‐term effects on sediment yield. The application of either land use changes (i.e. reforestation) or check‐dams to control sediment yield depends on the objective of the management and the specific environmental conditions of each area. Copyright © 2008 John Wiley & Sons, Ltd.     

A new soil freeze/thaw discriminant algorithm using AMSR‐E passive microwave imagery

The soil freeze–thaw controls the hydrological and carbon cycling and thus affects water and energy exchanges at land surface. This article reported a newly developed algorithm for distinguishing the freeze/thaw status of surface soil. The algorithm was based on information from Advanced Microwave Scanning Radiometer Enhanced (AMSR‐E) which records brightness temperature (Tb) in the afternoon and after midnight. The criteria and discriminant functions were obtained from both radiometer observations and model simulations. First of all, the microwave radiation from freeze–thaw soil was examined by carrying out experimental measurements at 18·7 and 36·5 GHz using a Truck‐mounted Multi‐frequency Microwave Radiometer (TMMR) in the Heihe River of China. The experimental results showed that the soil moisture is a key component that differentiates the microwave radiation behaviours during the freeze–thaw process, and the differences in soil temperature and emissivity between frozen and thawed soils were found to be the most important criteria. Secondly, a combined model was developed to consider the impacts of complex ground surface conditions on the discrimination. The model simulations quite followed the trend of in situ observations with an overall relation coefficient (R) of approximately 0·88. Finally, the ratio of Tb18·7H (horizontally polarized Tb at 18·7 GHz) to Tb36·5V was considered primarily as the quasi‐emissivity, which is more reasonable and explicit in measuring the microwave radiation changes in soil freezing and thawing than the spectral gradient. By combining Tb36·5V to indicate the soil temperature variety, a Fisher linear discrimination analysis was used to establish the discriminant functions. After being corrected by TMMR measurements, the new discriminant algorithm had an overall accuracy of 86% when validated by 4‐cm soil temperature. The multi‐year discriminant results also provided a good agreement with the classification map of frozen ground in China. Copyright © 2011 John Wiley & Sons, Ltd.     

Development of a river delta: a case study of Cimanuk river mouth,Indonesia

This paper describes delta development processes with particular reference to Cimanuk Delta in Indonesia. Cimanuk river delta, the most rapidly growing river delta in Indonesia, is located on the northern coast of Java Island. The delta is subject to ocean waves of less than 1 m height due to its position in the semi‐enclosed Java Sea in the Indonesian archipelago. The study has been carried out using a hydrodynamic model that accounts for sediment movement through the rivers and estuaries. As an advanced approach to management of river deltas, a numerical model, namely MIKE‐21, is used as a tool in the management of Cimanuk river delta. From calibration and verification of hydrodynamic model, it was found that the best value of bed roughness was 0·1 m. For the sediment‐transport model, the calibration parameters were adjusted to obtain the most satisfactory results of suspended sediment concentration and volume of deposition. By comparing the computed and observed data in the calibration, the best values of critical bed shear stress for deposition, critical bed shear stress for erosion and erosion coefficient were 0·05 N m^?2, 0·15 N m^?2, and 0·00001 kg m^?2 s^?1, respectively. The calibrated model was then used to analyse sensitivity of model parameters and to simulate delta development during the periods 1945–1963 and 1981–1997. It was found that the sensitive model parameters were bed shear stresses for deposition and erosion, while the important model inputs were river suspended sediment concentration, sediment characteristics and hydrodynamic. The model result showed reasonable agreement with the observed data. As evidenced by field data, the mathematical model proves that the Cimanuk river delta is a river‐dominated delta because of its protrusion pattern and very high sediment loads from the Cimanuk river. It was concluded that 86% of sediment load from the Cimanuk river was deposited in the Cimanuk delta. Copyright © 2008 John Wiley & Sons, Ltd.