Wavelet analysis of inter‐annual variability in the runoff regimes of glacial and nival stream catchments,Bow Lake,Alberta

Continuous wavelet analyses of hourly time series of air temperature, stream discharge, and precipitation are used to compare the seasonal and inter‐annual variability in hydrological regimes of the two principal streams feeding Bow Lake, Banff National Park, Alberta: the glacial stream draining the Wapta Icefields, and the snowmelt‐fed Bow River. The goal is to understand how water sources and flow routing differ between the two catchments. Wavelet spectra and cross‐wavelet spectra were determined for air temperature and discharge from the two streams for summers (June–September) 1997–2000, and for rainfall and discharge for the summers of 1999 and 2000. The diurnal signal of the glacial runoff was orders of magnitude higher in 1998 than in other years, indicating that significant ice exposure and the development of channelized glacial drainage occurred as a result of the 1997–98 El Niño conditions. Early retreat of the snowpack in 1997 and 1998 led to a significant summer‐long input of melt runoff from a small area of ice cover in the Bow River catchment; but such inputs were not apparent in 1999 and 2000, when snow cover was more extensive. Rainfall had a stronger influence on runoff and followed quicker flow paths in the Bow River catchment than in the glacial catchment. Snowpack thickness and catchment size were the primary controls on the phase relationship between temperature and discharge at diurnal time scales. Wavelet analysis is a fast and effective means to characterize runoff, temperature, and precipitation regimes and their interrelationships and inter‐annual variability. The technique is effective at identifying inter‐annual and seasonal changes in the relative contributions of different water sources to runoff, and changes in the time required for routing of diurnal meltwater pulses through a catchment. However, it is less effective at identifying changes/differences in the type of the flow routing (e.g. overland flow versus through flow) between or within catchments. Copyright © 2003 John Wiley & Sons, Ltd.     

Hysteresis‐based analysis of overland metal transport

Introducing a concept of equivalent mass depth of flow, this study describes the phenomenon of non‐point source pollutant (metal) transport for pavement (or overland) flow in analogy with wave propagation in wide open channels. Hysteretic and normal mass rating curves are developed for runoff rate and mass of 12 dissolved and particulate‐bound metal elements (pollutants) using the rainfall‐runoff and water quality data of the 15 × 20 m² instrumented pavement in Cincinnati, USA. Normal mass rating curves developed for easy computation of pollutant load are found to be of a form similar to Manning's, which is valid for open channel flows. Based on the hysteresis analysis, wave types for dissolution and mixing of particulate‐bound metals are identified. The analysis finds that the second‐order partial‐differential equation normally used for metal transport does not have the efficacy to describe fully the strong non‐linear phenomena such as is described for various metal elements by dynamic waves. In addition, the proportionality concept of the popular SCS‐CN concept is extended for determining the potential maximum metal mass M_p of all the 12 elements transported by a rain storm and related to the antecedent dry period (ADP). For the primary metal zinc element, M_p is found to increase with the ADP and vice versa. Copyright © 2003 John Wiley & Sons, Ltd.     

Modelling the effect of low soil temperatures on transpiration by Scots pine

For ecosystem modelling of the Boreal forest it is important to include processes associated with low soil temperature during spring‐early summer, as these affect the tree water uptake. The COUP model, a physically based SVAT model, was tested with 2 years of soil and snow physical measurements and sap flow measurements in a 70‐year‐old Scots pine stand in the boreal zone of northern Sweden. During the first year the extent and duration of soil frost was manipulated in the field. The model was successful in reproducing the timing of the soil warming after the snowmelt and frost thaw. A delayed soil warming, into the growing season, severely reduced the transpiration. We demonstrated the potential for considerable overestimation of transpiration by the model if the reduction of the trees' capacity to transpire due to low soil temperatures is not taken into account. We also demonstrated that the accumulated effect of aboveground conditions could be included when simulating the relationship between soil temperature and tree water uptake. This improved the estimated transpiration for the control plot and when soil warming was delayed into the growing season. The study illustrates the need of including antecedent conditions on root growth in the model in order to catch these effects on transpiration. The COUP model is a promising tool for predicting transpiration in high‐latitude stands. Copyright © 2006 John Wiley & Sons, Ltd.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Post-solidification cooling and the age of Io's lava flows

The modeling of thermal emission from active lava flows must account for the cooling of the lava after solidification. Models of lava cooling applied to data collected by the Galileo spacecraft have, until now, not taken this into consideration. This is a flaw as lava flows on Io are thought to be relatively thin with a range in thickness from ∼1 to 13 m. Once a flow is completely solidified (a rapid process on a geological time scale), the surface cools faster than the surface of a partially molten flow. Cooling via the base of the lava flow is also important and accelerates the solidification of the flow compared to the rate for the ‘semi-infinite’ approximation (which is only valid for very deep lava bodies). We introduce a new model which incorporates the solidification and basal cooling features. This model gives a superior reproduction of the cooling of the 1997 Pillan lava flows on Io observed by the Galileo spacecraft. We also use the new model to determine what observations are necessary to constrain lava emplacement style at Loki Patera. Flows exhibit different cooling profiles from that expected from a lava lake. We model cooling with a finite-element code and make quantitative predictions for the behavior of lava flows and other lava bodies that can be tested against observations both on Io and Earth. For example, a 10-m-thick ultramafic flow, like those emplaced at Pillan Patera in 1997, solidifies in ∼450 days (at which point the surface temperature has cooled to ∼280 K) and takes another 390 days to cool to 249 K. Observations over a sufficient period of time reveal divergent cooling trends for different lava bodies [examples: lava flows and lava lakes have different cooling trends after the flow has solidified (flows cool faster)]. Thin flows solidify and cool faster than flows of greater thickness. The model can therefore be used as a diagnostic tool for constraining possible emplacement mechanisms and compositions of bodies of lava in remote-sensing data.     

Motive force of Qinghai-Tibet Plateau moving to east

A calculation formula on spherical pattern of Qinghai-Tibet plateau moving model is established. Tibet massif moves norward by east in speed of 28 mm/a, Ganshu-Qinghai massif moves to northeast in speed of 15 mm/a, Qomolangma Feng moves northward by a few east in speed of 35 ～42 mm/a. The low latitude perimeter is longer than the nigh latitude perimeter. When the Tibet massif moves northward, its latitude perimeter must be contracted and the Tibet massif must move eastward by Coriolis. Coriolis force is inertial in earth rotation. It makes the fall body turning to east and the rising block turning westward. In the Northern Hemisphere, it makes the northward body turning to east and the southward block turning to west.This is the reason why the tectonic zones of western Pacific are different from those of eastern Pacific.     

The influence of macrophyte growth,typical of eutrophic conditions,on river flow velocities and turbulence production

The influence of emergent and submerged macrophytes on flow velocity and turbulence production is demonstrated in a 140 m reach of the River Blackwater in Farnborough, Hampshire, UK. Macrophyte growth occurs in patches and is dominated by Sparganium erectum and Sparganium emersum. In May 2001, patches of S. erectum were already established and occupied 18% of the channel area. The flow adjusted to these (predominantly lateral) patches by being channelled through a narrower cross‐section. The measured velocity profiles showed a logarithmic form, with deviations attributable to topographic control. The channel bed was the main source of turbulence. In September 2001, in‐stream macrophytes occupied 27% of the channel, and overhanging bank vegetation affected 32% of the area. Overall flow resistance, described by Manning's n, showed a threefold increase that could be attributed to the growth of S. emersum in the middle of the channel. Velocity profiles showed different characteristic forms depending on their position relative to plant stems and leaves. The overall velocity field had a three‐dimensional structure. Turbulence intensities were generally higher and turbulence profiles tended to mirror the velocity profiles. Evidence for the generation of coherent eddies was provided by ratios of the root mean square velocities. Spectral analysis identified deviations from the Kolmogorov ?5/3 power law and provided statistical evidence for a spectral short‐cut, indicative of additional turbulence production. This was most marked for the submerged vegetation and, in some instances, the overhanging bank vegetation. The long strap‐like leaves of S. emersum being aligned approximately parallel to the flow and the highly variable velocity field created by the patch arrangement of macrophytes suggest that the dominant mechanism for turbulence production is vortex shedding along shear zones. Wake production around individual stems of S. emersum close to the bed may also be important locally. Copyright © 2006 John Wiley & Sons, Ltd.     

The Lower Permian Wasp Head Formation,Sydney Basin: high‐latitude,shallow marine sedimentation following the late Asselian to early Sakmarian glacial event in eastern Australia

The Lower Permian Wasp Head Formation (early to middle Sakmarian) is a ～95 m thick unit that was deposited during the transition to a non‐glacial period following the late Asselian to early Sakmarian glacial event in eastern Australia. This shallow marine, sandstone‐dominated unit can be subdivided into six facies associations. (i) The marine sediment gravity flow facies association consists of breccias and conglomerates deposited in upper shoreface water depths. (ii) Upper shoreface deposits consist of cross‐stratified, conglomeratic sandstones with an impoverished expression of the Skolithos Ichnofacies. (iii) Middle shoreface deposits consist of hummocky cross‐stratified sandstones with a trace fossil assemblage that represents the Skolithos Ichnofacies. (iv) Lower shoreface deposits are similar to middle shoreface deposits, but contain more pervasive bioturbation and a distal expression of the Skolithos Ichnofacies to a proximal expression of the Cruziana Ichnofacies. (v) Delta‐influenced, lower shoreface‐offshore transition deposits are distinguished by sparsely bioturbated carbonaceous mudstone drapes within a variety of shoreface and offshore deposits. Trace fossil assemblages represent distal expressions of the Skolithos Ichnofacies to stressed, proximal expressions of the Cruziana Ichnofacies. Impoverished trace fossil assemblages record variable and episodic environmental stresses possibly caused by fluctuations in sedimentation rates, substrate consistencies, salinity, oxygen levels, turbidity and other physio‐chemical stresses characteristic of deltaic conditions. (vi) The offshore transition‐offshore facies association consists of mudstone and admixed sandstone and mudstone with pervasive bioturbation and an archetypal to distal expression of the Cruziana Ichnofacies. The lowermost ～50 m of the formation consists of a single deepening upward cycle formed as the basin transitioned from glacioisostatic rebound following the Asselian to early Sakmarian glacial to a regime dominated by regional extensional subsidence without significant glacial influence. The upper ～45 m of the formation can be subdivided into three shallowing upward cycles (parasequences) that formed in the aftermath of rapid, possibly glacioeustatic, rises in relative sea‐level or due to autocyclic progradation patterns. The shift to a parasequence‐dominated architecture and progressive decrease in ice‐rafted debris upwards through the succession records the release from glacioisostatic rebound and amelioration of climate that accompanied the transition to broadly non‐glacial conditions.     

Carbonate aeolianites of western Saurashtra,India: experimental decipherment of the depositional mechanisms

Granular carbonate deposits of Late Pleistocene to Early Holocene age, commonly referred to as ‘miliolite limestone’, occur in a linear belt, parallel to the southern coast of Saurashtra, India. In the present study area these carbonate deposits are found in select valleys between ridges and mounds of pyroclastic material present in the Deccan trap plateau. Two different depositional histories have been proposed for these sediments. The presence of marine bioclasts led to the postulation of a marine origin for these deposits. The second school of thought propounded redeposition of the coastal sediments by aeolian processes. Although a few features could not be explained by the proposed aeolian model, critical comparison of these two views favoured the aeolian origin. The mode of occurrence, lithological and structural attributes, and microscopic evidence presented here, also support a possible aeolian origin for these deposits. Experimental observation indicates that these carbonate aeolianites represent backflow deposits, which accumulated because of the flow separation caused by the presence of topographic highs. The conspicuous concave‐up geometry of the deposit conformed to the shape of the separation bulb. In view of the inferred depositional mechanism, the disposition of the deposits and the signature of the palaeoflow direction suggest that the carbonate particles were derived from the north‐western coast of Saurashtra by strong south‐easterly winds. Massive granular carbonates with outsized basement clasts appear to be the product of avalanching of granular material from the higher contours because of oversteepening of the primary deposit.     

Late Quaternary contourites and glaciomarine sedimentation in the Fram Strait

Two sites in the eastern Fram Strait, the Vestnesa Ridge and the Yermak Plateau, have been surveyed and sampled providing a depositional record over the last glacial‐interglacial cycle. The Fram Strait is the only deep‐water connection from the Arctic Ocean to the North Atlantic and contains a marine sediment record of both high latitude thermohaline flow and ice sheet interaction. On the Vestnesa Ridge, the western Svalbard margin, a sediment drift was identified in 1226 m of water. Gravity and multicores from the crest of the drift recovered turbidites and contourites. ¹⁴C dating indicates an age range of 8287 to 26 900 years BP (Early Holocene to Late Weichselian). The Yermak Plateau is characterized by slope sediments in 961 m of water. Gravity and multicores recovered contourites and hemipelagites. ¹⁴C ages were between 8615 and 46 437 years BP (Early Holocene to mid‐Weichselian). Downcore dinoflagellate cyst analyses from both sites provide a record of changing surface water conditions since the mid‐Weichselian, suggesting variable sea ice extent, productivity and polynyas present even during the Last Glacial Maximum. Four layers of ice‐rafted debris were also identified and correlated within the cores. These events occurred ca at 9, 24 to 25, 26 to 27 and 43 ka, asynchronous with Heinrich layers in the wider north‐east Atlantic and here interpreted as reflecting instability in the Svalbard/Barents Ice sheet and the northward advection of warm Atlantic water during the Late Weichselian. The activity of the ancestral West Spitsbergen Current is interpreted using mean sortable silt records from the cores. On the Vestnesa Ridge drift the modern mass accumulation rate, calculated using excess ²¹⁰Pb, is 0·076 g cm^?2 year^?1. On the Yermak Plateau slope the modern mass accumulation rate is 0·053 g cm^?2 year^?1.