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
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.
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
σ 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.
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.
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:
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):
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).
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.
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.