The GPS Toolbox

The GPS Toolbox is dedicated to highlighting algorithms utilized by GPS engineers and scientists. If you have an interesting algorithm you would like to share with our readers or if you have a topic you would like to see covered in a future column, contact us at gps-toolbox@ngs.noaa.gov. To comment on the algorithms presented here, or to leave a request for an algorithm you may be looking for, visit our Web site (http://www.ngs.noaa.gov/gps-toolbox). ? 2000 John Wiley & Sons, Inc.     

Evidence for the occurrence of hydrogeochemical processes in the groundwater of Khoy plain,northwestern Iran,using ionic ratios and geochemical modeling

Rapid population growth, industrialization, and agricultural expansion in the Khoy area (northwestern Iran) have led to its dependence on groundwater and degradation of groundwater quality. This study attempts to decipher the major processes and factors that degrade the groundwater quality of the Khoy plain. For this purpose, 54 groundwater samples from unconfined and confined aquifers of the plain were collected in July 2017 and analyzed for major cations and anions (Na, K, Ca, Mg, HCO₃, SO_4, and Cl), minor ions (NO₃ and F), and Al. Magnesium and bicarbonate were identified as the dominant cation and anion, respectively. Several ionic ratios and geochemical modeling using PHREEQC indicated that the most important hydrogeochemical processes to affect groundwater quality in the plain were weathering and dissolution of evaporitic and silicate minerals, mixing, and ion exchange. There were smaller effects from evaporation and anthropogenic factors (e.g., industries). Results showed that the high salinity of the groundwater in the northeast area of the plain was due to the high solubility of the evaporitic minerals, e.g., halite and gypsum. Reverse ion exchange and the contribution of mineral dissolution were more significant than ion exchange in the northeastern part of the plain. Elevated salinity of the groundwater in the southeast was attributed mostly to reverse ion exchange and somewhat to evaporation.     

Distribution of Indus born mica along Gulf of Kachchh coast: Implications in understanding current dynamics

Indus is one of the major sources of sediments to the Gulf of Kachchh. Yet only its <63 micron fraction is studied in detail with regards to the offshore current dynamics. Hence here we present our study on characteristic signature of the Indus sediment load (i.e. mica minerals) in >63 micron size fraction along the coast of Gulf of Kachchh. The spatial distribution of mica minerals along the Gulf of Kachchh coast was studied which showed in general decreasing trend as we move along the northern and southern coast of the Gulf of Kachchh but, an increase in amount near the southern mouth at Okha. The study shows that the earlier proposed tidal barrier is ineffective in restricting movement of mica across the mouth of the gulf due to its characteristic transport mechanism. Also the presence of mudflats along the gulf of Kachchh coast plays a vital role as sediment receptors in the active sediment transport processes and mica minerals prove to be a promising simple tracer in studying the Indus born sediments in the region.     

Computing the Sensitivity Kernels for 2.5-D Seismic Waveform Inversion in Heterogeneous,Anisotropic Media

2.5-D modeling and inversion techniques are much closer to reality than the simple and traditional 2-D seismic wave modeling and inversion. The sensitivity kernels required in full waveform seismic tomographic inversion are the Fréchet derivatives of the displacement vector with respect to the independent anisotropic model parameters of the subsurface. They give the sensitivity of the seismograms to changes in the model parameters. This paper applies two methods, called ‘the perturbation method’ and ‘the matrix method’, to derive the sensitivity kernels for 2.5-D seismic waveform inversion. We show that the two methods yield the same explicit expressions for the Fréchet derivatives using a constant-block model parameterization, and are available for both the line-source (2-D) and the point-source (2.5-D) cases. The method involves two Green’s function vectors and their gradients, as well as the derivatives of the elastic modulus tensor with respect to the independent model parameters. The two Green’s function vectors are the responses of the displacement vector to the two directed unit vectors located at the source and geophone positions, respectively; they can be generally obtained by numerical methods. The gradients of the Green’s function vectors may be approximated in the same manner as the differential computations in the forward modeling. The derivatives of the elastic modulus tensor with respect to the independent model parameters can be obtained analytically, dependent on the class of medium anisotropy. Explicit expressions are given for two special cases—isotropic and tilted transversely isotropic (TTI) media. Numerical examples are given for the latter case, which involves five independent elastic moduli (or Thomsen parameters) plus one angle defining the symmetry axis.     

Classifying zones of suitability for manual drilling using textural and hydraulic parameters of shallow aquifers: a case study in northwestern Senegal

A method is proposed that uses analysis of borehole stratigraphic logs for the characterization of shallow aquifers and for the assessment of areas suitable for manual drilling. The model is based on available borehole-log parameters: depth to hard rock, depth to water, thickness of laterite and hydraulic transmissivity of the shallow aquifer. The model is applied to a study area in northwestern Senegal. A dataset of boreholes logs has been processed using a software package (TANGAFRIC) developed during the research. After a manual procedure to assign a standard category describing the lithological characteristics, the next step is the automated extraction of different textural parameters and the estimation of hydraulic conductivity using reference values available in the literature. The hydraulic conductivity values estimated from stratigraphic data have been partially validated, by comparing them with measured values from a series of pumping tests carried out in large-diameter wells. The results show that this method is able to produce a reliable interpretation of the shallow hydrogeological context using information generally available in the region. The research contributes to improving the identification of areas where conditions are suitable for manual drilling. This is achieved by applying the described method, based on a structured and semi-quantitative approach, to classify the zones of suitability for given manual drilling techniques using data available in most African countries. Ultimately, this work will support proposed international programs aimed at promoting low-cost water supply in Africa and enhancing access to safe drinking water for the population.     

Space Weather: Physics,Effects and Predictability

The time varying conditions in the near-Earth space environment that may affect space-borne or ground-based technological systems and may endanger human health or life are referred to as space weather. Space weather effects arise from the dynamic and highly variable conditions in the geospace environment starting from explosive events on the Sun (solar flares), Coronal Mass Ejections near the Sun in the interplanetary medium, and various energetic effects in the magnetosphere–ionosphere–atmosphere system. As the utilization of space has become part of our everyday lives, and as our lives have become increasingly dependent on technological systems vulnerable to the space weather influences, the understanding and prediction of hazards posed by these active solar events have grown in importance. In this paper, we review the processes of the Sun–Earth interactions, the dynamic conditions within the magnetosphere, and the predictability of space weather effects on radio waves, satellites and ground-based technological systems today.     

Deep inside a neoproterozoic intra-oceanic arc: growth,differentiation and exhumation of the Amalaoulaou complex (Gourma,Mali)

We show here that the Amalaoulaou complex, in the Pan-African belt of West Africa (Gourma, Mali), corresponds to the lower and middle sections of a Neoproterozoic intra-oceanic arc. This complex records a 90–130-Ma-long evolution of magmatic inputs and differentiation above a subducting oceanic slab. Early c. 793 Ma-old metagabbros crystallised at lower crustal or uppermost mantle depths (25–30 km) and have geochemical characteristic of high-alumina basalts extracted from a depleted mantle source slightly enriched by slab-derived sedimentary components ((La/Sm)_N < 1; ε_Nd: +5.4–6.2; ⁸⁷Sr/⁸⁶Sr: 0.7027–0.7029). In response to crustal thickening, these mafic rocks were recrystallised into garnet-granulites (850–1,000°C; 10–12 kbar) and subject to local dehydration–melting reactions, forming trondhjemititic leucosomes with garnet–clinopyroxene–rutile residues. Slightly after the granulitic event, the arc root was subject to strong HT shearing during partial exhumation (detachment faults/rifting or thrusting), coeval with the emplacement of spinel- and garnet-pyroxenite dykes crystallised from a high-Mg andesitic parental magma. Quartz and hornblende-gabbros (700–660 Ma) with composition typical of hydrous volcanic rocks from mature arcs ((La/Sm)_N: 0.9–1.8; ε_Nd: +4.6 to +5.2; ⁸⁷Sr/⁸⁶Sr: 0.7028–0.7031) were subsequently emplaced at mid-arc crust levels (~15 km). Trace element and isotopic data indicate that magmas tapped a depleted mantle source significantly more enriched in oceanic sedimentary components (0.2%). Exhumation occurred either in two stages (700–660 and 623 Ma) or in one stage (623 Ma) with a final exhumation of the arc root along cold P-T path (550°C, 6–9 kbar; epidote–amphibolite and greenschist facies conditions) during the main Pan-African collision event (620–580 Ma). The composition of magmas forming the Cryogenian Amalaoulaou arc and the processes leading to intra-arc differentiation are strikingly comparable to those observed in the deep section of exposed Mezosoic oceanic arcs, namely the Kohistan and Talkeetna complex. This evolution of the Amalaoulaou oceanic arc and its accretion towards the West African craton belong to the life and closure of the Pharusian Ocean that eventually led to the formation of the Greater Gondwana supercontinent, a similar story having occurred on the other side of the Sahara with the Mozambique Ocean.     

Can Plant Competition and Diversity Reduce the Growth and Survival of Exotic <Emphasis Type="Italic">Phragmites australis</Emphasis> Invading a Tidal Marsh?

The rapid proliferation of Phragmites australis in North America has challenged resource managers to curb its expansion and reduce the loss of functional tidal marsh. We investigated whether native plant competition could reduce the ability of Phragmites to invade a tidal marsh, and if plant diversity (species richness, evenness, and composition) altered the competitive outcome. Immature Phragmites shoots and four native halophytes were transplanted to small but dense field plots (~1,200 shoots m⁻²) comprising three community structure types (Phragmites alone, Phragmites + 1 native species, and Phragmites + 4 native species). Interspecific competition significantly reduced Phragmites aboveground biomass, shoot length production, density, and survival by approximately 60%. Additionally, plots planted with greater native diversity contained Phragmites with the lowest growth and survival, potentially indicating diversity-enhanced resource competition. Competition consistently reduced the growth of Phragmites even under favorable conditions: lack of strong tidal flooding stresses as well as elevated nutrient pools.     

The Proterozoic,albitite-hosted,Valhalla uranium deposit,Queensland, Australia: a description of the alteration assemblage associated with uranium mineralisation in diamond drill hole V39

The Valhalla uranium deposit, located 40 km north of Mount Isa, Queensland, Australia, is an albitite-hosted, Mesoproterozoic U deposit similar to albitite-hosted uranium deposits in the Ukraine, Sweden, Brazil and Guyana. Uranium mineralisation is hosted by a thick package of interbedded fine-grained sandstones, arkoses and gritty siltstones that are bound by metabasalts belonging to the ca. 1,780 Ma Eastern Creek Volcanics in the Western Succession of the Mount Isa basin. Alteration associated with U mineralisation can be divided into an early, main and late stage. The early stage is dominated by laminated and intensely altered rock comprising albite, reibeckite, calcite, (titano)magnetite ± brannerite. The main stage of mineralisation is dominated by brecciated and intensely altered rocks that comprise laminated and intensely altered rock cemented by brannerite, apatite, (uranoan)-zircon, uraninite, anatase, albite, reibeckite, calcite and hematite. The late stage of mineralisation comprises uraninite, red hematite, dolomite, calcite, chlorite, quartz and Pb-, Fe-, Cu-sulfides. Brannerite has U–Pb and Pb–Pb ages that indicate formation between 1,555 and 1,510 Ma, with significant Pb loss evident at ca. 1,200 Ma, coincident with the assemblage of Rodinia. The oldest ages of the brannerite overlap with ⁴⁰Ar/³⁹Ar ages of 1,533 ± 9 Ma and 1,551 ± 7 Ma from early and main-stage reibeckite and are interpreted to represent the timing of formation of the deposit. These ages coincide with the timing of peak metamorphism in the Mount Isa area during the Isan Orogeny. Lithogeochemical assessment of whole rock data that includes mineralised and unmineralised samples from the greater Mount Isa district reveals that mineralisation involved the removal of K, Ba and Si and the addition of Na, Ca, U, V, Zr, P, Sr, F and Y. U/Th ratios indicate that the ore-forming fluid was oxidised, whereas the crystal chemistry of apatite and reibeckite within the ore zone suggests that F⁻ and were important ore-transporting complexes. δ¹⁸O values of co-existing calcite and reibeckite indicate that mineralisation occurred between 340 and 380°C and involved a fluid having δ¹⁸O_fluid values between 6.5 and 8.6‰. Reibeckite δD values reveal that the ore fluid had a δD_fluid value between −98 and −54‰. The mineral assemblages associated with early and main stages of alteration, plus δ¹⁸O_fluid and δD_fluid values, and timing of the U mineralisation are all very similar to those associated with Na–Ca alteration in the Eastern Succession of the Mount Isa basin, where a magmatic fluid is favoured for this style of alteration. However, isotopic data from Valhalla is also consistent with that from the nearby Mount Isa Cu deposit where a basinal brine is proposed for the transport of metals to the deposit. Based on the evidence to hand, the source fluids could have been derived from either or both the metasediments that underlie the Eastern Creek Volcanics or magmatism that is manifest in the Mount Isa area as small pegmatite dykes that intruded during the Isan Orogeny.