Cokriging nonstationary data

Universal cokriging is used to obtain predictions when dealing with multivariate random functions. An important type of nonstationarity is defined in terms of multivariate random functions with increments which are stationary of orderk. The covariance between increments of different variables is modeled by means of the pseudo-cross-covariance function. Criteria are formulated to which the parameters of pseudo-cross-covariance functions must comply so as to ensure positive-definiteness. Cokriging equations and the induced cokriging equations are given. The study is illustrated by an example from soil science.     

An Example of Island-Arc Petrogenesis: Geochemistry and Petrology of the Southern Luzon Arc, Philippines

Volcanism throughout the Luzon arc is associated with eastwardsubduction of the South China Sea floor along the Manila Trench.The southern section of the arc, the focus of this study, extendsfrom the Lingayen-Dingalan fault to the small islands just southof Luzon. Two segments appear to exist along this section ofthe arc the northern Bataan and southern Mindoro segments whichare separated by the Macolod Corridor. The volcanic rocks have typical arc phenocryst mineralogies:olivine, clinopyroxene, plagioclase, and titanomagnetite inthe most mafic rocks and clinopyroxene, plagioclase, orthopyroxene,titanomagnetite, ? amphibole in the more felsic samples. Complexzoning, sieve textures, and decoupling of incompatible traceelements suggest that processes such as assimilation have takenplace. The rocks from the study area range from basalts to rhyolitesand show typical calc-alkaline features. The rocks of the MacolodCorridor and Mindoro segment are particularly enriched in largeion lithophile elements (LILE), light rare earth elements (LREE),and radiogenic Sr compared with the Bataan segment. The datafall within the mantle array on Sr-Nd isotopic diagrams andgrade toward higher Sr and lower Nd isotopic values from northto south. A likely source for the volcanics of this study is either amid-ocean ridge basalt (MORB)-type mantle that undergoes higherdegrees of partial melting than regions involved in MORB generationor a previously depleted source. We suggest that the high fieldstrength element (HFSE) anomalies have been derived throughdifferential element partitioning during fluid transport fromthe subducted lithosphere to the mantle wedge. Continental crustal material seems to play a significant roleparticularly in the Macolod Corridor and the Mindoro segment,based on the high LILE, La/Sm ratios, radiogenic Sr isotopes,and ¹⁸O values. The Macolod Corridor and the Mindoro segmenthave undergone source contamination by crustal material fromthe North Palawan-Mindoro crustal block either during the collisionof this block with the Manila Trench or by subduction of sedimentsrich in this crustal material. A similar component has alsobeen detected in the Bataan segment but in minor amounts. Thetrace element and isotopic differences between the northernand southern sections of the arc are interpreted in terms ofvariable composition (i. e., variable amounts of a crustal componentintroduced from the Palawan-Mindoro crustal terrain) of themetasomatic fluids released into the source.     

Archean high-Mg monzodiorite–syenite,epidote skarn,and biotite–sericite gold lodes in the Granny Smith–Wallaby district,Australia: U–Pb and Re–Os chronometry of two intrusion-related hydrothermal systems

The Granny Smith (37 t Au production) and Wallaby deposits (38 t out of a 180 t Au resource) are located northeast of Kalgoorlie, in 2.7 Ga greenstones of the Eastern Goldfields Province, the youngest orogenic belt of the Yilgarn craton, Western Australia. At Granny Smith, a zoned monzodiorite–granodiorite stock, dated by a concordant titanite–zircon U–Pb age of 2,665 ± 3 Ma, cuts across east-dipping thrust faults. The stock is fractured but not displaced and sets a minimum age for large-scale (1 km) thrust faulting (D2), regional folding (D1), and dynamothermal metamorphism in the mining district. The local gold–pyrite mineralization, controlled by fractured fault zones, is younger than 2,665 ± 3 Ma. In augite–hornblende monzodiorite, alteration progressed from a hematite-stained alkali feldspar–quartz–calcite assemblage and quartz–molybdenite–pyrite veins to a late reduced sericite–dolomite–albite assemblage. Gold-related monazite and xenotime define a U–Pb age of 2,660 ± 5 Ma, and molybdenite from veins a Re–Os isochron age of 2,661 ± 6 Ma, indicating that mineralization took place shortly after the emplacement of the main stock, perhaps coincident with the intrusion of late alkali granite dikes. At Wallaby, a NE-trending swarm of porphyry dikes comprising augite monzonite, monzodiorite, and minor kersantite intrudes folded and thrust-faulted molasse. The conglomerate and the dikes are overprinted by barren (<0.01 g/t Au) anhydrite-bearing epidote–actinolite–calcite skarn, forming a 600-m-wide and >1,600-m-long replacement pipe, which is intruded by a younger ring dike of syenite porphyry pervasively altered to muscovite + calcite + pyrite. Skarn and syenite are cut by pink biotite–calcite veins, containing magnetite + pyrite and subeconomic gold–silver mineralization (Au/Ag = 0.2). The veins are associated with red biotite–sericite–calcite–albite alteration in adjacent monzonite dikes. Structural relations and the concordant titanite U–Pb age of the skarn constrain intrusion-related mineralization to 2,662 ± 3 Ma. The main-stage gold–pyrite ore (Au/Ag >10) forms hematite-stained sericite–dolomite–albite lodes in stacked D2 reverse faults, which offset skarn, syenite, and the biotite–calcite veins by up to 25 m. The molybdenite Re–Os age (2,661 ± 10 Ma) of the ore suggests a genetic link to intrusive activity but is in apparent conflict with a monazite–xenotime U–Pb age (2,651 ± 6 Ma), which differs from that of the skarn at the 95% confidence level. The time relationships at both gold deposits are inconsistent with orogenic models invoking a principal role for metamorphic fluids released during the main phase of compression in the fold belt. Instead, mineralization is related in space and time to late-orogenic, magnetite-series, high-Mg monzodiorite–syenite intrusions of mantle origin, characterized by Mg/(Mg + Fe_TOTAL) = 0.31–0.57, high Cr (34–96 ppm), Ni (22–63 ppm), Ba (1,056–2,321 ppm), Sr (1,268–2,457 ppm), Th (15–36 ppm), and rare earth elements (total REE: 343–523 ppm). At Wallaby, shared Ca–K–CO₂ metasomatism and Th-REE enrichment (in allanite) link Au–Ag mineralization in biotite–calcite veins to the formation of the giant epidote skarn, implicating a Th + REE-rich syenite pluton at depth as the source of the oxidized hydrothermal fluid. At Granny Smith, lead isotope data and the Rb–Th–U signature of early hematite-bearing wall-rock alteration point to fluid released by the source pluton of the differentiated alkali granite dikes.     

Modeling of water–rock interactions in an aquitard of sandy clay in the coastal area near Beihai,China

Unconsolidated sand, gravel and clay deposits near Beihai and in the Leizhou Peninsula in southern China form an unconfined aquifer, aquitard and a confined aquifer. Water and soil samples were collected from the two aquifers in the coastal Beihai area for the determination of chemical compositions, minerals and soluble ions. Hydrogeochemical modeling of three flow paths through the aquitard are carried out using PHREEQC to determine water–rock interactions along the flow paths. The results indicate that the dissolution of anorthite, fluorite, halite, rhodochrosite and CO₂, and precipitation of potash feldspar and kaolinite may be occurring when groundwater leaks through the aquitard from the unconfined aquifer to the confined aquifer. Cation exchanges between Na and Ca can also happen along the flow paths.     

Diel Patterns in Chlorophyll <Emphasis Type="Italic">a</Emphasis> Specific Absorption Coefficient and Absorption Efficiency Factor of Picoplankton

Diel patterns in the chlorophyll a specific absorption coefficient of surface picoplankton, a*_pico (γ) (m²·[mg chlorophyll a]⁻¹), were studied at 7 stations under daily cycle of in situ light condition in the western subarctic Pacific and Japan Sea. All the data were normalized by dividing the anomaly with daily averaged a*_pico (γ). Opposite diel patterns were observed for the normalized a*_pico (443) and a*_pico (675) with maximum toward dawn or dusk and minimum toward midday at 4 stations under low-irradiance (LI) conditions and vice versa at 3 stations under high-irradiance (HI) conditions. The absorption efficiency factors at red absorption peak, Q _a (675), were determined by reconstruction with intracellular chlorophyll a concentration and cell diameter. The normalized Q _a (675) also showed diel pattern with maximum toward midday and minimum toward dawn or dusk under LI. The diel pattern in a*_pico (675) and Q _a (675) were primarily caused by changes in intracellular chlorophyll a concentration due to photoadaptation under LI. The diel pattern in a*_pico (443) was influenced by pigmentation, as recognized by blue to red ratio [a*_pico (443)/a*_pico (675)] under HI. This study proposed that the opposite diel pattern in a*_pico (γ) might occur for a wide range of algal species. The results presented here have important consequences for the interpretation of diel variations in optical properties observed in the open ocean.     

CO<Subscript>2</Subscript> Flux Measurements over the Sea Surface by Eddy Correlation and Aerodynamic Techniques

The measurements of the vertical transport of CO₂ were carried out over the Sea of Japan using the specially designed pier of Kyoto University on September 20 to 22, 2000. CO₂ fluxes were measured by the eddy correlation and aerodynamic techniques. Both techniques showed comparable CO₂ fluxes during sea breeze conditions: −0.001 to −0.08 mg m⁻²s⁻¹ with the mean of −0.05 mg m⁻²s⁻¹. This means that the measuring site satisfies the fetch requirement for meteorological observations under sea breeze conditions. Moreover, the eddy diffusivity coefficient used in the aerodynamic technique is found to be consistent with the coefficient used in the eddy correlation technique. The present result leads us to conclude that the aerodynamic technique may be applicable to underway CO₂ flux measurements over the ocean and may be used in place of the bulk technique. The important point is the need to maintain a measuring accuracy of CO₂ concentration difference of the order of 0.1 ppmv on the research vessels or the buoys.     

Roles of Biogeochemical Processes in the Oceanic Carbon Cycle Described with a Simple Coupled Physical-Biogeochemical Model

To evaluate the contribution of biogeochemical processes to the oceanic carbon cycle and to calculate the ratio of calcium carbonate to organic carbon downward export, we have incorporated biological and alkalinity pumps in the yoked high-latitude exchange/interior diffusion-advection (YOLDA) model. The biogeochemical processes are represented by four parameters. The values of the parameters are tuned so that the model can reproduce the observed phosphate and alkalinity distributions in each oceanic region. The sensitivity of the model to the biogeochemical parameters shows that biological production rates in the euphotic zone and decomposition depths of particulate matters significantly influence horizontal and vertical distributions of biogeochemical substances. The modeled vertical fluxes of particulate organic phosphorus and calcium carbonate are converted to vertical carbon fluxes by the biological pump and the alkalinity pump, respectively. The downward carbon flux from the surface layer to the deep layer in the entire region is estimated to be 3.36 PgC/yr, which consists of 2.93 PgC/yr from the biological pump and 0.43 PgC/yr from the alkalinity pump, which is consistent with previous studies. The modeled rain ratio is higher with depth and higher in the Pacific and Indian Oceans than in the Atlantic Ocean. The global rain ratio at the surface layer is calculated to be 0.14 to 0.15. This value lies between the lower and higher ends of the previous estimates, which range widely from 0.05 to 0.25. This study indicates that the rain ratio is unlikely to be higher than 0.15, at least in the surface waters.