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.     

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.     

Detection of high-energy solar neutrons and protons by ground level detectors on April 15, 2001

In association with the large solar flare of April 15, 2001, the Chacaltaya neutron monitor observed a 3.6σ enhancement of the counting rate between 13:51 and 14:15 UT. Since the enhancement was observed beginning 11 min before the GLE, solar neutrons must be involved in this enhancement. The integral energy spectrum of solar neutrons can be expressed by a simple power law in energy with the index γ=-3.0±1.0. On the other hand, an integral energy spectrum of solar protons has been obtained in the energy range between 650 MeV and 12 GeV. The spectrum can also be expressed by a power law with the power index γ=-2.75±0.15. The flux of solar protons observed at Chacaltaya (at 12 GeV) was already one order less than the flux of the galactic cosmic rays. It may be the first simultaneous observation of the energy spectra of both high-energy protons and neutrons. Comparing the Yohkoh soft X-ray telescope images with the observed particle time profiles, an interesting picture of the particle acceleration mechanism has been deduced.     

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.     

X-ray observation of the Jovian impacts of comet Shoemaker - Levy 9

An ASCA observation of the Jovian impact of the comet Shoemaker-Levy 9 is reported. Four impacts of H, L, Q1 and R were observed and four impacts of B, C, G, and Q2 were observed within 60 minutes after their impacts. No significant flaring of X-ray emission was observed. Upper limit X-ray fluxes of 90 % confidence level, averaged 5 minutes just after the impacts, were 2.4 × 10^–13 erg sec^–1 cm^–2, 3.5 × 10^–13 erg sec^–1 cm^–2, 1.6 × 10^–13 erg sec^–1 cm^–2 and 2.9 × 10^–13 erg sec^–1 cm^–2 for the impacts of H, L, Q1 and R, respectively, in the 0.5(0.7 for H and Q1)–10 keV energy range. However, a hint of X-ray enhancement around Jupiter from July 17 to July 19 was detected with about 2 6 × 10^–14 erg sec^–1 cm^–2 in the 0.5–10 keV energy range.     

Impact of intra-daily SST variability on ENSO characteristics in a coupled model

This paper explores the impact of intra-daily Sea Surface Temperature (SST) variability on the tropical large-scale climate variability and differentiates it from the response of the system to the forcing of the solar diurnal cycle. Our methodology is based on a set of numerical experiments based on a fully global coupled ocean–atmosphere general circulation in which we alter (1) the frequency at which the atmosphere sees the SST variations and (2) the amplitude of the SST diurnal cycle. Our results highlight the complexity of the scale interactions existing between the intra-daily and inter-annual variability of the tropical climate system. Neglecting the SST intra-daily variability results, in our CGCM, to a systematic decrease of 15% of El Ni?o—Southern Oscillation (ENSO) amplitude. Furthermore, ENSO frequency and skewness are also significantly modified and are in better agreement with observations when SST intra-daily variability is directly taken into account in the coupling interface of our CGCM. These significant modifications of the SST interannual variability are not associated with any remarkable changes in the mean state or the seasonal variability. They can therefore not be explained by a rectification of the mean state as usually advocated in recent studies focusing on the diurnal cycle and its impact. Furthermore, we demonstrate that SST high frequency coupling is systematically associated with a strengthening of the air-sea feedbacks involved in ENSO physics: SST/sea level pressure (or Bjerknes) feedback, zonal wind/heat content (or Wyrtki) feedback, but also negative surface heat flux feedbacks. In our model, nearly all these results (excepted for SST skewness) are independent of the amplitude of the SST diurnal cycle suggesting that the systematic deterioration of the air-sea coupling by a daily exchange of SST information is cascading toward the major mode of tropical variability, i.e. ENSO.     

Decadal seesaw of the Central and Subtropical Mode Water formation associated with the Kuroshio Extension variability

Available Argo profiling float data from 2002 to 2011 were analyzed to examine the effect of the Kuroshio Extension (KE) current system variability on the formation of the Central Mode Water. Just north of the upstream portion of the KE at 140–152°E, formation of a lighter variety of the Central Mode Water in winter was active during the unstable period of the upstream KE in 2006–2009 and was reduced when the upstream KE was in the stable period of 2002–2005 and 2010–2011. This decadal formation variability is out of phase with that of the Subtropical Mode Water just south of the KE.