Mapping and compositional analysis of mare basalts in the Aristarchus region of the Moon using Clementine data

The process of accurately defining and outlining mare basalt units is nec- essary for constraining the stratigraphy and ages of basalt units, which are used to determine the duration and the flux of lunar volcanism. We use a combination of Clementine＇s five-band ultraviolet/visible data and Ti02 and FeO abundance distri- bution maps to define homogenous mare basalt units and characterize their composi- tional variations （with maturity） in the Aristarchus region. With 20 groups of distinct mare basaltic soils identified using the method in this paper, six additional spectrally defined areas and five basaltic units are constructed, and their mineralogic quanfiza- tion values provide new constraints on their temporal and spatial evolution. Our results indicate that the Aristarchus region has diverse basalt units and a complex history of volcanic evolution. We also demonstrate that the techniques, from which spectrally distinct mare basalts can be mapped, performed well in this study and can be confi- dently expanded to other mare regions of the Moon.

Experiment on diffuse reflection laser ranging to space debris and data analysis

Space debris poses a serious threat to human space activities and needs to be measured and cataloged. As a new technology for space target surveillance, the measurement accuracy of diffuse reflection laser ranging （DRLR） is much higher than that of microwave radar and optoelectronic measurement. Based on the laser ranging data of space debris from the DRLR system at Shanghai Astronomical Observatory acquired in March-April, 2013, the characteristics and precision of the laser ranging data are analyzed and their applications in orbit determination of space debris are discussed, which is implemented for the first time in China. The experiment indicates that the precision of laser ranging data can reach 39 cm-228 cm. When the data are sufficient enough （four arcs measured over three days）, the orbital accuracy of space debris can be up to 50 m.

Data processing and preliminary results of the Chang'e-3 VIS/NIR Imaging Spectrometer in-situ analysis

The Chang＇e-3 Visible and Near-infrared Imaging Spectrometer （VNIS） is one of the four payloads on the Yutu rover. After traversing the landing site during the first two lunar days, four different areas are detected, and Level 2A and 2B ra- diance data have been released to the scientific community. The released data have been processed by dark current subtraction, correction for the effect of temperature, radiometric calibration and geometric calibration. We emphasize approaches for re- flectance analysis and mineral identification for in-situ analysis with VNIS. Then the preliminary spectral and mineralogical results from the landing site are derived. After comparing spectral data from VNIS with data collected by the Ma instrument and samples of mare that were returned from the Apollo program, all the reflectance data have been found to have similar absorption features near 1000 nm except lunar sample 71061. In addition, there is also a weak absorption feature between 1750-2400nm on VNIS, but the slopes of VNIS and Ma reflectance at longer wavelengths are lower than data taken from samples of lunar mare. Spectral parameters such as Band Centers and Integrated Band Depth Ratios are used to analyze mineralogical features. The results show that detection points E and N205 are mixtures of high-Ca pyroxene and olivine, and the composition of olivineat point N205 is higher than that at point E, but the compositions of detection points S3 and N203 are mainly olivine-rich. Since there are no obvious absorption features near 1250 nm, plagioclase is not directly identified at the landing site.

Data processing and error analysis for the CE-1 Lunar microwave radiometer

The microwave radiometer （MRM） onboard the Chang＇ E-1 （CE-I） lu- nar orbiter is a 4-frequency microwave radiometer, and it is mainly used to obtain the brightness temperature （TB） of the lunar surface, from which the thickness, temperature, dielectric constant and other related properties of the lunar regolith can be derived. The working mode of the CE-1 MRM, the ground calibration （including the official calibration coefficients）, as well as the acquisition and processing of the raw data are introduced. Our data analysis shows that TB increases with increasing frequency, decreases towards the lunar poles and is significantly affected by solar illumination. Our analysis also reveals that the main uncertainty in TB comes from ground calibration.

Analysis of the geomorphology surrounding the Chang'e-3 landing site

Chang＇e-3 （CE-3） landed on the Mare Imbrium basin in the east part of Sinus Iridum （19.51°W, 44.12°N）, which was China＇s first soft landing on the Moon and it started collecting data on the lunar surface environment. To better understand the environment of this region, this paper utilizes the available high-resolution topography data, image data and geological data to carry out a detailed analysis and research on the area surrounding the landing site （Sinus Iridum and 45 km×70 km of the landing area） as well as on the topography, landform, geology and lunar dust of the area surrounding the landing site. A general topographic analysis of the surrounding area is based on a digital elevation model and digital elevation model data acquired by Chang＇e-2 that have high resolution; the geology analysis is based on lunar geological data published by USGS; the study on topographic factors and distribution of craters and rocks in the surrounding area covering 4km^4km or even smaller is based on images from the CE-3 landing camera and images from the topographic camera; an analysis is done of the effect of the CE-3 engine plume on the lunar surface by comparing images before and after the landing using data from the landing camera. A comprehensive analysis of the results shows that the landing site and its surrounding area are identified as typical lunar mare with flat topography. They are suitable for maneuvers by the rover, and are rich in geological phenomena and scientific targets, making it an ideal site for exploration.

Phase analysis of sunspot group numbers on both solar hemispheres

Cross-correlation analysis and wavelet transform methods are proposed to investigate the phase relationship between the monthly sunspot group numbers in the solar northern and southern hemispheres. It is found that （1） the monthly sunspot group numbers in the northern hemisphere begin two months earlier than those in the southern one, which should lead to phase asynchrony between them but with a slight effect; （2） the Schwabe cycle length for the monthly sunspot group numbers in the two hemispheres obviously differs from each other, and the mean Schwabe cycle length of the monthly sunspot group numbers in the northern hemisphere is slightly larger than that in the southern one; （3） the monthly sunspot group numbers in the northern hemisphere precede those in the southern hemisphere during the years of about 1874- 1927, after which, the southern hemisphere leads the northern hemisphere in the years 1928-1964, and then the northern hemisphere leads in time till the present.

A Periodicity Analysis of the Light Curve of 3C 454.3

We analyzed the radio light curves of 3C 454.3 at frequencies 22 and 37 GHz taken from the database of Metsaeovi Radio Observatory, and found evidence of quasi-periodic activity. The light curves show great activity with very complicated non-sinusoidal variations. Two possible periods, a very weak one of 1.57 ± 0.12 yr and a very strong one of 6.15 ±0.50 yr were consistently identified by two methods, the Jurkevich method and power specmun estimation. The period of 6.15 ± 0.50 yr is consistent with results previously reported by Ciaramella et al. and Webb et al. Applying the binary black hole model to the central structure we found black hole masses of 1.53 × 10^9M⊙ and 1.86 × 10^8M⊙, and predicted that the next radio outburst is to take place in 2006 March and April.

Performance analysis of parallel gravitational N-body codes on large GPU clusters

We compare the performance of two very different parallel gravitational N-body codes for astrophysical simulations on large Graphics Processing Unit（GPU） clusters, both of which are pioneers in their own fields as well as on certain mutual scales- NBODY6＋＋ and Bonsai. We carry out benchmarks of the two codes by analyzing their performance, accuracy and efficiency through the modeling of structure decomposition and timing measurements. We find that both codes are heavily optimized to leverage the computational potential of GPUs as their performance has approached half of the maximum single precision performance of the underlying GPU cards. With such performance we predict that a speed-up of200- 300 can be achieved when up to 1k processors and GPUs are employed simultaneously. We discuss the quantitative information about comparisons of the two codes, finding that in the same cases Bonsai adopts larger time steps as well as larger relative energy errors than NBODY6＋＋, typically ranging from10- 50 times larger, depending on the chosen parameters of the codes. Although the two codes are built for different astrophysical applications, in specified conditions they may overlap in performance at certain physical scales, thus allowing the user to choose either one by fine-tuning parameters accordingly.

Statistical analysis of geomagnetic field variations during the partial solar eclipse on 2011 January 4 in Turkey

Some geophysical parameters, such as those related to gravitation and the geomagnetic field, could change during solar eclipses. In order to observe geomagnetic fluctuations, geomagnetic measurements were carded out in a limited time frame during the partial solar eclipse that occurred on 2011 January 4 and was observed in Canakkale and Ankara, Turkey. Additionally, records of the geomagnetic field spanning 24 hours, obtained from another observatory （in Iznik, Turkey）, were also analyzed to check for any peculiar variations. In the data processing stage, a polynomial fit, following the application of a running average routine, was applied to the geomagnetic field data sets. Geomagnetic field data sets indicated there was a characteristic decrease at the beginning of the solar eclipse and this decrease can be well-correlated with previous geomagnetic field measurements that were taken during the total solar eclipse that was observed in Turkey on 2006 March 29. The behavior of the geomagnetic field is also consistent with previous observations in the literature. As a result of these analyses, it can be suggested that eclipses can cause a shielding effect on the geomagnetic field of the Earth.

Example-based super-resolution for single-image analysis from the Chang'e-1 Mission

Due to the low spatial resolution of images taken from the Chang＇e-1 （CE-I） orbiter, the details of the lunar surface are blurred and lost. Considering the limited spatial resolution of image data obtained by a CCD camera on CE-1, an example-based super-resolution （SR） algorithm is employed to obtain high- resolution （HR） images. SR reconstruction is important for the application of image data to increase the resolution of images. In this article, a novel example-based algorithm is proposed to implement SR reconstruction by single-image analysis, and the computational cost is reduced compared to other example-based SR methods. The results show that this method can enhance the resolution of images using SR and recover detailed information about the lunar surface. Thus it can be used for surveying HR terrain and geological features. Moreover, the algorithm is significant for the HR processing of remotely sensed images obtained by other imaging systems.

Photometric solution and period analysis of the contact binary system AH Cnc

Photometric observations of AH Cnc, a W UMa-type system in the open cluster M67, were car- fled out by using the 50BIN telescope. About 100h of time-series/3- and V-band data were taken, based on which eight new times of light minima were determined. By applying the Wilson-Devinney method, the light curves were modeled and a revised photometric solution of the binary system was derived. We con- firmed that AH Cnc is a deep contact （f = 51%）, low mass-ratio （q - 0.156） system. Adopting the distance modulus derived from study of the host cluster, we have re-calculated the physical parameters of the binary system, namely the masses and radii. The masses and radii of the two components were estimated to be respectively 1.188（4-0.061） Me, 1.332（4-0.063） RQ for the primary component and 0.185（4-0.032） Me, 0.592（4-0.051） Re for the secondary. By adding the newly derived minimum timings to all the available data, the period variations of AH Cnc were studied. This shows that the orbital period of the binary is con- tinuously increasing at a rate of dp/dt = 4.29 x 10-10 d yr-1. In addition to the long-term period increase, a cyclic variation with a period of 35.26 yr was determined, which could be attributed to an unresolved tertiary component of the system.

Spectral Synthesis via Mean Field approach to Independent Component Analysis

We apply a new statistical analysis technique, the Mean Field approach to Independent Component Analysis（MF-ICA） in a Bayseian framework, to galaxy spectral analysis. This algorithm can compress a stellar spectral library into a few Independent Components（ICs）, and the galaxy spectrum can be reconstructed by these ICs. Compared to other algorithms which decompose a galaxy spectrum into a combination of several simple stellar populations, the MF-ICA approach offers a large improvement in efficiency. To check the reliability of this spectral analysis method, three different methods are used：（1） parameter recovery for simulated galaxies,（2） comparison with parameters estimated by other methods, and（3） consistency test of parameters derived with galaxies from the Sloan Digital Sky Survey. We find that our MF-ICA method can not only fit the observed galaxy spectra efficiently, but can also accurately recover the physical parameters of galaxies. We also apply our spectral analysis method to the DEEP2 spectroscopic data, and find it can provide excellent fitting results for low signal-to-noise spectra.

The use of laser altimetry data in Chang'E-1 precision orbit determination

Accurate altimetric measurement not only can be applied to the calculation of a topography model but also can be used to improve the quality of the orbit reconstruction in the form of crossovers. Altimetry data from the Chang＇E-1 （CE-1） laser altimeter are analyzed in this paper. The differences between the crossover constraint equation in the form of height discrepancies and in the form of minimum distances are mainly discussed. The results demonstrate that the crossover constraint equation in the form of minimum distances improves the CE- l orbit precision. The overlap orbit performance has increased - 30% compared to the orbit using only tracking data. External assessment using the topography model also shows orbit improvement. The results will be helpful for recomputing ephemeris and improving the CE-1 topography model.

Random forest algorithm for classification of multiwavelength data

We introduced a decision tree method called Random Forests for multiwavelength data classification. The data were adopted from different databases, including the Sloan Digital Sky Survey （SDSS） Data Release five, USNO, FIRST and ROSAT. We then studied the discrimination of quasars from stars and the classification of quasars, stars and galaxies with the sample from optical and radio bands and with that from optical and X-ray bands. Moreover, feature selection and feature weighting based on Random Forests were investigated. The performances based on different input patterns were compared. The experimental results show that the random forest method is an effective method for astronomical object classification and can be applied to other classification problems faced in astronomy. In addition, Random Forests will show its superiorities due to its own merits, e.g. classification, feature selection, feature weighting as well as outlier detection.

Retrieving volume FeO and TiO_2 abundances of lunar regolith with CE-2 CELMS data using BPNN method

The volume FeO and TiO_2 abundances(FTAs) of lunar regolith can be more important for understanding the geological evolution of the Moon compared to the optical and gamma-ray results. In this paper, the volume FTAs are retrieved with microwave sounder(CELMS) data from the Chang'E-2 satellite using the back propagation neural network(BPNN) method. Firstly, a three-layered BPNN network with five-dimensional input is constructed by taking nonlinearity into account. Then, the brightness temperature(TB) and surface slope are set as the inputs and the volume FTAs are set as the outputs of the BPNN network.Thereafter, the BPNN network is trained with the corresponding parameters collected from Apollo, Luna,and Surveyor missions. Finally, the volume FTAs are retrieved with the trained BPNN network using the four-channel TBderived from the CELMS data and the surface slope estimated from Lunar Orbiter Laser Altimeter(LOLA) data. The rationality of the retrieved FTAs is verified by comparing with the Clementine UV-VIS results and Lunar Prospector(LP) GRS results. The retrieved volume FTAs enable us to re-evaluate the geological features of the lunar surface. Several important results are as follows. Firstly, very-low-Ti(<1.5 wt.%) basalts are the most spatially abundant, and the surfaces with TiO_2> 5 wt.% constitute less than 10% of the maria. Also, two linear relationships occur between the FeO abundance(FA) and the TiO_2 abundance before and after the threshold, 16 wt.% for FA. Secondly, a new perspective on mare volcanism is derived with the volume FTAs in several important mare basins, although this conclusion should be verified with more sources of data. Thirdly, FTAs in the lunar regolith change with depth to the uppermost surface,and the change is complex over the lunar surface. Finally, the distribution of volume FTAs hints that the highlands crust is probably homogeneous, at least in terms of the microwave thermophysical parameters.

Space weathering of the Moon from in situ detection

Space weathering is an important surface process that occurs on the Moon and other airless bodies, especially those that have no magnetic field. The optical effects of the Moon's space weathering have largely been investigated in the laboratory for lunar samples and lunar analogues. However, duplication of pristine regolith on Earth is not possible. Here we report on space weathering from the unique perspective of the "Yutu" rover, which was part of the Chang'e-3(CE-3) mission, building on our previous work.Measurement of the visually undisturbed uppermost regolith as well as locations that have been affected by rocket exhaust from the spacecraft by the Visible-Near Infrared Spectrometer(VNIS) revealed that the returned samples provide biased information about the pristine lunar regolith. The uppermost surficial regolith is much more weathered than the regolith immediately below, and the finest fraction is rich in space weathered products. These materials are very dark and attenuated throughout the visible and near-infrared(VNIR) wavelengths, hence reducing the reflectance and masking the absorption features. The effects on the spectral slope caused by space weathering are wavelength-dependent: the visible and near-infrared continuum slope(VNCS) increases while the visible slope(VS) decreases. In the visible wavelengths, the optical effects of space weathering and Ti O_2 are identical: both reduce albedo and blue the spectra. This suggests that a new Ti O_2 abundance algorithm is needed. Optical maturity indices are related to composition and hence only locally meaningful. Since optical remote sensing can only sense the uppermost few microns of regolith and since this surface tends to be very weathered, the interpretation of surface composition using optical remote sensing data needs to be carefully evaluated. Sampling the uppermost surface is suggested.

Data processing and initial results of Chang'e-3 lunar penetrating radar

To improve our understanding of the formation and evolution of the Moon, one of the payloads onboard the Chang＇e-3 （CE-3） rover is Lunar Penetrating Radar （LPR）. This investigation is the first attempt to explore the lunar subsurface structure by using ground penetrating radar with high resolution. We have probed the subsur- face to a depth of several hundred meters using LPR. In-orbit testing, data processing and the preliminary results are presented. These observations have revealed the con- figuration of regolith where the thickness of regolith varies from about 4 m to 6 m. In addition, one layer of lunar rock, which is about 330 m deep and might have been accumulated during the depositional hiatus of mare basalts, was detected.

Ancient subsurface structure beneath crater Clavius:constraint by recent high-precision gravity and topography data

With the increasing precision of the GRAIL gravity field models and topography from LOLA, it is possible to investigate the substructure beneath crater Clavius. An admittance between gravity and topography data is commonly used to estimate selenophysical parameters, including load ratio, crustal thickness and density, and elastic thickness. Not only a surface load, but also a subsurface load is considered in estimation. The algorithm of particle swarm optimization(PSO) with a swarm size of 400 is employed as well.Results indicate that the observed admittance is best-fitted by the modeled admittance based on a spherical shell model, which was proved to be unsatisfactory in the previous study. The best-fitted load ratio f is around-0.194. Such a small load ratio conforms to the direct proportion between the nearly uncompensated topography and its corresponding negative gravity anomaly. It also indicates that a surface load dominates all the loads. Constrained within 2σSTD, a small crustal thickness(～30 km) and a crustal density of ～2587 kg m-3are found, quite close to the results from previous GRAIL research. Considering the well constrained crustal thickness and density, the best-fitted elastic thickness(～7 km) is rational. This result is slightly smaller than the previous study(～12 km). Such difference can be attributed to the difference in crustal density used and the precision of gravity and topography data. Considering that the small difference between the modeled gravity anomaly and observations is quite small, a parameter inversed here could be an indicator of the subsurface structure beneath Clavius.

Data processing for the Active Particle-induced X-ray Spectrometer and initial scientific results from Chang'e-3 mission

The Active Particle-induced X-ray Spectrometer （APXS） is an important payload mounted on the Yutu rover, which is part of the Chang＇e-3 mission. The sci- entific objective of APXS is to perform in-situ analysis of the chemical composition of lunar soil and rock samples. The radioactive sources, 55Fe and 109Cd, decay and produce a-particles and X-rays. When X-rays and a-particles interact with atoms in the surface material, they knock electrons out of their orbits, which release energy by emitting X-rays that can be measured by a silicon drift detector （SDD）. The elements and their concentrations can be determined by analyzing their peak energies and in- tensities. APXS has analyzed both the calibration target and lunar soil once during the first lunar day and again during the second lunar day. The total detection time lasted about 266 min and more than 2000 frames of data records have been acquired. APXS has three operating modes： calibration mode, distance sensing mode and detection mode. In detection mode, work distance can be calculated from the X-ray counting rate collected by SDD. Correction for the effect of temperature has been performed to convert the channel number for each spectrum to X-ray energy. Dead time correction is used to eliminate the systematic error in quantifying the activity of an X-ray pulse in a sample and derive the real count rate. We report APXS data and initial results during the first and second lunar days for the Yutu rover. In this study, we analyze the data from the calibration target and lunar soil on the first lunar day. Seven major elements, including Mg, A1, Si, K, Ca, Ti and Fe, have been identified. Comparing the peak areas and ratios of calibration basalt and lunar soil the landing site was found to be depleted in K, and have lower Mg and A1 but higher Ca, Ti, and Fe. In the future, we will obtain the elemental concentrations of lunar soil at the Chang＇e-3 landing site using APXS data.

Scientific data products and the data pre-processing subsystem of the Chang'e-3 mission

The Chang＇e-3 （CE-3） mission is China＇s first exploration mission on the surface of the Moon that uses a lander and a rover. Eight instruments that form the scientific payloads have the following objectives： （1） investigate the morphological features and geological structures at the landing site; （2） integrated in-situ analysis of minerals and chemical compositions; （3） integrated exploration of the structure of the lunar interior; （4） exploration of the lunar-terrestrial space environment, lunar sur- face environment and acquire Moon-based ultraviolet astronomical observations. The Ground Research and Application System （GRAS） is in charge of data acquisition and pre-processing, management of the payload in orbit, and managing the data products and their applications. The Data Pre-processing Subsystem （DPS） is a part of GRAS. The task of DPS is the pre-processing of raw data from the eight instruments that are part of CE-3, including channel processing, unpacking, package sorting, calibration and correction, identification of geographical location, calculation of probe azimuth angle, probe zenith angle, solar azimuth angle, and solar zenith angle and so on, and conducting quality checks. These processes produce Level 0, Level 1 and Level 2 data. The computing platform of this subsystem is comprised of a high-performance computing cluster, including a real-time subsystem used for processing Level 0 data and a post-time subsystem for generating Level 1 and Level 2 data. This paper de- scribes the CE-3 data pre-processing method, the data pre-processing subsystem, data classification, data validity and data products that are used for scientific studies.