CO2-induced dissolution of low permeability carbonates. Part II: Numerical modeling of experiments

We used the 3D continuum-scale reactive transport models to simulate eight core flood experiments for two different carbonate rocks. In these experiments the core samples were reacted with brines equilibrated with pCO₂ = 3, 2, 1, 0.5 MPa (Smith et al., 2013 [27]). The carbonate rocks were from specific Marly dolostone and Vuggy limestone flow units at the IEAGHG Weyburn-Midale CO₂ Monitoring and Storage Project in south-eastern Saskatchewan, Canada. Initial model porosity, permeability, mineral, and surface area distributions were constructed from micro tomography and microscopy characterization data. We constrained model reaction kinetics and porosity–permeability equations with the experimental data. The experimental data included time-dependent solution chemistry and differential pressure measured across the core, and the initial and final pore space and mineral distribution. Calibration of the model with the experimental data allowed investigation of effects of carbonate reactivity, flow velocity, effective permeability, and time on the development and consequences of stable and unstable dissolution fronts.The continuum scale model captured the evolution of distinct dissolution fronts that developed as a consequence of carbonate mineral dissolution and pore scale transport properties. The results show that initial heterogeneity and porosity contrast control the development of the dissolution fronts in these highly reactive systems. This finding is consistent with linear stability analysis and the known positive feedback between mineral dissolution and fluid flow in carbonate formations. Differences in the carbonate kinetic drivers resulting from the range of pCO₂ used in the experiments and the different proportions of more reactive calcite and less reactive dolomite contributed to the development of new pore space, but not to the type of dissolution fronts observed for the two different rock types. The development of the dissolution front was much more dependent on the physical heterogeneity of the carbonate rock. The observed stable dissolution fronts with small but visible dissolution fingers were a consequence of the clustering of a small percentage of larger pores in an otherwise homogeneous Marly dolostone. The observed wormholes in the heterogeneous Vuggy limestone initiated and developed in areas of greater porosity and permeability contrast, following pre-existing preferential flow paths.Model calibration of core flood experiments is one way to specifically constrain parameter input used for specific sites for larger scale simulations. Calibration of the governing rate equations and constants for Vuggy limestones showed that dissolution rate constants reasonably agree with published values. However the calcite dissolution rate constants fitted to the Marly dolostone experiments are much lower than those suggested by literature. The differences in fitted calcite rate constants between the two rock types reflect uncertainty associated with measured reactive surface area and appropriately scaling heterogeneous distribution of less abundant reactive minerals. Calibration of the power-law based porosity–permeability equations was sensitive to the overall heterogeneity of the cores. Stable dissolution fronts of the more homogeneous Marly dolostone could be fit with the exponent n = 3 consistent with the traditional Kozeny–Carman equation developed for porous sandstones. More impermeable and heterogeneous cores required larger n values (n = 6–8).     

Co-sequestration of SO2 with supercritical CO2 in carbonates: An experimental study of capillary trapping,relative permeability,and capillary pressure

In this study we performed three categories of steady- and unsteady-state core-flooding experiments to investigate capillary trapping, relative permeability, and capillary pressure, in a scCO₂ + SO₂/brine/limestone system at elevated temperature and pressure conditions, i.e., 60 °C and 19.16 MPa. We used a Madison limestone core sample acquired from the Rock Springs Uplift in southwest Wyoming. We carried out two sets of steady-state drainage-imbibition relative permeability experiments with different initial brine saturations to study hysteresis. We found that the final scCO₂ + SO₂ drainage relative permeability was very low, i.e., 0.04. We also observed a rapid reduction in the scCO₂-rich phase imbibition relative permeability curve, which resulted in a high residual trapping. The results showed that between 62.8% and more than 76% of the initial scCO₂ + SO₂ at the end of drainage was trapped by capillary trapping mechanism (trapping efficiency). We found that at higher initial brine saturations, the trapping efficiency was higher. The maximum initial and residual scCO₂-rich phase saturations at the end of primary drainage and imbibition were 0.525 and 0.329, respectively. Each drainage-imbibition cycle was followed by a dissolution process to re-establish S_w = 1. The dissolution brine relative permeabilities for both cycles were also obtained. We characterized the scCO₂ + SO₂/brine capillary pressure hysteresis behavior through unsteady-state primary drainage, imbibition, and secondary drainage experiments. We observed negative imbibition capillary pressure curve indicative of possible wettability alteration throughout the experiments due to contact with scCO₂ + SO₂/brine fluid system. The trapping results were compared to those reported in literature for other carbonate core samples. We noticed slightly more residual trapping in our sample, which might be attributed to heterogeneity, different viscosity ratio, and pore-space topologies. The impact of dynamic effects, i.e., high brine flow rate imbibition tests, on trapping of the scCO₂-rich phase was also explored. We performed two imbibition experiments with relatively high brine flow rates. The residual scCO₂ saturation dropped to 0.291 and 0.262 at the end of the first and second imbibition tests, i.e., 11.5% and 20.4%, respectively, compared to 0.329 under capillary-dominated regime.     

Effect of supercritical CO2 on carbonates: Savonnières sample case study

CO₂ geosequestration is an efficient way to reduce greenhouse gas emissions into the atmosphere. Carbonate rock formations are one of the possible targets for CO₂ sequestration due to their relative abundance and ability to serve as a natural trapping reservoir. The injected supercritical CO₂ can change properties of the reservoir rocks such as porosity, permeability, tortuosity, and specific surface area due to dissolution and precipitation processes. This, in turn, affects the reservoir characteristics, i.e., their elastic properties, storage capacity, stability, etc. The tremendous progresses made recently in both microcomputed X‐ray tomography and high‐performance computing make numerical simulation of physical processes on actual rock microstructures feasible. However, carbonate rocks with their extremely complex microstructure and the presence of microporosity that is below the resolution of microcomputed X‐ray tomography scanners require novel, quite specific image processing and numerical simulation approaches. In the current work, we studied the effects of supercritical CO₂ injection on microstructure and elastic properties of a Savonnières limestone. We used microtomographic images of two Savonnières samples, i.e., one in its natural state and one after injection and residence of supercritical CO₂. A statistical analysis of the microtomographic images showed that the injection of supercritical CO₂ led to an increase in porosity and changes of the microstructure, i.e., increase of the average volume of individual pores and decrease in the total number of pores. The CO₂ injection/residence also led to an increase in the mean radii of pore throats, an increase in the length of pore network segments, and made the orientation distribution of mesopores more isotropic. Numerical simulations showed that elastic moduli for the sample subjected to supercritical CO₂ injection/residence are lower than those for the intact sample.     

Mechanism for calcite dissolution and its contribution to development of reservoir porosity and permeability in the Kela 2 gas field,Tarim Basin,China

This study is undertaken to understand how calcite precipitation and dissolution contributes to depth-related changes in porosity and permeability of gas-bearing sandstone reservoirs in the Kela 2 gas field of the Tarim Basin, Northwestern China. Sandstone samples and pore water samples are col-lected from well KL201 in the Tarim Basin. Vertical profiles of porosity, permeability, pore water chem-istry, and the relative volume abundance of calcite/dolomite are constructed from 3600 to 4000 m below the ground surface within major oil and gas reservoir rocks. Porosity and permeability values are in-versely correlated with the calcite abundance, indicating that calcite dissolution and precipitation may be controlling porosity and permeability of the reservoir rocks. Pore water chemistry exhibits a sys-tematic variation from the Na2SO4 type at the shallow depth (3600-3630 m), to the NaHCO3 type at the intermediate depth (3630―3695 m),and to the CaCl2 type at the greater depth (3728―3938 m). The geochemical factors that control the calcite solubility include pH, temperature, pressure, Ca2 concen-tration, the total inorganic carbon concentration (ΣCO2), and the type of pore water. Thermodynamic phase equilibrium and mass conservation laws are applied to calculate the calcite saturation state as a function of a few key parameters. The model calculation illustrates that the calcite solubility is strongly dependent on the chemical composition of pore water, mainly the concentration difference between the total dissolved inorganic carbon and dissolved calcium concentration (i.e., [ΣCO2] -[Ca2 ]). In the Na2SO4 water at the shallow depth, this index is close to 0, pore water is near the calcite solubility. Calcite does not dissolve or precipitate in significant quantities. In the NaHCO3 water at the intermedi-ate depth, this index is greater than 0, and pore water is supersaturated with respect to calcite. Massive calcite precipitation was observed at this depth interval and this intensive cementation is responsible for decreased porosity and permeability. In the CaCl2 water at the greater depth, pore water is un-der-saturated with respect to calcite, resulting in dissolution of calcite cements, as consistent with microscopic dissolution features of the samples from this depth interval. Calcite dissolution results in formation of high secondary porosity and permeability, and is responsible for the superior quality of the reservoir rocks at this depth interval. These results illustrate the importance of pore water chemis-try in controlling carbonate precipitation/dissolution, which in turn controls porosity and permeability of oil and gas reservoir rocks in major sedimentary basins.     

In situ experiments and seismic analysis of existing buildings. Part I: experimental investigations

Recent results of in situ measurements and their interest for a seismic assessment of existing buildings are presented and analysed. The present paper (Part I) is devoted to the experimental programme. The response to ambient vibrations, harmonic excitation and shock loading is recorded on intact buildings but also after their structure or their vicinity was modified. These tests aim to identify the dynamic behaviour of ordinary intact buildings built in a conventional practise. Moreover, taking advantage of their demolition, it was possible (through these tests) to determine the actual influence of the light work elements, full precast facade panels, bearing masonry walls, and the presence of neighbouring joined buildings. These experiments realized on real buildings show that information gathered from ambient measurements provide reliable and efficient data of real interest for a clear understanding of the actual building behaviour. The advantage of integrating these data in the vulnerability assessment is presented and discussed in the next paper (Part II). Copyright © 2005 John Wiley & Sons, Ltd.     

Global atmospheric response to specific linear combinations of the main SST modes. Part I: numerical experiments and preliminary results

This article investigates through numerical experiments the controversial question of the impact of El NinCo-Southern Oscillation (ENSO) phenomena on climate according to large-scale and regional-scale interhemispheric thermal contrast. Eight experiments (two considering only inversed Atlantic thermal anomalies and six combining ENSO warm phase with large-scale interhemispheric contrast and Atlantic anomaly patterns) were performed with the Météo-France atmospheric general circulation model. The definition of boundary conditions from observed composites and principal components is presented and preliminary results concerning the month of August, especially over West Africa and the equatorial Atlantic are discussed. Results are coherent with observations and show that interhemispheric and regional scale sea-surface-temperature anomaly (SST) patterns could significantly modulate the impact of ENSO phenomena: the impact of warm-phase ENSO, relative to the atmospheric model intercomparison project (AMIP) climatology, seems stronger when embedded in global and regional SSTA patterns representative of the post-1970 conditions [i.e. with temperatures warmer (colder) than the long-term mean in the southern hemisphere (northern hemisphere)]. Atlantic SSTAs may also play a significant role.     

Analysis and interpretation of anomalous conductivity and magnetic permeability effects in time domain electromagnetic data: Part I: Numerical modeling

Time domain electromagnetic (TDEM) response is usually associated with eddy currents in conductive bodies, since this is the dominant effect. However, other effects, such as displacement currents from dielectric processes and magnetic fields associated with rock magnetization, can contribute to TDEM response. In this paper we analyze the effect of magnetization on TDEM data. We use a 3-D code based on finite-difference method, developed by Wang and Hohmann [Geophysics 58 (1993) 797], to study transient electromagnetic field propagation through a medium containing bodies with both anomalous conductivity and anomalous magnetic permeability. The remarkable result is that the combination of anomalous conductivity and permeability within the same body could increase significantly the anomalous TDEM response in comparison with purely conductive or purely magnetic anomalies. This effect has to be taken into account in interpretation of TDEM data over electrical inhomogeneous structures with potentially anomalous magnetic permeability.     

基于数字岩心的碳酸盐岩孔隙结构对弹性性质的影响研究(下篇):储层孔隙结构因子表征与反演

碳酸盐岩复杂的孔隙结构如何影响其弹性性质一直是地球物理研究的难点问题,在此基础上如何半定量甚至是定量地对碳酸盐岩储层预测,特别是如何有效地获取孔隙结构参数相关的地震属性体一直是油气工业界追求的目标,本研究从数字岩心角度入手,联合测井以及地震数据尝试探究这一问题的解决方案.首先针对代表不同孔隙结构类型的有限数目的碳酸盐岩样品获得其对应的高精度数字岩心数据体,为了获得更加可靠的具有地球物理含义的弹性性质随孔隙度变化的统计规律,我们通过子网格的技术,在有限数目的碳酸盐岩数字岩心数据体上获得了大量的数字岩心子网格样本,对于每个子网格样本可以分别获得其对应的数字岩心图像孔隙度、表征孔隙软硬程度的孔隙结构参数（γ）、以及基于有限元法模拟的弹性性质,由此基于数字岩心的研究思路,我们最终获得了基于孔隙结构因子表征与分类下的弹性性质与孔隙度的定量化解释量版.与此同时,在地震尺度上通过叠前地震资料获取的纵横波及密度属性体后,基于如上获得的定量化解释量版,我们最终获得了针对碳酸盐岩储层的新的属性体——孔隙结构参数（γ）属性体,这使得在地震尺度上预测碳酸盐岩储层的孔隙结构类型成为可能,也使利用地震数据在孔隙结构参数表征与分类下的碳酸盐岩储层反演精度的提高成为可能.

     

In situ experiments and seismic analysis of existing buildings. Part II: Seismic integrity threshold

The interest of in situ measurements (presented in Part I paper) for a seismic assessment of existing buildings is analysed in this paper. It is shown that the experimental modal characteristics obtained on regular concrete structures are described successfully by suited Timoshenko beam modelling. For a given structure, taking into account the experimental data, the corresponding beam model, and choosing the maximum tensile strain of concrete as damage criterion for key structural elements, a maximum level of the ground acceleration can be determined. This so‐called seismic integrity threshold is directly related to the onset of structural damages. This new approach is illustrated on one of the studied buildings. The advantages of using ambient vibrations survey for the vulnerability assessment of existing buildings are discussed. Copyright © 2005 John Wiley & Sons, Ltd.     

Identifying outliers and assessing the accuracy of amino acid racemization measurements for geochronology: I. Age calibration curves

Numerical ages derived from amino acid racemization (AAR) geochronology are typically based on calibration curves that relate the extent of AAR to the age of independently dated specimens. Here, we compare options for developing calibration curves and quantifying age uncertainties using AAR data from 481 late Holocene shells, and AMS ¹⁴C analyses of 36 shells of four molluscan taxa (Ethalia, Natica, Tellina, and Turbo) collected from shallow sediment cores from a back-reef lagoon of the central Great Barrier Reef. The four taxa differed substantially in the quality of their geochronogical results. Explicitly including data from specimens alive at the time of collection improves calibration curves, but weighting numerical ages based on their uncertainty has no effect. Calibration curve statistics do not adequately assess calibration uncertainty. The relation between ages inferred from different amino acids is recommended for identifying aberrant specimens and quantifying the uncertainty of inferred ages. For this study, the AAR ages based on two amino acids (aspartic acid and glutamic acid) exceed 200 yr or 20% of their mean inferred age in 15% of the specimens. Once these were removed, the mean age error (1σ) for individual specimens based on two amino acids analyzed in duplicate subsamples ranged from 53 to 142 yr for Tellina and Turbo, respectively, or about a 30% age error for these relatively young shells. This compares favorably with analytical errors estimated at 50 yr or 5%. The presence of notable outliers undetectable using data from single amino acids emphasizes the importance of analyzing multiple amino acids.     

Intercomparison of oceanic and atmospheric forced and coupled mesoscale simulations. Part I: Surface fluxes

A mesoscale non-hydrostatic atmospheric model has been coupled with a mesoscale oceanic model. The case study is a four-day simulation of a strong storm event observed during the SEMAPHORE experiment over a 500 × 500 km² domain. This domain encompasses a thermohaline front associated with the Azores current. In order to analyze the effect of mesoscale coupling, three simulations are compared: the first one with the atmospheric model forced by realistic sea surface temperature analyses; the second one with the ocean model forced by atmospheric fields, derived from weather forecast re-analyses; the third one with the models being coupled. For these three simulations the surface fluxes were computed with the same bulk parametrization. All three simulations succeed well in representing the main oceanic or atmospheric features observed during the storm. Comparison of surface fields with in situ observations reveals that the winds of the fine mesh atmospheric model are more realistic than those of the weather forecast re-analyses. The low-level winds simulated with the atmospheric model in the forced and coupled simulations are appreciably stronger than the re-analyzed winds. They also generate stronger fluxes. The coupled simulation has the strongest surface heat fluxes: the difference in the net heat budget with the oceanic forced simulation reaches on average 50 Wm⁻² over the simulation period. Sea surface-temperature cooling is too weak in both simulations, but is improved in the coupled run and matches better the cooling observed with drifters. The spatial distributions of sea surface-temperature cooling and surface fluxes are strongly inhomogeneous over the simulation domain. The amplitude of the flux variation is maximum in the coupled run. Moreover the weak correlation between the cooling and heat flux patterns indicates that the surface fluxes are not responsible for the whole cooling and suggests that the response of the ocean mixed layer to the atmosphere is highly non-local and enhanced in the coupled simulation.     

The role and significance of extracellular polymers in activated sludge. Part I: Literature review

In biological wastewater treatment, bacteria tend to aggregate by forming flocs, biofilms and even granules. A large part of the floc structure in activated sludge is composed of extracellular polymers. The extracellular polymeric substances (EPS) are mainly responsible for the structural and functional integrity of biofilms/flocs and they are considered to be important for the physico‐chemical properties of activated sludge. This review considers the composition of EPS and proportion between the constituents. Moreover it includes the relationship between the EPS constituents and sludge properties such as sludge retention time, hydrophobicity and surface charge.     

北京城区二氧化碳浓度和通量的梯度变化特征——I浓度与虚温

根据北京塔7层涡动系统2012年5月至2013年12月的湍流观测数据,分析了北京城区二氧化碳浓度在不同高度层次的日变化和月变化特征,并初步给出不同季节和日变化时间段内二氧化碳的浓度垂直廓线.结果表明:二氧化碳浓度整体随高度而下降;各观测层均有浓度的明显日变化,夏季最为明显,冬季相对平缓;近地层浓度直接受城市供暖、地表植被、交通运输等碳源影响,更高观测层浓度则受对流输送和天气过程影响较大;垂直方向上,冬季浓度变化范围最大,夏季层间浓度变化最明显;在一天中的任何时刻,近地面层二氧化碳浓度的日变化最低值一般出现在夏季,50m以上则出现在春季,浓度最高值总是出现在冬季;根据对二氧化碳浓度四季垂直廓线变化的分析可以看出,边界层二氧化碳浓度强烈受到碳源、下垫面植被、大气稳定度、环境温度和天气过程等因素的影响.     

Geochemistry of meteoric diagenesis in carbonate islands of the northern Bahamas: 2. Geochemical modelling and budgeting of diagenesis

Geochemical characterization and numerical modelling of surface water and ground water, combined with hydrological observations, provide quantitative estimates of meteoric diagenesis in Pleistocene carbonates of the northern Bahamas. Meteoric waters equilibrate with aragonite, but water‐ rather than mineral‐controlled reactions dominate. Dissolutional lowering of the undifferentiated bedrock surface is an order of magnitude slower than that within soil‐filled topographic hollows, generating small‐scale relief at a rate of 65–140 mm ka^?1 and a distinctive pocketed topography. Oxidation of organic matter within the subsoil and vadose zones generates an average P of 4·0 × 10^?3 atm, which drives dissolution during vadose percolation and/or at the water table. However, these dissolution processes together account for <60% of the average rock‐derived calcium in groundwaters pumped from the freshwater lens. The additional calcium may derive from oxidation of organic carbon within the lens, accounting for the high P of the lens waters. Mixing between meteoric waters of differing chemistry is diagenetically insignificant, but evapotranspiration from the shallow water table is an important drive for subsurface cementation. Porosity generation in the shallow vadose zone averages 1·6–3·2% ka^?1. Phreatic meteoric diagenesis is focused near the water table, where dissolution generates porosity at 1·4–2·8% ka^?1. Maximum dissolution rates, however, are similar to those of evaporation‐driven precipitation, which occludes porosity of 4·0 ± 0·6% ka^?1. This drives porosity inversion, from primary interparticle to secondary mouldic, vug and channel porosity. In the deeper freshwater lens, oxidation of residual organic carbon and reoxidation of reduced sulphur species from deeper anaerobic oxidation of organic carbon may generate porosity up to 0·06% ka^?1. Meteoric diagenesis relies critically on hydrological routing and vadose thickness (controlled by sea level), as well as the geochemical processes active. A thin vadose zone permits direct evaporation from the water table and drives precipitation of meteoric phreatic cements even where mineral stabilization is complete. Copyright © 2007 John Wiley & Sons, Ltd.     

Experimental study of identification and control of structures using neural network. Part 2: control

Experimental verifications of a recently developed active structural control method using neural networks are presented in this paper. The experiments were performed on the earthquake simulator at the University of Illinois at Urbana—Champaign. The test specimen was a 1/4 scale model of a three-storey building. The control system consisted of a tendon/pulley system controlled by a single hydraulic actuator at the base. The control mechanism was implemented through four active pre-tensioned tendons connected to the hydraulic actuator at the first floor. The structure modelling and system identification has been presented in a companion paper. (Earthquake Engng. Struct. Dyn. 28 , 995–1018 (1999)). This paper presents the controller design and implementation. Three controllers were developed and designed: two neurocontrollers, one with a single sensor feedback and the other with three sensor feedback, and one optimal controller with acceleration feedback. The experimental design of the neurocontrollers is accomplished in three steps: system identification, multiple emulator neural networks training and finally the neurocontrollers training with the aid of multiple emulator neural networks. The effectiveness of both neurocontrollers are demonstrated from experimental results. The robustness and the relative stability are presented and discussed. The experimental results of the optimal controller performance is presented and assessed. Comparison between the optimal controller and neurocontrollers is presented and discussed. Copyright © 1999 John Wiley & Sons, Ltd.     

Bridge pier failure probabilities under combined hazard effects of scour, truck and earthquake. Part I: occurrence probabilities

In many regions of the world, a bridge will experience multiple extreme hazards during its expected service life. The current American Association of State Highway and Transportation Officials (AASHTO) load and resistance factor design (LRFD) specifications are formulated based on failure probabilities, which are fully calibrated for dead load and nonextreme live loads. Design against earthquake loads is established separately. Design against scour effect is also formulated separately by using the concept of capacity reduction (or increased scour depth). Furthermore, scour effect cannot be linked directly to an LRFD limit state equation, because the latter is formulated using force-based analysis. This paper (in two parts) presents a probability-based procedure to estimate the combined hazard effects on bridges due to truck, earthquake and scour, by treating the effect of scour as an equivalent load effect so that it can be included in reliability-based bridge failure calculations. In Part I of this series, the general principle of treating the scour depth as an equivalent load effect is presented. The individual and combined partial failure probabilities due to truck, earthquake and scour effects are described. To explain the method of including non-force-based natural hazards effects, two types of common scour failures are considered. In Part II, the corresponding bridge failure probability, the occurrence of scour as well as simultaneously having both truck load and equivalent scour load are quantitatively discussed.     

Sensitivity of runoff and soil erosion to climate change in two Mediterranean watersheds. Part I: model parameterization and evaluation

Climate change is expected to effect storm runoff and erosion processes in Mediterranean watersheds at multiple spatial scales. Models are typically applied to estimate these impacts; however, the scarcity of spatially distributed data for parameterization, calibration and validation often prevents application of these models, particularly for larger catchments. This report, the first part of a two‐part article, presents an application and evaluation of the MEFIDIS model for two Mediterranean meso‐scale watersheds (115 and 290 km²) in a data‐scarce environment. A multi‐scale assessment method was used that combines quantitative validation and qualitative evaluation, consisting of three steps: (1) calibration at the small (field) scale using results from rainfall simulation experiments; (2) calibration and validation for catchment‐scale results while changing catchment‐scale parameters only (channel roughness and a parameter controlling the distribution of saturated areas); and (3) qualitative evaluation of within‐watershed erosion processes using empirical estimates of sediment delivery ratio and gully location. The results indicate that calibrating MEFIDIS at the field scale can provide reasonable results for catchment runoff and sediment export and for within‐watershed erosion processes. Copyright © 2009 John Wiley & Sons, Ltd.     

Paired‐basin comparison of hydrological response in harvested and undisturbed hardwood forests during snowmelt in central Ontario: I. Streamflow,groundwater and flowpath behaviour

Impacts of forest harvesting on groundwater properties, water flowpaths and streamflow response were examined 4 years after the harvest using a paired‐basin approach during the 2001 snowmelt in a northern hardwood landscape in central Ontario. The ability of two metrics of basin topography (Beven and Kirkby's ln(a/tan β) topographic index (TI) and distance to stream channel) to explain intra‐basin variations in groundwater dynamics was also evaluated. Significant relationships between TI and depth to potentiometric surface for shallow groundwater emerged, although the occurrence of these relationships during the melt differed between harvested and control basins, possibly as a result of interbasin differences in upslope area contributing to piezometers used to monitor groundwater behaviour. Transmissivity feedback (rapid streamflow increases as the water table approaches the soil surface) governed streamflow generation in both basins, and the mean threshold depths at which rapid streamflow increases corresponded to small rises in water level were similar for harvested (0·41 ± 0·05 m) and forested (0·38 ± 0·04 m) basins. However, topographic properties provided inconsistent explanations of spatial variations in the relationship between streamflow and depth to water at a given piezometer for both basins. Streamflow from the harvested basin exceeded that from the forested basin during the 2001 melt, and hydrometric and geochemical tracer results indicated greater runoff from the harvested basin via surface and near‐surface pathways. These differences are not solely attributable to harvesting, since the difference in spring runoff from the harvested basin relative to the forested control was not consistently larger than under pre‐harvest conditions. Nevertheless, greater melt rates following harvesting appear to have increased the proportion of water delivery to the stream channel via surface and near‐surface pathways. Copyright © 2005 John Wiley & Sons, Ltd.