Shuzo Takemoto

Simulation of the Indonesian land gravity data using a digital terrain model data

The Indonesian gravity field is neither accurately nor comprehensively determined, especially due to inadequacy of land gravity data. This study deals with determination of Indonesian land gravity and proposes the solution to data unavailability by means of a simulation technique. The simulation was carried out by combining short wavelength topographic effects from GTOPO30 and long wavelength information from EGM96. The simulated result was then compared with the observed gravity data. Over Java, Sumatra and Sulawesi islands, using three methods commonly used on the computation of topographic effect; topography, isostatic and RTM (Residual Terrain Model), it was estimated that error propagation by the GTOPO30 into the simulated gravity is about 4.5 to 11.7 mgal, with the RTM method was affected less than others. It was also shown that the simulated gravity from the RTM method gave the best agreement with STD (Standard Deviation) differences of 17 to 42 mgal compared to the observed data. This result was achieved after applying optimal RTM parameters over the Indonesian area: a reference field of 25′–27.5′ and density of 2–2.2 gr/cm3. Compared to STD differences between EGM96 and observed data, that between the simulated gravity and observed data improved by 2.5–7 mgal, and gave more detailed gravity features, especially over areas of high topography.

Crustal deformations associated with the great Sumatra-Andaman earthquake deduced from continuous GPS observation

We analyzed continuous GPS data from more than 20 sites in Asia, Australia and islands in Indian Ocean in order to detect crustal deformations associated with the Sumatra-Andaman earthquake of December 26, 2004. Coseismic steps can be recognized at sites about 3,000 km away from the epicenter such as Kunming in south China, Quezon in Philippines, and Diego Garcia Island in central Indian Ocean. The largest displacement of about 26 cm is found at Phuket in Thailand about 600 km away from the epicenter, about twice as large as that at Sampari, the nearest site in northern Sumatra. These observations suggest that as large slip as 14 m occurred beneath the Nicobar Islands. Large postseismic displacements are observed at Phuket and Sampari after the mainshock, but the former is three times larger than the latter. This suggests that the spatial distribution of afterslip is different from the coseismic slip distribution. The temporal variation of postseismic displacements can be explained by a logarithmic function derived from rate-state dependent friction law with short characteristic time. The area where coseismic displacements from the Nias earthquake of March 28, 2005 are detected is much smaller than that from the December mainshock, but displacement at Sampari is larger than that during the mainshock. These displacements suggest less than 4 m slip on a shallow dipping thrust fault and resultant moment release is smaller than that estimated from seismological data. Finally, total moment released by afterslip amounts to 3.83 × 1022 Nm which is equivalent to Mw 8.99 for about five months, including the afterslip for the Nias earthquake.

The CLIO project

The CLIO project including a 100 m baseline cryogenic gravitational wave laser interferometer and a 100 m baseline geophysical strain meter was conducted in the Kamioka mine in Japan to investigate the technical feasibility of the large-scale cryogenic gravitational wave telescope (LCGT), which is planned to be constructed in the same Kamioka mine with 30 times longer baseline than CLIO, and to demonstrate the collaborative operation between these instruments about long-term continuous operation and gravitational wave signal veto analysis. About the cryogenic gravitational wave interferometer, the whole vacuum system and four cryostats, which house and cool sapphire mirrors, were constructed, and the required vacuum level of 10 −6 mbar and the temperature of 8 K at the inner radiation shield in the cryostat were achieved. About the geophysical strain meter, the obtained geophysical strain in the Kamioka mine was successfully simulated with a finite element model with a good agreement with less than 5% error. The strain meter also verified a permanent ground step change of micrometre order due to some earthquakes. We present the recent progress about the CLIO project.

Accurate estimation of atmospheric effects on gravity observations made with a superconducting gravity meter at Kyoto

Abstract Atmospheric effects on gravity observations at Kyoto were estimated by using meteorological data sets at an interval of 12 h during a 4 month period from July to October in 1993. The effects owing to the air mass near the gravity station were evaluated by numerical integrals, and those distant from the station were calculated by using spherical harmonic expansions of meteorological data. The error in the calculated atmospheric effects was of the order of 0.1 μgal at most, except for the error related to the response of the oceans near the station based on the inverted barometric loading mode. About 90% of the atmospheric effects were attributed to local atmospheric variations within 50 km of the station. The remaining effects owing to the air mass outside this zone were of the order of 1 μgal, in which different features were recognized as compared with the effects owing to the regional air mass around Kyoto. The atmospheric effects thus estimated were compared with gravity data obtained by a superconducting gravity meter at Kyoto. The residuals showed gravity changes of a few microgals, a part of which might be caused by sources such as variations of the ground water level around the station.

Restoration of the depression structure at the eastern part of central Kyushu, Japan by means of dislocation modeling

Abstract We have attempted to restore a subsurface structure in the eastern part of central Kyushu, Japan, by combining fault motions which were modeled as dislocation planes embedded in an elastic isotropic half space. The simulated crustal deformation pattern was compared with the subsurface structures estimated from gravity anomalies and/or seismic prospecting. The modeling procedure successfully restored all tectonic basins in the area without any need for motive forces for uplift or subsidence. The results also suggest that two major tectonic events have occurred in this region. Those are (1) the formation of half-graben caused by north–south extension, and (2) the formation of the pull-apart basin caused by east–west compression.

Status of the CLIO project

The CLIO project involves the Cryogenic Laser Interferometer Observatory (CLIO) detector complex for gravitational wave detection and the Kamioka Laser Interferometric Strainmeter for the acquisition of geophysical data. CLIO has been constructed to demonstrate the feasibility of a future project, the Large-scale Cryogenic Gravitational wave Telescope (LCGT). It will utilize the low seismic and stable environment of the Kamioka mine as well as sapphire mirrors and suspension fibres at low temperature to reduce thermal noise. We designed CLIO to have a noise level limited by the thermal noise of sapphire mirrors and sapphire suspension fibres, which vary from 3 × 10−19 m Hz−1/2 at 300 K to 2 × 10−20 m Hz−1/2 at 20 K around 100 Hz. The strainmeter has already succeeded in monitoring the Earths tidal motion with a strain sensitivity of 2 × 10−12. The seismic noise veto between these same-scale interferometers is expected to provide an effective means of data selection for the gravitational wave signal analysis, and the ground motion data obtained by the strainmeter will help to maintain the stable operation of CLIO.

Design and construction status of CLIO

Construction of CLIO (cryogenic laser interferometer observatory) with 100 m baseline length has begun in the Kamioka mine. The tunnel for CLIO has been dug and infrastructure work is now in progress. CLIO is the final step to LCGT (large scale cryogenic gravitational wave telescope) and the first practical construction of a cryogenic interferometer in the world. The objective of CLIO is to demonstrate two of three features of LCGT, which are to utilize the quietness and stable environment of the underground site and to adopt cryogenic sapphire mirrors for thermal noise reduction. Also, it is a joint project by gravitational wave and geophysics researchers. CLIO has a locked Fabry–Perot configuration equipped with ring mode cleaners and cryocoolers to cool the sapphire mirrors to 20 K. The noise level of CLIO is designed to trace the thermoelastic noise of sapphire mirrors which varies from 10−18 m Hz−1/2 at 300 K to 10−19 m Hz−1/2 at 20 K around 100 Hz. A 7 m single-arm cryogenic test facility has been built at ICRR (Institute for Cosmic Ray Research), while the 20 m room temperature interferometer is in operation at Kamioka. Technical knowledge developed by these prototypes will be leveraged to realize CLIO.

Current status of the CLIO project

CLIO (Cryogenic Laser Interferometer Observatory) is a Japanese gravitational wave detector project. One of the main purposes of CLIO is to demonstrate thermal-noise suppression by cooling mirrors for a future Japanese project, LCGT (Large-scale Cryogenic Gravitational Telescope). The CLIO site is in Kamioka mine, as is LCGT. The progress of CLIO between 2005 and 2007 (room- and cryogenic-temperature experiments) is introduced in this article. In a room-temperature experiment, we made efforts to improve the sensitivity. The current best sensitivity at 300 K is about 6 × 10-21/√Hz around 400 Hz. Below 20 Hz, the strain (not displacement) sensitivity is comparable to that of LIGO, although the baselines of CLIO are 40-times shorter (CLIO: 100m, LIGO: 4km). This is because seismic noise is extremely small in Kamioka mine. We operated the interferometer at room temperature for gravitational wave observations. We obtained 86 hours of data. In the cryogenic experiment, it was confirmed that the mirrors were sufficiently cooled (14 K). However, we found that the radiation shield ducts transferred 300K radiation into the cryostat more effectively than we had expected. We observed that noise caused by pure aluminum wires to suspend a mirror was suppressed by cooling the mirror.

Direct polynomial approach to nonlinear distance (ranging) problems

In GPS atmospheric sounding, geodetic positioning, robotics and photogrammetric (perspective center and intersection) problems, distances (ranges) as observables play a key role in determining the unknown parameters. The measured distances (ranges) are however normally related to the desired parameters via nonlinear equations or nonlinear system of equations that require explicit or exact solutions. Procedures for solving such equations are either normally iterative, and thus require linearization or the existing analytical procedures require laborious forward and backward substitutions. We present in the present contribution direct procedures for solving distance nonlinear system of equations without linearization, iteration, forward and backward substitution. In particular, we exploit the advantage of faster computers with large storage capacities and the computer algebraic softwares of Mathematica, Maple and Matlab to test polynomial based approaches. These polynomial (algebraic based) approaches turn out to be the key to solving distance nonlinear system of equations. The algebraic techniques discussed here does not however solve all general types of nonlinear equations but only those nonlinear system of equations that can be converted into algebraic (polynomial) form.

Understanding the concept of outlier and its relevance to the assessment of data quality: Probabilistic background theory

In recent years an increasing interest in the studies on outlier can be observed, however, for the time being there exists no general definition of outlier. In the present paper we introduced a generic descriptive definition of outlier. We observed that the outlier problems had so far been treated in statistical way without paying proper attention to probabilistic-theoretic backgrounds. In view of this gap, we made an attempt to establish a probabilistic background theory. Within this framework, the large deviations are considered as probabilistic-theoretic model of outlier, and the interrelationship of the laws of large numbers, the central limit theorems and the large deviations are clarified. These considerations are specialized for the case of statistical sample, which is important from the point of view of the assessment of data quality. Some methodological and historical aspects of geodesy, geophysics and astronomy are mentioned, too. We revealed that the data analysis carried out by Kepler in the process of discovery of his famous elliptic law of planetary motion has relevance to the outlier problem. This methodologically interesting fact is a new result in the history of geosciences. We established that the accuracy of Chebyshev inequality increases as the deviation of the random variable involved from its expectation, increases. The possibility of application of Chebyshev inequality to the outlier problem is pointed out.