Introduction

Environmental conditions at the time of ice formation largely determine its structure. Ice crystals are uniaxial and the optic axis corresponds to the c-axis. Random c-axis orientation is commonly observed near the top surface of newly formed sea ice in the Arctic. Once a cover has formed, the ice structure is characterized by long vertical columns that extend downward in the growth direction of the ice sheet, a result of quiescent, unidirectional growth. Under these conditions, a selective growth process occurs and the c-axes of the crystals become primarily oriented in the horizontal plane of the ice sheet (Reference Weeks and AckleyWeeks and Ackley, 1982). In the presence of a predominant current direction, strong c-axis alignment develops in the direction of the current, with generally decreasing scatter as ice thickness increases (Reference Weeks and GowWeeks and Gow, 1978), and this ice structure causes anisotropic material behavior in all directions. Uniaxial compression data on first-year sea ice (Reference WangWang, 1979; Reference Richter-Menge, Cox and PerronRichter-Menge and others, 1987) have indicated a strong dependence of peak compressive strength on c-axis alignment and on the angle between an applied load and the dominant c-axis direction. Therefore, to interpret data from mechanical property tests, we must define the relative orientation and alignment of the ice fabric.

The techniques used in the analysis of ice fabrics were originally developed in structural petrology (see, for example, Reference Knopf and IngersonKnopf and Ingerson, 1938; Reference FairbairnFairbairn, 1949; Reference Turner and WeissTurner and Weiss, 1963). Crystal-orientation measurements usually involve optical measurements of the c-axis orientations. Reference LangwayLangway (1958) described techniques for obtaining ice-crystal c-axis orientation data using a Rigsby universal stage. One orientation measurement is made for each ice crystal in a sample, and these data are plotted on a Schmidt equal-area net that represents the surface of a hemisphere of unit radius. The points on the net comprise an orientation diagram that depicts the relative spatial concentration of the data. An orientation diagram may indicate a random or a patterned structure, depending on the dominant features of the diagram. A random fabric refers to a homogeneous distribution of the data and represents an isotropic material. In this configuration, there is an equal probability of finding points in equal-area elements anywhere on the net. In contrast, the most significant feature of an anisotropic material is the preferred orientation indicated by the grouping of points on the net. A “girdle” corresponds to data that are distributed along a great circle of the net, indicating a preferred planar orientation of the c-axes. The pole of this great circle is termed the girdle axis. An area of highly concentrated data points (point maximum) indicates a linear-preferred orientation of the crystals in the fabric. The statistical significance of the orientation diagram increases if the main features are reproducible in different comparable samples from the same homogeneous body.

Reference PearsonPearson (1901) used statistical arguments to develop the equation of a line or plane that provided the closest fit to points in space when all variables contain error. The solution depends on knowledge of the means, standard deviation and correlations of the variables. A significant result was that the plane of best fit contains the line of best fit. Reference WatsonWatson (1966) presented a matrix of sums of direction cosines of vectors representing crystal orientation in a Cartesian coordinate system. The reasoning presented was that the greatest moment of inertia of the points would be about the eigenvector corresponding to the minimum eigenvalue of the matrix. Reference MardiaMardia (1972) used the same reasoning to obtain this matrix and interpret the results. Both Watson and Mardia proposed distributions for the data on the sphere, and developed statistical analyses based on these assumptions. Reference Diggle and FisherDiggle and Fisher (1985) described a program that computed these eigenvalues and eigenvectors and quantitatively contoured spherical data. The analysis of ice-fabric diagrams has been largely visual, frequently based on approximate data-concentration contours drawn on the net. However, Reference Herron and LangwayHerron and Langway (1982) applied the eigenvalue/eigenvector method of Reference MardiaMardia (1972) to study various fabric types, including small-circle girdles and multi-maxima patterns. The results were interpreted qualitatively and it was not clear whether any fabrics could be quantitatively assessed by this method.

In this paper, we will seek a plane passing through the origin that minimizes the sum of the squared normal distances from the data, and obtain the dominant c-axis orientation in this plane. A detailed derivation from simple geometric arguments is developed, yielding least-squares equations that minimize the orthogonal distance between the data and the best line and plane. The eigenvalue problem that results is the same as that obtained by Pearson, Watson and Mardia. We identify an implicit assumption in this method of equal-measurement uncertainty in each coordinate at all points on the sphere. Normalized eigenvalues provide quantitative measures of physical distance of the data from the plane and line, specifying the directional characteristics of the c-axes of crystals in a sample. Mean angular measures of variability are also developed. This analysis clearly indicates the fabrics that are well-described by the eigenvalue/eigenvector method, and provides a framework for developing related quantitative methods for other fabric types. We demonstrate the capabilities of the analysis on data sets representing samples of first-year sea ice. The results are viewed on special Schmidt nets that represent data hemispheres defined by the best plane and the predominant basal plane, in addition to the xy-plane.

Measurement Errors and Data Fitting

Measurement of the c-axis orientation of an ice crystal with a Rigsby universal stage provides an azimuth angle, the direction of inclination as right or left and the type of measurement as polar or equatorial. From these data, we obtain a pair of angles (θ,φ) that define the orientation (Reference Ferrick and ClaffeyFerrick and Claffey, 1992). The line representing each crystal in the sample is plotted through the origin of a unit sphere. Each of these lines intersects the surface of the sphere at one point in each hemisphere and the lower hemisphere is traditionally studied. Our analysis represents the c-axis of each ice crystal as a three-dimensional unit vector from the origin in the half-space below the xy-plane, yielding an array of points on the surface of a hemisphere of unit radius p. The three-dimensional Cartesian coordinates of each point on this hemisphere (Fig. 1) are obtained directly from the spherical cordinates (ρ, θ, φ) as

(1)

Fig. 1. Sketch of unit vectors 1, 2 ..., m,... Ν representing the c-axis orientations of ice crystals in a sample. The hemisphere, Z ≤ 0, is shown by convention.

The relative uncertainties in x, y and z caused by measurement error must be understood to identify an appropriate method for fitting the ice-fabric data. Reference LangwayLangway (1958) listed several sources of error in the measurement of c-axis orientation, and indicated that the errors in azimuth and inclination should each be less than 5 °. We observe that the same statement holds for errors in θ and φ. Taking these errors as random and each of magnitude β, we have from Equation (1) that

(2)

Expanding Equation (2) and grouping terms according to the error, we obtain

(3)

where Α, Β, C, D and Ε are composed of sines and cosines of φ and θ. The relationships given in Equation (1) are contained in Equation (3), but modified by errors of ß. An error in θ only affects x and y. The resulting uncertainty in each coordinate is periodic, depending on θ, and out of phase with the other. The amplitude of these errors approaches zero near the pole and a maximum at the equator. Errors in φ affect all three Cartesian coordinates. The uncertainty in z is larger near the equator and smaller near the pole than the larger of x and y. Uncertainties in x and y resulting from errors in φ again vary individually with θ, displaying maximum amplitudes near the pole and approaching zero near the equator.

Classical least-squares methods require a dependent variable but the coordinates of each c-axis are all independent. Reference ReedReed (1989) presented a method for fitting a line to points in the plane when both coordinates of these points are independent and uncertain due to measurement error. The method allows the errors in the fit to be weighted according to the relative uncertainties in the measurement of x and y. The measurement errors discussed above could be considered by extending this method to three dimensions. However, we note for small β that cos β ≈ cos² β ≈ 1 and sin β ≈ 0, indicating that errors in x,y and z resulting from measurement errors are relatively small. As a first approximation, we will choose equal weighting in each direction at every point and the best-fit line and plane will minimize the perpendicular distances from the data.

Determination of the Best Plane by Orthogonal Least Squares

The unit vectors representing the c-axis orientation of each crystal in a sample have a common point at the origin of the unit sphere. The problem we consider here is to find the plane of best fit to fabric data that contains the origin, and to provide quantitative measures of the quality of the fit. The form of the equation of a plane through the origin is

(4)

We choose the function F(x, y,z; c ₁, c ₂, c ₃) with this same form and depending linearly on parameters c ₁, c ₂, c ₃ as

(5)

where φ₁ = x, φ₂ = y, φ₃ = z are a specified set of mutually orthogonal functions and the c _i are unknowns to be determined. The pole of a plane is the point of intersection Ρ with the hemisphere of a line through the origin normal to the plane. The unit normal η to the plane of best fit has the form

(6)

where i ₁, i ₂ and i ₃ are unit vectors in the x, y and z directions, respectively, and repeated indices indicate summation.

A unit vector in three-dimensional space represents each crystal, and the total number of crystals in a sample Ν used to fit the plane should be much larger than 3. The unit vector representing the mth crystal (Fig. 2) is

(7)

This vector intersects the unit sphere at (x_m, y_m, z_m) and we can evaluate F(x_m, y_m, z_m; c₁, c₂, c₃) = F_m with Equation (5) as

(8)

Fig. 2. Sketch of the unit vector u_m representing the mth crystal in an ice sample and its projections on to the best-fit plane or the basal plane, and on to η or c the unit normal vector to the corresponding plane. The angles between u_m and the best (plane, line) are (α_m, α_m), respectively.

The vector is the projection of u_m on to n, representing the normal vector from the plane to the point (x_m, y_m, Z_m)

(9)

where (u _m. n) is a scalar product between unit vectors. Throughout this development, the notation (~) over a vector indicates that it does not have unit length. The normal distance d _m from point m to the plane is

(10)

and we observe that

(11)

The sign of d _m distinguishes distances on opposite sides of the plane.

The orthogonal projection of u_m in the best-fit plane is represented by the vector , and

(12)

This vector joins the origin and the point x_m, y_m, z_m, the orthogonal projection of (x_m, y_m, z_m) in the plane. The vector is unique and the best approximation of u_m in the plane. An equivalent statement is that the closest point in the plane to (x_m, y_m, z_m) is (x_m, y_m, z_m). From Equation (12), we can determine as

(13)

where represent elements of TV-dimensional vectors respectively.

We want to choose the (¾ values that specify the plane through the origin with normal distances | d_m between the points representing the JV crystals in a sample and the plane that are as small as possible. The Nvector d = (d ₁, d ₂, …, d_N ) ^T represents the individual normal distances from the data to the plane, where Τ indicates the transpose. Minimum normal distance to the plane is equivalent to maximum length of , the projection of u_m on to the plane. We will seek a least-squares fit and define £ as a function of the unknown coefficients

(14)

Because we consider data projected on to a plane, the hemisphere in which the data appears is arbitrary. The sum of the squares of the lengths of will be maximized and the sum of the squared distances from the plane d _m will be minimized where the gradient of Ε vanishes.

(15)

with indicating coefficients of the plane of best fit.

We now perform the differentiation indicated in Equation (15) with respect to c_j, where each choice of j yields a scalar equation of the form

(16)

From Equation (13) we obtain

(17)

] where δ_ij is the Kronecker delta. We observe from Equation (10) that d_m = c_kϕ_km , and

(18)

Inserting Equation (13), (17) and (18) into Equation (16), we obtain the normal equations as Equation (2)

or expanding and re-arranging as

(19)

Finally, because n is a unit vector and i is a summation index, and Equation (19) becomes

(20)

In Equation (20) we observe that a constant, yielding an eigenvalue problem,

(21)

where each term a_ij of the 3×3 symmetric matrix A is

(22)

The eigenvalue λ is the sum of the squared normal distances of the data from the corresponding plane. Equation (21) in homogeneous form indicates that non-trivial solutions exist if and only if

(23)

The determinant given in Equation (23) yields a cubic equation called the characteristic polynomial

(24)

Following Reference BeyerBeyer (1987), we obtain the solution as

(25)

where

As p, q and r are real, the eigenvalues will also be real. Because A is a real symmetric matrix, it is similar to a diagonal matrix Β composed of the eigenvalues of A, and therefore the eigenvalues of A are real. Similar matrices have the same trace and the same determinant. The trace of a matrix is the sum of the elements on the principal diagonal. The elements a₁₁, a₂₂ and a₃₃ of matrix A represent the sum of the squares of the distances between the data and the three planes defined by the coordinate axes. As each crystal is represented by a unit vector,

(26)

indicating that the sum of the eigenvalues is N, the total number of crystals in the sample. The determinants of A and Β are the product of the eigenvalues. The eigenvalues of A are non-negative if A is positive semi-definite; that is, x^T Ax ≥ 0 for any x. Since the eigenvalues represent sums of squared distances, λ_i ≥ 0 and matrix A is positive semi-definite. We will designate the eigenvalues in increasing order according to magnitude as λ ₁ ≤ λ ₂ ≤ λ ₃, and define normalized eigenvalues as

(27)

The normalized eigenvalues give the mean squared normal distance between the points on the unit sphere defined by the unit vectors and the plane normal to the corresponding eigenvector. These eigenvalues provide a measure of the fit that is equivalent to the variance in classical dependent variable least-squares methods.

Eigenvectors of a real symmetric matrix corresponding to different eigenvalues are orthogonal and, because the eigenvalues are real, the eigenvectors can be taken to be real. The vector ν is an eigenvector for A belonging to the eigenvalue λ if

(28)

The lengths of these eigenvectors are arbitrary, and we normalize them to unit length to obtain an orthonormal basis in three-dimensional space. Each eigenvector represents the unit normal to a plane, and the corresponding eigenvalue gives the sum of the squared normal deviations of the data from that plane. The minimum eigenvalue defines the plane of best least-squares fit to the data, and the higher eigenvalues are associated with the remaining mutually orthogonal planes through the origin. With the origin fixed, the eigenvector basis represents a coordinate system that is rotated relative to the coordinate axes.

The eigenvectors written in columns form the matrix P. The elements of Ρ are the direction cosines between each eigenvector and the coordinate axes. The angle α_ij between the eigenvector v_j, and the i_i is

(29)

Because the columns of Ρ are orthonormal, Ρ is an orthogonal matrix, P ^–1 = P^T is also orthogonal, and det P = ±1. Matrices A and Β are related through Ρ as

(30)

representing a singular value decomposition of A. The diagonal elements of Β are the singular values c/j as well as the eigenvalues of A. With singular values ordered by their magnitude in the same way as the eigenvalues, the condition number of matrix A of full rank is

(31)

Condition numbers ≫ 1 indicate that A is nearly singular.

The sum of the squared normal distances between the data and the best plane, given by λ_x, provides a measure of the planar structure of an ice sample. Values of λ₁ or λ′₁ approaching zero indicate an increasingly planar ice structure. If λ₁ = 0, then cond(A) = ⋡ and the data are perfectly planar. A visual representation of the error is obtained from the angle a_m between the c-axis of an individual crystal and its projection in the best plane (Fig. 2)

(32)

The absolute value in Equation (32) is needed if the angle between u_m and n is greater than π/2. The average angular deviation α between the data and the best plane is a parameter we term the planar spread that can be readily determined and understood:

(33)

A small planar spread indicates a small mean angle between the c-axes of crystals in a sample and the plane.

A development that parallels the above to find the best plane by classical dependent variable least squares has been given by Reference Ferrick and ClaffeyFerrick and Claffey (1992). This standard method seeks a minimum error in the dependent variable and is more commonly available than the orthogonal method, but the results are sensitive to the choice of dependent variable. The normal distances between the data and the plane taking z as the dependent variable can be determined as

(34)

The sum of the squared normal distances given in Equation (34) is the same measure as the minimum eigenvalue. Similar expressions result when x or y is taken as the dependent variable.

Alignment of The C-Axes

An area of highly concentrated c-axis data on the Schmidt net indicates a linear preferred orientation of the ice crystals of a sample, and suggests the need to determine the predominant optic-axis orientation of the fabric. We will locate this linear orientation by following a development parallel to that used to determine the best plane. The unit vector c represents the unknown preferred c-axis orientation of the crystals in a sample, and is expressed as

(35)

where the primes distinguish these coefficients from those of the unit vector n given in Equation (6). The plane through the origin that is normal to c represents the predominant basal plane orientation of the crystals in the sample and is described by

(36)

As before, the unit vector u_m represents the mth crystal c-axis orientation and intersects the unit hemisphere at (x_m, y_m, z_m). In order to take advantage of the detail given in the previous development, the vector projections of u_m on to the plane and its normal vector are again respectively. Then, with c replacing n and a different plane, Figure 2 represents our present condition.

When searching for the plane of best fit, we sought to minimize the squared normal distance. However, the closest representation of u_m by the unit vector c requires that we maximize the sum of the squares of the lengths of . The predominant basal plane orientation is then the plane of maximum squared normal error with the data. Each step in the previous development applies except that c_i is replaced by . Note that the matrix A given in Equation (22) is unchanged because it depends only on the coefficients of the intersections of the individual unit vectors u_m with the unit hemisphere. Therefore, we are solving the same eigenvalue problem as before. The eigenvectors obtained in Equation (28) are orthogonal, and c is contained in the best plane. Reference KassKass (1989) has shown for the general case that spaces of closest fit are nested. Every p-dimensional sub-space of closest fit lies in one of dimension p + 1. The third eigenvector s together with c form the plane of best fit to the c-axes, and s and n form the predominant basal plane (Fig. 3). The direction cosines of the angles between the preferred linear orientation and the coordinate axes are given in matrix P, and the angles can be obtained from Equation (29).

Fig. 3. Unit vector c representing the linear prefmed c-axis orientations of the ice fabric, and the long axis of the columnar crystals of an ice sample are represented by n. The unit vectors n and s provide the predominant orientation of the basal plane of the sample. A unit vector σ represents the direction of load application on the sample, σ is the projection of this vector on to the basal plane, σ_c is the angle between the applied load and c, and σ_z is the angle between the load and n.

The mean squared normal distance between the data and the vector c is . Again, using angles to visualize error, the angle δ_m between an individual unit vector and its projection on to c (Fig.2), always ≤90 °, is obtained as

(37)

If 8_m obtained with Equation (37) is > 90 °, δ_m is replaced by its supplement. The average angle δ between the data and c is

(38)

The parameter δ is termed the linear spread, and a small δ indicates that the angles between the preferred orientation and the data are also small.

As a group, the eigenvalues provide measures of the structure of the ice sample. An eigenvalue of zero occurs if the data are two-dimensional, and a pair of zero eigenvalues represent one-dimensional data. More generally, small together with large occurs with aligned data, and small with significantly larger and indicates planar data. If, in addition to small the crystal orientations distribute symmetrically about a great circle of the sphere. With every set of orthonormal vectors will serve as the eigenvectors. The data are maximally dispersed with respect to lines and planes, and the material is isotropic.

The eigenvalue/eigenvector analysis has been applied to fabrics displaying other patterns, including multi-maxima and small-circle girdles. Fitting these data with a best line or plane yields results that can only be qualitatively interpreted. However, when a fabric displays several point maxima, it may be possible to objectively subdivide the data. The best line for each sub-set would quantify the orientation and alignment of the corresponding point maximum. The intersection of cones with orientation diagrams was drawn by Reference Kohnen and GowKohnen and Gow (1979) to bound small-circle girdle fabrics. A small-circle girdle fabric could be assessed quantitatively by finding a circular cone with its axis defined by the best line and its apex at the origin. The surface of the cone would be located to minimize the sum of the squared distances with the data, and variability measures could again be defined in terms of distance and angle.

If the c-axes in a given sample are sufficiently aligned, the mechanical properties of the ice will be affected. If σ is a unit vector in the direction of an applied force, the angle σ_c between the load and the dominant c-axis direction is

(39)

and the complement of σ_c is the angle between the load and the basal plane (Fig.3). For a columnar ice sample, n gives the predominant direction of crystal elongation and growth. The angle σ_z between the load and the vector n is obtained as

(40)

Schmidt-Net Representations

The details necessary for automated plotting of the data, the plane of best fit and its pole, and vector of best fit on the Schmidt net in the xy-plane have been given by Reference Ferrick and ClaffeyFerrick and Claffey (1992). The intercept of c must fall on the great circle of best planar fit.

During field ice coring and thin-section preparation, a sample intended as horizontal may deviate by an angle of several degrees. Sea ice with horizontal c-axis alignment will then be represented by a best plane at this angle. The linear dimensions on the periphery of the Schmidt net are distorted, and it is difficult to judge normal distances between the data in this region and the trace of the best plane. Points that appear near the perimeter and directly across the net from each other represent crystals with close planar alignment. For these reasons, viewing the data on Schmidt nets drawn on alternative planes is frequently advantageous. Viewing the data on a Schmidt net drawn in the best-fit plane with n vertical eliminates sample-preparation error for horizontally aligned sea ice. The great circle of the best plane falls on the perimeter of the net and the pole of this plane appears at the origin. With c vertical, the plane of the Schmidt net is the predominant basal plane of the sample, and the data are transformed from the perimeter to the middle of the net. The great circle of the plane of best fit must then pass through the origin. Accurate visual assessments of the linear and planar preferred orientations of the sample are possible on this net because of minimal distortion of linear distance near the origin. With these planes for mapping, the importance of net distortion, sample preparation and measurement accuracy are minimized by providing optimal views of the crystal fabric.

Unit vectors in the Cartesian coordinate directions were used in Equation (7) to obtain a unit vector representing each crystal in a sample. These unit vectors are related to an orthogonal coordinate system of eigenvectors by the matrix Ρ as

(41)

This relationship is depicted in Figure 4 with eigenvectors n, c and s as the (primed) unit eigenvectors. The relationship between a general vector in the Cartesian system and its transform in the eigenvector system is (Reference HildebrandHildebrand, 1965)

(42)

and Ρ is termed the transformation matrix. Because Ρ is orthogonal, transformations using Ρ are orthogonal. Orthogonal transformations maintain length and preserve angle, and can be interpreted as a combination of rotations and reflections. The order of the eigenvector placement in Ρ determines the transformation, and n or c is made vertical.

Fig. 4. Relationship between the Cartesian and eigenvector coordinate systems for visualizing alternate Schmidt-net representations of ice-fabric data.

Applications

We analyzed the ice fabric of many samples of first-year sea ice taken from the Beaufort Sea. Both orthogonal and dependent-variable least-squares analyses were used and compared. These cases are presented in Figure 5 in an arbitrary order from highest to lowest mean normal error of the best plane obtained from the orthogonal analysis. In all cases, this approach provided a lower bound to the normal error of the dependent-variable methods. The individual samples are designated by the dependent-variable analysis with mean normal error closest to that of the orthogonal analysis. The dashed lines connecting the mean normal errors for the dependent-variable solutions are erratic. Different dependent-variable choices produced dramatically different planes with widely varying normal errors. The coordinate axis most nearly normal to the plane of the data provided the lowest mean error. Diminishing mean normal error generally corresponds to diminishing planar spread. However, the minimizations of these two parameters are not equivalent, explaining the lack of perfect agreement between the trends plotted in Figure 5. The mean normal errors for most of these cases are less than 0.1, corresponding to planar spreads of less than 12 °.

Fig. 5. Mean squared normal error for several sea-ice samples from orthogonal and x-,y- and z-dependent variable least-squares analyses. Planar spread from the best plane is given for these same cases. The cases are arranged arbitrarily according to of the orthogonal analysis and named in sequence according to the best dependent-variable solution.

The same cases were also analyzed for preferred alignment of the c-axes. Mean squared normal distance between the data and the best vector c is The maximum normalized eigenvalue and the linear spread are displayed in Figure 6. Values of that exceed 0.9 correspond to linear spreads of less than 15 °. Linear spread is larger than planar spread because it represents angles with a particular line in the plane, while planar spread represents angles with the plane itself. Comparing Figures 5 and 6, we observe that increasingly planar sea-ice fabrics do not necessarily correspond to the degree of preferred alignment.

Fig. 6. Maximum normalized eigenvalue and linear spread in degrees for several sea-ice samples. The order and labeling of the cases is the same as in Figure 5.

Schmidt-net plots of the data representing cases z ₁, z ₁₀ ^T and z ₁₁ are presented in Figure 7 and parameters developed from these data are given in Table 1, Included on each net are the intersections with the hemisphere of the plane of best fit, the pole Ρ of this plane, and the vector c of best linear fit. The best line is always contained in the best plane. In case z ₁ the normalized eigenvalues are approximately equal, the planar and linear spreads are large, and the matrix A is non-singular. These conditions indicate that the data are uniformly distributed over the surface of the sphere. Cases z ₁₀ ^T and z ₁₁ are very different from z ₁ but similar to each other. For this group of samples λ₁ is smallest for z ₁₁ indicating that these data are the most closely represented by a plane. Sample z ₁₀ ^T has the smallest value together with a small and is the most linear case in this group. The distribution of data indicated by the eigenvalues is confirmed by the angles that quantify the planar and linear spreads α and δ , respectively, for each sample. The matrix A is nearly singular in the latter two cases as and the condition number is large.

Table 1. Normalized eigenvalues and condition number of matrix A, and planar and linear spread of the data for selected cases. The type of the distribution is interpreted from these parameters

Fig. 7. Schmidt-net plots of cases z ₁, z ₁₀ ^T and z ₁₁ including the intersections with the hemisphere of the best plane,the pole Ρ of this plane, and the best vector c.

A group of samples taken in close proximity and from the same vertical position in the ice sheet should have similar structure. Together, the individual and collective analyses of these samples provide quantitative measures of comparison that indicate the sample size needed to represent the ice fabric at that location. Three samples designated z ₁₀ ^T , z ₁₀ ^M , z ₁₀ ^B were taken from an ice core at 1.3 m from the surface. The Schmidt-net plots of both the data and the computed fits are presented in Figure 7 for z ₁₀ ^T and in Figure 8 for the other samples and the composite data. Each of these cases has nearly the same planar structure, pole and linear structure. This similarity is quantified by the eigenvalues, planar spreads and linear spreads given in Table 1.

Fig. 8. Schmidt-net plots of cases z ₁₀ ^M , z ₁₀ ^B and composite case z₁₀ including the intersections with the hemisphere of the best plane, the pole Ρ of this plane, and the best vector c.

All the Schmidt nets presented above depict the lower hemisphere on the xj-plane. Alternative Schmidt nets are obtained for the hemispheres below the best fit and predominant basal planes using Equation (42) and a Ρ matrix with different eigenvector placement in each case. As a result of these transformations, the position of the data on the net shifts, correcting for sample-preparation error and allowing a visual assessment of the distances from the points to the best line and plane. Figure 9 gives the xy and the alternative nets for sample z ₁₀ ^T . The data are near the perimeter of the net in the standard xy-plot.

Fig. 9. Schmidt-net plots of case z ₁₀ ^T in the standard xy-plane with the z-axis vertical, in the best plane with n vertical, and in the predominant basal plane with cvertical.

The net in the best plane displays a balanced distribution of data on opposite sides of the net. The net in the basal plane depicts the points as a single group near the center of the net where linear distance is accurately represented and the linear and planar fits can be readily evaluated and compared to other samples.

Conclusions

A quantitative characterization of ice fabrics is critical for understanding the mechanical properties of sea ice, but was not previously available. An orthogonal least-squares analysis of uniaxial crystal-orientation data was developed from geometric arguments with unit vectors representing individual crystal orientations. Minimization of the perpendicular distances with a best line or plane provided an eigenvalue problem that was identical to that obtained by other investigators using different methods. Normalized eigenvalues give the mean squared normal distance of the data from the line or plane, and corresponding eigenvectors provide the dominant c-axis, planar and basal plane orientations, and the direction of columnar-crystal elongation. The preferred c-axis orientation is always contained in the plane of best fit. The method is the basis of a relatively simple algorithm for computer analysis of large volumes of orientation data.

The formulation of a least-squares method greatly influences the results. This observation was demonstrated for many samples of first-year sea ice by comparing the mean squared normal distance of the data with planes obtained using the classical dependent-variable least-squares approaches and the orthogonal method. The dependent-variable solutions produced dramatically different planes of best fit with erratic and widely varying normal errors. The error approached the minimum given by the orthogonal method when the dependent-variable direction was almost normal to the best plane. Orthogonal least-squares and other analogous methods producing the eigenvalue problem all rely on the implicit assumption of equal measurement error in all coordinate directions at all points on the unit sphere. We find that this assumption is a first approximation for optical data obtained with the universal stage.

Normalized eigenvalues give the mean squared normal distance between the data on the unit sphere and the plane through the origin normal to the corresponding eigenvector. This measure of the planar fit to data is equivalent to the variance in classical dependent-variable least-squares methods. The majority of the sea-ice samples studied had planar fabrics and several had aligned fabrics, characterized by mean squared normal distances of less than 0.1 with the data. These relative distances are well represented by the angular measures of linear and planar spread; however, the distance and angular error measures are not equivalent. For these samples, increasingly planar orientations of the c-axes do not correspond to increasingly linear fabrics.

The normalized eigenvalues allow quantitative comparisons between samples, and of composite data representing a collection of several samples. Proximate samples from the same vertical position in the ice sheet had nearly identical structures. This similarity, quantified by the computed eigenvalues and eigenvectors, was displayed on Schmidt nets for the individual samples and the composite. The capability to view the data and the analytical results on Schmidt nets in the planes defined by the eigenvectors was developed, providing information to improve the interpretation of the data and the fits by minimizing the importance of sample preparation, net distortion and measurement accuracy limitations.