Unraveling driving factors for large rock–ice avalanche mobility

Large rock–ice avalanches have attracted attention from scientists for decades and some of these events have caused high numbers of fatalities. A relation between rock slope instabilities in cold high mountain areas and climate change is currently becoming more evident and questions about possible consequences and hazard scenarios in densely populated high mountain regions leading beyond historical precedence are rising. To improve hazard assessment of potential rock–ice avalanches, their mobility is a critical factor. This contribution is an attempt to unravel driving factors for the mobility of large rock–ice avalanches by synthesizing results from physical laboratory experiments and empirical data from 64 rock–ice avalanches with volumes >1x10⁶ m³ from glacierized high mountain regions around the world. The influence of avalanche volume, water and ice content, low‐friction surfaces, and topography on the apparent coefficient of friction (as a measure of mobility) is assessed. In laboratory experiments granular ice in the moving mass was found to reduce bulk friction up to 20% while water led to a reduction around 50% for completely saturated material compared with dry flows. Evidence for the effects of water as a key driving factor to enhance mobility was also found in the empirical data, while the influence of the ice content could not be confirmed to be of much relevance in nature. Besides liquefaction, it was confirmed that mobility increases with volumes and that frictional surface characteristics such as flow paths over glaciers are also dominant variables determining mass movement mobility. Effects of the topography along the flow path as well as channeling are assumed to be other critical factors. The results provide an empirical basis to roughly account for different path and flow characteristics of large rock–ice avalanches and to find appropriate ranges for friction parameters for scenario modeling and hazard assessments. Copyright © 2011 John Wiley & Sons, Ltd.     

A large landslide event in a post‐glacial landscape: rethinking glacial legacy

Threlkeld Knotts (c. 500 m above sea level) in the English Lake District has hitherto been considered to be a glacially‐modified intrusion of microgranite. However, its surface features are incompatible with glacial modification; neither can these nor the subsurface structures revealed by ground‐penetrating radar (GPR) be explained by post‐glacial subaerial processes acting on a glacially‐modified microgranite intrusion. Here we re‐interpret Threlkeld Knotts as a very large post‐glacial landslide involving the microgranite, with an estimated volume of about 4 × 10⁷ m³. This interpretation is tested against published and recent information on the geology of the site, the glacial geomorphic history of the area and newly‐acquired GPR data. More than 60 large post‐Last Glacial Maximum (LGM) rock–slope failures have significantly modified the glaciated landscape of the Lake District; this is one of the largest. Recognition of this major landslide deposit in such a well‐studied environment highlights the need to continuously re‐examine landscapes in the light of increasing knowledge of geomorphic processes and with available technology in currently active or de‐glaciating environments. Copyright © 2012 John Wiley & Sons, Ltd.     

Coupling patterns between para‐glacial and permafrost degradation responses in Antarctica

The recently deglaciated environments in maritime permafrost regions are usually affected by very active paraglacial processes. Elephant Point is an ice‐free area of 1.16 km2 located in the SW of Livingston island (South Shetland Islands, Antarctica). Between 1956–2010 the retreat of the ice cap covering most part of this island has exposed 17.3% of the land surface in this peninsula. Two geomorphological units were identified in this new ice‐free area: a moraine extending from the western to the eastern coastlines and a relatively flat proglacial surface. The glacier in 1956 sat in contact with the northern slope of the moraine, but its accelerated retreat ‐ in parallel to the warming trend recorded in the Antarctic Peninsula ‐ left these areas free of glacier ice. Subsequently, the postglacial evolution was controlled by the relaxation phase typical of paraglacial systems. The typology and intensity of geomorphological processes show a significantly different dynamics between the southern and northern slopes of the moraine. This pattern is related to the different stage of paraglacial adjustment in both slopes. In the southern side, on coarser sediments, pronival ramparts, debris flows and alluvial fans are distributed, with a low to moderate activity of slope processes. In the northern side, mass wasting processes are extremely active on fine‐grained unconsolidated sediments. Ice‐rich permafrost is being degraded by thermokarst processes. Landslides and mudflows transfer large amounts of sediments down‐slope. The surface affected by retrogressive‐thaw slumps in the moraine has been quantified in 24,172 m2, which accounts for 9.6% of its surface. The abundance of kettle‐lakes is also indicative of the degradation of the ground ice. Paraglacial processes are expected to continue in the moraine and proglacial area in the near future, although their intensity and duration will depend on the magnitude and rate of future climate trends in the northern Antarctic Peninsula. Copyright © 2015 John Wiley & Sons, Ltd.     

Hazardous processes and events from glacier and permafrost areas: lessons from the Chilean and Argentinean Andes

Glacier and permafrost hazards such as glacial‐lake outburst floods and rock–ice avalanches cause significant socio‐economic damages worldwide, and these processes may increase in frequency and magnitude if the atmospheric temperature rises. In the extratropical Andes nearly 200 human deaths were linked to these processes during the twentieth century. We analysed bibliographical sources and satellite images to document the glacier and permafrost dynamics that have caused socio‐economic damages in this region in historic time (including glacial lake outburst floods, ice and rock–ice avalanches and lahars) to unravel their causes and geomorphological impacts. In the extratropical Andes, at least 15 ice‐dammed lakes and 16 moraine‐dammed lakes have failed since the eighteenth century, causing dozens of floods. Some floods rank amongst the largest events ever recorded (5000 × 10⁶ m³ and 229 × 10⁶ m³, respectively). Outburst flood frequency has increased in the last three decades, partially as a consequence of long‐term (decades to centuries) climatic changes, glaciers shrinkage, and lake growth. Short‐term (days to weeks) meteorological conditions (i.e. intense and/or prolonged rainfall and high temperature that increased meltwater production) have also triggered outburst floods and mass movements. Enormous mass failures of glaciers and permafrost (> 10 × 10⁶ m³) have impacted lakes, glaciers, and snow‐covered valleys, initiating chain reactions that have ultimately resulted in lake tsunamis and far‐reaching (> 50 km) flows. The eruption of ice‐covered volcanoes has also caused dozens of damaging lahars with volumes up to 45 × 10⁶ m³. Despite the importance of these events, basic information about their occurrence (e.g. date, causes, and geomorphological impact), which is well established in other mountain ranges, is absent in the extratropical Andes. A better knowledge of the processes involved can help to forecast and mitigate these events. Copyright © 2014 John Wiley & Sons, Ltd.     

Enhanced rock‐slope failure following ice‐sheet deglaciation: timing and causes

The temporal pattern of rock‐slope failures (RSFs) following Late Pleistocene deglaciation on tectonically stable terrains is controversial: previous studies variously suggest (1) a rapid response due to removal of supporting ice (‘debuttressing’), (2) a progressive decline in RSF frequency, and (3) a millennial‐scale delay before peak RSF activity. We test these competing models through beryllium‐10 (¹⁰Be) exposure dating of five closely‐spaced quartzite RSFs on the Isle of Jura, Scotland, to establish the relationship between timing of failure and those of deglaciation, episodes of rapid warming and periods of rapid glacio‐isostatic uplift. All five dated RSFs occurred at least 720–2240 years after deglaciation, with the probability of failure peaking ~2 ka after deglaciation, consistent with millennial‐scale delay model (3). This excludes debuttressing as an immediate cause of failure, though it is likely that time‐dependent stress release due to deglacial unloading resulted in progressive development of failure planes within the rock. Thaw of permafrost ice in joints is unlikely to have been a prime trigger of failure as some RSFs occurred several centuries after the onset of interstadial warming. Conversely, the timespan of the RSFs coincides with the period of maximum glacio‐isostatic crustal uplift, suggesting that failure was triggered by uplift‐driven seismic events acting on fractured rock masses. Implications of this and related research are: (1) that retreat of the last Pleistocene ice sheets across tectonically‐stable mountainous terrains was succeeded by a period of enhanced rock‐slope failure due to deglacial unloading and probably uplift‐driven seismicity; (2) that the great majority of RSFs in the British Isles outside the limits of Loch Lomond Stadial (= Younger Dryas) glaciation are of Lateglacial (pre‐Holocene) age; and (3) numerous RSFs must also have occurred inside Loch Lomond Stadial (LLS) glacial limits, but that runout debris was removed by LLS glaciers. Copyright © 2013 John Wiley & Sons, Ltd.     

Effect of permafrost thaw on the dynamics of lakes recharged by ice‐jam floods: case study of Yukon Flats,Alaska

Large river floods are a key water source for many lakes in fluvial periglacial settings. Where permeable sediments occur, the distribution of permafrost may play an important role in the routing of floodwaters across a floodplain. This relationship is explored for lakes in the discontinuous permafrost of Yukon Flats, interior Alaska, using an analysis that integrates satellite‐derived gradients in water surface elevation, knowledge of hydrogeology, and hydrologic modelling. We observed gradients in water surface elevation between neighbouring lakes ranging from 0.001 to 0.004. These high gradients, despite a ubiquitous layer of continuous shallow gravel across the flats, are consistent with limited groundwater flow across lake basins resulting from the presence of permafrost. Permafrost impedes the propagation of floodwaters in the shallow subsurface and constrains transmission to ‘fill‐and‐spill’ over topographic depressions (surface sills), as we observed for the Twelvemile‐Buddy Lake pair following a May 2013 ice‐jam flood on the Yukon River. Model results indicate that permafrost table deepening of 1–11 m in gravel, depending on watershed geometry and subsurface properties, could shift important routing of floodwater to lakes from overland flow (fill‐and‐spill) to shallow groundwater flow (‘fill‐and‐seep’). Such a shift is possible in the next several hundred years of ground surface warming and may bring about more synchronous water level changes between neighbouring lakes following large flood events. This relationship offers a potentially useful tool, well suited to remote sensing, for identifying long‐term changes in shallow groundwater flow resulting from thawing of permafrost. Copyright © 2015 John Wiley & Sons, Ltd.     

Physical modelling of rainfall‐ and snowmelt‐induced erosion of stony slope underlain by permafrost

Physical modelling experiments have been carried out in a cold room to test on a small scale, the effects of water supply during the thaw of an experimental slope with permafrost. Permafrost was maintained at depth and a thin active layer was frozen and thawed from the surface. Data from the experiments relate to two different conditions, first with moderate rainfall, and second with heavy rainfall during the thaw period. When moderate rainfall is applied during thaw phases, the experimental slope is slightly degraded. At the scale of the experiment, erosion processes involve frost jacking of the coarse blocks, frost creep and gelifluction that induce slow and gradual down slope displacements of the active layer, but also small landslides leading to large but slow mass movements with short displacements. Changes in experimental slope morphology are marked by the initiation of a small‐scale drainage network and the development of a little crest line which shows a progressive upslope migration. With such boundary conditions, there is not enough water supply to evacuate downslope the whole of the eroded material and a topographic smoothing is observed. When heavy rainfall is applied during thaw periods, rapid mass wasting (small mud‐flows and debris flows) become prominent. Slope failures are largely controlled by the water saturation of the active layer and by the occurrence of steeper slopes. At the scale of the experiment, rates of erosion and maximum incision increase by about 100% leading to significant slope degradation with marked and specific scars comparable to gullying. These morphological changes are dependant on both the size and the frequency of catastrophic events. These experiments provide detailed data that could improve the knowledge of the physical parameters that control the initiation, at a small‐scale, of erosion processes on periglacial slopes with a thin active layer and/or with thin cover of mobilizable slope deposits. Copyright © 2010 John Wiley & Sons, Ltd.     

A geomorphology‐based integrated stream–aquifer interaction model for semi‐gauged catchments

Many of the existing stream–aquifer interaction models available in the literature are very complex with limited applicability in semi‐gauged and ungauged catchments. In this study, to estimate the influent and effluent subsurface water fluxes under limited geo‐hydrometeorological data availability conditions, a simple stream–aquifer interaction model, namely, the variable parameter McCarthy–Muskingum (VPMM) hillslope‐storage Boussinesq (hsB) model, has been developed. This novel model couples the VPMM streamflow transport with the hsB groundwater flow transport modules in online mode. In this integrated model, the surface water–groundwater flux exchange process is modelled by the Darcian approach with the variable hydraulic heads between the river stage and groundwater table accounting for the rainfall forcing. Considering the exchange fluxes in the hyporheic zone and lateral overland flow contribution, this approach is field tested in a typical 48‐km stretch of the Brahmani River in eastern India to simulate the streamflow and its depth with the minimum Nash–Sutcliffe efficiency of 94% and 88%; the maximum root mean square error of 134 m³/s and 0.35 m; and the minimum index of agreement of 98% and 97%, respectively. This modelling approach could be very well utilized in data‐scarce world‐river basins to estimate the stream–aquifer exchange flux due to rainfall forcings.     

Hydrochemical variations during flood pulses in the south‐west China peak cluster karst: impacts of CaCO3–H2O–CO2 interactions

High‐resolution measurements of rainfall, water level, pH, conductivity, temperature and carbonate chemistry parameters of groundwater at two adjacent locations within the peak cluster karst of the Guilin Karst Experimental Site in Guangxi Province, China, were made with different types of multiparameter sonde. The data were stored using data loggers recording with 2 min or 15 min resolution. Waters from a large, perennial spring represent the exit for the aquifer's conduit flow, and a nearby well measures water in the conduit‐adjacent, fractured media. During flood pulses, the pH of the conduit flow water rises as the conductivity falls. In contrast, and at the same time, the pH of groundwater in the fractures drops, as conductivity rises. As Ca²⁺ and HCO₃^? were the dominant (>90%) ions, we developed linear relationships (both r² > 0·91) between conductivity and those ions, respectively, and in turn calculated variations in the calcite saturation index (SI_C) and CO₂ partial pressure (P) of water during flood pulses. Results indicate that the P of fracture water during flood periods is higher than that at lower flows, and its SI_C is lower. Simultaneously, P of conduit water during the flood period is lower than that at lower flows, and its SI_C also is lower. From these results we conclude that at least two key processes are controlling hydrochemical variations during flood periods: (i) dilution by precipitation and (ii) water–rock–gas interactions. To explain hydrochemical variations in the fracture water, the water–rock–gas interactions may be more important. For example, during flood periods, soil gas with high CO₂ concentrations dissolves in water and enters the fracture system, the water, which in turn has become more highly undersaturated, dissolves more limestone, and the conductivity increases. Dilution of rainfall is more important in controlling hydrochemical variations of conduit water, because rainfall with higher pH (in this area apparently owing to interaction with limestone dust in the lower atmosphere) and low conductivity travels through the conduit system rapidly. These results illustrate that to understand the hydrochemical variations in karst systems, considering only water–rock interactions is not sufficient, and the variable effects of CO₂ on the system should be evaluated. Consideration of water–rock–gas interactions is thus a must in understanding variations in karst hydrochemistry. Copyright © 2004 John Wiley & Sons, Ltd.     

Permian back‐arc extension in central Inner Mongolia,NE China: Elemental and Sr–Nd–Pb–Hf–O isotopic constraints from the Linxi high‐MgO diabase dikes

Early Permian (272 ± 2 Ma) diabase dikes from the Linxi area in central Inner Mongolia of NE China have high MgO (10.4 – 12.3 wt%), Cr (301 – 448 ppm) and Ni (167 – 233 ppm) concentrations, and show enrichments in large ion lithophile element (LILE) and light rare earth elements (REE) but depletions in high field strength element (HFSE, e.g., Nb and Ta), with depleted mantle‐type Sr [⁸⁷Sr/⁸⁶Sr (i) = 0.70315 – 0.70362], Nd [ε_Nd (t) = +6.8 – +7.4], Pb [²⁰⁶Pb/²⁰⁴Pb (i) = 18.10 – 18.16] and zircon Hf [ε_Hf (t) = +14.7 – +19.1] isotopic compositions, but slightly higher zircon δ¹⁸O (5.2 – 6.0 ‰ with an average of 5.7 ‰) than normal mantle. The combined geochemical data indicate their derivation from a depleted mantle metasomatized by recycled crustal component. Elemental and isotopic modeling results suggest that the primary magma was produced through 5 % to 10 % melting of a depleted mantle, which contained approximately 1 % sediment fluid released from the subducted paleo‐Asian Ocean. Considering the widespread distribution of contemporaneous mafic rocks across the central Inner Mongolia, which show REE patterns from E‐MORBs to normal MORBs, we propose a petrogenetic link between the Early Permian mafic magmatism and a back‐arc extension in response to northward subduction of the paleo‐Asian Ocean. The Permian mafic magmatism and the new age constraints from the metamorphic and sedimentary records in this area tend to indicate the ultimate closure of the paleo‐Asian Ocean by the end of Paleozoic.     

Climate‐invariant area–slope relations in channel heads initiated by surface runoff

We use a dataset of 38 field‐mapped channel heads from a semi‐arid environment in western Colorado to examine relationships between contributing drainage area (A) and local hillslope gradient (θ) in relation to dominant process of initiation (surface runoff versus subsurface flow). Channel heads resulting primarily from subsurface flow have significantly greater values of A, longer basin lengths, and shallower local gradients than channel heads resulting primarily from surface runoff. We also compare the data from western Colorado to six analogous datasets from more humid regions in other portions of the United States and in southeast Australia. Comparison of the confidence intervals for the exponent values of A–θ regression lines reveals that the confidence intervals for the exponent of western Colorado channel heads with both surface and subsurface flow overlap with the confidence intervals for the exponent of all other datasets. This suggests that A–θ relationships do not differ significantly between diverse geographic locations. Copyright © 2017 John Wiley & Sons, Ltd.     

Zircon fission‐track and U–Pb double dating using femtosecond laser ablation–inductively coupled plasma–mass spectrometry: A technical note

We present a new LA–ICP–MS system for zircon fission‐track (FT) and U–Pb double dating, whereby a femtosecond laser combined with galvanometric optics simultaneously ablates multiple spots to measure average surface U contents. The U contents of zircon measured by LA–ICP–MS and standardized with the NIST SRM610 glass are comparable to those measured by the induced FT method, and have smaller analytical errors. LA–ICP–MS FT dating of seven zircon samples including three IUGS age standards is as accurate as the external detector method, but can give a higher‐precision age depending on the counting statistics of the U content measurement. Double dating of the IUGS age standards gives FT and U–Pb ages that are in agreement. A chip of the Nancy 91500 zircon has a homogeneous U content of 84 ppm, suggesting the possibility of using this zircon as a matrix‐matched U‐standard for FT dating. When using the Nancy 91500 zircon as a U‐standard, a zeta calibration value of 42–43 year cm² for LA–ICP–MS FT dating is obtained. While this value is strictly valid only for the particular session, it can serve as a reference for other studies.     

Groundwater–surface‐water interactions in a braided river: a tracer‐based assessment

Natural tracers (alkalinity and silica) were used to infer groundwater–surface‐water exchanges in the main braided reach of the River Feshie, Cairngorms, Scotland. Stream‐water samples were collected upstream and downstream of the braided section at fortnightly intervals throughout the 2001–2002 hydrological year and subsequently at finer resolution over two rainfall events. The braided reach was found to exert a significant downstream buffering effect on the alkalinity of these waters, particularly at moderate flows (4–8 m³ s^?1/?Q_30–70). Extensive hydrochemical surveys were undertaken to characterize the different source waters feeding the braids. Shallow groundwater flow systems at the edge of the braided floodplain, recharged by effluent streams and hillslope drainage, appeared to be of particular significance. Deeper groundwater was identified closer to the main channel, upwelling through the hyporheic zone. Both sources contributed to the significant groundwater–surface‐water interactions that promote the buffering effect observed through the braided reach. Their impact was less significant at higher flows (>15 m³ s^?1/>Q₁₀) when acidic storm runoff from the peat‐covered catchment headwaters dominated, as well as under baseflow conditions (<4 m³ s^?1/<Q₇₀), when upstream alkalinity was already buffered owing to headwater groundwater sources assuming dominance. The significant temporally and spatially dynamic influence of these groundwater–surface‐water interactions was therefore seen to have important implications for both catchment functioning and instream ecology. Copyright © 2004 John Wiley & Sons, Ltd.     

A combined rainfall infiltration model based on Green–Ampt and SCS‐curve number

The Green–Ampt infiltration equation is an incomplete governing equation for rainfall infiltration due to the absence of an inertia term. The estimation of the capillary pressure head at the wetting front is difficult to determine. Thus, a major limitation of the Green–Ampt model is the constant, non‐zero surface ponding depth. This paper proposes an integrated rainfall infiltration model based on the Green–Ampt model and the SCS‐CN model. It achieves a complete governing equation for rainfall infiltration by momentum balance and the water budget based on the Green–Ampt assumption, and uses the curve number from the SCS‐CN method to calculate the initial abstraction, which is used as a basic parameter for the governing equation of the intensity of rainfall loss during the runoff period. The integrated rainfall infiltration model resolves the dilemma for capillary pressure head estimation, overcomes the limitation of constant, non‐zero surface ponding depth, and facilitates the calculation of runoff for individual flood simulations. Copyright © 2014 John Wiley & Sons, Ltd.     

Implementation of a second‐order paraxial boundary condition for a water‐saturated layered half‐space in plane strain

A half‐space finite element and a transmitting boundary are developed for a water‐saturated layered half‐space using a paraxial boundary condition. The exact dynamic stiffness of a half‐space in plane strain is derived and a second‐order paraxial approximation of the stiffness is obtained. A half‐space finite element and a transmitting boundary are then formulated. The development is verified by comparison of the dynamic stiffness of impermeable and permeable rigid strip foundations with other published results. The advantage of using the paraxial boundary condition in comparison with the rigid boundary condition is examined. It is shown that the paraxial boundary condition offers significant gain and the resulting half‐space finite element and transmitting boundary can represent the effects of a water‐saturated layered half‐space with good accuracy and efficiency. In addition, the numerical method described herein maintains the strengths and advantages of the finite element method and can be easily applied to demanding problems of soil–structure interaction in a water‐saturated layered half‐space. Copyright © 2010 John Wiley & Sons, Ltd.     

Spatial and temporal variability of 0‐ to 5‐m soil–water storage at the watershed scale

Dynamic relationships among rainfall patterns, soil water distribution, and plant growth are crucial for sustainable conservation of soil and water resources in water‐limited ecosystems. Spatial and temporal variation in deep soil water content at a watershed scale have not yet been characterized adequately due to the lack of deep soil water data. Deep soil–water storage (SWS) up to a depth of 5 m (n = 73) was measured at 19 sampling occasions at the LaoYeManQu watershed on the Chinese Loess Plateau (CLP). At a depth of 0–1.5 m, the annual mean SWS was highly correlated with rain intensity, and the correlation decreased with depth, but within the layers at 1.5–5.0 m, the changes in SWS indicated a lag between precipitation and the replenishment of soil water. Geostatistical parameters of SWS were also highly dependent on depth, and the mean SWS presented similar spatial structures in two adjacent layers. Temporal stability of SWS as indicated by mean relative difference, standard deviation of the relative difference (SDRD), and mean absolute bias error (MABE) was significantly weaker at the shallow than at deeper layers. Soil separates and organic carbon content controlled the spatial pattern of SWS at the watershed scale. One representative location (Site 57) was identified to estimate the mean SWS in the 1‐ to 5‐m layer of the watershed. Semivariograms of the SDRD and MABE were best fitted by an isotropic spherical model, and their spatial distributions were depth‐dependent. Both temporal stability and spatial variability of SWS increased over depth. This study is helpful for deep SWS estimation and sustainable management of soil and water on the CLP, and for other similar regions around the world.     

Block‐and‐ash flow deposit of the Narcondam Volcano: Product of dacite–andesite dome collapse in the Burma–Java subduction complex

Narcondam Island in the Andaman Sea represents a dacite–andesite dome volcano in the volcanic chain of the Burma–Java subduction complex. The pyroclasts of andesitic composition are restricted to the periphery of the dome predominantly in the form of block‐and‐ash deposits and minor base surge deposits. Besides pyroclastic deposits, andesitic lava occurs dominantly at the basal part of the dome whereas dacitic lava occupies the central part of the dome. The pyroclasts are represented by non‐vesiculated to poorly vesiculated blocks of andesite, lapilli, and ash. The hot debris derived from dome collapse was deposited initially as massive to reversely‐graded beds with the grain support at the lower part and matrix support at the upper part. This sequence is overlain by repetitive beds of lapilli breccia to tuff breccia. These deposits are recognized as a basal avalanche rather than lahar deposit. This basal avalanche was punctuated by an ash‐cloud surge deposit representing a sequence of thinly bedded units of normal graded unit to parallel laminated beds.     

Mean and turbulent flow fields in a simulated ice‐covered channel with a gravel bed: some laboratory observations

Northern rivers experience freeze‐up over the winter, creating asymmetric under‐ice flows. Field and laboratory measurements of under‐ice flows typically exhibit flow asymmetry and its characteristics depend on the presence of roughness elements on the ice cover underside. In this study, flume experiments of flows under a simulated ice cover are presented. Open water conditions and simulated rough ice‐covered flows are discussed. Mean flow and turbulent flow statistics were obtained from an Acoustic Doppler Velocimeter (ADV) above a gravel‐bed surface. A central region of faster flow develops in the middle portion of the flow with the addition of a rough cover. The turbulent flow characteristics are unambiguously different when simulated ice covered conditions are used. Two distinct boundary layers (near the bed and in the vicinity of the ice cover, near the water surface) are clearly identified, each being characterized by high turbulent intensity levels. Detailed profile measurements of Reynolds stresses and turbulent kinetic energy indicate that the turbulence structure is strongly influenced by the presence of an ice cover and its roughness characteristics. In general, for y/d > 0·4 (where y is height above bed and d is local flow depth), the addition of cover and its roughening tends to generate higher turbulent kinetic energy values in comparison to open water flows and Reynolds stresses become increasingly negative due to increased turbulence levels in the vicinity of the rough ice cover. The high negative Reynolds stresses not only indicate high turbulence levels created by the rough ice cover but also coherent flow structures where quadrants one and three dominate. Copyright © 2012 John Wiley & Sons, Ltd.     

Rosenbrock‐based algorithms and subcycling strategies for real‐time nonlinear substructure testing

In this paper, Rosenbrock‐based algorithms originally developed for real‐time testing of linear systems with dynamic substructuring are extended for use on nonlinear systems. With this objective in mind and for minimal overhead, both two‐ and three‐stages linearly implicit real‐time compatible algorithms were endowed with the Jacobian matrices requiring only one evaluation at the beginning of each time step. Moreover, these algorithms were improved with subcycling strategies. In detail, the paper briefly introduces Rosenbrock‐based L‐Stable Real‐Time (LSRT) algorithms together with linearly implicit and explicit structural integrators, which are now commonly used to perform real‐time tests. Then, the LSRT algorithms are analysed in terms of linearized stability with reference to an emulated spring pendulum, which was chosen as a nonlinear test problem, because it is able to exhibit a large and relatively slow nonlinear circular motion coupled to an axial motion that can be set to be stiff. The accuracy analysis on this system was performed for all the algorithms described. Following this, a coupled spring‐pendulum example typical of real‐time testing is analysed with respect to both stability and accuracy issues. Finally, the results of representative numerical simulations and real‐time substructure tests, considering nonlinearities both in the numerical and the physical substructure, are explored. These tests were used to demonstrate how the LSRT algorithms can be used for substructuring tests with strongly nonlinear components. Copyright © 2010 John Wiley & Sons, Ltd.