State of stress within a bending spherical shell and its implications for subducting lithosphere

The state of stress within a bending spherical shell has some special features that are caused by sphericity. While most lithospheres are more like spherical shells than flat plates, our ideas of the state of stress have been dominated by flat-plate models. As a consequence, we might be missing some important aspects of the state of stress within subducting lithospheres. In order to examine this problem, we analyse spherical-shell bending problems from basic equations. We present two approaches to solve spherical-shell bending problems: one by the variational approach, which is suitable for global-scale problems, and the other by the asymptotic equation, which is valid to first order in h/R , where h is the thickness of the lithosphere and R is its curvature radius (i.e. under the assumption of small curvature). The form of the equation for displacement shows that wavelengths of deformation are determined by the spherical (elastic) effect and the gravitational buoyancy effect, for which only the latter effect is included in the usual flat-plate formulations. In the case of the Earth, the buoyancy force is dominant and, consequently, spherical effects are suppressed to a large extent; this explains why flat-plate models have been successful for Earth's lithospheric problems. On the other hand, the state of stress shows interesting spherical effects: while bending (fibre) stress along the subduction zone is always important, bending stress along the trench-strike direction can also be important, in particular when the subduction zone arc is small. Numerical results also indicate that compressive normal stress along the trench-strike direction is important when a subduction zone arc is large. These two stresses, the bending stress and the compressive normal stress, both along the trench-strike direction, may have important implications for intraplate earthquakes at subduction zones.     

The solution of the inverse gravity potential problem for cylindrical bodies

Summary. The inverse gravity potential problem consists in the determination of the form and the density of the body by its exterior gravity potential. We describe two similar classes of bodies for which this problem has a unique constructive solution.
(1) The first class contains the cylindrical bodies with finite length, arbitrary form of section and ρ( R , ø, z) =ρ₁( z )ρ₂( R , ø) density distribution, where z is the cylindrical coordinate; R , ø are the polar coordinates in a section plane. This class is important for prospecting geophysics in that it allows us to determine in a unique and constructive way, the function ρ₁( R , ø), the length, form and orientation of the cylinder if we know the function ρ₁( z ) and the exterior potential. The classical moment problem of functions is the basis for the solution of this problem.
(2) The analogous problem for the class of the spherical cylinders, or bodies bounded by arbitrary similar sections of two different concentric spheres and the radial lateral surface, appears when bodies of planetary size are studied. (An example of these bodies would be the Moon mascons.) The density distribution of these cylinders is ρ(τ, θ, ø) =ρ₁(τ)ρ₂(θ, ø) where τ, θ, ø are the spherical coordinates. The function ρ₁(θ, ø), length and form of spherical sections can be uniquely determined by exterior potential if we know the function ρ₁(τ). We propose a new constructive method for harmonic continuation of the gravity potential into the region containing the perturbing masses for the solution of the problem.     

Constraints on the frequency–magnitude relation and maximum magnitudes in the UK from observed seismicity and glacio-isostatic recovery rates

Earthquake populations have recently been shown to have many similarities with critical-point phenomena, with fractal scaling of source sizes (energy or seismic moment) corresponding to the observed Gutenberg–Richter (G–R) frequency–magnitude law holding at low magnitudes. At high magnitudes, the form of the distribution depends on the seismic moment release rate M˙ and the maximum magnitude m _max . The G–R law requires a sharp truncation at an absolute maximum magnitude for finite M˙ . In contrast, the gamma distribution has an exponential tail which allows a soft or 'credible' maximum to be determined by negligible contribution to the total seismic moment release. Here we apply both distributions to seismic hazard in the mainland UK and its immediate continental shelf, constrained by a mixture of instrumental, historical and neotectonic data. Tectonic moment release rates for the seismogenic part of the lithosphere are calculated from a flexural-plate model for glacio-isostatic recovery, constrained by vertical deformation rates from tide-gauge and geomorphological data. Earthquake focal mechanisms in the UK show near-vertical strike-slip faulting, with implied directions of maximum compressive stress approximately in the NNW–SSE direction, consistent with the tectonic model. Maximum magnitudes are found to be in the range 6.3–7.5 for the G–R law, or 7.0–8.2 m _L for the gamma distribution, which compare with a maximum observed in the time period of interest of 6.1 m _L . The upper bounds are conservative estimates, based on 100 per cent seismic release of the observed vertical neotectonic deformation. Glacio-isostatic recovery is predominantly an elastic rather than a seismic process, so the true value of m _max is likely to be nearer the lower end of the quoted range.     

Rapid computation of magnetic anomalies with demagnetization included, for arbitrarily shaped magnetic bodies

Summary. The potential function ø for a magnetic body of susceptibility μ in a medium of susceptibility μ* satisfies the integral equation
Here Φ* is the potential function for the region without the heterogeneity and R is the distance from the point of observation to the point on the surface, s , of the body. δΦ /δn is the normal derivative, in the direction of the outward normal. The equation allows for the effects of demagnetization. For numerical purposes the surfaces can be divided into N facets over which δΦ/δ n is a constant. The unknown quantities δΦ/δn_j can be found from the system of equations defined by:
The prime on the summation sign denotes that the summation does not include the i th element. The magnetic field in the direction of the unit vector P( P ₁, P ₂, P₃ ) is given by:      

Inversion of sonic waveforms with unknown source and receiver functions

A seismogram that is several times the length of the source-receiver wavelet is windowed into two parts—these may overlap—to obtain two seismograms with approximately the same source function but different Green's functions. A similarly windowed synthetic seismogram gives two corresponding synthetic seismograms. The spectral product of the window 1 data with the window 2 synthetic is equal to the spectral product of the window 1 synthetic with the window 2 data only if the correct earth model is used to compute the synthetic. This partition principle is applied to well-log sonic waveform data from Ocean Drilling Project hole 806B, a slow formation, and used there to estimate Poisson's ratio from a single seismogram whose transmitter and receiver functions are unknown. A multichannel extension of the algorithm gives even better results. The effective borehole radius R _b, was included in the inversion procedure, because of waveform sensitivity to R _b. Inversion results for R _b agreed with the sonic caliper, but not the mechanical caliper; thus if R _b is not included in the inversion its value should be taken from the sonic caliper.     

Effects of stress on the two-dimensional permeability tensor of natural fracture networks

The effects of stress on the 2-D permeability tensor of natural fracture networks were studied using a numerical method (Universal Distinct Element Code). On the basis of three natural fracture networks sampled around Dounreay, Scotland, numerical modelling was carried out to examine the fluid flow in relation to the variations in burial depth, differential stress and loading direction. It was demonstrated that the permeability of all the networks decreased with depth due to the closure of aperture. The permeability approached the minimum value at some depth below which little further variation occurred. Also, differential stress had a significant effect on both the magnitude and direction of permeability. The permeability generally decreased with increasing major horizontal stress for a fixed minor horizontal stress, but the various networks considered showed different behaviours. A factor, termed the average deviation angle of maximum permeability ( A _m), was defined to describe quantitatively the deviation degree of the direction of the major permeability component from the applied major stress direction. For networks whose behaviour is controlled by set(s) of systematic fractures, A _m is significantly greater than zero, whereas those comprised of non-systematic fractures have A _m close to zero. In general, fractured rock masses, especially those with one or more sets of systematic fractures, cannot be treated as equivalent porous media. Specification of the geometry of the network is a necessary, but not sufficient, condition for models of fluid flow. Knowledge of the in situ stress, and the deformation it induces, is necessary to predict the behaviour of the rock mass.     

Shear-wave anisotropy: spatial and temporal variations in time delays at Parkfield, Central California

Shear-wave splitting is analysed on data recorded by the High Resolution Seismic Network (HRSN) at Parkfield on the San Andreas fault, Central California, during the three-year period 1988-1990. Shear-wave polarizations either side of the fault are generally aligned in directions consistent with the regional horizontal maximum compressive stress, at some 70° to the fault strike, whereas at station MM in the immediate fault zone, shear-wave polarizations are aligned approximately parallel to the fault. Normalized time delays at this station are found to be about twice as large as those in the rock mass either side. This suggests that fluid-filled cracks and fractures within the fault zone are elastically or seismically different from those in the surrounding rocks, and that the alignment of fault-parallel shear-wave polarizations are associated with some fault-specific phenomenon.
Temporal variations in time delays between the two split shear-waves before and after a M_L = 4 earthquake can be identified at two stations with sufficient data: MM within the fault zone and VC outside the immediate fault zone. Time delays between faster and slower split shear waves increase before the M_L = 4 earthquake and decrease near the time of the event. The temporal variations are statistically significant at 68 per cent confidence levels. Earthquake doublets and multiplets also show similar temporal variations, consistent with those predicted by anisotropic poroelasticity theory for stress modifications to the microcrack geometry pervading the rock mass. This study is broadly consistent with the behaviour observed before three other earthquakes, suggesting that the build-up of stress before earthquakes may be monitored and interpreted by the analysis of shear-wave splitting.     

Magnetic induction fields (2–30 cpd) on Hawaii Island and their implications regarding electrical conductivity in the oceanic mantle

Summary. Geomagnetic time-variations observed at several sites on the island of Hawaii are analysed for the effects of island bathymetry as well as for the inductive response of the deeper mantle. The data are generally consistent with the deep conductivity profile derived using lower frequency, electromagnetic data from the Island of Oahu. Hawaii data fit better if that model is modified to give the upper 200 km of the mantle a lower conductivity of 0.02 S/m compared to 0.1 S/m for Oahu. The data are represented by a complex, frequency-dependent function of location, T _u, relating the vertical variation Z to a component U of the horizontal variation ( T_u = Z/U ). The direction of U is nearly frequency independent at each site but is different for each site. Below a frequency of about 30 cycles per day, the functions, T _u, at any two sites are found to be related by a real constant. This suggests that the deeper conductivity structure is the same beneath each site. This result is consistent with quasi-static induction in a non-uniformly conducting thin sheet above a stratified conductivity structure. The response of such a model can be written as T _u= Aq , where q is a quasi-uniform, complex, frequency-response function characterizing the effect of the deep conductivity and A is a spatially dependent parameter parameterizing the effect of variable conductivity in the thin sheet. The parameter A may be estimated by fitting observational estimates of T _u to models of deep conductivity structure.     

TRM deviations in anisotropic assemblages of multidomain magnetite

Anisotropic assemblages of multidomain magnetite particles develop an anisotropy of magnetic susceptibility (AMS), which in turn induces deviations of thermo-remanent magnetization (TRM) from the field direction. From the theories of multidomain TRM acquisition, it is shown that the TRM anisotropy tensor has its eigenvalue ratios ( T _i) related to the principal weak-field susceptibility ratios ( P _i) by the order of magnitude T _i≃ P ²_i. This relation has been experimentally verified on two sets of highly anisotropic rock samples. The exponent has been determined to be 1.94 in the samples from a Peruvian gabbro, and 1.81 in those from the granite of Flamanville (NW France). Accounting for experimental difficulties in determining the TRM anisotropy tensors, these exponents are judged to agree well with the expected one. It is therefore stressed that AMS measurements provide a good means of evaluating the magnetic field direction from deviated TRM directions, providing magnetic carriers are mainly multidomain magnetites.     

Calculating horizontal stress orientations with full or partial knowledge of the tectonic stress tensor

Earthquakes potentially serve as abundant and cost-effective gauges of tectonic stress provided that reliable means exist of extracting robust stress parameters. Several algorithms have been developed for this task, each of which typically provides information on the orientations of the three principal stresses and a single stress magnitude parameter. A convenient way of displaying tectonic stress results is to map the azimuth of maximum horizontal compressive stress, which is usually approximated using the azimuth of the larger subhorizontal principal stress. This approximation introduces avoidable errors that depend not only on the principal stress axes' plunges but also on the value of the stress magnitude parameter. Here we outline a method of computing the true direction of maximum horizontal compressive stress ( S _H) and show that this computation can be performed using only the four stress parameters obtained in routine focal mechanism stress estimation. Using theoretical examples and new stress inversion results obtained with focal mechanism data from the central Grímsey lineament, northern Iceland, we show that the S _H axis may differ by tens of degrees from its commonly adopted proxy. In order to most appropriately compare tectonic stress estimates with other geophysical parameters, such as seismic fast directions or geodetically measured strain rate tensors, or to investigate spatiotemporal variations in stress, we recommend that full use be made of the routinely estimated stress parameters and that a formal axis of maximum horizontal compression be calculated.     

Aseismic fault movement before the 1995 Kobe earthquake detected by a GPS survey: implication for preseismic stress localization?

By inversion analysis of the baseline changes and horizontal displacements observed with GPS (Global Positioning System) during 1990–1994, a high-angle reverse fault was detected in the Shikoku-Kinki region, southwest Japan. The active blind fault is characterized by reverse dip-slip (0.7±0.2  m yr⁻¹ within a layer 17–26  km deep) with a length of 208±5  km, a (down-dip) width of 9±2  km, a dip-angle of 51°±2° and a strike direction of 40°±2° (NE). Evidence from the geological investigation of subfaults close to the southwestern portion of the fault, two historical earthquakes ( M _L=7.0, 1789 and 6.4, 1955) near the centre of the fault, and an additional inversion analysis of the baseline changes recorded by the nationwide permanent GPS array from 18 January to 31 December 1995 partially demonstrates the existence of the fault, and suggests that it might be a reactivation of a pre-existing fault in this region. The fact that hardly any earthquakes ( M _L>2.0) occurred at depth on the inferred fault plane suggests that the fault activity was largely aseismic. Based on the parameters of the blind fault estimated in this study, we evaluated stress changes in this region. It is found that shear stress concentrated and increased by up to 2.1 bar yr⁻¹ at a depth of about 20  km around the epicentral area of the 1995 January 17  Kobe earthquake ( M _L=7.2, Japan), and that the earthquake hypocentre received a Coulomb failure stress of about 5.6 bar yr⁻¹ during 1990–1994. The results suggest that the 1995  Kobe earthquake could have been induced or triggered by aseismic fault movement.     

A study of geomagnetic storms

Summary. An attempt is made to find interplanetary magnetic field and solar-wind parameters which control the development of geomagnetic storms. For this purpose, the interplanetary energy flux is estimated in terms of the Poynting flux ( E × B /4π), and its time variations are compared with the rate of energy dissipation in terms of the ring-current particle injection u _i( t ), Joule dissipation in the ionosphere u_j ( t ) and auroral particle injection u_p ( t ) for 15 major geomagnetic storms.
It is shown that the growth of geomagnetic storms, namely the time variations of the rate of the total energy dissipation, u ( t ) = u _i( t ) + u _j( t ) + u _p( t ), is closely related to the Poynting flux by the following relation:
where l ₀≅ 7 R _E and θ' is a measure of the angle between the interplanetary magnetic field vector and the magnetospheric field vector at the front of the magnetosphere in the equatorial plane. Further, it is shown that within a factor of 2 for each storm period.
A large increase of u ( t ) is associated with substorm activity. Thus, the energy flux ɛ( t ) entering the magnetosphere is dissipated through magneto-spheric substorm processes within the magnetosphere, and their accumulated effects can be understood as geomagnetic storm phenomena.     

Dynamic stress field of a kinematic earthquake source model with k-squared slip distribution

Source models such as the k -squared stochastic source model with k -dependent rise time are able to reproduce source complexity commonly observed in earthquake slip inversions. An analysis of the dynamic stress field associated with the slip history prescribed in these kinematic models can indicate possible inconsistencies with physics of faulting. The static stress drop, the strength excess, the breakdown stress drop and critical slip weakening distance D _c distributions are determined in this study for the kinematic k -squared source model with k -dependent rise time. Several studied k -squared models are found to be consistent with the slip weakening friction law along a substantial part of the fault. A new quantity, the stress delay, is introduced to map areas where the yielding criterion of the slip weakening friction is violated. Hisada's slip velocity function is found to be more consistent with the source dynamics than Boxcar, Brune's and Dirac's slip velocity functions. Constant rupture velocities close to the Rayleigh velocity are inconsistent with the k -squared model, because they break the yielding criterion of the slip weakening friction law. The bimodal character of D _c / D _tot frequency–magnitude distribution was found. D _c approaches the final slip D _tot near the edge of both the fault and asperity. We emphasize that both filtering and smoothing routinely applied in slip inversions may have a strong effect on the space–time pattern of the inferred stress field, leading potentially to an oversimplified view of earthquake source dynamics.     

Source mechanism of the deep Colombian earthquake of 1970 July 31 from the free oscillation data

Summary. The method proposed by Mendiguren to determine the source parameters from free oscillation data is applied to the 1970 July 31 deep Colombian earthquake. The results indicate a source propagating horizontally for about 150 km along the lithosphere and cutting across its width. The slab behaves as a guide for source propagation. The horizontal propagation velocity is determined as 3.8 km/s. The intensity of the source grew proportionately to the second power of the propagation distance. This rate of source intensity growth may be interpreted either by a fan-shaped fault model or by a cone-shaped volume source. The average slip and stress drop are estimated as 360 cm and 300 bar for the fault model. For the volume source model the transformational shear strain and stress are estimated as 11 × 10⁻⁵ and 160 bar. There is no evidence of a double couple radiation preceding the P origin time. It is shown that the isotropic and deviatoric components of the moment tensor cannot be uniquely resolved when only observations of a single mode are available. It is observed that, from a statistical basis, the available ₀ S_n data for Colombian shock can be equally well explained by a pure deviatoric source model or by a source model including an isotropic component. Numerical experiments indicate that the inclusion of higher mode data does not change this situation. But, on the other hand, numerical experiments show that the available data and the scheme used for the inversion would not result in a solution including an artificial implosive component if the actual source were pure deviatoric. If the departure from a pure deviatoric source is produced by noise, it has to be non-random, as it could be produced by lateral heterogeneities not included in the inversion scheme.     

Surface motion of a fluid planet induced by impacts

In order to approximate the free-surface motion of an Earth-sized planet subjected to a giant impact, we have described the excitation of body and surface waves in a spherical compressible fluid planet without gravity or intrinsic material attenuation for a buried explosion source. Using the mode summation method, we obtained an analytical solution for the surface motion of the fluid planet in terms of an infinite series involving the products of spherical Bessel functions and Legendre polynomials. We established a closed form expression for the mode summation excitation coefficient for a spherical buried explosion source, and then calculated the surface motion for different spherical explosion source radii a (for cases of   a / R = 0.001  to 0.035, R is the radius of the Earth) We also studied the effect of placing the explosion source at different radii r ₀ (for cases of   r ₀/ R = 0.90  to 0.96) from the centre of the planet. The amplitude of the quasi-surface waves depends substantially on a / R , and slightly on   r ₀/ R   . For example, in our base-line case,   a / R = 0.03, r ₀/ R = 0.96  , the free-surface velocity above the source is 0.26 c , whereas antipodal to the source, the peak free surface velocity is 0.19 c . Here c is the acoustic velocity of the fluid planet. These results can then be applied to studies of atmosphere erosion via blow-off caused by asteroid impacts.     

Dynamic nucleation process of shallow earthquake faulting in a fault zone

Two distinct phases are commonly observed at the initial part of seismograms of large shallow earthquakes: low-frequency and low-amplitude waves following the onset of a P wave ( P ₁) are interrupted by the arrival of the second impulsive phase P₂ enriched with high-frequency components. This observation suggests that a large shallow earthquake involves two qualitatively different stages of rupture at its nucleation.
We propose a theoretical model that can naturally explain the above nucleation behaviour. The model is 2-D and the deformation is assumed to be anti-plane. A key clement in our model is the assumption of a zone in which numbers of pre-existing cracks are densely distributed; this cracked zone is a model for the fault zone. Dynamic crack growth nucleated in such a zone is intensely affected by the crack interactions, which exert two conflicting effects: one tends to accelerate the crack growth, and the other tends to decelerate it. The accelerating and decelerating effects are generally ascribable to coplanar and non-coplanar crack interactions, respectively. We rigorously treat the multiple interactions among the cracks, using the boundary integral equation method (BIEM), and assume the critical stress fracture criterion for the analysis of spontaneous crack propagation.
Our analysis shows that a dynamic rupture nucleated in the cracked zone begins to grow slowly due to the relative predominance of non-coplanar interactions. This process radiates the P₁ phase. If the crack continues to grow, coalescence with adjacent coplanar cracks occurs after a short time. Then, coplanar interactions suddenly begin to prevail and crack growth is accelerated; the P₂ phase is emitted in this process. It is interpreted that the two distinct phases appear in the process of the transition from non-coplanar to coplanar interaction predominance.     

Lithospheric flexure due to prograding sediment loads: implications for the origin of offlap/onlap patterns in sedimentary basins

Abstract Simple elastic plate models have been used to determine the stratigraphic patterns that result from prograding sediment loads. The predicted patterns, which include coastal offlap/onlap and downlap in a basinward direction, are generally similar to observations of stratal geometry from Cenozoic sequences of the U.S. Atlantic and Gulf Coast margins. Coastal offlap is a feature of all models in which the water depth and elastic thickness of the lithosphere, T _e (which is a measure of the long-term strength of the lithosphere), are held constant, and is caused by a seaward shift in the sediment load and its compensation as progradation proceeds. The coastal offlap pattern is reduced if sediments prograde into a subsiding basin, since subsidence causes an increase in the accommodation space and loading landward of a prograding wedge. The stratal geometry that results is complex, however, and depends on the sediment supply, the amount of subsidence, and T _e. If the sediment supply to a subsiding basin proceeds in distinct 'pulses' (due, say, to different tectonic events in a source region) then it is possible to determine the relationship between stratal geometry and T _e. Coastal offlap and downlap are features of most models where the lithosphere either has a constant T _e slowly increases T_e with time, or changes T _e laterally; however, in the case where sediments prograde onto lithosphere that rapidly increases T _e with rime, the offlap can be replaced by onlap. Lithospheric flexure due to prograding sediment loads is capable of producing a wide variety of stratal geometries and may therefore be an important factor to take into account when evaluating the relative role of tectonics and eustatic sea-level changes in controlling the stratigraphic record.     

Run-up from impact tsunami

We report on calculations of the on-shore run-up of waves that might be generated by the impact of subkilometre asteroids into the deep ocean. The calculations were done with the COULWAVE code, which models the propagation and shore-interaction of non-linear moderate- to long-wavelength waves  ( kh < π)  using the extended Boussinesq approximation. We carried out run-up calculations for several different situations: (1) laboratory-scale monochromatic wave trains onto simple slopes; (2) 10–100 m monochromatic wave trains onto simple slopes; (3) 10–100 m monochromatic wave trains onto a compound slope representing a typical bathymetric profile of the Pacific coast of North America; (4) time-variable scaled trains generated by the collapse of an impact cavity in deep water onto simple slopes and (5) full-amplitude trains onto the Pacific coast profile. For the last case, we also investigated the effects of bottom friction on the run-up. For all cases, we compare our results with the so-called 'Irribaren scaling': The relative run-up   R / H ₀=ξ= s ( H ₀/ L ₀)^−1/2  , where the run-up is   R , H ₀  is the deep-water waveheight, L ₀ is the deep-water wavelength, s is the slope and ξ is a dimensionless quantity known as the Irribaren number. Our results suggest that Irribaren scaling breaks down for shallow slopes   s ≤ 0.01  when  ξ < 0.1 − 0.2  , below which   R / H ₀  is approximately constant. This regime corresponds to steep waves and very shallow slopes, which are the most relevant for impact tsunami, but also the most difficult to access experimentally.     

Transient and impulse responses of a one-dimensional linearly attenuating medium—II. A parametric study

Summary. We investigate one-dimensional waves in a standard linear solid for geophysically relevant ranges of the parameters. The critical parameters are shown to be T*= t_u/Q_m where t _u is the travel time and Q_m the quality factor in the absorption band, and τ⁻¹ _m , the high-frequency cut-off of the relaxation spectrum. The visual onset time, rise time, peak time, and peak amplitude are studied as functions of T* and τ _m. For very small τ _m , this model is shown to be very similar to previously proposed attenuation models. As τ _m grows past a critical value which depends on T* , the character of the attenuated pulse changes. Seismological implications of this model may be inferred by comparing body wave travel times with a'one second'earth model derived from long-period observations and corrected for attenuation effects assuming a frequency independent Q over the seismic band. From such a comparison we speculate that there may be a gap in the relaxation spectrum of the Earth's mantle for relaxation times shorter than about one second. However, observational constraints from the attenuation of body waves suggest that such a gap might in fact occur at higher frequencies. Such a hypothesis would imply a frequency dependence of Q in the Earth's mantle for short-period body waves.     

The best estimate of Fisher's precision parameter k

Summary. Three estimates for Fisher's precision parameter, K , are derived and their relevance to inference testing in palaeomagnetic work is discussed. It is shown that the statistic K (= ( N - 2)/ N - R )), an unbiased estimate for K , is essentially irrelevant. The maximum likelihood estimate K (= N /( N — R )) is shown to be the appropriate estimate for the 'goodness of fit' test. The estimate k _F(=( N - 1)/( N - R )), commonly referred to as 'the best estimate', is shown to be the appropriate statistic for comparisons of sample dispersions and its derivation by Fisher is re-examined. Further, (1/ k _F) is shown to be an unbiased estimate for (1/ K ) and to be analogous to the unbiased estimate of variance in a normal distribution. The limits of applicability of these various results are examined.