Non-Linear Theory and Power-Law Models for Information Integration and Mineral Resources Quantitative Assessments

Singular physical or chemical processes may result in anomalous amounts of energy release or mass accumulation that, generally, are confined to narrow intervals in space or time. Singularity is a property of different types of non-linear natural processes including cloud formation, rainfall, hurricanes, flooding, landslides, earthquakes, wildfires, and mineralization. The end products of these non-linear processes can be modeled as fractals or multifractals. Hydrothermal processes in the Earth’s crust can result in ore deposits characterized by high concentrations of metals with fractal or multifractal properties. Here we show that the non-linear properties of the end products of singular mineralization processes can be applied for prediction of undiscovered mineral deposits and for quantitative mineral resource assessment, whether for mineral exploration or for regional, national and global planning for mineral resource utilization. In addition to the general theory and framework for the non-linear mineral resources assessment, this paper focuses on several power-law models proposed for characterizing non-linear properties of mineralization and for geoinformation extraction and integration. The theories, methods, and computer system discussed in this paper were validated using a case study dealing with hydrothermal Au mineral potential in southern Nova Scotia, Canada.     

Inversely-Mapped Analytical Solutions for Flow Patterns around and within Inclined Elliptic Inclusions in Fluid-Saturated Rocks

In this paper, an inverse mapping is used to transform the previously-derived analytical solutions from a local elliptical coordinate system into a conventional Cartesian coordinate system. This enables a complete set of exact analytical solutions to be derived rigorously for the pore-fluid velocity, stream function, and excess pore-fluid pressure around and within buried inclined elliptic inclusions in pore-fluid-saturated porous rocks. To maximize the application range of the derived analytical solutions, the focal distance of an ellipse is used to represent the size of the ellipse, while the length ratio of the long axis to the short one is used to represent the geometrical shape of the ellipse. Since the present analytical solutions are expressed in a conventional Cartesian coordinate system, it is convenient to investigate, both qualitatively and quantitatively, the distribution patterns of the pore-fluid flow and excess pressure around and within many different families of buried inclined elliptic inclusions. The major advantage in using the present analytical solution is that they can be conveniently computed in a global Cartesian coordinate system, which is widely used in many scientific and engineering computations. As an application example, the present analytical solutions have been used to investigate how the dip angle of an inclined elliptic inclusion affects the distribution patterns of the pore-fluid flow and excess pore-fluid pressure when the permeability ratio of the elliptic inclusion is of finite but nonzero values.     

On the use of gravity and magnetic anomalies for locating probable areas of metallic mineralization in South Sinai, Egypt

The present study deals with processing and interpreting the potential gravity and magnetic data in order to locate promising sites for metallic mineral occurrences in the basement rocks of South Sinai. Two promising sites were anticipated by combining the geophysical results and the available geological information. Each site was profiled and interpreted to the nearly approximated subsurface feature model, and the geometrical parameters of the expected subsurface structures may form mineralized bodies such as depth, thickness, width, dipping, density contrast, and magnetic susceptibility contrast could be delineated. Detailed ground Very Low Frequency (VLF) and Vertical Magnetic Gradient (VMG) methods were carried out in the site A2, eight VLF and VMG anomalies were determined, and their depths were calculated.     

Changing grading of soil: effect on critical states

Examples of situations are presented where the grading of a soil changes during its lifetime either by crushing of particles leading to an increase of fine material or by slow transport of fine particles with seepage leading to a decrease of fine material. Such grading changes influence the basic constitutive properties of the soil, in particular properties such as critical states which are dependent on the available range of densities of packing. Discrete element modelling is used to show the dependence of critical state conditions on grading and the way in which the particle assembly seeks out new critical state conditions as the grading changes.     

The gravitational field of topographic-isostatic masses and the hypothesis of mass condensation II-the topographic-isostatic geoid

In Part I we focussed on a convergent representation of the gravitational potential generated bytopographic masses on top of the equipotential surface atMean Sea Level, thegeoid, and by those masses which compensate topography. Topographic masses have also been condensated, namely represented by a single layer. Part II extends the computation of the gravitational field of topographic-isostatic masses by a detailed analysis of itsforce field in terms ofvector-spherical harmonic functions. In addition, the discontinuous mass-condensated topographic gravitational force vector (head force) is given. Once we identify theMoho discontinuity asone interface of isostatically compensated topographical masses, we have computed the topographic potential and the gravitational potential which is generated by isostatically compensated masses atMean Sea Level, the geoid, and illustrated by various figures of geoidal undulations. In comparison to a data oriented global geoid computation ofJ. Engels (1991) the conclusion can be made that the assumption of aconstant crustal mass density, the basic condition for isostatic modeling, does not apply. Insteaddensity variations in the crust, e.g. betweenoceanic and continental crust densities, have to be introduced in order to match the global real geoid and its topographic-isostatic model. The performed analysis documents that thestandard isostatic models based upon aconstant crustal density areunreal.     

LREE distribution in perovskite,apatite and titanite from South West Ugandan xenoliths and kamafugite lavas

Summary In-situ microprobe LREE analyses of perovskite and titanite (La, Ce, Nd), and apatite (La, Ce), from SW Ugandan clinopyroxenite xenoliths and kamafugite lavas indicate that LREE distribution in these minerals is determined by a number of factors related to their different parageneses: In particular LREE content is affected by whether the LREE-bearing minerals have crystallised from metasomatic carbonate or from silicate (i.e. metasomatic or magmatic) melts in the mantle. In this situation LREE partition favours carbonate over silicate melts. Distribution of LREE in perovskite and apatite crystallised from magmatic mantle melts or mantle-derived lavas is chiefly determined by preference of LREE for perovskite > apatite > titanite. LREE zoning in perovskite is influenced by changes in melt structure: increasing melt polymerisation enhancing mineral_LREE/melt_LREE partition into perovskite rims in magmatic xenoliths; decreasing melt polymerisation depleting LREE in lava perovskite rims. This zoning is reinforced by perovskite competition with apatite for LREE: perovskite (cores/rims) co-crystallising with apatite is reduced in LREE. There are 37 instances of perovskitewith Ce below detection while La and Nd levels are normal. These occur in both xenoliths and lavas; in grain zones or whole grains. Likewise Ce alone of the LREE is below detection in six out of ten titanite analyses. These observations are interpreted as evidence for increased f_{O 2}, Ce⁴ + being excluded from these mineral structures. Recognition of these various processes can elucidate the interpretation of bulk rock and bulk mineral LREE signatures in kamafugite volcanism.

---

LREE Verteilung in Perovskit, Apatit und Titanit aus Xenolithen und kamafugitischen Laven Südwest-Ugandas

Zusammenfassung In-situ LREE Analysen von Perovskit und Titanit (La, Ce, Nd) und Apatit (La, Ce) aus Klinopyroxenit-Xenolithen und kamafugitischen Laven Südwest-Ugandas zeigen, daß die LREE Verteilung in diesen Mineralen durch eine Vielzahl von Faktoren, die mit Unterschieden in den Paragenesen zusammenhängen, bestimmt wird: Der LREE-Gehalt wird im besonderen davon bestimmt, ob die LREE-führenden Minerale aus metasomatischen Karbonat- oder aus (metasomatischen oder magmatischen) Silikatschmelzen im Mantel auskristallisierten. Dabei erfolgt die LREE Fraktionierung zu Gunsten der Karbonatschmelzen. Die LREE-Verteilung von Perovskit und Apatit, die aus magmatischen Mantelschmelzen oder -laven kristallisierten, wird vorrangig durch den bevorzugten Einbau der LREE in Perovskit > Apatit > Titanit kontrolliert. Der LREE Zonarbau von Perovskit wird durch die Änderungen der Schmelzstruktur beinflußt: Verstärkte Schmelzpolymerisation führt zu verstärkter Mineral_LFEE/Schmelze_LREE Fraktionierung in den Perovskiträndern magmatischer Xenolithe, eine Abnahme der Schmelzpolymerisation hingegen resultiert in einer Abreicherung der LREE in den Perovskiträndern. Diese Art der Zonierung wird durch den Wettbewerb von Perovskit mit Apatit um die LREE verstärkt. Perovskit (Kerne/Ränder), der mit Apatit gemeinsam auskristallisierte, ist ärmer an LREE. 37 Fälle, in denenCe nicht nachweisbar war, La und Nd aber in normaler Konzentration auftreten, wurden sowohl in den Xenolithen als auch in den Laven gefunden; und zwar entweder in Kornbereichen oder in ganzen Körnern. Vergleichsweise liegt Ce nur in sechs von zehn Titanitproben unterhalb der Nachweisgrenze. Diese Beobachtungen werden als Hinweise auf erhöhte SauerstoffFugazitäten, bei denen Ce^4– aus der Mineralstruktur ausgeschlossen wird, angesehen.Ein Verständnis dieser verschiedenen Prozesse kann zur besseren Interpretation von LREE Gesamtgesteins- und Gesamtmineral-Signaturen in Kamafugiten beitragen.

---

---

With 3 Figures     

Detection of a new interstellar molecular ion, H2COH+ (protonated formaldehyde)

A new interstellar molecular ion, H2COH+ (protonated formaldehyde), has been detected toward Sgr B2, Orion KL, W51, and possibly in NGC 7538 and DR21(OH). Six transitions were detected in Sgr B2(M). The 1(1,0)-1(0,1) transition was detected in all sources listed above. Searches were also made toward the cold, dark clouds TMC-1 and L134N, Orion (3N, 1E), and a red giant, IRC + 10216, without success. The excitation temperatures of H2COH+ are calculated to be 60-110 K, and the column densities are on the order of 10(12)-10(14) cm-2 in Sgr B2, Orion KL, and W51. The fractional abundance of H2COH+ is on the order of 10(-11) to 10-(9), and the ratio of H2COH+ to H2CO is in the range 0.001-0.5 in these objects. The values in Orion KL seem to be consistent with the "early time" values of recent model calculations by Lee, Bettens, & Herbst, but they appear to be higher than the model values in Sgr B2 and W51 even if we take the large uncertainties of column densities of H2CO into account. We suggest production routes starting from CH3OH may play an important role in the formation of H2COH+.     

Towards a nonperturbation transport theory in hetereogeneous aquifers

A large body of existing theories of flow and contaminant transport in aquifers ignore the presence of recharge, eliminate the boundary conditions, neglect transient conditions in groundwater flow, conceive hydraulic gradients as linear, and require parameter variability to be stationary and Gaussian. The most outstanding and difficult to justify assumption is the subjective small size of the stochastic terms (i.e., small perturbation methods), which usually is forced by considering the logarithm of the hydraulic conductivity. Several problems in flow and contaminant subsurface hydrology, such as the enhanced dispersion parameters with plume size or time after injection, remain to be observed in the light of a stochastic theory that allows a more realistic consideration of physical and hydrologic properties. In this article, an attempt is made to reformulate a contaminant transport equation (the variable dispersion equation, VDE) with transport parameters in terms of regional hydrologic and aquifer hydraulic properties, such as recharge rate, spatially random transmissivity, hydraulic gradient, aquifer thickness, and soil porosity. Subsequently, a general analytic procedure, the method of decomposition, is used to derive a solution to the VDE. This procedure does not require small perturbation, logarithmic transformations, or specific probability law assumptions. Comparison tests with existing theoretical and field results are given. The tests illustrate the enhanced dispersion and shifting concentration effects produced by the variable dispersion equation. Finally a generalization of the method to nonstationary dispersion in three-dimensional domains is proposed.     

Trace Element and Isotope Geochemistry of the Volcanic Rocks of Bequia, Grenadine Islands, Lesser Antilles Arc: a Study of Subduction Enrichment and Intra-crustal Contamination

Pliocene volcanics on the island of Bequia comprise two interbeddedsuites of basalts and andesites. The isotopically homogeneoussuite (IHS) has a limited range of Sr—Nd—Pb isotopes(⁸⁷Sr/⁸⁶Sr 0.7040–0.7046, ¹⁴³ Nd/¹⁴⁴ Nd 0.5130 and ²⁰⁶Pb/²⁰⁴Pb 19.36–19.51), and mantle-like ¹⁸O values (5.5in clinopyroxene). The isotopically diverse suite (IDS) is characterizedby much wider ranges of radiogenic isotopes (⁸⁷ Sr/⁸⁶Sr 0.7048–0.7077,¹⁴³ Nd/¹⁴⁴ Nd 0.5128–0.5123 and ²⁰⁶ Pb/²⁰⁴ Pb 19.7–20.2),in which all of the Sr and Pb ratios are higher and Nd ratiosare lower than those of the IHS. The IDS is also characterizedby high ¹⁸ O values, up to 7.6 in clinopyroxene. The Sr andPb isotope ratios are too high, and the Nd isotope ratios aretoo low in the IDS for any of these lavas to be derived fromunmodified depleted mantle. Both suites are petrologically very similar and their majorelement compositions and phenocryst contents suggest that theywere formed largely by fractional crystallization of a hydroustholeiitic melt at pressures <3 kbar. The isotopic ratiosand enrichments in large ion lithophile elements (LILE), andto some extent light rare earth elements (LREE), as comparedwith mid-ocean ridge basalts (MORB), of the IHS lavas suggestthat they were derived from a depleted mantle source which hadbeen re-enriched by the addition of 1–4% of a subductioncomponent. This component probably comprised a mixture of dehydrationfluids, and perhaps minor siliceous melts, released from subductingsediments and mafic crust. The extreme isotopic ranges, largeenrichments in incompatible elements, more fractionated LREEpatterns and higher ¹⁸ O values of the IDS lavas are interpretedas resulting from 10–55% assimilation—fractionalcrystallization of sediments, derived from the Guyana Shield,which are present in the arc crust, by IHS type melts. KEY WORDS: trace elements; radiogenic isotopes; arc lavas; Lesser Antilles ^*Corresponding author.     

Total ozone variations in the tropical belt: An application for quality of ground based measurements

Summary The study of the regime of ozone variations in the huge tropical belt (25° S to 25° N), which are, in general, very small and zonally nearly symmetric, permits to establish a statistical model for estimating the ozone deviations using Total Ozone Mapping Spectrometer (TOMS) data. The equatorial stratospheric winds at 25 and 50hPa and the solar flux at 10.7 cm are used as major predictors and the linear trend was also estimated. The 10m/sec stratospheric wind change is related to1.2% ozone change at the equator, to practically no change in the 8–15° belts and up to 1.4% change with opposite phase over the tropics in spring but nearly zero change in fall. The solar cycle related amplitude is about 1.4% per 100 units of 10.7 cm solar flux. The ozone trends are negative: not significant over the equator and about –2% per decade (significant at 95% level) over the tropics. The latter could have been enforced by the 2 to 4% lower ozone values during 1991–1993, part of which might be related to the effects of the Mt. Pinatubo eruption, but might also be due to the strong QBO. The estimated deviations are verified versus reliable observations and the very good agreement permits applying the model for quantitative quality control of the reported ozone data from previous years. The standard deviation of the difference between observed ozone deviations and those estimated from the model is only 0.9–1.6% for yearly mean, that means instruments used for total ozone observations in the tropical belt should have systematic error of less than 1%. Cases when the discrepancies between the model and reported observations at a given station exceed 2–3% for time interval of 2 or more years should be verified.With 17 Figures