Hydrocode modelling of hypervelocity impacts on ice

Experimental data are now widely available for the size of craters resulting from hypervelocity impacts of millimetre scale projectiles onto water ice targets. At such size scales the bowl shaped crater formed in ductile materials, or in larger scale impacts, is here surrounded by a large spallation zone due to the brittle nature of the ice. Modelling of these impacts therefore has to take account of this spallation. Here we used the iSALE2 hydrocode to simulate such impacts and compared the results to experimental data. We found that it was possible to reproduce the experimental data over a range of speeds (1–7 km s⁻¹) for aluminium and copper projectiles. Initially, to reproduce the large spallation regions around the craters it was assumed that above a certain degree of damage, material was removed by spallation. However this simple one-parameter model failed to model the crater depth adequately. Accordingly, to obtain the best agreement of the simulations with the experimental data, a two-step ice strength was introduced, whereby above a critical amount of damage (0.95), the yield strength reduced from 1 MPa (intact) to 70 kPa (damaged). As a result, experimental data for crater depth and diameter and the results of the simulations agree to within 6% for diameter and 5% for depth over the impact energy range used in the experiments (1–240 J).     

Integrated method for crater detection from topography and optical images and the new PH9224GT catalogue of Phobos impact craters

In a large majority of lunar and planetary surface images, impact craters are the most abundant geological features. Therefore, it is not surprising that crater detection algorithms (CDAs) are one of the most studied subjects of image processing and analysis in lunar and planetary science. In this work we are proposing an Integrated CDA, consisting of: (1) utilization of DEM (digital elevation map)-based CDA; (2) utilization of an optical-based CDA; (3) re-projection of used datasets and crater coordinates from normal to rotated view and back; (4) correction of the brightness and contrast of a used optical image; and (5) tile generation for the optical-based CDA and an assembling of results with an elimination of multiple detections, in combination with a pyramid approach down to the resolution of the available DEM image; and (6) a final integration of the results of DEM-based and optical-based CDAs, including a removal of duplicates. The proposed CDA is applied to one specific asteroid-like body, the small Martian moon Phobos. The experimental evaluation of the proposed CDA is done by a manual verification of crater-candidates and a search for uncatalogued craters. The evaluation has shown that the proposed CDA was used successfully for cataloging Phobos craters. The major result of this paper is the PH9224GT – currently the most complete global catalogue of the 9224 Phobos craters. The possible applications of the new catalogue are: (1) age estimations for any selected location; and (2) comparison/evaluation of the different chronology and production functions for Phobos. This confirms the practical applicability of the new Integrated CDA – an additional result of this paper, which can be used in order to considerably extend the current crater catalogues.     

The first new application of the Mathematical Theory of Stochastic Processes to lunar and planetary science: Topography-Profile Diagrams of Mars

The Mathematical Statistics Theory (MST) and the Mathematical Theory of Stochastic Processes (MTSP) are different branches of the more general Mathematical Probability Theory (MPT) that can be used to investigate physical processes through mathematics. Each model of a stochastic process, according to MTSP, can provide one or more interpretations in the MST domain. A large body of work on impact crater statistics according to MST exists, showing cumulative crater frequency (N km⁻²) as a function of age (years) for some particular crater diameter. However, this is only one possible representation in the MST domain of the bombardment of the planetary surface modeled as a stochastic process according to MTSP. The idea that other representations are possible in the MST domain of the same stochastic process from MTSP has been recently presented. The importance of the approach is that each such mathematical-based interpretation can provide a large amount of new information. Coupled with MOLA data, Topography-Profile Diagrams (TPD) are one of the many examples that can provide a large amount of new information regarding the history of Mars. TPD consists of: (1) Topography-Profile Curve (TPC), which is a representation of the planet’s topography, (2) Density-of-Craters Curve (DCC), which represents density of craters, (3) Filtered-DCC (FDCC), which represents DCC filtered by a low-pass filter, included with the purpose of reducing the noise, and (4) Level-of-Substance-Over-Time Curve (LSOTC), which represents interpretation of the influence on the distribution of craters shown by FDCC. TPC uniquely corresponds to the computation of TPD, whereas DCC depends on algorithms for computing the elevation of each crater according to the topography, center coordinates, and radius of impact crater, and FDCC relies on the architecture of the custom designed low-pass filter for filtering DCC. However, all variations of DCC and FDCC, which includes the various impact crater data sets, showed a correlation among the density of craters and elevation over 70–80% of the planet surface. Additionally, if we assume that the ocean primarily caused the noted correlation, LSOTC offers a mathematical approach for estimating topographic change of the ocean’s extent over time. Accordingly, TPD is the first new practical application of MTSP to lunar and planetary sciences, showing correlation of topography to a physical process.     

一种基于多尺度边缘提取的陨石坑检测算法

针对星体表面的陨石坑可用于探测器的自主导航、障碍识别等任务,提出一种基于多尺度边缘提取的陨石坑检测算法。该算法首先利用高斯金字塔得到不同尺度的陨石坑图像;其次,针对不同尺度的陨石坑图像,利用EDPF边缘提取算法对陨石坑进行边缘提取,并连接关键边缘像素点为直线段来近似表示图像边缘;然后将具有相同偏转方向的边缘直线段连接成圆弧,并将有相似半径和中心的圆弧拟合成候选圆和椭圆;最后对候选圆、椭圆进行验证。该算法的优点在于,能够准确地检测出陨石坑,有较高的检测率,且对存在较多陨石坑的图像有较好的检测结果。     

月球撞击坑辐射纹的形成和消失机制及可能的年代学应用

辐射纹是年轻撞击坑周缘呈辐射状分布的明亮细条纹,是月表最显著的地貌特征之一,也是月球科学领域热点研究之一.根据目前对于月球撞击坑辐射纹的形成和类型的认识,分析辐射纹消失相关的空间风化、撞击导致的物质混合等地质过程;比较不同形成年龄撞击坑辐射纹的光学成熟度(OMAT)剖面,发现溅射物逐渐成熟过程中OMAT剖面的演变过程,...     

基于深度学习的大型陨石坑识别方法研究

陨石坑是天体表面最为显著的地形特征，传统陨石坑识别方法主要是对小型陨石坑正负样本的二分类问题研究，且效率和精度均不高。以星体宏观视角下的大型陨石坑作为研究对象，结合图像处理和神经网络等方面的知识，创建了来自不同数据源的陨石坑样本数据库，研究了数据源对网络模型泛化能力的影响，提出了一种效率更高的陨石坑多分类识别方法。在非极大值抑制（NMS）算法基础上，提出了一种精度更高的陨石坑检测算法。经过参数优化和实验验证，构建的基于深度学习的多尺度多分类陨石坑自动识别网络框架取得了较高的准确率，在同源验证集上识别率可达0.985，在异源验证集上识别率可达0.863，并且有效改善了目标检测时检测框冗余及误检测的问题。      

Detection of sub-kilometer craters in high resolution planetary images using shape and texture features

Counting craters is a paramount tool of planetary analysis because it provides relative dating of planetary surfaces. Dating surfaces with high spatial resolution requires counting a very large number of small, sub-kilometer size craters. Exhaustive manual surveys of such craters over extensive regions are impractical, sparking interest in designing crater detection algorithms (CDAs). As a part of our effort to design a CDA, which is robust and practical for planetary research analysis, we propose a crater detection approach that utilizes both shape and texture features to identify efficiently sub-kilometer craters in high resolution panchromatic images. First, a mathematical morphology-based shape analysis is used to identify regions in an image that may contain craters; only those regions - crater candidates - are the subject of further processing. Second, image texture features in combination with the boosting ensemble supervised learning algorithm are used to accurately classify previously identified candidates into craters and non-craters. The design of the proposed CDA is described and its performance is evaluated using a high resolution image of Mars for which sub-kilometer craters have been manually identified. The overall detection rate of the proposed CDA is 81%, the branching factor is 0.14, and the overall quality factor is 72%. This performance is a significant improvement over the previous CDA based exclusively on the shape features. The combination of performance level and computational efficiency offered by this CDA makes it attractive for practical application.     

Experimental study of impact-cratering damage on brittle cylindrical column model as a fundamental component of space architecture

The cylindrical column of brittle material processed from soil and rock is a fundamental component of architectures on the surface of solid bodies in the solar system. One of the most hazardous events for the structure is damaging by hypervelocity impacts by meteoroids and debris. In such a background, cylindrical columns made of plaster of Paris and glass-bead-sintered ceramic were impacted by spherical projectiles of nylon, glass, and steel at velocity of about 1–4.5 km/s. Measured crater radii, depth, and excavated mass expressed by a function of the cylinder radius are similar irrespective of the target material, if those parameters are normalized by appropriate parameters of the crater produced on the flat-surface target. The empirical scaling relations of the normalized crater radii and depth are provided. Using them, crater dimensions and excavated mass of crater on cylindrical surface of any radius can be predicted from the existing knowledge of those for flat surface. Recommendation for the minimum diameter of a cylinder so as to resist against a given impact is provided.     

Single crater-aided inertial navigation for autonomous asteroid landing

In this paper, a novel crater-aided inertial navigation approach for autonomous asteroid landing mission is developed. It overcomes the major deficiencies of existing approaches in the literature, which mainly focuses on the case where craters are abundant in the camera field of view. As a result, traditional crater based methods require at least three craters to achieve crater matching, which limits their application in final landing phase where craters are scarce in the camera’s field of view. In contrast, the proposed algorithm enables single crater based crater matching based on a novel 2D-3D crater re-projection model. The re-projection model adopts inertial measurements as a reference, and re-projects the 3D crater model onto descent images to achieve the matching to its counterpart. An asteroid landing simulation toolbox is developed to validate the performance of the proposed approach. Through comparison with the state-of-the-art local image feature and crater based navigation algorithms, the proposed approach is validated to achieve a competitive performance in terms of feature matching and pose estimation accuracy with a much lighter computational cost.     

Bayesian network-based extraction of lunar impact craters from optical images and DEM data

Impact craters are among the most noticeable geomorphological features on the planetary surface and yield significant information about terrain evolution and the history of the solar system. Thus, the recognition of impact craters is an important branch of modern planetary studies. Aiming at addressing problems associated with the insufficient and inaccurate detection of lunar impact craters, a decision fusion method within the Bayesian network (BN) framework is developed in this paper to handle multi-source information from both optical images and associated digital elevation model (DEM) data. First, we implement the edge-based method for efficiently searching crater candidates which are the image patches that can potentially contain impact craters. Secondly, the multi-source representations of an impact crater derived from both optical images and DEM data are proposed and constructed to quantitatively describe the two-dimensional (2D) and three-dimensional (3D) morphology, consisting of Histogram of Oriented Gradient (HOG), Histogram of Multi-scale Slope (HMS) and Histogram of Multi-scale Aspect (HMA). Finally, a BN-based framework integrates the multi-source representations of impact craters, which can provide reductant and complementary information, for distinguishing craters from non-craters. To evaluate the effectiveness and robustness of the proposed method, experiments were conducted on three lunar scenes using both orthoimages from the Lunar Reconnaissance Orbiter (LRO) and DEM data acquired by the Lunar Orbiter Laser Altimeter (LOLA). Experimental results demonstrate that integrating optical images with DEM data significantly decreases the number of false positives compared with using optical images alone, with F₁-score of 84.8% on average. Moreover, compared with other existing fusion methods, our proposed method was quite advantageous especially for the detection of small-scale craters with diameters less than 1000 m.     

Lunar apex–antapex cratering asymmetry as an impactor recorder in the Earth–Moon system

The degree of apex–antapex cratering asymmetry of a synchronously rotating satellite primarily depends on the mean encounter velocity of impactors with respect to the planetary system and the orbital velocity of the satellite. This means that we can estimate the mean encounter velocity of impactors by observing the apex–antapex cratering asymmetry, if the relationship between these is known. To apply this technique to the Moon, we attempt to derive the relationship between the mean encounter velocity of impactors and the degree of the lunar cratering asymmetry as a function of time, considering the temporal variation in the lunar orbital velocity during the last 4.0 Gyr. We used the cratering asymmetry of Zahnle et al. [Zahnle, K., Schenk, P., Sobieszczyk, S. et al. Differential cratering of synchronously rotating satellites by ecliptic comets. Icarus 153, 111–129, 2001] to obtain the relationship. Applying this relationship enables us to estimate the impactor’s velocity of the Earth–Moon system from an investigation of the spatial distribution of lunar craters. Furthermore, we re-evaluate the cratering asymmetry’s influence on lunar cratering chronology.     

A new approach based on crater detection and matching for visual navigation in planetary landing

This paper provides an approach of crater detection and matching to visual navigation in planetary landing missions. The approach aims to detect craters on the planetary surface and match them to a landmark database during the descent phase of a planetary landing mission. Firstly an image region pairing method is proposed to detect the crater by using an image region feature detector. Then a WTA-rule is adopted to match the detected crater to the crater in database. To further reduce the false matching rate, an efficient method for reducing false matches using parameters of crater in 3-D database is proposed. Real images of planetary terrain and a semi-physical planetary landing simulation platform are utilized to test the performance of the approach, simulation results show the proposed approach is able to match the required number of craters to the database for pin-point planetary landing with a low rate of false detection and false matching, which will lead to an improved planetary landing precision.     

Midday bottomside electron density profiles during moderate solar activity and comparison with IRI-2001

Bottomside electron density (Ne-h) profiles during midday (10–14 h) are analyzed using modern digital ionosonde observations at a low-middle latitude station, New Delhi (28.6N, 77.2E, dip 42.4N), for the period from January 2003 to December 2003, pertaining to moderate solar activity (MSA). Each individual profile is normalized with respect to the peak height and density (hmF2, NmF2) of the F2-region. These profiles are compared with those obtained from the International Reference Ionosphere (IRI-2001) model. Bilitza [Bilitza, D. International Reference Ionosphere 2000. Radio Sci. 36 (2), 261–275, 2001] using both the options namely: Gulyaeva’s model [Gulyaeva, T.L. Progress in ionospheric informatics based on electron density profile analysis of ionograms. Adv. Space Res. 7 (6) 39–48, 1987] and B0 Tab. option [Bilitza, D., Radicella, S.M., Reinisch, B.W., Adeniyi, J.O., Mosert Gonzalez, M.E., Zhang, S.R., Obrou, O. New B0 and B1 models for IRI. Adv. Space Res. 25 (1), 89–95, 2000]. The study reveals that during summer and equinox, the IRI model with B0 Tab. option in general, produces better agreement with the observed median profiles, while the IRI predictions using Gulyaeva’s option, overestimate the electron density distribution at all the heights below the F2-peak. However, during winter, in general, the IRI model, using both the options, reveals shows fairly good agreement with the observations.     

Application of hydrocode modelling to the study of hypervelocity impact crater morphology

In order to obtain a better understanding and model of the natural and artificial particulate environment from measurements of impact damage features on returned spacecraft materials, it is necessary to be able to determine how the size and shape of an impact feature are related to the parameters of the impacting particle. The AUTODYN-3D hydrocode has been used to study the effects of projectile density, velocity and impact angle on the depth, diameter and ellipticity of the impact craters. The results are used to determine the distributions of crater depth to crater diameter ratios and of crater ellipticities to be expected on an aluminium surface exposed to an isotropic distribution of incident particles of given densities and velocities. Comparison of these calculated distributions with those observed for craters on aluminium clamps on various faces of the Long Duration Exposure Facility shows that particles with a wide range of densities, including significant proportions both greater and smaller than that of aluminium, were responsible for these craters.     

A machine learning approach to crater detection from topographic data

Craters are distinctive features on the surfaces of most terrestrial planets. Craters reveal the relative ages of surface units and provide information on surface geology. Extracting craters is one of the fundamental tasks in planetary research. Although many automated crater detection algorithms have been developed to exact craters from image or topographic data, most of them are applicable only in particular regions, and only a few can be widely used, especially in complex surface settings. In this study, we present a machine learning approach to crater detection from topographic data. This approach includes two steps: detecting square regions which contain one crater with the use of a boosting algorithm and delineating the rims of the crater in each square region by local terrain analysis and circular Hough transform. A new variant of Haar-like features (scaled Haar-like features) is proposed and combined with traditional Haar-like features and local binary pattern features to enhance the performance of the classifier. Experimental results with the use of Mars topographic data demonstrate that the developed approach can significantly decrease the false positive detection rate while maintaining a relatively high true positive detection rate even in challenging sites.     

Hybrid method for crater detection based on topography reconstruction from optical images and the new LU78287GT catalogue of Lunar impact craters

Impact craters are ubiquitous and well-studied structures of high geological relevance on the surfaces of the Earth’s Moon, the terrestrial planets, the asteroids and the satellites of the outer planets. Therefore, it is not surprising that crater detection algorithms (CDAs) are one of the most studied subjects of image processing and analysis in lunar and planetary science. In this paper we are proposing a Hybrid CDA: a modified DEM (digital elevation map) reconstruction method used as a step in an existing CDA based on Hough transform. The new Hybrid CDA consists of: (1) reconstruction of topography from optical images using a shape from shading approach; (2) utilization of the DEM-based CDA; (3) correction of brightness and contrast of optical images used in order to be more suitable for evaluation of detections. An additional result of this work is a new method for evaluation of topography reconstruction algorithms, using a DEM-based CDA and an earlier approach for evaluation of CDAs. The new Hybrid CDA was tested using two Chandrayaan-1 Moon Mineralogy Mapper (M³) images and two excerpts of the Lunar Reconnaissance Orbiter (LRO) Wide Angle Camera (WAC) global optical image mosaic. As a result, the number of craters inside these four regions increased considerably from 1754 (as available in the previous LU60645GT catalogue) to 19 396 craters (as available in the resulting new LU78287GT catalogue). This confirmed the practical applicability of the new Hybrid CDA, which can be used in order to considerably extend current crater catalogues.     

Effects of geomagnetic secular variations on cosmic ray access to the terrestrial environment

Taking advantage of the cutoff computations performed for more than a hundred locations from 1955 to 1995 [every 5 years; Shea, M.A., Smart, D.F. Vertical cutoff rigidities for cosmic ray stations since 1955, in: Proceedings of the ICRC 2001, Hamburg, Copernicus Gesellschaft, vol. 10, pp. 4063–4066, 2001], we carefully checked the relationship between the vertical cutoff rigidity and the McIlwain parameter introduced by Shea et al. [Shea, M.A., Smart, D.F., Gentile, L.C. Estimating cosmic ray vertical cutoff rigidities as a function of the McIlwain L-parameter for different epochs of the geomagnetic field. Phys. Earth Planet. Int., 48, 200–205, 1987]. We derived an updated algorithm that can be used outside the polar and equatorial regions, avoiding time consuming computations. Results for the European area and 1990 epoch suggest that the fast evaluation is accurate within 0.1 GV in 26 out of the 30 considered locations.     

Laboratory investigations of the survivability of bacteria in hypervelocity impacts.

It is now well established that material naturally moves around the Solar System, even from planetary surface to planetary surface. Accordingly, the idea that life is distributed throughout space and did not necessarily originate on the Earth but migrated here from elsewhere (Panspermia) is increasingly deemed worthy of consideration. If life arrived at the Earth from space, its relative speed will typically be of order many km s-1, and the resulting collision with the Earth and its atmosphere will be in the hypervelocity regime. A mechanism for the bacteria to survive such an impact is required. Therefore a programme of hypervelocity impacts in the laboratory at (4.5 +/- 0.6) km s-1 was carried out using bacteria (Rhodococcus) laden projectiles. After impacts on a variety of target materials (rock, glass and metal) attempts were made to culture Rhodococcus from the surface of the resulting craters and also from the target material ejected during crater formation. Control shots with clean projectiles yielded no evidence for Rhodococcus growth from any crater surface or ejecta. When projectiles doped with Rhodococcus were used no impact crater surface yielded colonies of Rhodococcus. However, for four shots of bacteria into rock (two on chalk and two on granite) the ejecta was afterwards found to give colonies of Rhodococcus. This was not true for shots onto glass. In addition, shots into aerogel (density 96 kg m-3) were also carried out (two with clean projectiles and two with projectiles with Rhodococcus). This crudely simulated aero-capture in a planetary atmosphere. No evidence for Rhodococcus growth was found from the projectiles captured in the aerogel from any of the four shots.     

Electric field and ion density anomalies in the mid latitude ionosphere: Possible connection with earthquakes?

The problem of earthquake prediction has stimulated the search for a correlation between seismic activity and ionospherical anomalies. We found observational evidence of possible earthquake effects in the near-equatorial and low latitude ionosphere; these ionospheric anomalies have been proposed by Gousheva et al. [Gousheva, M., Glavcheva, R., Danov, D., Angelov P., Hristov, P., Influence of earthquakes on the electric field disturbances in the ionosphere on board of the Intercosmos-Bulgaria-1300 satellite. Compt. Rend. Acad. Bulg. Sci. 58 (8) 911–916, 2005a; Gousheva, M., Glavcheva, R., Danov, D., Angelov, P., Hristov, P., Kirov, B., Georgieva, K., Observation from the Intercosmos-Bulgaria-1300 satellite of anomalies associated with seismic activity. In: Poster Proceeding of 2nd International Conference on Recent Advances in Space Technologies: Space in the Service of Society, RAST ‘2005, June 9–11, Istanbul, Turkey, pp. 119–123, 2005b; Gousheva, M., Glavcheva, R., Danov, D., Angelov, P., Hristov, P., Kirov, B., Georgieva, K., Satellite monitoring of anomalous effects in the ionosphere probably related to strong earthquakes. Adv. Space Res. 37 (4), 660–665, 2006]. This paper presents new results from observations of the quasi-static electric field and ion density on board INTERCOSMOS-BULGARIA-1300 satellite in the mid latitude ionosphere above sources of moderate earthquakes. Data from INTERCOSMOS-BULGARIA-1300 satellite and seismic data (World Data Center, Denver, Colorado, USA) for magnetically quiet and medium quiet days are juxtaposed in time-space domain. For satellite’s orbits in the time period 15.09–01.10.1981 an increase in the horizontal and vertical components of the quasi-static electric field and fluctuations of the ion density are observed over zones of forthcoming seismic events. Some similar post effects are observed too. The emphasis of this paper is put on the anomalies which specify the mid latitude ionosphere. The obtained results contain important information because they confirm our previous results for near-equatorial and low latitude regions.     

Investigation of precursory phenomena in the ionosphere,atmosphere and groundwater before large earthquakes of M > 6.5

Precursory phenomena in the ionosphere, atmosphere and groundwater before large earthquakes (M > 6.5) are extensively investigated toward the earthquake prediction. Upward tornado type seismic clouds occurred near the epicenter associated with strong LF-VLF radio noises from lightning discharges in the evening of January 9, 1995 [Yamada, T., Oike, K. On the increase of electromagnetic noises before and after the 1995 Hyogo-Ken Nanbu earthquake. In: Hayakawa M. (Ed.), Atmospheric and Ionospheric Electromagnetic Phenomena Associated with Earthquakes. TERRAPUB, Tokyo, pp. 417–427, 1999] and anomalous foEs increases up to 10 MHz were detected at Shigaraki, 90 km of the epicenter and at Kokubunji, 500 km east of the epicenter [Ondoh, T. Anomalous sporadic-E layers observed before M7.2 Hyogo-ken Nanbu earthquake; Terrestrial gas emanation model. Adv. Polar Upper Atmos. Res. 17, 96–108, 2003; Ondoh, T. Anomalous sporadic-E ionization before a great earthquake, Adv. Space Research 34, 1830–1835, 2004] associated with strong ELF noises from lightning discharges in the daytime on January 15, 1995 [Hata, M., Fujii, T., Takumi, I. EM precursor of large-scale earthquakes in Japan, in: Abstracts of International Workshop on Seismo Electromagnetics (IWSE 2005), Univ. Electro-Communications, Chofu, Tokyo, Japan, March 15–17, pp. 182–186, 2005] before the M7.2 Hyogoken–Nanbu earthquake of January 17, 1995. The anomalous foEs increases occurred at epicentral distances within 500 km that are the same as those of the terrestrial gas emanations along active faults before large earthquakes [King, C.-Y. Gas geochemistry applied to earthquake prediction: An overview. J. Geophys. Res. 91(B12), 12269–12281, 1986]. The anomalous foEs increases seem to be a seismic precursor because geomagnetic and solar conditions were very quiet all day on January 15,1995 and the normal foEs in Japanese winter is below 6 MHz. No significant pre-seismic geomagnetic field variation was detected at epicentral distance of 100 km before this earthquake [Ondoh, T., Hayakawa, M. Anomalous occurrence of sporadic-E layers before the Hyogoken–Nanbu earthquake, M7.2 of January 17, 1995. In: Hayakawa, M. (Ed.), Atmospheric and Ionospheric Electromagnetic Phenomena Associated with Earthquakes, TERRAPUB, Tokyo, pp. 629–639, 1999; Ondoh, T., Hayakawa, M. Seismo discharge model of anomalous sporadic E ionization before great earthquakes. In: Hayakawa, M., O.A. Molchanov, (Eds.), Seismo Electromagnetics: Lithosphere–Atmosphere–Ionosphere Couplings, TERRAPUB, Tokyo, pp. 385–390, 2002; Ondoh. T., Hayakawa, M. Synthetic study of precursory phenomena of the M7.2 Hyogo-ken Nanbu earthquake. Phys. Chem. Earth 31, 378–388, 2006]. The foF2 decrease and h’F increase occurred before the M7.8 Hokkaido Nansei-Oki earthquake of July 12,1993 in a geomagnetic quiet period [Ondoh, T. Ionospheric disturbances associated with great earthquake of Hokkaido southwest coast, Japan of July 12, 1993. Phys. Earth Planet. Interiors. 105, 261–269, 1998; Ondoh, T. Seismo ionospheric phenomena. Adv. Space Res. 26, 8, 1267–1272, 2000]. Characteristic phase changes at terminator times of Omega 10.2 kHz waves passing 70 km of the epicenter extended toward darker local times by 1 h for 3 days before this earthquake due to lowering of the wave reflection height or ion density increases in the D region [Hayakawa, M., Molchanov, O. A., Ondoh, T., Kawai, E. The precursory signature effect of the Kobe earthquake on VLF subionospheric signals. J. Commun. Res. La., 43, 00. 169–180, 1996]. The radon concentration in the atmosphere over Ashiya fault, Kobe [Yasuoka, Y., Shinogi, M. Anomaly in atmospheric radon concentration: a possible precursor of the 1995 Kobe, Japan, earthquake. Health Phys. 72(5), 759–761, 1997] and in the groundwater at 17 m well in Nishinomiya, Japan [Igarashi, G., Saeki, S., Takahata, N., Sumikawa, K., Tasaki, S., Sasaki, Y., Takahashi, M., Sano, Y. Ground-water radon anomaly before the Kobe earthquake in Japan. Science 269, 60–61, 1995] had gradually increased since 2 months before the M7.2 earthquake, increased suddenly in December 1994, and rapidly returned to the normal low level of October, 1994 [Yasuoka, Y., Shinogi, M. 1997. Anomaly in atmospheric radon concentration: a possible precursor of the 1995 Kobe. Japan, earthquake. Health Phys. 72(5), 759–761.]. Radon concentration changes in the groundwater before the M 7.0 Izu-Oshima-kinkai earthquake, Japan on January 14, 1978 [Wakita, H., Nakamura, Y., Notsu, K., Noguchi, M., Asada, T. 1980. Radon anomaly: a possible precursor of the 1978 Izu-Oshima-kinkai earthquake. Science 207, 882–883] and the M6.8 Chengkung earthquake, Taiwan on December 10, 2003 [Kuo, T., Fan, K., Chen, W., Kuochen, H., Han, Y., Wang, C., Chang, T., Lee, Y. Radon anomaly at the Antung Hot Spring before the Taiwan M6.8 Chengkung earthquake. Proceedings, Thirty-First Workshop on Geothermal Reservoir Engineering, Stanford University, Stanford, California, January 30–February 1, 2006, SGP-TR-179, 2006] are also investigated to find common features of the groundwater radon concentration changes before large earthquakes (M > 6.5) in comparison with those before the M7.2 Hyogoken–Nanbu earthquake. Groundwater radon concentrations before the 3 large earthquakes had shown common characteristic changes of gradually initial ones from the normal level since about 2 months before the earthquake onsets, rapid decreases down to the minimum, and quick increases up to the maximum at 7–20 days before the earthquake onsets, respectively. These are very useful characteristics of pre-seismic radon anomaly for the earthquake prediction or warning. Promising observations toward the earthquake prediction are also discussed.