Joint inversion of long‐offset and central‐loop transient electromagnetic data: Application to a mud volcano exploration in Perekishkul,Azerbaijan

Mud volcanism is commonly observed in Azerbaijan and the surrounding South Caspian Basin. This natural phenomenon is very similar to magmatic volcanoes but differs in one considerable aspect: Magmatic volcanoes are generally the result of ascending molten rock within the Earth's crust, whereas mud volcanoes are characterised by expelling mixtures of water, mud, and gas. The majority of mud volcanoes have been observed on ocean floors or in deep sedimentary basins, such as those found in Azerbaijan. Furthermore, their occurrences in Azerbaijan are generally closely associated with hydrocarbon reservoirs and are therefore of immense economic and geological interest. The broadside long‐offset transient electromagnetic method and the central‐loop transient electromagnetic method were applied to study the inner structure of such mud volcanoes and to determine the depth of a resistive geological formation that is predicted to contain the majority of the hydrocarbon reservoirs in the survey area. One‐dimensional joint inversion of central‐loop and long‐offset transient electromagnetic data was performed using the inversion schemes of Occam and Marquardt. By using the joint inversion models, a subsurface resistivity structure ranging from the surface to a depth of approximately 7 km was determined. Along a profile running perpendicular to the assumed strike direction, lateral resistivity variations could only be determined in the shallow depth range using the transient electromagnetic data. An attempt to resolve further two‐dimensional/three‐dimensional resistivity structures, representing possible mud migration paths at large depths using the long‐offset transient electromagnetic data, failed. Moreover, the joint inversion models led to ambiguous results regarding the depth and resistivity of the hydrocarbon target formation due to poor resolution at great depths (>5 km). Thus, 1D/2D modelling studies were subsequently performed to investigate the influence of the resistive terminating half‐space on the measured long‐offset transient electromagnetic data. The 1D joint inversion models were utilised as starting models for both the 1D and 2D modelling studies. The results tend to show that a resistive terminating half‐space, implying the presence of the target formation, is the favourable geological setting. Furthermore, the 2D modelling study aimed to fit all measured long‐offset transient electromagnetic Ex transients along the profile simultaneously. Consequently, 3125 2D forward calculations were necessary to determine the best‐fit resistivity model. The results are consistent with the 1D inversion, indicating that the data are best described by a resistive terminating half‐space, although the resistivity and depth cannot be determined clearly.     

The influence of induced polarization on long-offset transient electromagnetic data

The diffusion of electromagnetic fields is dependent not only on conductivity, but also on magnetic permeability, dielectric permittivity and polarizability, i.e. dispersive conductivity. The long‐offset transient electromagnetic (LOTEM) method is mainly used to determine the spatial distribution of conductivity in the subsurface. However, earlier work on loop‐loop TEM suggests that transient EM methods can also be affected by induced polarization (IP). Numerous 1D forward calculations were carried out to study the IP effect on LOTEM data, using the Cole‐Cole relaxation model to simulate the polarizability of the ground. Besides the polarizability of each layer, the IP effect depends on the LOTEM field set‐up and the spatial distribution of conductivity in the ground. In particular, near‐surface layers with high chargeabilities can significantly distort the late time transients of the electric field components in the vicinity of the transmitter. The influence of polarizable layers on the magnetic field components can be neglected under normal circumstances. In 1997 and 1999, LOTEM measurements were carried out at Mt. Vesuvius in Italy to explore the geological structure of the volcano. Sensitivity studies on the effect of polarizable layers suggest that high chargeabilities in connection with conductive layers at greater depths would result in a detectable distortion of the electric field transients. Although the simultaneous IP measurements revealed high chargeabilities in a near‐surface layer, no evidence of IP effects could be found in the measured LOTEM data. We conclude that the observed chargeabilities are local and that 3D effects are probably present in the data. Another aspect is the measurement of the system response, which is usually measured by placing a receiver very close to the transmitter. Therefore, large distortions can be expected if near‐surface polarizable layers exist. This was verified in practice by field measurements in an area with high chargeabilities in Longerich, Cologne.     

Multi‐transient electromagnetic repeatability experiment over the North Sea Harding field‡

We present results of synthetic time‐lapse and real repeatability multi‐transient electromagnetic surveys over the North Sea Harding field. Using Archie's law to convert porosity and fluid saturation to resistivity we created 3D isotropic models of the reservoir resistivity at different stages of production from the initial state in 1996 through to complete hydrocarbon production by 2016 and, for each stage, we simulated an east‐west transient electromagnetic survey line across Harding. Unconstrained 1D full‐waveform Occam inversions of these synthetic data show that Harding should be detectable and its lateral extent reasonably well‐defined. Resistivity changes caused by hydrocarbon production from initial pre‐production state to production of the oil rim in 2011 are discernible as are significant changes from 2011–2016 during the modelled gas blowdown phase. The 2D repeatability surveys of 2007 and 2008 tied two wells: one on and the other off the structure. Between the two surveys the segment of the field under investigation produced 3.9 million barrels of oil – not enough to generate an observable time‐lapse electromagnetic anomaly with a signal‐to‐noise ratio of 40 dB. Processing of the 2007 and 2008 data included deconvolution for the measured source current and removal of spatially‐correlated noise, which increased the signal‐to‐noise ratio of the recovered impulse responses by about 20 dB and resulted in a normalized root‐mean‐square difference of 3.9% between the data sets. 1D full‐waveform Occam inversions of the real data showed that Harding was detectable and its lateral extent was also reasonably well‐defined. The results indicate that the multi‐transient electromagnetic method is suitable for exploration, appraisal and monitoring hydrocarbon production.     

Analysing two‐dimensional effects in central loop transient electromagnetic sounding data using a semi‐synthetic tipper approach

We present a simple and feasible approach to analyse and identify two‐dimensional effects in central loop transient electromagnetic sounding data and the correspondingly derived quasi two‐dimensional conductivity models. The proposed strategy is particularly useful in minimising interpretation errors. It is based on the calculation of a semi‐synthetic transient electromagnetic tipper at each sounding and for each observational transient time point. The semi‐synthetic transient electromagnetic tipper is derived from the measured vertical component of the induced voltage and the synthetically calculated horizontal component. The approach is computationally inexpensive and involves one two‐dimensional forward calculation of an obtained quasi two‐dimensional conductivity section. Based on a synthetic example, we demonstrate that the transient electromagnetic tipper approach is applicable in identifying which transient data points and which corresponding zones in a derived quasi two‐dimensional subsurface model are affected by two‐dimensional inhomogeneities. The one‐dimensional inversion of such data leads to false models. An application of the semi‐synthetic transient electromagnetic tipper to field data from the Azraq basin in Jordan reveals that, in total, eight of 80 investigated soundings are affected by two‐dimensional structures although the field data can be fitted optimally using one‐dimensional inversion techniques. The largest semi‐synthetic tipper response occurs in a 300 m‐wide region around a strong lateral resistivity contrast. The approach is useful for analysing structural features in derived quasi two‐dimensional sections and for qualitatively investigating how these features affect the transient response. To avoid misinterpretation, these identified zones corresponding to large tipper values are excluded from the interpretation of a quasi two‐dimensional conductivity model. Based on the semi‐synthetic study, we also demonstrate that a quantitative interpretation of the horizontal voltage response (e.g. by inversion) is usually not feasible as it requires the exact sensor position to be known. Although a tipper derived purely from field data is useful as a qualitative tool for identifying two‐dimensional distortion effects, it is only feasible if the sensor setup is sufficiently accurate. Our proposed semi‐synthetic transient electromagnetic tipper approach is particularly feasible as an a posteriori approach if no horizontal components are recorded or if the sensor setup in the field is not sufficiently accurate.     

Three-dimensional numerical modeling of full-space transient electromagnetic responses of water in goaf

The full-space transient electromagnetic response of water-filled goaves in coal mines were numerically modeled. Traditional numerical modeling methods cannot be used to simulate the underground full-space transient electromagnetic field. We used multiple transmitting loops instead of the traditional single transmitting loop to load the transmitting loop into Cartesian grids. We improved the method for calculating the z-component of the magnetic field based on the characteristics of full space. Then, we established the fullspace 3D geoelectrical model using geological data for coalmines. In addition, the transient electromagnetic responses of water-filled goaves of variable shape at different locations were simulated by using the finite-difference time-domain (FDTD) method. Moreover, we evaluated the apparent resistivity results. The numerical modeling results suggested that the resistivity differences between the coal seam and its roof and floor greatly affect the distribution of apparent resistivity, resulting in nearly circular contours with the roadway head at the center. The actual distribution of apparent resistivity for different geoelectrical models of water in goaves was consistent with the models. However, when the goaf water was located in one side, a false low-resistivity anomaly would appear on the other side owing to the full-space effect but the response was much weaker. Finally, the modeling results were subsequently confirmed by drilling, suggesting that the proposed method was effective.     

A conceptual model of supra‐glacial lake formation on debris‐covered glaciers based on GPR facies analysis

Supra‐glacial lakes and ponds can create hotspots of mass loss on debris‐covered glaciers. While much research has been directed at understanding lateral lake expansion, little is known about the rates or processes governing lake deepening. To a large degree, this knowledge gap persists due to sparse observations of lake beds. Here we report on the novel use of ground penetrating radar (GPR) surveys to simultaneously collect supra‐glacial lake bathymetry and bottom composition data from Spillway Lake (surface area of 2.4 × 10⁵ m²; volume of 9.5 × 10⁴ m³), which is located in the terminus region of the Ngozumpa Glacier in the Khumbu region of the Nepal Himalaya. We identified two GPR bottom signals corresponding to two sedimentary facies of (1) sub‐horizontal layered fine sediment drape and (2) coarse blocky diamict. We provide an understanding of the changes in subaqueous debris distribution that occur through stages of lake expansion by combining the GPR results with in situ observations of shoreline deposits matching the interpreted facies. From this, we present an updated conceptual model of supra‐glacial lake evolution, with the addition of data on the evolving debris environment, showing how dominant depositional processes can change as lakes evolve from perched lakes to multi‐basin base‐level lakes and finally onto large moraine‐dammed lakes. Throughout lake evolution, processes such as shoreline steepening, lakebed collapse into voids and conduit interception, subaerial and subaqueous calving and rapid areal expansion alter the spatial distribution and makeup of lakebed debris and sediments forcing a number of positive and negative feedbacks on lake expansion. Copyright © 2016 John Wiley & Sons, Ltd.     

Field test of sub‐basalt hydrocarbon exploration with marine controlled source electromagnetic and magnetotelluric data

The recent use of marine electromagnetic technology for exploration geophysics has primarily focused on applying the controlled source electromagnetic method for hydrocarbon mapping. However, this technology also has potential for structural mapping applications, particularly when the relative higher frequency controlled source electromagnetic data are combined with the lower frequencies of naturally occurring magnetotelluric data. This paper reports on an extensive test using data from 84 marine controlled source electromagnetic and magnetotelluric stations for imaging volcanic sections and underlying sediments on a 128‐km‐long profile. The profile extends across the trough between the Faroe and Shetland Islands in the North Sea. Here, we focus on how 2.5D inversion can best recover the volcanic and sedimentary sections. A synthetic test carried out with 3D anisotropic model responses shows that vertically transverse isotropy 2.5D inversion using controlled source electromagnetic and magnetotelluric data provides the most accurate prediction of the resistivity in both volcanic and sedimentary sections. We find the 2.5D inversion works well despite moderate 3D structure in the synthetic model. Triaxial inversion using the combination of controlled source electromagnetic and magnetotelluric data provided a constant resistivity contour that most closely matched the true base of the volcanic flows. For the field survey data, triaxial inversion of controlled source electromagnetic and magnetotelluric data provides the best overall tie to well logs with vertically transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data a close second. Vertical transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data provided the best interpreted base of the volcanic horizon when compared with our best seismic interpretation. The structural boundaries estimated by the 20‐Ω·m contour of the vertical resistivity obtained by vertical transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data gives a maximum geometric location error of 11% with a mean error of 1.2% compared with the interpreted base of the volcanic horizon. Both the model study and field data interpretation indicate that marine electromagnetic technology has the potential to discriminate between low‐resistivity prospective siliciclastic sediments and higher resistivity non‐prospective volcaniclastic sediments beneath the volcanic section.     

瞬变电磁法正演计算进展

详细介绍了瞬变电磁法正演计算的方法、现状和发展趋势.瞬变电磁法一维正演计算需要将电磁场从频率域转换至时间域,转换方法有三种,分别是Gaver-Stehfest算法、余弦变换和Guptasarma算法.在这三种方法中,使用较多的是Gaver-Stehfest算法和余弦变换,Gaver-Stehfest算法速度较快,但精度不及余弦变换.瞬变电磁法的数值模拟主要集中于2.5维和三维,使用的数值计算方法有积分方程法、有限差分法、有限单元法和SLDM法.积分方程法主要在三维数值模拟中使用,现已很少使用;有限差分法和有限单元法是目前瞬变电磁法2.5维和三维数值模拟的主要方法;SLDM法主要应用于三维数值模拟.我国瞬变电磁法正演计算成果主要集中在回线源激发的瞬变电磁场一维数值计算和利用有限单元法进行2.5维和三维数值模拟.瞬变电磁法正演计算的发展趋势有:数值算法的改进、提高计算效率和研究地形对瞬变电磁场的影响规律.     

Estimating daily inflows of large lakes using a water‐balance‐based runoff coefficient scaling approach

Daily river inflow time series are highly valuable for water resources and water environment management of large lakes. However, the availability of continuous inflow data for large lakes is still relatively limited, especially for large lakes situated within humid plain regions with tens or even hundreds of tributaries. In this study, we choose the fifth largest freshwater Lake Chaohu in China as our study area to introduce a new approach to reconstruct historical daily inflows at ungauged subcatchments of large lakes. This approach makes use of water level, lake surface rainfall, evaporation from the lake, and catchment rainfall observations. Rainfall–runoff relationship at a reference catchment was analysed to select rainfall input and estimate run‐off coefficient firstly, and the run‐off coefficient was then transferred to ungauged subcatchments to initially estimate daily inflows. Run‐off coefficient was scaled to adjust daily inflows at ungauged subcatchments according to water balance of the lake. This approach was evaluated using sparsely measured inflows at eight subcatchments of Lake Chaohu and compared with the commonly used drainage area ratio method. Results suggest that the inflow time series reconstructed from this approach consistent well to corresponding observations, with mean R² and Nash–Sutcliffe efficiency values of 0.69 and 0.6, respectively. This approach outperforms drainage area ratio method in terms of mean R² and Nash–Sutcliffe efficiency values. Accuracy of this approach holds well when the number of water‐level station being used decreased from four to one.     

瞬变电磁法中心回线装置资料解释方法的改进

中心回线装置是瞬变电磁勘探中最常用的装置之一.中心回线装置的视电阻率一般从回线中心点场的公式出发导出.在现场实际施工中,为提高工作效率,把发射回线中部三分之一的区域作为观测范围.通过对场分布特征的分析研究表明,与中心点相比,观测区边缘处的感应电动势数值偏离达15%~25%,这与广泛存在于华北型煤系中赋存深度为400~1000 m的陷落柱、导水小断层等引起的异常相比,已经不可忽略.大定源回线公式可以准确地表示任意场点的感应电动势,由此导出视电阻率无边缘效应影响;观测按照中心回线方式并保持在近区进行,可使感应电动势和视电阻率之间的转换简单而直接.理论和实际应用结果表明,这些措施进一步提高了瞬变电磁探测中心回线方法的探测精度.在有上覆低阻屏蔽层的情况下,对埋深500~700 m的陷落柱给出了明显的异常反映,并被井下掘进所验证.     

Trapping of aeolian sediments and build‐up of the ice cover of a dry‐based Antarctic lake

Among the perennially frozen lakes of the Dry Valleys of South Victoria Land (Antarctica), some are dry‐based, i.e. frozen to the bottom. One of these is studied here by a multiparametric investigation (isotopic composition in δD and δ¹⁸O, ions, gas and ice texture analyses). A sediment layer about 10 cm thick appearing at a depth of 3·5 m is also studied by grain size, X‐ray diffraction and scanning electron microscope analyses. The information retrieved indicates that this ice‐block lake results from a build‐up in two steps and explains how aeolian sediments were included as a layer into the ice. Copyright © 2002 John Wiley & Sons, Ltd.     

Small fixed-loop transient electromagnetic in tunnel forward geological prediction

Tunnel forward geological prediction is indispensable to tunnel construction. The transient electromagnetic method is widely used in the field of tunnel forward geological prediction. A small fixed-loop transient electromagnetic method adapted to the limited space of a tunnel face is proposed for tunnel forward geological prediction for the first time. This method makes use of a large number of non-central points inside the loop. The obtained numerical simulation and field data show opposite trends owing to mutual inductance and self-inductance. A correction method for the field value of non-central points is proposed. The influences of different parameters are verified and the correction function is established by a test method. The validity of the correction function is verified through controlled field tests. The small fixed-loop transient electromagnetic has been successfully applied to a tunnel under the Weiwu expressway. This method can effectively improve the detection depth and better detect unfavourable geological bodies, such as karsts in front of the tunnel face.     

Coil response inversion for very early time modelling of helicopter‐borne time‐domain electromagnetic data and mapping of near‐surface geological layers

Very early times in the order of 2–3 μs from the end of the turn‐off ramp for time‐domain electromagnetic systems are crucial for obtaining a detailed resolution of the near‐surface geology in the depth interval 0–20 m. For transient electromagnetic systems working in the off time, an electric current is abruptly turned off in a large transmitter loop causing a secondary electromagnetic field to be generated by the eddy currents induced in the ground. Often, however, there is still a residual primary field generated by remaining slowly decaying currents in the transmitter loop. The decay disturbs or biases the earth response data at the very early times. These biased data must be culled, or some specific processing must be applied in order to compensate or remove the residual primary field. As the bias response can be attributed to decaying currents with its time constantly controlled by the geometry of the transmitter loop, we denote it the ‘Coil Response’. The modelling of a helicopter‐borne time‐domain system by an equivalent electronic circuit shows that the time decay of the coil response remains identical whatever the position of the receiver loop, which is confirmed by field measurements. The modelling also shows that the coil response has a theoretical zero location and positioning the receiver coil at the zero location eliminates the coil response completely. However, spatial variations of the coil response around the zero location are not insignificant and even a few cm deformation of the carrier frame will introduce a small coil response. Here we present an approach for subtracting the coil response from the data by measuring it at high altitudes and then including an extra shift factor into the inversion scheme. The scheme is successfully applied to data from the SkyTEM system and enables the use of very early time gates, as early as 2–3 μs from the end of the ramp, or 5–6 μs from the beginning of the ramp. Applied to a large‐scale airborne electromagnetic survey, the coil response compensation provides airborne electromagnetic methods with a hitherto unseen good resolution of shallow geological layers in the depth interval 0–20 m. This is proved by comparing results from the airborne electromagnetic survey to more than 100 km of Electrical Resistivity Tomography measured with 5 m electrode spacing.     

Pressure‐driven,shoreline currents in a perennially ice‐covered,pro‐glacial lake in Antarctica,identified from a LiCl tracer injected into a pro‐glacial stream

The distribution of streamwater within ice‐covered lakes influences sub‐ice currents, biological activity and shoreline morphology. Perennially ice‐covered lakes in the McMurdo Dry Valleys, Antarctica, provide an excellent natural laboratory to study hydrologic–limnologic interactions under ice cover. For a 2 h period on 17 December 2012, we injected a lithium chloride tracer into Andersen Creek, a pro‐glacial stream flowing into Lake Hoare. Over 4 h, we collected 182 water samples from five stream sites and 15 ice boreholes. Geochemical data showed that interflow travelled West of the stream mouth along the shoreline and did not flow towards the lake interior. The chemistry of water from Andersen Creek was similar to the chemistry of water below shoreline ice. Additional evidence for Westward flow included the morphology of channels on the ice surface, the orientation of ripple marks in lake sediments at the stream mouth and equivalent temperatures between Andersen Creek and water below shoreline ice. Streamwater deflected to the right of the mouth of the stream, in the opposite direction predicted by the Coriolis force. Deflection of interflow was probably caused by the diurnal addition of glacial runoff and stream discharge to the Eastern edge of the lake, which created a strong pressure gradient sloping to the West. This flow directed stream momentum away from the lake interior, minimizing the impact of stream momentum on sub‐ice currents. It also transported dissolved nutrients and suspended sediments to the shoreline region instead of the lake interior, potentially affecting biological productivity and bedform development. Copyright © 2014 John Wiley & Sons, Ltd.     

Investigation into the use of 2D elastic waveform inversion from look‐ahead walk‐away VSP surveys

Pore pressure in sediments beneath salt in the Gulf of Mexico varies widely creating a potentially dangerous and difficult drilling challenge. Estimating elastic parameters of sediments beneath salt is key to the prediction of pore pressure and reducing the drilling risk in exiting the base of the salt. In this paper we investigate the ability of 2D waveform inversion to recover the elastic parameters in the sedimentary layer beneath the salt from a walk‐away VSP (vertical seismic profile) carried out with the receivers in the salt, with the objective of estimating pore pressure at the base of the salt (to be estimated using traditional methods). We propose an effective method for performing the inversion and apply this method to a blind test of a large and realistic synthetic dataset. To facilitate the design of a VSP survey suitable to this type of inversion we also present an analysis into the effects of the receiver array receiver spacing. It is shown that the resulting velocity estimates are sufficiently accurate to predict the pore pressure within the limits required for drilling.     

Multi‐method virtual electromagnetic experiments for developing suitable monitoring designs: A fictitious CO2 sequestration scenario in Northern Germany

We present a numerical study for 3D time‐lapse electromagnetic monitoring of a fictitious CO₂ sequestration using the geometry of a real geological site and a suite of suitable electromagnetic methods with different source/receiver configurations and different sensitivity patterns. All available geological information is processed and directly implemented into the computational domain, which is discretized by unstructured tetrahedral grids. We thus demonstrate the performance capability of our numerical simulation techniques. The scenario considers a CO₂ injection in approximately 1100 m depth. The expected changes in conductivity were inferred from preceding laboratory measurements. A resistive anomaly is caused within the conductive brines of the undisturbed reservoir horizon. The resistive nature of the anomaly is enhanced by the CO₂ dissolution regime, which prevails in the high‐salinity environment. Due to the physicochemical properties of CO₂, the affected portion of the subsurface is laterally widespread but very thin. We combine controlled‐source electromagnetics, borehole transient electromagnetics, and the direct‐current resistivity method to perform a virtual experiment with the aim of scrutinizing a set of source/receiver configurations with respect to coverage, resolution, and detectability of the anomalous CO₂ plume prior to the field survey. Our simulation studies are carried out using the 3D codes developed in our working group. They are all based on linear and higher order Lagrange and Nédélec finite‐element formulations on unstructured grids, providing the necessary flexibility with respect to the complex real‐world geometry. We provide different strategies for addressing the accuracy of numerical simulations in the case of arbitrary structures. The presented computations demonstrate the expected great advantage of positioning transmitters or receivers close to the target. For direct‐current geoelectrics, 50% change in electric potential may be detected even at the Earth's surface. Monitoring with inductive methods is also promising. For a well‐positioned surface transmitter, more than 10% difference in the vertical electric field is predicted for a receiver located 200 m above the target. Our borehole transient electromagnetics results demonstrate that traditional transient electromagnetics with a vertical magnetic dipole source is not well suited for monitoring a thin horizontal resistive target. This is due to the mainly horizontal current system, which is induced by a vertical magnetic dipole.     

瞬变电磁法的探测深度问题

用解析分析、时域有限差分、时-频分析的方法,以地面中心回线装置和阶跃脉冲激励源为例,分析讨论了瞬变电磁测深法的勘探深度问题,以便为野外勘探工作设计提供依据,达到预期的探测目的.解析计算证实了瞬变场在地下以有限速度传播,数值模拟表示出了准静态条件下瞬变场的反射.研究结果表明,由于时间域电磁场遵循因果律,瞬变电磁法的探测深度主要由观测时间决定. 瞬变电磁场的初始传播速度与大地电阻率无关,继后在大地色散作用下,阶跃脉冲前沿逐渐变得平缓,各频率分量的传播速度与电阻率有关,在低阻地层中探测同样的深度需要较长的观测时间. 最大探测深度是在给定时间内电磁波往返地下某一深度的单程距离,最小探测深度受仪器性能的限制,但是埋藏较浅的异常体也有可能在晚时段被观测到.从时-频密度谱中可得到瞬变电磁场信号时间与频率的关系.     

Spatial patterns of groundwater‐lake exchange – implications for acid neutralization processes in an acid mine lake

Exchange of groundwater and lake water with typically quite different chemical composition is an important driver for biogeochemical processes at the groundwater‐lake interface, which can affect the water quality of lakes. This is of particular relevance in mine lakes where anoxic and slightly acidic groundwater mixes with oxic and acidic lake water (pH < 3). To identify links between groundwater‐lake exchange rates and acid neutralization processes in the sediments, exchange rates were quantified and related to pore‐water pH, sulfate and iron concentrations as well as sulfate reduction rates within the sediment. Seepage rates measured with seepage meters (?2.5 to 5.8 L m^‐2 d^‐1) were in reasonable agreement with rates inverted from modeled chloride profiles (?1.8 to 8.1 L m^‐2 d^‐1). Large‐scale exchange patterns were defined by the (hydro)geologic setting but superimposed by smaller scale variations caused by variability in sediment texture. Sites characterized by groundwater upwelling (flow into the lake) and sites where flow alternated between upwelling and downwelling were identified. Observed chloride profiles at the alternating sites reflected the transient flow regime. Seepage direction, as well as seepage rate, were found to influence pH, sulfate and iron profiles and the associated sulfate reduction rates. Under alternating conditions proton‐consuming processes, for example, sulfate reduction, were slowed. In the uppermost layer of the sediment (max. 5 cm), sulfate reduction rates were significantly higher at upwelling (>330 nmol g^‐1 d^‐1) compared to alternating sites (<220 nmol g^‐1 d^‐1). Although differences in sulfate reduction rates could not be explained solely by different flux rates, they were clearly related to the prevailing groundwater‐lake exchange patterns and the associated pH conditions. Our findings strongly suggest that groundwater‐lake exchange has significant effects on the biogeochemical processes that are coupled to sulfate reduction such as acidity retention and precipitation of iron sulfides. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydrologic characteristics of lake‐ and stream‐side riparian wetted margins in the McMurdo Dry Valleys,Antarctica

Water is a limiting factor for life in the McMurdo Dry Valleys (MDV), Antarctica. The active layer (seasonally thawed soil overlying permafrost) accommodates dynamic hydrological and biological processes for 10–16 weeks per year. Wetted margins (visually wetted areas with high moisture content) adjacent to lakes and streams are potential locations of great importance in the MDV because of the regular presence of liquid water, compared with the rest of the landscape where liquid water is rare. At 11 plots (four adjacent to lakes, seven adjacent to streams), soil particle size distribution, soil electrical conductivity, soil water content and isotopic signature, width of the wetted margin, and active layer thaw depth were characterised to determine how these gradients influence physicochemical properties that determine microbial habitat and biogeochemical cycling. Sediments were generally coarse‐grained in wetted margins adjacent to both lakes and streams. Wetted margins ranged from 1·04 to 11·01 m in average length and were found to be longer at lakeside sites than streamside. Average thaw depths ranged from 0·12 to 0·85 m, and were found to be deepest under lake margins. Lake margins also had much higher soil electrical conductivity, steeper topographic gradients, but more gradual soil moisture gradients than stream margins. Patterns of soil water δ¹⁸O and δD distribution indicate capillary action and evaporation from wetted margins; margin pore waters generally demonstrated isotopic enrichment with distance from the shore, indicating evaporation of soil water. Lake margin pore waters were significantly more negative in D_XS (D_XS = δD‐8δ¹⁸O) than streamside pore waters, indicating a longer history of evaporation there. Differences between lake and stream margins can be explained by the more consistent availability of water to lake margins than stream margins. Differences in margin characteristics between lakes and streams have important consequences for the microbial habitat of these margins and their functional role in biogeochemical cycling at these terrestrial–aquatic interfaces. Copyright © 2009 John Wiley & Sons, Ltd.     

Three‐dimensional inversion of towed streamer electromagnetic data

A towed streamer electromagnetic system capable of simultaneous seismic and electromagnetic data acquisition has recently been developed and tested in the North Sea. We introduce a 3D inversion methodology for towed streamer electromagnetic data that includes a moving sensitivity domain. Our implementation is based on the 3D integral equation method for computing responses and Fréchet derivatives and uses the re‐weighted regularized conjugate gradient method for minimizing the objective functional with focusing regularization. We present two model studies relevant to hydrocarbon exploration in the North Sea. First, we demonstrate the ability of a towed electromagnetic system to detect and characterize the Harding field, a medium‐sized North Sea hydrocarbon target. We compare our 3D inversion of towed streamer electromagnetic data with 3D inversion of conventional marine controlled‐source electromagnetic data and observe few differences between the recovered models. Second, we demonstrate the ability of a towed streamer electromagnetic system to detect and characterize the Peon discovery, which is representative of an infrastructure‐led shallow gas play in the North Sea. We also present an actual case study for the 3D inversion of towed streamer electromagnetic data from the Troll field in the North Sea and demonstrate our ability to image all the Troll West Oil and Gas Provinces and the Troll East Gas Province. We conclude that 3D inversion of data from the current generation of towed streamer electromagnetic systems can adequately recover hydrocarbon‐bearing formations to depths of approximately 2 km. We note that by obviating the need for ocean‐bottom receivers, the towed streamer electromagnetic system enables electromagnetic data to be acquired over very large areas in frontier and mature basins for higher acquisition rates and relatively lower cost than conventional marine controlled‐source electromagnetic methods.