Age model, physical properties and paleoceanographic implications of the middle Pleistocene core sediments in the Choshi area, central Japan

Abstract   A continuous, well-preserved core was obtained from the Choshi area, on the Pacific side of Japan, to investigate paleoceanographic and paleoclimatic changes around the northwestern Pacific region during the middle Pleistocene. Siliciclastic sequences in the core are divided into five formations – the Obama, Yokone, Kurahashi and Toyosato Formations in the Inubo Group and the Katori Formation, in ascending order. Examination of calcareous nannofossils and magnetic polarities detected four datums in the core sediments of the Inubo Group: the top of Reticulofenestra asanoi , the base of Helicosphaera inversa , the top of Pseudoemiliania lacunosa and the Brunhes–Matuyama boundary. Fourteen marine isotope stages (MIS24–MIS11) were identified in the δ¹⁸O and δ¹³C records based on detected datums and the graphic correlation with the standard stack oxygen isotope curve. Magnetic susceptibility and gamma-ray attenuation porosity evaluator density were also measured and low values characterize the glacial intervals. Biogenic sedimentation by primary production may be larger during the glacial periods because of invasions of nutrient-rich northern surface-waters related to the southward shift of the Kuroshio front in the Choshi area.     

Geology, Mineralogy, and Formation Environment of the Disseminated Gold-Silver Telluride Bulawan Deposit, Negros Occidental, Philippines

Abstract: The Bulawan deposit is located in the porphyry copper belt of southwest Negros island, Philippines. Propylitic, K–feldspar, sericitic, and carbonate alteration types can be distinguished in the deposit. Propylite alteration occurs mainly in Cretaceous-Eocene andesitic lavas and agglomerates while K–feldspar, sericite and carbonate alteration types occur mostly in the Middle Miocene dacite porphyry breccia pipes and stocks which were intruded into the andesites. K-feldspar zones occur in the inner parts of the sericitized zone. Sericite alteration overprinted the propylitized and K-feldspar alteration zones, at lower temperature than epidote and chlorite in the propylitized zone. Carbonate alteration is associated with the mineralization in the center of the breccia pipes and along faults. Mineralization consists of gold-silver telluride ores that are hosted by the carbonate– and sericite-altered dacite porphyry breccia pipes. The Bulawan ores occur mainly as disseminations, but unlike many epithermal gold deposits, lack classical epithermal colloform and crustiform quartz veins. The ore minerals are sphalerite, galena, chalcopyrite, pyrite and tetrahedite-tennantite with minor amounts of electrum, calaverite, petzite, sylvanite, hessite, tellurobismuthite, coloradoite, altaite, and rucklidgeite. Electrum and telluride minerals are associated mostly with calcite and dolomite-ankerite minerals. Fluid inclusions in quartz and calcite in clasts of propylitized andesite in the breccia pipes homogenize from about 300° to 400°C while fluid inclusions in quartz, calcite and sphalerite within the dacite porphyry breccia pipes homogenize between 300° to 310°C. The ores were formed around 300°C from hydrothermal solutions with salinity of about 6. 6 wt % NaCl equivalent. The presence of sylvanite and calaverite as intergrowths with each other, and the Ag content of calaverite are consistent with the above temperature estimate. Based on paragenesis, the Bulawan deposit formed in a pyrite-stable environment, with pH between 3. 4 and 5. 5, f_O2 between 10^-32 to 10^-30 atm, f_S2 between 10^-9.8 to 10^-7.8 atm, f_Te2 between 10^-8.9 to 10^-6.5 atm, and total sulfur content about 10^-2.8 molal. The dominant reduced sulfur species in the ore solutions may have been H₂S_(aq), and the likely aqueous tellurium species were H₂Te_(aq) and H₂TeO_3(aq). The ore minerals in the Bulawan deposit were probably formed by mixing of slightly saline and low salinity fluids.     

Chaotic Itinerancy and Thermalization in a One-Dimensional Self-Gravitating System

Long time behavior of one-dimensional self-gravitating system is numerically investigated. Itinerant behavior among quasiequilibria is discovered, and its origin is clarified. The macroscopic distributions of the quasiquilibria are different from the isothermal distribution, but the transient states, which appears between two successive quasiequilibria, resemble the isothermal distribution. Averaging over the is itinerancy gives approximately the thermal equilibrium. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

Shock metamorphism on the moon

Shock metamorphism of the lunar samples is discussed. All types of lunar glasses formed by various-size collision-type impact are found as impact glass, ropy glass and agglutinates. The agglutinates bonded by crystal and glassy materials contain hydrogen and helium from the solar wind components. Lunar shocked minerals of plagioclase and silica show anomalous compositions and densities. There are typical two formation processes on planetary materials formed by shock events; that is (1) shocked quartz formed by silica-rich target rocks (esp. on evolved planets of the Earth and Mars), and (2) shocked silica with minor Al contents formed from plagioclase-rich primordial crusts of the Moon. The both shocked silica grows to coarse-grain normal crystals after high-temperature metamorphism which cannot distinguish the original main formation event of impact process.     

IOD and ENSO impacts on the extreme stream-flows of Citarum river in Indonesia

Extreme stream-flow events of Citarum River are derived from the daily stream-flows at the Nanjung gauge station. Those events are identified based on their persistently extreme flows for 6 or more days during boreal fall when the seasonal mean stream-flow starts peaking-up from the lowest seasonal flows of June–August. Most of the extreme events of high-streamflows were related to La Ni?a conditions of tropical Pacific. A few of them were also associated with the negative phases of IOD and the newly identified El Ni?o Modoki. Unlike the cases of extreme high streamflows, extreme low streamflow events are seen to be associated with the positive IODs. Nevertheless, it was also found that the low-stream-flow events related to positive IOD events were also associated with El Ni?o events except for one independent event of 1977. Because the occurrence season coincides the peak season of IOD, not only the picked extreme events are seen to fall under the IOD seasons but also there exists a statistically significant correlation of 0.51 between the seasonal IOD index and the seasonal streamflows. There also exists a significant lag correlation when IOD of June–August season leads the streamflows of September–November. A significant but lower correlation coefficient (0.39) is also found between the seasonal streamflow and El Ni?o for September–November season only.     

Probabilistic prediction of Indian summer monsoon rainfall using global climate models

Probabilistic seasonal predictions of rainfall that incorporate proper uncertainties are essential for climate risk management. In this study, three different multi-model ensemble (MME) approaches are used to generate probabilistic seasonal hindcasts of the Indian summer monsoon rainfall based on a set of eight global climate models for the 1982–2009 period. The three MME approaches differ in their calculation of spread of the forecast distribution, treated as a Gaussian, while all three use the simple multi-model subdivision average to define the mean of the forecast distribution. The first two approaches use the within-ensemble spread and error residuals of ensemble mean hindcasts, respectively, to compute the variance of the forecast distribution. The third approach makes use of the correlation between the ensemble mean hindcasts and the observations to define the spread using a signal-to-noise ratio. Hindcasts are verified against high-resolution gridded rainfall data from India Meteorological Department in terms of meteorological subdivision spatial averages. The use of correlation for calculating the spread provides better skill than the other two methods in terms of rank probability skill score. In order to further improve the skill, an additional method has been used to generate multi-model probabilistic predictions based on simple averaging of tercile category probabilities from individual models. It is also noted that when such a method is used, skill of probabilistic forecasts is improved as compared with using the multi-model ensemble mean to define the mean of the forecast distribution and then probabilities are estimated. However, skill of the probabilistic predictions of the Indian monsoon rainfall is too low.     

The role of the intra-daily SST variability in the Indian monsoon variability and monsoon-ENSO–IOD relationships in a global coupled model

The impact of diurnal SST coupling and vertical oceanic resolution on the simulation of the Indian Summer Monsoon (ISM) and its relationships with El Ni?o-Southern Oscillation (ENSO) and Indian Ocean Dipole (IOD) events are studied through the analysis of four integrations of a high resolution Coupled General Circulation Model (CGCM), but with different configurations. The only differences between the four integrations are the frequency of coupling between the ocean and atmosphere for the Sea Surface Temperature (SST) parameter (2 vs. 24?h coupling) and/or the vertical oceanic resolution (31 vs. 301 levels) in the CGCM. Although the summer mean tropical climate is reasonably well captured with all the configurations of the CGCM and is not significantly modified by changing the frequency of SST coupling from once to twelve per day, the ISM–ENSO teleconnections are rather poorly simulated in the two simulations in which SST is exchanged only once per day, independently of the vertical oceanic resolution used in the CGCM. Surprisingly, when 2?h SST coupling is implemented in the CGCM, the ISM–ENSO teleconnection is better simulated, particularly, the complex lead-lag relationships between the two phenomena, in which a weak ISM occurs during the developing phase of an El Ni?o event in the Pacific, are closely resembling the observed ones. Evidence is presented to show that these improvements are related to changes in the characteristics of the model’s El Ni?o which has a more realistic evolution in its developing and decaying phases, a stronger amplitude and a shift to lower frequencies when a 2-hourly SST coupling strategy is implemented without any significant changes in the basic state of the CGCM. As a consequence of these improvements in ENSO variability, the lead relationships between Indo-Pacific SSTs and ISM rainfall resemble the observed patterns more closely, the ISM–ENSO teleconnection is strengthened during boreal summer and ISM rainfall power spectrum is in better agreement with observations. On the other hand, the ISM–IOD teleconnection is sensitive to both SST coupling frequency and the vertical oceanic resolution, but increasing the vertical oceanic resolution is degrading the ISM–IOD teleconnection in the CGCM. These results highlight the need of a proper assessment of both temporal scale interactions and coupling strategies in order to improve current CGCMs. These results, which must be confirmed with other CGCMs, have also important implications for dynamical seasonal prediction systems or climate change projections of the monsoon.