# Introduction

This thesis is concerned primarily with the determination of the scattering of ice-coupled waves by imperfections in large expanses of sea ice, with a possible future application being the development of a technique that is able to sense ice thickness remotely. Such ice sheets are found mainly in the Arctic Ocean but they can also be found close to the Antarctic continent, although the sheets there are of a more transient nature as they are regularly broken up by incoming waves from the surrounding Southern Ocean. They may also be dispersed by strong katabatic winds from the continent.

The sea ice of both the Arctic and the Antarctic plays a very important role in the earth's climate. Section 1.1 below presents a general overview of this. As well as describing the specific regions we intend to model, that section also describes how work involving wave interactions with sea ice can be generalized to provide results for certain situations that occur in marine engineering.

The following sections describe and justify our choice of model for polar sea ice (Section 1.2.1), review and discuss mathematical techniques used previously and currently in use in the wave scattering literature (Section 1.2.2), and provide the general outline that this thesis will follow (Section 1.3).

# Arctic and Antarctic Sea Ice

This section begins by giving an overview of the important role that polar sea ice plays in the earth's environment (Section 1.1.1). This includes the reflection of sunlight and the production of bottom water which helps to drive global currents. Section 1.1.2 then proceeds to describe the specific regions we intend to model, while Section 1.1.3 discusses how work involving wave interactions with sea ice can be generalized to provide results for certain situations that occur in marine engineering.

Given the importance of sea ice in the environment, a convenient, remote ice-thickness-sensing technique for monitoring the condition of the polar ice sheets would be a valuable asset in monitoring the overall health of the planet. Techniques currently in use are submarine sonar, most recently coupled with LIDAR---the former is able to detect the under-sea profile of the ice, while the latter uses a laser operated from an aeroplane to detect the ice's upper profile---and EMS, electromagnetic sensing from a specially-fitted helicopter. If ocean wave scattering could be used to deduce average ice thickness, it might be possible to measure it using suitably placed strain gauges or tilt meters, which would eliminate the need for aeroplanes or submarines to be in operation for extended periods of time.

Recall also that by breaking up sea ice, waves affect the way in which the ice interacts with the rest of the climate system; this process should be included in climate models. Consequently, any understanding that wave-ice research can contribute to such models is extremely valuable.

## Environmental Role of Sea Ice

The ice of both the Arctic and Antarctic regions is crucial in the maintenance of the earth's climate. Arctic ice particularly also plays a role in the lives of many human and animal populations---for example, the Sami, the indigenous people of Lapland, use the sea ice for reindeer herding, while polar bears use it as a hunting ground. Animals such as seals and penguins also use Antarctic sea ice as a base from which to search for fish.

Sea ice, the type of ice that this thesis is most concerned with, covers about 7% of the area of the world's ocean, providing a solid barrier between the ocean and atmosphere that hinders the free exchange of heat and moisture between the two. It is also involved in helping to drive currents, most notably the massive thermohaline conveyor which links the Atlantic and Indian with the Pacific and Southern Oceans, and, aided also by the freshwater terrestrial ice shelves of Antarctica and the Greenland ice cap, reflects a large amount of incoming solar radiation.

By reflecting so much solar energy, the ice provides some resistance to increasing global temperatures, and by melting it can absorb still more. However, the melting of the ice reduces the amount of energy that can be reflected, and so the rate of warming increases. In particular, this will allow the amount of heat entering the ocean to increase, contributing to sea level rise due to thermal expansion, a major source of sea level rise. The melting of terrestrial ice also increases sea levels, and decreases the salinity of polar waters.

Sea ice also has a role in driving the world's currents. Its ice formation in the Antarctic causes the production of dense salty water, which then sinks to become Antarctic bottom water (AABW). Similarly North Atlantic deep water (NADW) is partly produced by sea ice formation in the Greenland Sea. Deep water produced then spreads out across the ocean floor (i.e. away from the poles), rising to the surface because of mixing due to roughness in the sea bed and wind, and returning to the poles in wind-driven currents such as the Gulf Stream. Along with the thermal forcing of warm water from the tropics travelling to cooler waters, cooling and sinking, this is the origin of the massive thermohaline conveyor current, which travels between both poles.

Before the existence of plumes was discovered, it was thought that the above process was due entirely to convection, which is impossible in fresh water. (For convection to occur, one needs density to increase as temperature decreases---this is not the case with fresh water, which has a density maximum at $4^\circ$C.) Therefore it was also thought that if the salinity of the polar oceans decreased enough, convection would shut off, with a disastrous effect on currents like the Gulf Stream which bring warm water to North Atlantic countries such as the United Kingdom and the United States of America.

However, deep water can also be formed by smaller scale events such as the refreezing of a newly formed lead. This process produces a plume---a "packet" of dense salty water that sinks to the sea floor. With the extra bottom water produced by these entities, global currents would not be as affected as was first thought by a decrease in polar salinity. Nonetheless, a significant decrease in the formation of sea ice, and thus bottom water, would inevitably affect global currents to some extent.

## Regions to be Modelled

We are interested in modelling the scattering of waves by imperfections in large sea ice sheets, particularly keeping in mind the possibility of remotely estimating ice thickness. In the northern hemisphere, which we are most interested in, such sheets exist mainly in the central Arctic Ocean, due to its sheltered location, although there are also some shore fast ice sheets skirting the Arctic coastline.

The sea ice which forms in places like the Bering Sea or the Greenland Sea does not normally grow into very large sheets due to the presence of high amplitude waves originating in the Pacific and Atlantic Oceans respectively. The ice in both seas is broken up into pack ice, although much of the ice in the Greenland sea is older, thicker ice that has come from the Arctic Ocean. This ice serves to attenuate waves from the Atlantic, and filter out those with shorter periods. As a result, only waves with periods greater than about 6s penetrate into the Arctic Ocean from the Atlantic Hunkins 1962. The size of the Atlantic Ocean also limits the possible lengths of waves originating inside it, with a result that only periods less than about 20s are possible.

The ice in the Bering Sea also plays a similar role. However, the main barriers to waves from the Pacific Ocean reaching the Arctic Ocean are the Aleutian Islands and the Bering Strait. These only allow very long waves (those with periods greater than about 22s) to reach the Arctic from the Pacific (the greater fetches of the Pacific Ocean permits much longer waves to develop within it).

As well as this relative lack of short wave energy, another factor which aids in the growth of large ice sheets in the Arctic is the Beaufort Gyre current. This current circulates around the North Pole, trapping the ice in the centre of the ocean. The result is that the same ice is sometimes able to survive for several years, enabling it to grow to thicknesses of up to 10m, although most of the undeformed sheets are around 1m-3m Wadhams 1995.

The Arctic ice sheets are criss-crossed with sometimes quite large pressure ridges, cracks and open and refrozen leads. Cracks and leads form as adjacent ice sheets respond differently to forces such as those exerted by winds and currents, causing them to drift apart. As a lead opens up, the water exposed to the air freezes; this new, thinner ice then may contribute to the sail and keel of a pressure ridge as winds and currents change and the two ice sheets are pushed back together again. Pressure ridges may also be produced by rafting as ice sheets of similar thicknesses collide and buckle, or in shear as adjacent floes drift past one another.

Antarctic sea ice is not afforded as much protection from incoming waves and is thus generally younger and thinner than Arctic sea ice. (Except in very sheltered bays and inlets, and to the east of the Antarctic peninsula, Antarctic sea ice usually ranges from 0.5m to 1m in thickness.) In winter it is usually formed from roughly circular plates of pancake ice that freeze together to give larger and larger floes the further they are from the ice edge. Like the ice described above that occurs in the Bering Sea, the region near the margins of the ice cover (called the Marginal Ice Zone) is made up of pack ice which reflects and attenuates incoming waves, with the result that the sea ice closer to the continent usually only experiences low amplitude waves of periods in the range 10--20s. The floes are thus able to reach extremely large sizes, even to the point of being able to be described as a continuous ice sheet (although vestiges of their pancake ice origins are left in the form of abundant slightly raised circular edges, giving the ice a "stony field" appearance (Lange et al. 1989). Cracks and leads can be found between neighbouring ice sheets.

The Arctic process of cracks/leads forming then contributing to ridges also occur in the Antarctic. However, much fewer ridges are observed in the Antarctic because the ice is both thinner and less constrained. Substantial ridges are only observed in the ice of the Weddell Sea, where they are formed by ice buckling against the Antarctic Peninsula.

## Engineering Applications of Ice Research

Although this thesis is primarily concerned with providing solutions for situations involving sea ice, some of its results, as well as similar results obtained by other researchers, can also be applied in marine engineering.

A simple application could be to scattering due to a series of floating breakwaters---this could be solved by the method given by Meylan 1993, which describes the effect of a series of ice floes modelled by floating elastic plates.

Another more topical example is provided by the 1000m$\times$121m$\times$3m floating "Mega-Float" platform that was constructed in 2000 in Tokyo Bay, Japan, with the intention of undertaking experiments to determine the suitability of such a structure as an airport (Kashiwagi 2000; Watanabe et al. 2004]]). (Such a structure is also referred to as a Very Large Floating Structure, or VLFS.) With property in Tokyo at such a premium, such an airport would avoid using valuable space on land.

Verification experiments performed included actual take-offs and landings, with pilots reporting experiencing no mechanical differences from landing at land airports (Kanzawa et al. 2002). Landing distances and other factors in safe landing were also found to be satisfactory, as were other considerations such as under-sea aircraft noise, which was found to be at a level that would not be detrimental to the Mega-Float's marine environment.

Given the apparent success of this first trial, the design of a 5km long Mega-Float has been completed, including trial landings with pilots using flight simulators. In order to do this, the hydroelastic response of the airport to waves must first be calculated. By approximating the Mega-Float by a thin elastic plate, work by Meylan and Squire 1994, which models a single ice floe, or alternatively by Tkacheva 2002 could be used to estimate this response. Such models would approximate the Mega-Float as having the side perpendicular to the incident wave of infinite length. Such a solution would be most accurate closer to its centre. Three-dimensional solutions also exist (Namba and Ohkusu 1999, who give a shallow water solution; Wang and Meylan 2003; Peter and Meylan 2004).

Work in this thesis could be used as well---there could easily be the possibility of a hinge joint or a similar type of joint in such a large structure, and strains near such a joint could be estimated by varying the working done in Chapter 5 (also solved in different situations by Squire and Dixon 2000; Evans and Porter 2003a).

Of course, when dealing with specific structures, one would need to refine our approaches somewhat. Inoue et al. 2002 use a thin plate approximation to first obtain the water pressure and instantaneous accelerations at the surface, which are then incorporated into a finite element structural analysis. Clearly techniques that provide the hydrodynamic results quickly would be of great use in such a situation.

# Mathematical Model and Existing Solution Techniques

In Section 2.1, this section begins by presenting and justifying the mathematical model that is generally used for ice sheets, namely the thin elastic plate model. There is now considerable evidence that this model is an accurate one.

In Section 2.2, we then describe the various methods that have been and are being used in the solution of scattering problems involving thin plates floating on fluids. The geometries of the plates have become successively harder, often with several different approaches being attempted for each geometry before the most efficient method is settled upon.

Before proceeding to do this, however, it is worth defining some terms. The first is dispersion. Its name follows from the fact that waves of different periods travel at different speeds. The dispersion relation is the relationship that quantifies this, and forces the wavenumber, $\alpha_0$, to satisfy an equation that depends on the wave's radial frequency, $\omega$, and the properties of the ice (such as its Young's modulus, thickness and density). In general, $\alpha_0$ must be found numerically---an explicit formula in terms of the other quantities rarely exists. As a consequence of the dispersion relation, the phase velocity of the wave ($\omega/\alpha_0$) and also the wavelength ($2\pi/\alpha_0$) change with period ($\tau=2\pi/\omega$) and ice properties. In particular, for constant ice properties the wavelength increases with period, and for constant period it increases with ice thickness.

The last term for the moment, which appears in the title of this thesis, is flexural gravity wave. Gravity waves are waves on a fluid surface that are driven by gravity and a restoring buoyancy force. Flexural gravity waves are waves where this interaction is modified by the inertia and the rigidity of an elastic material floating on the surface.

## An Ice Sheet as a Thin Elastic Plate

An extensive amount of literature exists to support the modelling of an ice-sheet as a linear thin elastic plate, for which the governing equations are presented and discussed in Appendix A (also cf. Mindlin 1951; Fung 1965). Greenhill 1887 was the first to propose modelling a floating ice sheet by a thin elastic beam on a fluid foundation, suggesting a dispersion relation based on the Euler-Bernoulli beam theory. He also wrote a further paper, humourously entitled "Skating on thin ice" Greenhill 1916, in which he elaborated further on the dispersion of waves in ice, as well as touching on other related problems. (Another noteworthy comment that was made in that paper is that ice was the first material for which an experimental value of Young's modulus was obtained.)

This elastic model for ice has since been corroborated experimentally by numerous researchers, with slight improvements to Greenhill's theory being made with the passage of time. For example, during a series of seismic experiments carried out by Ewing et al. 1934, the generation of flexural waves in thin lake/canal ice was observed in addition to the longitudinal and transverse waves being investigated at the time. Extending the theoretical component of Greenhill's work to include the compressibility of the water beneath the ice sheet and to model the ice as a thin plate rather than as a beam (i.e. as a three-dimensional structure rather than as a two-dimensional one), the measured group velocities of the flexural waves corresponded quite well to those predicted by their theory (Ewing and Crary 1934).

A more extensive set of seismic experiments on lake ice to create ice-coupled waves artificially was done by Press et al. 1951, who also observed air-coupled flexural waves travelling at the speed of sound (discussed in further detail by Press and Ewing 1951a). Again, the measured flexural wave dispersion agreed quite well with the thin plate dispersion theory.

Press and Ewing 1951b developed the elastic model of ice further by allowing for the horizontal and vertical displacements to vary arbitrarily in the vertical direction, although they neglect the effect of gravity. In the thin plate model, the horizontal displacements are neglected, while the vertical displacements are taken to be linear in the vertical coordinate. The latter approximations are valid when the thickness of the ice is small in comparison to the wavelength. Indeed, in the large wavelength limit, the dispersion relation of Press and Ewing 1951b explicitly converges to the thin plate relation (allowing for the fact that the gravity term also becomes negligible for large wavelengths).

In the short wave limit, where gravity does have more of an effect, their dispersion relation reproduces the relation for an infinitely thick ice sheet with unattenuated Rayleigh waves travelling along the ice-air surface and with attenuated Rayleigh waves travelling along the ice-water interface. The disadvantage of the theory of Press and Ewing 1951b is that they were only able to solve their dispersion relation in the above long and short wave limits.

Oliver et al. 1954 did a similar seismic investigation to Press et al. 1951 in the Beaufort Sea and the Arctic Ocean amongst Arctic sea ice, and applied the theory of Press and Ewing 1951b in their data analysis. In particular, studies of flexural wave velocities in shore-fast ice off Barter Island showed that for smaller wavelengths dispersion results were in very good agreement with theoretical predictions. Experimental results for longer waves were still quite good, although not as good as for shorter waves. The deviation was attributed to sea-bottom effects due to the relative shallowness of the water (3m). Oliver et al. 1954 also compared measured ice thicknesses in different areas with those deduced from the group velocities of flexural waves or from the frequency at which the air-coupled waves occurred. Probably due to ice inhomogeneity, the flexural wave dispersion was not as reliable a guide to thickness as the frequency of the air-coupled wave. In general, however, the air-coupled frequency still significantly underpredicted the ice thickness (by around 11% to 15%).

Clearly, in seismological studies, the effect of water compression will become more important, especially if the charge is exploded in the liquid. However, in less destructive studies of ice-ocean interaction, such as the one carried out in Notre Dame Bay, Newfoundland by Squire and Allan 1980 involving strain gauges, neglecting that effect is less consequential. Moreover, in the analysis of Press et al. 1951, carried out before the work of Press and Ewing 1951b, the success of the thin plate model in describing the observed dispersion seems to suggest that the thin plate assumptions about the horizontal and vertical motion are not too significant. In any event, Squire and Allan 1980 were able to verify that the thin plate dispersion relation predicted wavelengths that agreed with measured ones to within experimental error. In addition, by considering the problem of a semi-infinite floating ice sheet, they were able to predict the variation of strain in the ice with distance from the ice edge, given the spectrum of the incident waves. Although no such data were available, they were able to estimate the incident wave spectrum by calculating the ratio of the measured strains to the theoretical strains given a uniform spectrum of unit amplitude. If the calculated incident wave spectrum was independent of the position of the strain-meter producing the experimental results, it would serve as a good check of the theory, although it would not confirm it absolutely. Results for three different meters were qualitatively quite similar for three different locations, and differences could possibly have been able to be attributed to the theoretical strains being based on an incomplete solution Wadhams 1973.

Mindlin 1951 also develops a thick plate model which allows for rotational and shear effects inside the plate. For a semi-infinite ice sheet, theoretical predictions (e.g. Fox and Squire 1990; Balmforth and Craster 1999) show that those effects make little difference, although they might be more important near a structure of smaller width such as a pressure ridge.

## Existing Solution Techniques

In summarising the development of mathematical solution techniques, we will concentrate on two particular problems that have been the focus of a lot of previous sea ice research, before showing how these methods could be applied and extended to the treatment of irregularities in sea ice sheets.

### Shore-fast Sea Ice

One particular problem that illustrates the progression of solution techniques for floating plate problems is that of the scattering of incoming ocean waves by shore-fast sea ice, modelled by a semi-infinite sheet of ice. The first attempt at a solution was by Weitz and Keller 1950 who solved a simplified problem by using the Wiener-Hopf technique to obtain an extremely simple expression for the modulus of the reflection coefficient, $R$. This formula is given by

$R=\frac{\alpha_0-k_0}{\alpha_0+k_0}$

where $\alpha_0$ and $k_0$ are the wavenumbers of the propagating waves in the water and in the ice respectively (both are real), and was pointed out by both Shapiro and Simpson 1953 and Keller and Weitz 1953. (Note that although Wadhams 1986 showed that there is an error in the paper by Shapiro and Simpson 1953, their formula for $R$ is still valid as that mistake only pertained to an energy conservation theorem and an expression for the modulus of the transmission coefficient.) Their solution was obtained by ignoring the elastic properties of the ice, and simply modelling the ice as a collection of independent point masses. This model is referred to as the "mass loading model", which, although somewhat inaccurate for ice sheets (for example it predicts that wavelengths should shorten on moving from water into ice, Squire et al. 1995), has proved reasonably successful for wave propagation through frazil or pancake ice (Wadhams and Holt 1991). It effectively reduces the order of the thin plate boundary condition from five to one, eliminating the need to apply additional conditions at the ice-edge.

Evans and Davies 1968 subsequently used the Wiener-Hopf method (Noble 1958; Roos 1969) to solve the full problem. At the time the infinite products that formed a large part of the solution proved to be too difficult to compute and so they were not able to present any results, except for certain limiting cases (such as the shallow water limit). Incomplete mode-matching approaches were then taken (e.g. Hendrickson and Webb 1963, Wadhams 1973, Squire 1978, Squire 1984, Wadhams 1986), before Fox and Squire 1990 computed the solution for the full set of eigenfunctions using a conjugate gradient technique. Fox and Squire 1994 used the same method to complete further studies on this problem, investigating the strain in the ice, and also the effect of shore fast ice on an incoming directional wave spectrum of specified structure.

Gol'dshtein and Marchenko 1989 presented a later Wiener-Hopf solution for infinite depth, but again only certain limiting situations were discussed. However, in the late 1990's and early $21^st$ century improvements in computing power enabled other authors to compute results using the Wiener-Hopf solution of Evans and Davies 1968. Balmforth and Craster 1999 turned the required infinite products into integrals which were evaluated by quadrature, while Chung and Fox 2002 showed that the products themselves could be evaluated directly with relatively little effort. Ironically Tkacheva 2001 finally showed that if the inertia term in the thin plate equation was neglected for normally incident waves (as it can for most wavelengths), then $|R|$ could be calculated by simply using the correct value for $k_0$, the wavenumber in the ice, in the formula (1.1) given by Keller and Weitz 1953.

Other authors have returned to the mode-matching approach of Fox and Squire 1990. Sahoo et al. 2001 defined an inner product enabling the solution to be found by using $N$ eigenfunctions and inverting an $N\times N$ matrix, while Linton and Chung 2003 effectively showed that the equations Sahoo et al. 2001 had set up could be solved analytically using residue calculus. In the process they also demonstrated its equivalence to the Wiener-Hopf solution, and confirmed the formula of Tkacheva 2001. Earlier Chakrabarti 2000 had used an infinite depth mode-matching scheme to set up a singular integral equation, the equivalent problem to a residue calculus problem when the eigenvalue spectrum is continuous. This singular integral equation was solved by transforming it into a Hilbert problem Roos 1969, a generalization of the Wiener-Hopf problem (although in practice both are solved in exactly the same way).

### An Ice Strip of Finite Width in Open Water

A second problem that has generated a lot of the mathematics used herein is that of the scattering of water waves by an ice strip of finite width surrounded by open water. By restricting the incoming waves to normal incidence, it was initially intended as a two-dimensional model for a single ice floe such as might occur in the MIZ (Meylan 1993; Meylan and Squire 1994). Although modelling of ice floes has since moved away from this approach towards a more three-dimensional one (e.g. Meylan et al. 1997; [[Peter and Meylan 2004),the results and techniques are still quite applicable to situations where the length of the floating object is large in comparison to other characteristic lengths of the problem. Examples include the Mega-Float discussed in Section 1.3 (Kanzawa et al. 2002).

The solution in the shallow water limit is relatively straightforward (Stoker 1957) and simply involves the inversion of an $8\times8$ matrix. In his Ph.D. thesis, Meylan 1993 solved the problem for both finite and infinite depth by using Green's theorem and the two-dimensional open water Green's function to set up an integral equation over the floe (the infinite depth solution was also published in Meylan and Squire 1994. It is also possible to extend the mode-matching technique of Sahoo et al. 2001 to treat this problem, and the author is aware of several researchers who have done this, although such a solution is still unpublished. This technique could be improved by using the residue calculus method in the same way that is it is used in Appendix F of this thesis, which treats the problem of the scattering by three adjacent ice sheets of different thicknesses. With minor adjustments, either or both of the outer ice sheets may be replaced with open water (cf. Section 1.3 of Scattering by a Single Lead; published papers on the residue calculus method include those by Linton and Chung 2003, Chung and Linton 2005).

Tkacheva 2002 solved a pair of coupled Wiener-Hopf equations to produce another approach to this problem, using a similar method to the one that is used in Scattering by a Single Lead of this thesis, and Meylan 2002 solved the time dependent problem by using a spectral approach.

### Irregularities in a Continuous Ice Sheet

The first type of irregularity in a continuous ice sheet considered was an abrupt change in ice properties, i.e. the scattering by a joint/edge separating two semi-infinite ice sheets, each with different properties. Barrett and Squire 1996 solved this problem numerically by extending the method of Fox and Squire 1990 used in dealing with the shore fast ice problem, while Chung and Fox 2002b solved it by extending their own Wiener-Hopf solution (Chung and Fox 2002a).

A special case of the above problem is that of the scattering by a crack, i.e. by a free edge separating two identical semi-infinite ice sheets. Squire and Dixon 2000 extended the work of Meylan 1993 to derive the infinite depth Green's function for an ice sheet; this was then used to obtain an analytical expression for the reflection and transmission coefficients for a normally incident wave upon a crack (in ice floating on water of infinite depth). Williams and Squire 2002 later used the same method to derive similar expressions for the same coefficients when the incident wave arrived at an oblique angle. The finite depth problem was solved by Evans and Porter 2003a who used two different methods to calculate the finite depth solution---the first involved using a mode-matching technique and the other used the finite depth Green's function for an ice sheet.

Chung and Linton 2005 extended this mode-matching approach to treat a crack of finite width (also known as an open lead), using the residue calculus technique to improve their solution. Williams and Squire 2004b also presented results using a mode-matching technique for an open lead (although without detailing their method). Williams and Squire 2004b could also have used the variable thickness method outlined in the same paper, and first presented in their Williams and Squire 2004a paper, by simply setting the thickness of the variable region to be constantly zero.

Marchenko 1997 also solved the problem of cracks in ice for infinite depth, by taking a Fourier transform directly (as opposed to finding the transform of the Green's function), allowing the ice to have arbitrarily many (parallel) cracks. The scattering by many cracks was also treated by Dixon and Squire 2001, who considered the effect of randomising the separation of cracks, while Evans and Porter 2003a applied Floquet's theorem to investigate the scattering by an infinite series of identically spaced parallel cracks.

An additional problem treated by Marchenko 1997 was the problem of scattering by a pressure ridge, in which he modelled the ridge as an elastic beam joined to the surrounding ice cover by an elastic hinge joint. Approaching the problem in this manner enables the solution to be found in much the same way as the crack solution was found. Williams and Squire 2004a treated the ridge as a thin elastic plate of finite width, but with a varying thickness, and used Green's theorem to reduce the problem to an integral equation over the ridge's width (also presented in Scattering by a Single Ridge of this thesis). Dixon and Squire 2001 had previously solved this problem in the special case where the thickness of the ridge was constant by replacing the open water Green's function in the method of Meylan 1993 (also reported in Meylan and Squire 1994) with the thin plate Green's function. Both Williams and Squire 2004a and Dixon and Squire 2001 neglected the submergence of the ridge/berg, although this effect could be allowed for in the latter paper in the same way that Meylan 1993 solved the problem of a submerged floe.

A different approach to the ridge problem was taken by Porter and Porter 2004, who used a variational approach and a mild slope approximation to reduce it to a sixth order ordinary differential equation. Their method allowed a pressure ridge of arbitrary thickness and keel depth to be considered by assuming that a negligible amount of evanescent wave action was produced as the incident wave was scattered. In addition, an arbitrary sea floor topography could be treated simultaneously. (Evanescent waves are waves with non-real wavenumbers whose amplitudes decay exponentially as the distance from their point of generation increases. Their wavenumbers are generally pure imaginary, although two complex wavenumbers are possible in ice; the waves corresponding to those wavenumbers are sometimes distinguished from the others by being called damped traveling waves, or simply complex modes. Traveling or propagating waves have purely real wavenumbers.)

A surprising implication of the work of Porter and Porter 2004 is that for a large enough water depth, removing the submergence would have negligible impact on the scattering by a ridge as long as its overall rigidity is preserved. This means that, assuming that the gradient of the ridge's rigidity is not too large, the integral equation method of Williams and Squire 2004a will produce accurate results for such depths, as long as the total ice thickness is increased to allow for the keel.

Williams and Squire 2004a presented results for two ridges as well, comparing the exact solutions to a wide-spacing approximation which becomes more accurate as the separation between the ridges is increased. This approximation assumes that the evanescent waves generated at each ridge have decayed to a negligible amplitude by the time they have reached the next. For this reason it is called the No Evanescent Waves (NEW) approximation, although it is actually equivalent to the wide spacing approximation used, for example, by Newman 1965 or Evans 1990. Equivalent results were presented by Meylan 1993 for two floes.

The NEW approximation was extended by Williams and Squire 2004a to allow for arbitrarily many ridges. In addition, when the effect of the ridge separations was averaged out, the result was extremely well approximated by the so-called serial approximation. The serial approximation is obtained by simply multiplying the transmission coefficients for each ridge to obtain an overall transmission coefficient. It was introduced to allow the general scattering properties of a given ice field to be studied without the complicated interference produced in the presence of many ridges.

Williams and Squire 2004b extended the work of Williams and Squire 2004a to allow the ice field to contain leads and cracks, and the serial approximation also proved to describe the average scattering behaviour of fields containing those features accurately.

This result was taken by Williams 2004 to suggest that the potential for using flexural waves in remote ice-thickness-sensing could initially be considered by simply using the serial approximation. In that paper, reflection curves for a given set of ridges were compared for different ice thicknesses and were found to be fairly different over a relatively wide range of periods (from about 2s to 11s). Thus there seemed to be some hope that scattering results for different thicknesses would be sufficiently different to use those results in determining the ice thickness.