Geochemical study of arsenic and other trace elements in groundwater and sediments of the Old Brahmaputra River Plain,Bangladesh

The geochemical study of groundwaters and core sediments from the Old Brahmaputra plain of Bangladesh was conducted to investigate the distribution of arsenic and related trace elements. Groundwaters from tube wells are characterized by pH of 6.4–7.4, dissolved oxygen (DO) of 0.8–1.8 mg/l, Ca contents of 5–50 mg/l, and Fe contents of 0.2–12.9 mg/l. Arsenic concentrations ranged from 8 to 251 μg/l, with an average value of 63 μg/l. A strong positive correlation exists between As and Fe (r ² = 0.802; p = 0.001) concentrations in groundwater. The stratigraphic sequences in the cores consist of yellowish silty clays at top, passing downward into grayish to yellowish clays and sands. The uppermost 3 m and lower parts (from 13 to 31 m) of the core sediments are oxidized (average oxidation reduction potential (ORP) +170 and +220 mV, respectively), and the ORP values gradually become negative from 3 to 13 m depths (−35 to −180 mV), indicating that anoxic conditions prevail in the shallow aquifers of the Brahmaputra plain. Age determinations suggest that clay horizons at ~10 m depth were deposited at around 2,000 and 5,000 years BP (¹⁴C ages) during the transgressive phase of sea-level change. Elevated concentrations of As, Pb, Zn, Cu, Ni, Cr, and V are present in the silts and clays, probably due to adsorption onto clay particles. Significant concentrations of As occur in black peat and peaty sediments at depths between 9 and 13 m. A strong positive correlation between As and Fe was found in the sediments, indicating As may be adsorbed onto Fe oxides in aquifer sediments.     

Transformation between surface spherical harmonic expansion of arbitrary high degree and order and double Fourier series on sphere

In order to accelerate the spherical harmonic synthesis and/or analysis of arbitrary function on the unit sphere, we developed a pair of procedures to transform between a truncated spherical harmonic expansion and the corresponding two-dimensional Fourier series. First, we obtained an analytic expression of the sine/cosine series coefficient of the \(4 \pi \) fully normalized associated Legendre function in terms of the rectangle values of the Wigner d function. Then, we elaborated the existing method to transform the coefficients of the surface spherical harmonic expansion to those of the double Fourier series so as to be capable with arbitrary high degree and order. Next, we created a new method to transform inversely a given double Fourier series to the corresponding surface spherical harmonic expansion. The key of the new method is a couple of new recurrence formulas to compute the inverse transformation coefficients: a decreasing-order, fixed-degree, and fixed-wavenumber three-term formula for general terms, and an increasing-degree-and-order and fixed-wavenumber two-term formula for diagonal terms. Meanwhile, the two seed values are analytically prepared. Both of the forward and inverse transformation procedures are confirmed to be sufficiently accurate and applicable to an extremely high degree/order/wavenumber as \(2^{30}\,{\approx }\,10^9\). The developed procedures will be useful not only in the synthesis and analysis of the spherical harmonic expansion of arbitrary high degree and order, but also in the evaluation of the derivatives and integrals of the spherical harmonic expansion.     

A method solving kepler's equation without transcendental function evaluations

We developed two approximations of the Newton-Raphson method. The one is a sort of discretization, namely to search an approximate solution on pre-specified grid points. The other is a Taylor series expansion. A combination of these was applied to solving Kepler's equation for the elliptic case. The resulting method requires no evaluation of transcendental functions. Numerical measurements showed that, in the case of Intel Pentium processor, the new method is three times as fast as the original Newton-Raphson method. Also it is more than 2.5 times as fast as Halley's method, Nijenhuis's method, and others.     

Polar magnetic substorms

From the world distribution of geomagnetic disturbance, the connection between the electric current in the ionosphere, the field-aligned current and asymmetric equatorial ringcurrent in the magnetosphere is discussed. The partial ring-current in the afternoon-evening region, whose intensity is closely correlated with the AE-index, usually develops and decays earlier than the symmetric ring-current in the course of magnetic storms. The partial ringcurrent seems to have a direct connection with the positive geomagnetic bay in high latitudes in the evening hours through the ionizing effect of the particles leaking from the partial ringcurrent. The dawn-to-dusk electric field in the magnetospheric tail is transferred to the polar ionosphere, producing there the twin vortex Hall current responsible for polar cap geomagnetic variation. The magnetic effect of the associated Pedersen current in the ionosphere is shown to be small but still worth considering. The electrojet near midnight along the auroral oval is thought to appear when the electric conductivity of the ionosphere is locally increased under the presence of large scale dawn-to-dusk electric field. The occasional appearance of a localized abnormal geomagnetic disturbance with reversed direction near the geomagnetic pole seems to suggest the occasional reversal of electric field near the outer surface of the magnetospheric tail, especially when the interplanetary magnetic field is northward.     

Calibration and validation of the ocean color version-3 product from ADEOS OCTS

We present calibration and validation results of the OCTS’s ocean color version-3 product, which mainly consists of the chlorophyll-a concentration (Chl-a) and the normalized water-leaving radiance (_nL_w). First, OCTS was calibrated for the inter-detector sensitivity difference, offset, and absolute sensitivity using external calibration source. It was also vicariously calibrated using in-situ measurements for water (Chl-a and_nL_w) and atmosphere (optical thickness), which were acquired synchronously with OCTS under cloud-free conditions. Second, the product was validated using selected 17 in-situ Chl-a and 11 in-situ_nL_w measurements. We confirmed that Chl-a was estimated with an accuracy of 68% for Chl-a less than 2 mg/m³, and_nL_w from 94% (band 2) to 128% (band 4). Geometric accuracy was improved to 1.3 km. Stripes were significantly reduced by modifying the detector normalization factor as a function of input radiance.     

Correction of atmospheric effect on ADEOS/OCTS ocean color data: Algorithm description and evaluation of its performance

This paper first describes the atmospheric correction algorithm for OCTS visible band data used at NASDA/EOC. Sharing a basic structure with Gordon and Wang’s Sea WiFS algorithm, it uses 10 candidate aerosol models including the “Asian dust model” introduced in consideration of the unique feature of aerosols over the east Asian waters. Based on the observations at 670 and 865 nm bands, the algorithm selects a pair of aerosol models that account best for the observed spectral reflectances, and synthesizes the aerosol reflectance used for the atmospheric correction. Two different schemes for determining the value of the parameter for the aerosol model selection are presented and their anticipated estimation error is analyzed in terms of retrieved water reflectance at 443 nm. The results of our numerical simulation show that the standard deviation of the estimation error of the “weighted average” scheme is mostly within the permissible level of ±0.002, reducing the error by 18% on average compared to the “simple average” scheme. The paper further discusses the expected error under the old CZCS-type atmospheric correction, which assumes constant aerosol optical properties throughout the given image. Although our algorithm has a better performance than the CZCS algorithm, further analysis shows that the error induced by the assumption taken in the algorithm that the water-leaving radiance at 670 nm band is negligibly small may be large in high pigment concentration waters, indicating the necessity for future improvements.     

Genesis of the Kamioka Skarn Deposits: An Important Role of Clinopyroxene Skarn and Graphite-bearing Limestone in Precipitating Sulfide Ore

Abstract: A genetical relationship between skarn formation and mineralization is investigated for the Kamioka skarn deposits which are the largest Zn-Pb producer in Japan. In the Mozumi deposit, one of main deposits in the Kamioka mining area as well as Tochibora and Maruyama, clinopyroxene skarn was generally subjected to later replacement by garnet or magnetite–calcite–quartz during the Zn-Pb mineralization. The replacement of hedenbergitic clinopyroxene by andraditic garnet resulted in the formation of diopsidic clinopyroxene relicts. With the progress of replacement, the S/So value (So: an estimated area occupied by an original clinopyroxene grain in a thin section, S: a total area of relict clinopyroxene fragments) which is an index of the degree of replacement decreases from 0. 7 to 0. 1, and the hedenbergite mole percent of relict clinopyroxene decreases drastically from about 65 to less than 40. A close association of andraditic garnet and sphalerite suggests that heden-bergitic clinopyroxene skarn played an important role to reduce the relatively oxic ore-forming fluid enriched in Zn²⁺ and SO4^2– and to precipitate sphalerite from the fluid. Ferrous iron in the hedenbergitic clinopyroxene skarn was oxidized to form andraditic garnet. Besides this garnet formation, the mineral assemblage of magnetite–calcite–quartz replaced the clinopyroxene skarn at the time of mineralization. In both cases, the reduction of relatively oxic ore-forming fluid by hedenbergitic clinopy-roxene skarn at the later stage brought about the precipitation of sulfide minerals. In contrast, these types of later replacement are not found in the Tochibora deposit. Instead, graphite-bearing crystalline limestone and relatively fresh clinopyroxene skarn are common. Mineralized clinopyroxene skarn has high graphite carbon contents relative to barren one, suggesting that the amount of graphite in the skarn was an important controlling factor for mineralization. It is very likely that the graphite played a role of reducing agent during the mineralization in the Tochibora deposit.     

Application of common reflection angle migration for imaging deformation structures in an inner accretionary wedge,Nankai Trough,Japan

To better image deformation structures within the inner accretionary wedge of the Nankai Trough, Japan, we apply common reflection angle migration to a legacy two-dimensional seismic data set acquired with a 6 km streamer cable. In this region, many seismic surveys have been conducted to study the seismogenic zone related to plate subduction. However, the details of the accreted sediments beneath the Kumano forearc basin are still unclear due to the poor quality of seismic images caused by multiple reflections, highly attenuated signals, and possibly complex geological structures. Generating common image gathers in the subsurface local angle domain rather than the surface offset domain is more advantageous for imaging geological structures that involve complex wave paths and poor illumination. By applying this method, previously unseen structures are revealed in the thick accreted sediments. The newly imaged geometric features of reflectors, such as the folds in the shallow part of the section and the deep reflectors with stepwise discontinuities, imply deformation structures with multiple thrust faults. The reflections within the deep accreted sediments (approximately 5 km) are mainly mapped to far angles (30°–50°) in the common reflection angles, which correspond to the recorded offset distances greater than 4.5 km. This result indicates that the far offset/angle information is critical to image the deformation structures at depth. The new depth image from the common reflection angle migration provides seismic evidence of multiple thrust faults and their relationship with the megathrust fault that is essential for understanding the structure and evolution of the Nankai Trough seismogenic zone.     

Simulation analyses of buildings damaged in the 1995 Kobe,Japan, earthquake considering soil–structure interaction

A series of studies was conducted on three buildings of steel reinforced concrete structures with RC shear walls damaged in the 1995 Hyogo-ken Nanbu earthquake. These buildings are located in an area where structural damage centred around. Two of these buildings suffered severe damage, while the third was not structurally damaged. Our studies deal with site inspections, including micro-tremor measurement of buildings, the evaluation of input motions, and the response analyses considering soil–structure interaction. The results of simulation analyses of the two severely damaged buildings correspond to their actual damage state. From the response analyses of the one slender building with no structural damage, it was concluded that uplifting is the main reason it did not suffer any structural damage. Through these studies, the importance of soil–structure interaction and effective input motion is fully understood. Copyright © 1999 John Wiley & Sons, Ltd.