Hostname: page-component-586b7cd67f-t7fkt Total loading time: 0 Render date: 2024-11-23T03:15:30.077Z Has data issue: false hasContentIssue false

In situ measurements of the direct-current conductivity of Antarctic Sea ice: implications for airborne electromagnetic Sounding of Sea-ice thickness

Published online by Cambridge University Press:  14 September 2017

J.E. Reid
Affiliation:
School of Earth Sciences, University of Tasmania, Private Bag 79, Hobart, Tasmania 7001, Australia, E-mail: james.reid@utas.edu.au
A. Pfaffling
Affiliation:
Alfred Wegener Institute for Polar and Marine Research, PO Box 120161, D-27515 Bremerhaven, Germany
A.P. Worby
Affiliation:
Antarctic Climate and Ecosystems Cooperative Research Centre and Australian Antarctic Division, Private Bag 80, Hobart, Tasmania 7001, Australia
J.R. Bishop
Affiliation:
Mitre Geophysics Pty Ltd, PO Box 974, Sandy Bay, Tasmania 7006, Australia
Rights & Permissions [Opens in a new window]

Abstract

Airborne, Ship-borne and Surface low-frequency electromagnetic (EM) methods have become widely applied to measure Sea-ice thickness. EM responses measured over Sea ice depend mainly on the Sea-water conductivity and on the height of the Sensor above the Sea-ice–sea-water interface, but may be Sensitive to the Sea-ice conductivity at high excitation frequencies. We have conducted in Situ measurements of direct-current conductivity of Sea ice using Standard geophysical geoelectrical methods. Sea-ice thickness estimated from the geoelectrical Sounding data was found to be consistently underestimated due to the pronounced vertical-to-horizontal conductivity anisotropy present in level Sea ice. At five Sites, it was possible to determine the approximate horizontal and vertical conductivities from the Sounding data. The average horizontal conductivity was found to be 0.017 Sm–1, and that in the vertical direction to be 9–12 times higher. EM measurements over level Sea ice are Sensitive only to the horizontal conductivity. Numerical modelling has Shown that the assumption of zero Sea-ice conductivity in interpretation of airborne EM data results in a negligible error in interpreted thickness for typical level Antarctic Sea ice.

Type
Research Article
Copyright
Copyright © The Author(s) [year] 2006 

Introduction

Low-frequency geophysical airborne electromagnetic (em) methods have become an important technique for remotely estimating Sea-ice thickness. airborne em methods have been used extensively in the arctic (Reference Kovacs, Valleau and HolladayKovacs and others, 1987; Reference Kovacs and HolladayKovacs and holladay, 1990; Reference Liu and BeckerLiu and becker, 1990; Reference Haas, Edeler, Schürmann, Lobach and SengpielHaas and others, 2002; Reference Prinsenberg, Holladay, Lee, Chung, Sayed, Kashiwagi, Setoguchi and HongPrinsenberg and others, 2002), and Surveys have recently also been conducted in antarctica (Reference Pfaffling, Haas and ReidPfaffling and others, 2004). airborne em methods provide reliable estimates of the thickness of level Sea ice: the alfred wegener institute (awi) helicopter em (hem-bird) System considered in this Study is able to determine modal thickness to within 10 cm for Sea ice up to 3 m thick. Sea-ice thickness may, however, be Significantly over- or underestimated in areas of deformed Sea ice (e.g. pressure ridges), both as a result of the relatively poor lateral resolution of low-frequency em methods (e.g. Reference Kovacs, Holladay and BergeronKovacs and others, 1995), and because data are routinely interpreted using Simple, rapid ‘layered-earth’ interpretation algorithms.

Low-frequency em techniques for measuring Sea-ice thickness rely on the fact that the electrical conductivity of Sea ice is typically one to two orders of magnitude Smaller than that of polar Sea water. the em transmitter therefore mainly induces eddy currents in the conductive Sea water, and Secondary magnetic fields measured at the em receiver depend mainly on the altitude of the em System above the Sea-ice–sea-water interface, and on the Sea-water conductivity. Sea-ice thickness (or Snow plus Sea-ice thickness) is determined by Subtracting the height of the em System above the Surface of the Sea ice or Snow, measured using a laser altimeter, from the depth to Sea water estimated from the em data. em Systems employed for Sea-ice thickness measurements utilize transmitter frequencies of a few hundred hz to 200 khz. in most interpretations of airborne em Sea-ice thickness data, it is assumed that the conductivity of the Sea ice is So low that it has no influence on the measured responses (e.g. Reference Pfaffling, Haas and ReidPfaffling and others, 2004). however, at higher transmitter frequencies (>50 khz), the measured airborne em response may also be affected by the finite conductivity of the Sea ice. the possibility exists that, in addition to thickness, Sea-ice conductivity can also be recovered from high-frequency airborne em data. Sea-ice conductivity data could potentially be used to derive information on the bulk brine volume of Sea ice, and hence to estimate its mechanical properties, as discussed by Reference KovacsKovacs (1996).

The ability of airborne em methods to determine Sea-ice conductivity depends on the conductivity and thickness of the Sea ice, and on the noise levels of the em System. in this Study, we use in Situ direct-current (d.c.) geoelectrical methods to characterize the electrical conductivity of east antarctic pack ice, in order to constrain interpretation of airborne em Sea-ice thickness data acquired during September–october 2003 as part of the antarctic remote ice Sensing experiment (arise). airborne em data were acquired using the awi hem-bird System (Reference Haas, Edeler, Schürmann, Lobach and SengpielHaas and others, 2002; Reference Pfaffling, Haas and ReidPfaffling and others, 2004). theoretical modelling has been conducted in order to investigate the possibility of recovering Sea-ice conductivity from practical awi hem-bird em data.

D.c. electrical conductivity of Sea ice

Sea ice is an inhomogeneous composite of pure ice, brine, air and, at low temperatures, precipitated Salts. a number of authors have reported results of in Situ conductivity (or resistivity) measurements on arctic and antarctic Sea ice. Reference Thyssen, Kohnen, Cowan and TimcoThyssen and others (1974), Reference KohnenKohnen (1976), Reference TimcoTimco (1979) and buckley and others (1986) report results of traditional d.c. resistivity measurements on Sea ice. Reference Morey, Kovacs and CoxMorey and others (1984) conducted time-domain reflectometry (tdr) measurements using ‘ladders’ of transmission lines frozen into growing arctic Sea ice, from which they were able to determine the d.c. conductivity. Reference Becker, Echert and G.B.Becker and others (1992) determined in Situ Sea-ice conductivity using high-frequency em methods (8–16 mhz) via measurements of the radiation resistance of a circular loop antenna placed on the Surface of the ice.

Columnar Sea ice exhibits a preferred fabric, with the orientation of the c axes of the ice crystals predominantly horizontal, perpendicular to the direction of Sea-ice growth (Reference TimcoTimco, 1979; Reference Weeks, Ackley and UntersteinerWeeks and ackley, 1986). brine inclusions within the Sea ice occur as thread-like cells arranged in thin, vertically oriented layers. the Shape of the brine cells is temperature-dependent: the cells are needle-shaped at low temperatures, but become disc-shaped at very warm temperatures (Reference TimcoTimco, 1979). as the Sea ice grows thicker, the entrapped brine drains out under gravity, forming vertically aligned drainage tubes. drainage tube diameters are typically 0.1–1 cm, although diameters of up to 10cm have been reported (Reference Weeks, Ackley and UntersteinerWeeks and ackley, 1986). under certain temperature conditions, the drainage tubes may extend through the entire thickness of the Sea ice (Reference GoldenGolden, 2001).

The vertical orientation of both the brine cells and drainage tubes imparts a Strong conductivity anisotropy to the Sea ice, with the bulk conductivity in the vertical direction being higher than that in the horizontal direction (Reference Thyssen, Kohnen, Cowan and TimcoThyssen and others, 1974). Reference Watson and BarkerThyssen and others (1974) verified the presence of this conductivity anisotropy in Situ, by making direct measurements of the horizontal and vertical conductivities in the walls of a pit excavated in thick, undeformed arctic Sea ice. buckley and others (1986) detected a Similar anisotropy in first-year antarctic Sea ice.

Results

Direct-current measurements of Sea-ice conductivity were made using a wenner electrode array to conduct geoelectrical Soundings (fig. 1a). Sixteen measurements were made on first-year pack ice within the region bounded by latitudes 63˚56.2' S and 65˚14.3' S and longitudes 109˚27.3' e and 117˚44.5' e. at each measurement Site, the electrically insulating Snow cover was removed, and resistance measurements were made at a Series of electrode Separations ranging from a = 0.1 m to a = 10 m(fig. 1a) with the midpoint of the array fixed in position. for the purposes of data interpretation, the measured resistances were then transformed to apparent resistivity (ρ a) according to the equation

Fig. 1. (a) Wenner electrode array. Current I is injected through current electrodes C1 and C2. The potential difference, Δv, resulting from the Subsurface current flow is measured between potential electrodes P1 and P2. The ratio Δv/I is the ground resistance R in ohms. (b) Comparison of observed and calculated data for a typical Wenner-array geoelectrical Sounding on Sea ice (site 1b, Table 1). Resistances measured at each electrode Spacing have been transformed to apparent conductivity, as described in the text. The calculated curve is the apparent-conductivity response of the best-fit three-layered model obtained by least-squares inversion of the observed data.

(1)

Where a is the electrode Spacing (m) and R is the measured resistance (Ω).

The apparent-resistivity vs electrode-separation data were then interpreted assuming a one-dimensional (layered) model, in which conductivity was assumed to be isotropic and to vary only with depth. the layered conductivity model was obtained using an iterative least-squares inversion process, in which the parameters of a Starting model were automatically adjusted until the apparent-resistivity curve computed for the model matched the field data. figure 1b Shows a typical comparison of observed data (symbols) with the calculated response of the best-fit three-layered model (dashed line) obtained via one-dimensional (1-d) inversion. note that observed and calculated data have been plotted as apparent conductivity (σ a=1/ρ a).

For the majority of the d.c. wenner-array measurements, the interpreted three-layered model comprised a thin, conductive near-surface layer, a thicker, less conductive middle layer and a highly conductive lower layer (sea water). this conductivity distribution is typical of that for young (first-year) Sea ice, and has been previously observed in both the arctic (Reference Thyssen, Kohnen, Cowan and TimcoThyssen and others, 1974) and antarctic (Reference Buckley, Staines and RobinsonBuckley and others, 1986). if the Sea-ice conductivity were isotropic, the Sum of the interpreted thicknesses of the upper two layers would yield the depth to Sea water below the Surface, i.e. Sea-ice thickness. figure 2 compares the Sea-ice thicknesses determined from the d.c. geoelectrical Sounding data with the actual thickness obtained by drilling. the thickness interpreted from the d.c. Sounding data consistently underestimates the drilled thickness. this is mainly the result of the high vertical conductivity of the Sea ice, in comparison to that in the horizontal direction (Reference MailletMaillet, 1947; Reference Thyssen, Kohnen, Cowan and TimcoThyssen and others, 1974). another factor which may contribute to underestimation of the thickness is where the depth to Sea water below the electrode array is highly variable: this renders the assumed 1-d interpretation model invalid and may result in erroneous interpreted Sea-ice thicknesses.

Fig. 2. Plot of Sea-ice thickness interpreted from East Antarctic d.c. geoelectric Soundings vs actual drilled thickness. Thicknesses were obtained from the d.c. Sounding data via inversion of field data assuming a three-layered model, as described in the main text. The Solid line indicates a 1 : 1 correspondence between interpreted and drilled thicknesses. The d.c. conductivity method consistently underestimates the actual Sea-ice thickness mainly as a result of the pronounced vertical to horizontal conductivity anisotropy in undeformed Sea ice.

In addition to the d.c. Soundings, 13 azimuthal d.c. conductivity measurements were carried out at the field Sites using an offset wenner array (Reference Watson and BarkerWatson and barker, 1999), in order to identify any anisotropy of conductivity in the horizontal plane. the offset wenner array is a variant of the traditional wenner array which enables the effect of azimuthal variations in conductivity resulting from anisotropy to be distinguished from those which arise due to other effects, Such as changes in Sea-ice thickness beneath the electrode array. our offset wenner data did not conclusively identify any anisotropy of conductivity in the horizontal plane. this result was expected, given that azimuthal conductivity variations develop as a result of a preferred horizontal orientation of the c axes of the ice crystals, which can occur when Sea ice forms in a fixed location relative to a prevailing current direction, Such as for landfast ice (Reference Weeks, Ackley and UntersteinerWeeks and ackley, 1986). as noted by Reference Weeks, Ackley and UntersteinerWeeks and ackley (1986), the dynamic nature of the antarctic ice pack means that floes are free to rotate during Sea-ice growth, and that no preferred horizontal c-axis orientation can develop.

Resolution of horizontal and vertical conductivities

Our wenner Sounding and offset wenner azimuthal data Show that the Sea-ice conductivity is horizontally isotropic, but that the vertical conductivity is higher than the horizontal conductivity. this type of electrical Structure is referred to as ‘transversely isotropic’ in the geophysical literature. Reference MailletMaillet (1947) has Shown that a transversely isotropic layer of thickness t a, with conductivities of σ h and σ v in the horizontal and vertical directions respectively, yields an identical d.c. Sounding response to an isotropic layer of thickness t = (σ h/σ v)1/2 t a and conductivity σ m=(σ h σ v)1/2. this equivalence is illustrated in figure 3. the quantity f=(σ h/σ v)1/2 = t/t a is the coefficient of anisotropy. in typical Situations, it is not possible to identify the presence of this type of anisotropy on the basis of Surface geoelectrical Sounding data. in the Sea-ice case, the difference between the interpreted thickness (t) and drilled thickness (t a) reveals the presence of the anisotropy: the ratio t/t a allows us to directly determine f and hence to resolve σ h and σ v independently, Since

Fig. 3. Equivalent transversely isotropic (left) and isotropic (right) 1-D d.c. geoelectrical models. For the anisotropic case, σ v and σ h are the conductivities in the vertical and horizontal directions respectively, and ta is the layer thickness. In the isotropic model, σ m and t are the equivalent upper-layer conductivity and thickness respectively. In both models, the conductivity of the bottom layer (sea water) is isotropic.

(2)

σ v is easily obtained once σ h has been determined as above. Reference Thyssen, Kohnen, Cowan and TimcoThyssen and others (1974) calculated horizontal and vertical conductivities using the Same approach.

We have performed a detailed analysis of five of our geoelectrical Soundings for which drilling revealed the Sea-ice thickness to be essentially constant below the electrode array, and where our assumed 1-d interpretation model was appropriate. at the other field Sites, drilled Sea-ice thickness below the electrode array varied by tens of centimetres to metres, and it was not possible to assume that the underestimation of Sea-ice thickness was Solely due to conductivity anisotropy. for each of the Soundings on level Sea ice, an adequate fit to the observed data was obtained using a three-layered model containing a near-surface conductive layer (Reference Thyssen, Kohnen, Cowan and TimcoThyssen and others, 1974; Reference Buckley, Staines and RobinsonBuckley and others, 1986). the Smallest electrode Spacing used for the field measurements (a = 0.1 m; fig. 1a) was too large to allow the conductivity and thickness of the uppermost layer to be determined (buckley and others, 1986). however, an equivalence analysis of the data indicates that the upper layer is at most a few centimetres thick. all of the equivalent models for the upper conductive layer have an identical longitudinal conductance (conductivity–thickness product); the calculated conductance of the layer, S 1, is given in table 1.

Table 1. Results of in Situ d.c. conductivity measurements on level East Antarctic pack ice. Measurements were made during the period 26 September–20 October 2003. Temperature (T) was measured at the Snow–sea-ice interface; Sea-water temperature was ≈–1.8˚C. t a is the mean drilled thickness of the Sea ice. t and σ m are the Sea-ice thickness and Second-layer conductivity interpreted from d.c. geoelectrical data assuming a layered isotropic model, and S 1 is the conductance (conductivity–thickness product) of the first layer. f is the coefficient of anisotropy, and σ h and σ v are the approximate horizontal and vertical conductivities, calculated as described in the text

The Small thickness of the upper conductive layer means that the Sea-ice thickness interpreted from the geoelectrical Sounding data is essentially equal to the thickness of the Second layer of the model. Subject to this assumption, it is possible to calculate f and hence to determine approximate horizontal and vertical conductivities for the Second model layer, which represents the bulk of the Sea ice. the results of these calculations are Summarized in table 1 and plotted in figure 4.

Fig. 4. Graph of in Situ Sea-ice conductivity and coefficient of anisotropy vs drilled thickness for five d.c. geoelectrical Soundings on level Sea ice in East Antarctica. σ h and σ v are the interpreted Sea-ice conductivities in the horizontal and vertical directions respectively. σ m is the conductivity of the equivalent isotropically conductive layer and is equal to the geometric mean of the horizontal and vertical conductivities.

The conductivities listed in table 1 are consistent with those from other published in Situ data. Reference Thyssen, Kohnen, Cowan and TimcoThyssen and others (1974) measured horizontal and vertical conductivities in the wall of a pit excavated in arctic Sea ice, and obtained σ h = 0.01 Sm–1 and σ v = 0.063 Sm–1. Reference Buckley, Staines and RobinsonBuckley and others (1986) noted the presence of conductivity anisotropy in antarctic Sea ice of thickness 1.35–1.75 m, but did not determine σ h and σ v. based on their published data, we have calculated values of σ h between 0.002 and 0.01 Sm–1 and σ v between 0.01 and 0.034 Sm–1. the in Situ tdr measurements of Reference Morey, Kovacs and CoxMorey and others (1984) measured horizontal conductivities for 1.4m thick first-year arctic Sea ice of between 0.009 Sm–1 near the top of the Sea ice and 0.13 Sm–1 near the more porous base. it Should be noted that conductivities interpreted from d.c. geoelectrical data assign bulk horizontal and vertical conductivities to the Sea ice and cannot be directly compared with the results of Reference Morey, Kovacs and CoxMorey and others (1984) which were made at a much Smaller horizontal Scale (∽1 cm) at a range of depths within the Sea ice.

Comparison with Archie’s law

In the absence of in Situ data, a first approximation to the Sea-ice conductivity can be estimated from core Samples using the empirical Reference ArchieArchie’s law (archie, 1942; Reference Haas, Gerland, Eicken and MillerHaas and others, 1997) which relates the bulk conductivity σ of the Sea ice to its porosity (assumed equivalent to the brine volume, v b) and the conductivity of the brine, σ b:

(3)

Values for the constant m, appropriate to Sea ice, have been determined by thyssen and others (1974) and Reference Morey, Kovacs and CoxMorey and others (1984). thyssen and others (1974) obtained m = 2.2, while morey and others (1984) obtained m = 1.55 near the top of the Sea ice, and 1.75 near the base. the change in the value of m with depth was attributed to a change in Structure of the Sea ice from a vertical c-axis orientation near the top of the ice to a predominantly horizontal c-axis orientation near the base. Reference Haas, Gerland, Eicken and MillerHaas and others (1997), Reference Worby, Griffin, Lytle and MassomWorby and others (1999) and Reference Tateyama, Uto, Shirasawa, Enomoto, Chung, Izumiyama, Sayed and HongTateyama and others (2004) all assumed a value of m = 1.75 for their archie’s law calculations.

Calculation of v b requires Salinity and temperature measurements along the length of the core. the Salinity measurement requires that the core Sections be melted, which has the important consequence of destroying the Structural fabric of the Sea ice, and hence any anisotropy. the archie’s law conductivities can also be affected by loss of brine when the core is extracted from the Sea ice.

We have calculated v b using the method of Reference Cox and WeeksCox and weeks (1983). the brine conductivity σ b for each core interval was calculated from the measured temperature using an equation given by Reference Stogryn and DesargantStogryn and desargant (1985):

(4)

Where T is the temperature in degrees celsius and σ b is in Sm–1.

figures 5 and 6 Show Salinity and temperature measurements for the cores from Sites 1a and 2, as well as bulk conductivities determined using archie’s law for two published values of m. these cores correspond to the Sites with the Smallest (site 1a) and largest (site 2) coefficients of anisotropy determined from the wenner d.c. Soundings. horizontal and vertical conductivities determined from the wenner Sounding data are indicated by the horizontal dashed lines in figures 5b and 6b. the archie’s law conductivities calculated for m = 1.75 are generally a better approximation to σ v than σ h. we have determined bulk core conductivities by Summing the conductances (conductivity times length) of each core Segment and dividing by the total length of the core. assuming m = 1.75, bulk conductivities calculated for cores 1a and 2 are 0.063 and 0.098 Sm–1. for m = 2.2, bulk conductivities are 0.023 Sm–1 (core 1a) and 0.041 Sm–1 (core 2).

Fig. 5. (a) Salinity and temperature profiles measured on a Sea-ice core from Site 1a (Table 1). (b) Bulk Sea-ice conductivity computed from the data in (a) using Archie’s law, for different values of the exponent m. The horizontal dashed lines denote the horizontal (σ h) and vertical (σ v) conductivities determined from interpretation of the in Situ d.c. geoelectrical Sounding data, as described in the text.

Fig. 6. (a) Salinity and temperature profiles measured on a Sea-ice core from Site 2 (Table 1). (b) Bulk Sea-ice conductivity computed from the data in (a) using Archie’s law, for different values of the exponent m. The horizontal dashed lines denote the horizontal (σ h) and vertical (σ v) conductivities determined from interpretation of the in Situ d.c. geoelectrical Sounding data.

The data Shown in table 1 indicate coefficients of anisotropy (f) for east antarctic pack ice of 0.27–0.7. thyssen and others (1974) obtained f = 0.26 for undeformed Sea ice, and found deformed Sea ice to be isotropic (f = 1) as a result of destruction of the preferred fabric in undeformed Sea ice. Reference Buckley, Staines and RobinsonBuckley and others (1986) obtained a coefficient of anisotropy of 0.5 for undeformed first-year antarctic Sea ice.

figure 4 Shows a consistent increase in the coefficient of anisotropy with decreasing Sea-ice thickness, although it is difficult to Say whether this is a real trend, given the limited number of measurements. congelation (columnar) Sea ice could be expected to Show Strong electrical anisotropy in comparison with frazil ice. although no Structural description was available for core 1a, three cores taken from the Same floe, and of almost identical thickness, contained an average of 75% congelation ice, with the uppermost Sections of all cores being composed of frazil and Snow ice. core 2 was composed of 58% congelation ice, with the uppermost Section of the core being composed of frazil and Snow ice. Site 1 Shows the Strongest anisotropy, consistent with its higher content of congelation ice, although the difference in congelation ice contents Seems too minor to explain the large difference in the coefficient of anisotropy at the two Sites. the weaker anisotropy at Site 2 could also be due to the relatively high temperature of the Sea ice at this Site. table 1 Suggests a general positive correlation between Sea-ice temperature and the coefficient of anisotropy. this may imply an increase in the lateral interconnection between the brine cells with increasing temperature.

Discussion

Implications for airborne EM Soundings of Sea-ice thickness

An arbitrarily oriented em transmitter located on or above the Surface of a layered earth will induce current flow only in the horizontal plane. em measurements over undeformed, transversely isotropic Sea ice are therefore Sensitive only to the horizontal conductivity. this has been formally demonstrated for the case of a transversely isotropic half-space by Reference SinhaSinha (1968), but the result generalizes to the case of a multilayered earth (Reference Yin and FraserYin and fraser, 2004).

We have conducted a theoretical investigation of the Sensitivity of the awi hem-bird System to Sea-ice conductivity. the awi hem-bird is a two-frequency System which employs a horizontal coplanar transmitter–receiver geometry. data are acquired at frequencies of 3680 and 112 khz, with transmitter–receiver Separations of 2.77 and 2.05m respectively. at each frequency, the receiver measures the components of the Secondary magnetic field both in-phase and out-of-phase with the primary magnetic field of the transmitter. the Secondary fields are expressed as parts per million (ppm) of the primary field. the altitude of the System above the Sea ice is monitored using a laser altimeter. a full description of the System, and the practical techniques used to recover Sea-ice thickness from the measured data, are given in Reference Haas, Edeler, Schürmann, Lobach and SengpielHaas and others (2002) and Reference Pfaffling, Haas and ReidPfaffling and others (2004).

We have calculated the theoretical awi hem-bird response of a Suite of two-layer 1-d models, in which the Sea-ice conductivity was varied while the Sea-ice thickness and Sea-water conductivity were held constant. the results of one Such modelling exercise are Shown in figure 7. in this case, the two-layered model consisted of a 2 m thick layer of Sea ice over Sea water of conductivity 2.77 Sm–1. the Sea-ice conductivity used for the model calculations ranged from 0.001 to 0.2 Sm–1 in increments of 0.001 Sm–1. random noise was added to the in-phase and out-of-phase components of the model response at each frequency. noise levels were derived from data recorded during the 2003 field program. noise was assumed to be gaussian, with Standard deviations of 6.38 and 5.81 ppm for the in-phase and out-of-phase components at 3680 hz, and 9.24 and 13.42 ppm for the in-phase and out-of-phase components at 112 khz. the noise-contaminated data were then interpreted using a 1-d least-squares inversion program in order to attempt to Simultaneously recover both Sea-ice conductivity and thickness. the inversion process also assumed a two-layered model, in which the Sea-ice thickness and conductivity were allowed to vary, while the Sea-water conductivity was held fixed at 2.77 Sm–1. this is Similar to field Situations, where the Sea-water conductivity is usually known from independent measurements Such as conductivity–temperature–depth Soundings. the Starting model for the inversion assumed a Sea-ice thickness of 4m, and conductivity of 0.1 Sm–1. the Sea-ice conductivity and thickness determined by inversion of the noisy Synthetic data are Shown in figure 7. the true Sea-ice thickness of 2.0 m is recovered quite well, with a maximum error of 12%. however, the Sea-ice conductivity is very poorly constrained, particularly for model conductivities of <0.05 Sm–1. the recovered conductivity exhibits variations of up to two orders of magnitude even for model Sea-ice conductivities above 0.05 Sm–1. figure 7 Shows that, given practical noise levels, the awi hem-bird System is insensitive to the Sea-ice conductivity, even over Sea ice of thickness 2 m.

Fig. 7. Sea-ice conductivity and thickness recovered by inversion of noise-contaminated Synthetic airborne EM data, for a Suite of models in which the Sea-ice thickness and Sea-water conductivities were held fixed at 2m and 2.77 Sm–1 respectively, while the Sea-ice conductivity varied between 0.001 and 0.2 Sm–1. The heavy diagonal line denotes a 1 : 1 relationship between the true and inverted Sea-ice conductivities. The height of the airborne EM System was assumed to be 15 m for all model calculations. A full description of the model parameters and noise levels is given in the text.

We have found the Sensitivity of the airborne em response to Sea-ice conductivity to be even lower over Sea ice <2m thick. our model calculations Suggest that, given the practical noise limits of the awi hem-bird System, the conductance of the Sea ice (σ h times drilled thickness) must be around 0.25–0.4 S before the Sea-ice conductivity can be recovered from the em data with any degree of confidence. this is Supported by further model calculations (not Shown), which indicate that level Sea ice of conductivity 0.05 Sm–1 must be in excess of 5 m thick before the conductivity can be reliably recovered (<10% error) from inversion of practical awi hem-bird data.

The maximum Sea-ice conductance determined from our in Situ measurements was approximately 0.014 S, including the contribution from the conductive near-surface layer (table 1). the horizontal conductivities determined from our in Situ d.c. Sounding data are So low that the effects of Sea-ice conductivity on awi hem-bird data can generally be ignored over level Sea ice in the east antarctic, where the modal Spring thickness is typically around 0.5 m (Reference Worby, Massom, I., Lytle, Heil and JeffriesWorby and others, 1998). this may not be the case during the late Summer, where high temperatures result in high brine volumes and a Significantly increased bulk Sea-ice conductivity (Reference HaasHaas, 1998).

Conclusions

Airborne em measurements over transversely isotropic Sea ice are Sensitive only to the conductivity in the horizontal direction. over level Sea ice, it is possible to determine approximate horizontal and vertical conductivities from d.c. geoelectrical Sounding data. five wenner-array geoelectrical Soundings on level east antarctic pack ice during Spring 2003 have yielded average horizontal and vertical conductivities of 0.017 and 0.073 Sm–1 respectively. these conductivities compare well with those from other in Situ measurements reported in the literature. the widely used form of archie’s law (with exponent m = 1.75) appears to Significantly overestimate the horizontal conductivity of the Sea ice.

High-frequency airborne em data from the awi hem-bird System have been Shown to be insensitive to the electrical conductivity of undeformed east antarctic Sea ice, given the practical noise levels measured during the September 2003 field Surveys. a much earlier field-based Study of airborne em measurements of Sea-ice thickness (Reference Kovacs and HolladayKovacs and holladay, 1990) reached a Similar conclusion for arctic Sea ice. the insensitivity of the System to Sea-ice conductivity was attributed to temperature-related drift affecting both the System electronics and transmitter frequency. although the awi hem-bird contains internal calibration coils which allow in-flight calibration factor adjustments during high-altitude drift correction of the data, further reductions in the noise levels (possibly coupled with an increased maximum transmitter frequency) are required before Sea-ice conductivity measurements will be possible over the relatively thin Sea ice typical of the east antarctic.

Acknowledgements

This research was Supported by australian antarctic Science advisory committee grant 2381. we thank the officers and crew of rsv Aurora Australis for field Support. k. tateyama kindly provided the Sea-ice core temperature data used to produce figures 5 and 6. we thank a. munro for assistance with data acquisition and processing.

References

Archie, G.E. 1942. The electrical resistivity log as an aid in determining Some reservoir characteristics. Trans. Am. Inst. Min. Metal. Petrol. Eng., 146, 54–64.Google Scholar
Becker, A., Echert, D.C. and G.B., White. 1992. Measurement of Sea ice conductivity by electromagnetic induction. In Oceanic Engineering Society (US), ed. Oceans 92: Mastering the Oceans through Technology. Proceedings, 26–29 October 1992, Newport, Rhode Island, Vol. 2. Piscataway, NJ, Institute of Electrical and Electronic Engineers, 748–752.Google Scholar
Buckley, R.G., Staines, M.P. and Robinson, W.H.. 1986. In Situ measurements of the resistivity of Antarctic Sea ice. Cold Reg. Sci. Technol., 12(3), 285–290.Google Scholar
Cox, G.F.N. and Weeks, W.F.. 1983. Equations for determining the gas and brine volumes in Sea-ice Samples. J. Glaciol., 29(102), 306–316.Google Scholar
Golden, K.M. 2001. Brine percolation and the transport properties of Sea ice. Ann. Glaciol., 33, 28–36.Google Scholar
Haas, C. 1998. Evaluation of Ship-based electromagnetic-inductive thickness measurements of Summer Sea-ice in the Bellingshausen and Amundsen Seas, Antarctica. Cold Reg. Sci. Technol., 27(1), 1–16.Google Scholar
Haas, C., Gerland, S., Eicken, H. and Miller, H.. 1997. Comparison of Sea-ice thickness measurements under Summer and winter conditions in the Arctic using a Small electromagnetic induction device. Geophysics, 62(3), 749–757.Google Scholar
Haas, C., Edeler, H., Schürmann, M., Lobach, J. and Sengpiel, J.-P.. 2002. First operation of AWI HEM-bird for Sea-ice thickness Sounding. Proc. Jahrestag. Deut. Geophys. Gesell., 62, 36–58.Google Scholar
Kohnen, H. 1976. On the DC-resistivity of Sea ice. Z. Gletscherkd. Glazialgeol., 11(2), 143–154.Google Scholar
Kovacs, A. 1996. Sea ice. Part 2. Estimating the full-scale tensile, flexural, and compressive Strength of first-year Sea ice. CRREL Rep. 96-11.Google Scholar
Kovacs, A. and Holladay, J.S.. 1990. Sea-ice thickness measurement using a Small airborne electromagnetic Sounding System. Geophysics, 55(10), 1327–1337.Google Scholar
Kovacs, A., Valleau, N.D. and Holladay, J.S.. 1987. Airborne electromagnetic Sounding of Sea ice thickness and Sub-ice bathymetry. Cold Reg. Sci. Technol., 14(3), 289–311.Google Scholar
Kovacs, A., Holladay, J.S. and Bergeron, C.J. Jr., 1995. Footprint/altitude ratio for helicopter electromagnetic Sounding of Sea-ice thickness: comparison of theoretical and field estimates. Geophysics, 60(2), 374–380.Google Scholar
Liu, G. and Becker, A.. 1990. Two-dimensional mapping of Sea-ice keels with airborne electromagnetics. Geophysics, 55(2), 239–248.Google Scholar
Maillet, R. 1947. The fundamental equations of electrical prospecting. Geophysics, 12(4), 529–556.Google Scholar
Morey, R.M., Kovacs, A. and Cox, G.F.N.. 1984. Electromagnetic properties of Sea ice. Cold Reg. Sci. Technol., 9(1), 53–75.Google Scholar
Pfaffling, A., Haas, C. and Reid, J.E.. 2004. Empirical inversion of HEM data for Sea ice thickness mapping. In Extended abstracts: 10th European meeting of environmental and engineering geophysics. Bremerhaven, Alfred Wegener Institute for Polar and Marine Research.Google Scholar
Prinsenberg, S.J., Holladay, S. and Lee, J.. 2002. Measuring ice thickness with EISFlowTM, a fixed-mounted helicopter electromagnetic-laser System. In Chung, J.S., Sayed, M., Kashiwagi, M., Setoguchi, T. and Hong, S.W., eds. Proceedings of the 12th International Offshore and Polar Engineering Conference, 26–31 May 2002, Kitakyushu, Japan, Vol. 1. Cupertino, CA, International Society of Offshore and Polar Engineers, 737–740.Google Scholar
Sinha, A.K. 1968. Electromagnetic fields of an oscillating magnetic dipole over an anisotropic earth. Geophysics, 33(2), 346–353.Google Scholar
Stogryn, A. and Desargant, G.J.. 1985. The dielectric properties of brine in Sea ice at microwave frequencies. IEEE Trans. Antennas Propag., 33(5), 523–532.Google Scholar
Tateyama, K., Uto, S., Shirasawa, K. and Enomoto, H.. 2004. Electromagnetic–inductive measurements for the undeformed and deformed Sea-ice and Snow in the East Antarctic. In Chung, J.S., Izumiyama, K., Sayed, M. and Hong, S.W., eds. Proceedings of the 14th International Offshore and Polar Engineering Conference, 23–28 May 2004, Toulon, France, Vol. I. Cupertino, CA, International Society of Offshore and Polar Engineers, 806–812.Google Scholar
Thyssen, F., Kohnen, H., Cowan, M.V. and Timco, G.W.. 1974. DC resistivity measurements on the Sea ice near Point Inlet, N.W.T. (Baffin Island). Polarforschung, 44(2), 117–126.Google Scholar
Timco, G.W. 1979. An analysis of the in Situ resistivity of Sea ice in terms of its microstructure. J. Glaciol., 22(88), 461–471.Google Scholar
Watson, K.A. and Barker, R.D.. 1999. Differentiating anisotropy and lateral effects using azimuthal resistivity offset Wenner Soundings. Geophysics, 64(3), 739–745.Google Scholar
Weeks, W.F. and Ackley, S.F.. 1986. The growth, Structure, and properties of Sea ice. In Untersteiner, N., ed. The geophysics of Sea ice. New York, Plenum Press, 9–164.Google Scholar
Worby, A.P., Massom, R.A., I., Allison, Lytle, V.I. and Heil, P.. 1998. East Antarctic Sea ice: a review of its Structure, properties and drift. In Jeffries, M.O., ed. Antarctic Sea ice: physical processes, interactions and variability. Washington, DC, American Geophysical Union, 41–67. (Antarctic Research Series 74.)Google Scholar
Worby, A.P., Griffin, P.W., Lytle, V.I. and Massom, R.A.. 1999. On the use of electromagnetic induction Sounding to determine winter and Spring Sea ice thickness in the Antarctic. Cold Reg. Sci. Technol., 29(1), 49–58.Google Scholar
Yin, C. and Fraser, D.C.. 2004. The effect of the electrical anisotropy on the response of helicopter-borne frequency-domain electromagnetic Systems. Geophys. Prospect., 52(5), 399–416.Google Scholar
Figure 0

Fig. 1. (a) Wenner electrode array. Current I is injected through current electrodes C1 and C2. The potential difference, Δv, resulting from the Subsurface current flow is measured between potential electrodes P1 and P2. The ratio Δv/I is the ground resistance R in ohms. (b) Comparison of observed and calculated data for a typical Wenner-array geoelectrical Sounding on Sea ice (site 1b, Table 1). Resistances measured at each electrode Spacing have been transformed to apparent conductivity, as described in the text. The calculated curve is the apparent-conductivity response of the best-fit three-layered model obtained by least-squares inversion of the observed data.

Figure 1

Fig. 2. Plot of Sea-ice thickness interpreted from East Antarctic d.c. geoelectric Soundings vs actual drilled thickness. Thicknesses were obtained from the d.c. Sounding data via inversion of field data assuming a three-layered model, as described in the main text. The Solid line indicates a 1 : 1 correspondence between interpreted and drilled thicknesses. The d.c. conductivity method consistently underestimates the actual Sea-ice thickness mainly as a result of the pronounced vertical to horizontal conductivity anisotropy in undeformed Sea ice.

Figure 2

Fig. 3. Equivalent transversely isotropic (left) and isotropic (right) 1-D d.c. geoelectrical models. For the anisotropic case, σv and σh are the conductivities in the vertical and horizontal directions respectively, and ta is the layer thickness. In the isotropic model, σm and t are the equivalent upper-layer conductivity and thickness respectively. In both models, the conductivity of the bottom layer (sea water) is isotropic.

Figure 3

Table 1. Results of in Situ d.c. conductivity measurements on level East Antarctic pack ice. Measurements were made during the period 26 September–20 October 2003. Temperature (T) was measured at the Snow–sea-ice interface; Sea-water temperature was ≈–1.8˚C. ta is the mean drilled thickness of the Sea ice. t and σm are the Sea-ice thickness and Second-layer conductivity interpreted from d.c. geoelectrical data assuming a layered isotropic model, and S1 is the conductance (conductivity–thickness product) of the first layer. f is the coefficient of anisotropy, and σh and σv are the approximate horizontal and vertical conductivities, calculated as described in the text

Figure 4

Fig. 4. Graph of in Situ Sea-ice conductivity and coefficient of anisotropy vs drilled thickness for five d.c. geoelectrical Soundings on level Sea ice in East Antarctica. σh and σv are the interpreted Sea-ice conductivities in the horizontal and vertical directions respectively. σm is the conductivity of the equivalent isotropically conductive layer and is equal to the geometric mean of the horizontal and vertical conductivities.

Figure 5

Fig. 5. (a) Salinity and temperature profiles measured on a Sea-ice core from Site 1a (Table 1). (b) Bulk Sea-ice conductivity computed from the data in (a) using Archie’s law, for different values of the exponent m. The horizontal dashed lines denote the horizontal (σh) and vertical (σv) conductivities determined from interpretation of the in Situ d.c. geoelectrical Sounding data, as described in the text.

Figure 6

Fig. 6. (a) Salinity and temperature profiles measured on a Sea-ice core from Site 2 (Table 1). (b) Bulk Sea-ice conductivity computed from the data in (a) using Archie’s law, for different values of the exponent m. The horizontal dashed lines denote the horizontal (σh) and vertical (σv) conductivities determined from interpretation of the in Situ d.c. geoelectrical Sounding data.

Figure 7

Fig. 7. Sea-ice conductivity and thickness recovered by inversion of noise-contaminated Synthetic airborne EM data, for a Suite of models in which the Sea-ice thickness and Sea-water conductivities were held fixed at 2m and 2.77 Sm–1 respectively, while the Sea-ice conductivity varied between 0.001 and 0.2 Sm–1. The heavy diagonal line denotes a 1 : 1 relationship between the true and inverted Sea-ice conductivities. The height of the airborne EM System was assumed to be 15 m for all model calculations. A full description of the model parameters and noise levels is given in the text.