MinIdent: An application to the identification and classification of amphiboles

Summary Amphibole data in the MinIdent database (Smith andLeibovitz, 1986) were initially entered using species names quoted in the original source. The database has been updated by reclassifying these early data using the program AMPHTAB supplied by N. M. S. Rock and by adding supplemental data from the more recent literature, with the species names again checked using AMPHTAB. Associated MinIdent mineral identification software was utilized to determine which minerals in the database most closely resemble a series of unknown specimens chemically, as expressed in the Chemical Matching Index, CM, a relative figure-of-merit. Chemical data fromMogessie and Tessadri (1982) and Hawthorne (1983) were used to check the agreement between MinIdent and AMPHTAB for the classification of 221 unknown amphiboles.With 450 amphibole analyses entered and compiled in MinIdent, the name assigned by AMPHTAB showed the highest value of CM in MinIdent for 127 of the 221 unknown amphiboles (57.5°/x) and the second highest value for another 32 (14.5%). A chemically adjacent amphibole field had the highest value of CM for 59 of the 221 unknowns (26.7%), where chemically adjacent refers to a change in one chemical parameter. The greatest discrepancy between the two programs occurred in the hornblendes, with an agreement of just 20%, although for 58% of the unknowns the species with the highest CM in MinIdent was in a chemical field adjacent to the species name assigned by AMPHTAB. In many cases the disagreement between MinIdent and AMPHTAB could be ascribed to a lack of data in MinIdent.A comparison of the two programs suggests that the assignment of a single name to an unknown amphibole by AMPHTAB with no direct indication of its reliability may be; misleading. Standard analytical errors are frequently sufficient to overlap the arbitrary boundaries between amphibole species fields. In such cases it may be preferable to use a program such as MinIdent which, rather than assigning an arbitrary amphibole name, presents a list of 20 amphiboles with the degree of similarity between them and the unknown amphibole indicated. MinIdent offers the additional benefit of allowing input of other than chemical data and bases the match between unknown and standard data upon all input data. This will become more of an advantage as instruments such as automated refractometers become available for routine use.

---

Zusammenfassung Ausgangspunkt war das Amphibol-Datenmaterial (Smith und Leibovitz, 1986) mit den dort verwendeten Artnamen. Diese Basisdaten wurden vervollständigt und erneuert durch Reklassifizierung mittels des AMPHTAB Programms, ergänzt durch N. M. S. Rock, und durch Hinzufügung weiterer Daten aus der neuesten Literatur, deren Speciesnamen wiederum mit AMPHTAB überprüft wurden. Außerdem wurde eine MinIdent Mineralidentifizierungs-Software verwendet, um die Minerale zu bestimmen, die in ihrem Chemismus am ehesten einer Serie von unbekannten Amphibol-Species entsprechen, wie sie im Chemical Matching Index (CM) aufscheinen. Zur Klassifikation von 221 unbekannten Amphibolen wurden chemische Daten von Mogessie und Tessadri (1982) verwendet um die Übereinstimmung zwischen MinIdent und AMPHTAB zu überprüfen.Unter den 450 in MinIdent zusammengestellten und eingegebenen Amphibolanalysen zeigen die bei AMPHTAB angegebenen die höchsten CM Werte, nämlich 127 von 221 unbekannten Amphibolen (57,5%) und weitere 32 (14,5%) die zweithöchsten Werte. Innerhalb eines chemisch benachbarten Amphibolfeldes hatten 59 der 221 unbekannten Amphibole (26,7%) die höchsten CM Werte, wobei unter achemisch benachbart die Änderung eines chemischen Parameters zu verstehen ist. Die größten Unterschiede zwischen den beiden Programmen traten bei den Hornblenden auf. Die Übereinstimmung lag bei nur 20%, obwohl bei 58% der unbekannten Amphibole die Species mit dem höchsten CM Wert in MinIdent in ein chemisches Feld zu liegen kamen, welches zu den bei AMPHTAB angegebenen Speciesnamen eine benachbarte Position einnimmt. Die Unterschiede zwischen MinIdent und AMPHTAB könnten in vielen Fällen auf ein Fehlen von Daten in MinIdent zurükzuführen sein.Ein Vergleich beider Programme deutet an, daß die Angabe eines Einzelnamens für ein unbekanntes Amphibol im AMPHTAB Programm ohne Angaben über die Zuverlässigkeit zu Mißverständnissen führen kann. Normale analytische Fehler können bereits dazu führen, daß die Grenzen zweier willkürlicher Amphibolfelder überlappen. In derartigen Fällen emphiehlt sich die Anwendung des MinIdent Programmes, welches eben nicht einen willkürlichen Amphibolnamen angibt, sondern eine Liste von 20 Amphibolen mit dem Grad ihrer Ähnlichkeit, und einem Hinweis auf den unbekannten Amphibol. MinIdent bietet den zusätzlichen Vorteil, daß man außer chemischen auch andere Daten eingeben kann, und stellt dann sämtliche Daten des unbekannten Amphibols den Standard Daten gegenüber. Dieser Klassifizierungsvorgang wird mit der zunehmenden Routineanwendung von automatischen Refraktometern verstärkte Anwendung finden.

     

Generation of correlated properties in heterogeneous porous media

The spatial distribution of rock properties in porous media, such as permeability and porosity, often is strongly variable. Therefore, these properties usefully may be considered as a random field. However, this variability is correlated frequently on length scales comparable to geological lengths (for example, scales of sand bodies or facies). To solve various engineering problems (for example, in the oil recovery process) numerical models of a porous medium often are used. A need exists then to understand correlated random fields and to generate them over discretized numerical grids. The paper describes the general mathematical methods required to do this, with one particular method (the nearest neighbor model) described in detail. How parameters of the mathematical model may be related to rock property statistics for the nearest neighbor model is shown. The method is described in detail in one, two, and three dimensions. Examples are given of how model parameters may be determined from real data.     

Subduction and back-arc activity at the Hikurangi convergent Margin,New Zealand

The Hikurangi Margin is a region of oblique subduction with northwest-dipping intermediate depth seismicity extending southwest from the Kermadec system to about 42°S. The current episode of subduction is at least 16–20 Ma old. The plate convergence rate varies along the margin from about 60 mm/a at the south end of the Kermadec Trench to about 45 mm/a at 42°S. The age of the Pacific lithosphere adjacent to the Hikurangi Trench is not known.The margin divides at about latitude 39°S into two quite dissimilar parts. The northern part has experienced andesitic volcanism for about 18 Ma, and back-arc extension in the last 4 Ma that has produced a back-arc basin onshore with high heaflow, thin crust and low upper-mantle seismic velocities. The extension appears to have arisen from a seawards migration of the Hikurangi Trench north of 39°S. Here the plate interface is thought to be currently uncoupled, as geodetic data indicate extension of the fore-arc basin, and historic earthquakes have not exceededM _s=7.South of 39°S there is no volcanism and a back-arc basin has been produced by downward flexure of the lithosphere due to strong coupling with the subducting plate. Heatflow in the basin is normal. Evidence for strong coupling comes from historic earthquakes of up to aboutM _s=8 and high rates of uplift on the southeast coast of the North Island.The reason for this division of the margin is not known but may be related to an inferred increase, from northeast to southwest, in the buoyancy of the Pacific lithosphere.     

New field and laboratory evidence for the origin of hyaloclastite flows on seamount summits

New field observations with the submersible ALVIN and photographic evidence from a study of the summits of seamounts near the East Pacific Rise show that hyaloclastite deposits occur commonly. Hyaloclastite outcrops were found on six volcanoes at depths from 1240 to 2500 m. These new observations plus laboratory study of new hyaloclastite specimens extend the results of previous studies. Most of the hyaloclastite samples are of hydrovolcanic eruptive origin, but a few show evidence of a predominantly sedimentary origin. Primarily from morphology, we identify several vent areas from which hyaloclastite presumably erupted. The surface appearance of the hyaloclastite deposits varies with distance to these vents, leading us to propose a facies model for deep-sea hyaloclastites on seamount summits. Hyaloclastites of hydromagmatic origin exhibit weak normal grading and bedding-parallel alignment of platy shards. They consist of blocky, sliver and fluidal basalt glass shards and lithics in a matrix that contains pelagic sediment. The shards themselves are remarkably free of even the tiniest crystals and are usually chemically homogeneous. We propose that the shards form mainly by cooling-contraction granulation, but cannot rule out the possibility of limited steam explosivity. Hyaloclastites are closely associated with submarine pahoehoe and we propose that a very rapid eruption rate, promoting clastic-dominated versus flow-dominated eruptive behavior, is the dominant control on hyaloclastite formation. We propose that shard formation occurs during submarine lava fountaining. Gravitational instability of the resulting slurry of shards, sea water and possibly steam causes gravity flow that carries the shards outward from the vent. Further field and modelling studies are needed to test these ideas and more quantitatively constrain the ascent mechanism, eruption dynamics and deposition of deep-sea hyaloclastites.     

Shape analysis of Pacific seamounts

Shape statistics have been compiled from 85 profiles of well-surveyed Pacific seamounts in the height range 140–3800 m. A flat-topped cone was fit to each seamount's cross-sectional profile maintaining the slopes of the sides as closely as possible. On each profile a basal widthd_b, a summit widthd_t, and a maximum heighth, were measured. The height-to-basal-radius ratio isξ_h is estimated by the ratio2hd_b and flatnessf by the ratiod_td_b. Slope angleφ = arctan(ε) is estimated fromε =2h(d_b − d_t). Summit height and basal radius are found to be highly correlated (r = 0.93). The 85-point sample mean of the height-to-basal-radius ratio isξ_h = 0.21 ± 0.08 implying that a seamount's summit height is typically one fifth its basal radius. Despite the high correlation, individual points show some scatter, and there may be groupings into different morphological types. For example, all but one of the seamounts with summit heights above 1000 m have values ofξ_h that are larger than the sample mean. The 85-point sample mean of flatness isf = 0.31 ± 0.18. Data points show a large scatter with values off varying between 0 (a pointy cone) and 0.69 (a flat-topped cone). A histogram representation of flatness, however, indicates that certain values off may be more common than others: the histogram shows a bimodal distribution with maxima occurring at values off in the ranges 0.10–0.20 and 0.35–0.50. Moreover, there is some evidence that the mean flatness decreases with summit height so that the preferred shape of a large-sized seamount may be a pointy cone. Slope angle has an 85-point sample mean ofφ = 18 ± 6°; individual values ofφ vary between 5° and 36°. In addition to having a lower than average mean flatness seamounts with heights above 2600 m also have a lower than average mean slope angle (15°). To determine which variables account for most of the observed variation in the seamount shapes, a multivariate principal component analysis was performed on the data using five shape variables (summit height, basal radius, summit radius, flatness, and slope). The analysis indicates that most of the variation is described by two variables: flatness and summit height.     

On the modelling of sand bedforms using the semivariogram

This study shows the usefulness of the semivariogram for modelling sand ripples created by water flows of varied flow intensity. A combination of two mathematical functions is fitted to each sample semivariogram, that is an exponential (or stochastic) component and a periodic component. The parameters of each of these components have direct physical meaning. A non-dimensional ratio combining the two parameters of the exponential model is interpreted as a regularity index (which increases with the degree of regularity of bedform arrangement). This regularity index is inversely related to the Froude number of the flow. The non-dimensional wavelength, estimated from the dominant periodic function, is also inversely and closely related to the Froude number. The wave height, accurately estimated from properties of the two fitted components, is a direct function of flow velocity and is also proportional to the standard deviation of bed elevations. The bedform shape introduces a considerable discrepancy between the generally assumed normal frequency distribution and the empirical distributions of bed height. The series of bed elevations are generally characterized by a mixture of normal distributions having the same variance but different means. The calculation of a covariance assuming a constant and single mean (as in spectral analysis) can therefore be misleading and the problem may be avoided by using the semivariogram.     

Magnotospheric imaging of high latitude ion outflows

High latitude ion outflows mostly consist of upward streaming O⁺ and He⁺ emanating from the ionosphere. At heights above 1000 km, these flows consist of cold and hot components which resonantly scatter solar extreme ultraviolet (EUV) light, however, the ion populations respond differently to Doppler shifting resulting from the large relative velocities between the ions and the Sun. The possibility of optical detection of the Doppler effect on the scattering rate will be discussed for the O⁺ (83.4 nm) ions. We have contrasted the EUV solar resonance images of these outflows by simulations of the 30.4 nm He⁺ and 83.4 nm O⁺ emissions for both quiet and disturbed geomagnetic conditions. Input data for the 1000 km level has been obtained from the EICS instrument aboard the Dynamics Explorer satellite. Our results show emission rates of 50 and 56 milli-Rayleighs at 30.4 nm for quiet and disturbed conditions and 65 and 75 milli-Rayleighs at 83.4 nm for quiet and disturbed conditions, respectively, obtained for a polar orbiting satellite and viewing radially outward. We also find that an imager at an equatorial distance of 9 R_E or more is in a favorable position for detecting ion outflows, particularly when the plasmapause is depressed in latitude. However, an occultation disk is necessary to obscure the bright plasmaspheric emissions.     

Within-storm variations in runoff and sediment export from a rapidly eroding coal-refuse deposit

Measurements of rainfall, runoff and sediment export from a barren deposit of coal mine refuse in south-western Indiana were collected during three storms in the summer and autumn of 1990. Interfluve sheetwash, sediment mass flux, sediment concentration and, to a lesser extent, trunk gully discharge all responded quickly to changes in rainfall intensity. Grain-size distributions varied considerably during storms, containing exclusively fine-grained sediment at low sediment discharges but very large quantities of coarse (> 2mm) sediment at peak sediment discharges. Although data from a fairly long, multipulsed storm indicate that sediment production is limited by supply, the imbricated layer of flat chips that exists at the surface of the deposit is apparently mobilized during most high-intensity pulses of rainfall, thereby producing large volumes of coarse sediment during summer thunderstorms.