Quantitative Analysis of Seeing with Height and Time at Muztagh-Ata Site Based on ERA5 Database

Seeing is an important index to evaluate the quality of an astronomical site.To estimate seeing at the Muztagh-Ata site with height and time quantitatively,the European Centre for Medium-Range Weather Forecasts reanalysis database(ERA5)is used.Seeing calculated from ERA5 is compared consistently with the Differential Image Motion Monitor seeing at the height of 12 m.Results show that seeing decays exponentially with height at the Muztagh-Ata site.Seeing decays the fastest in fall in 2021 and most slowly with height in summer.The seeing condition is better in fall than in summer.The median value of seeing at 12 m is 0.89 arcsec,the maximum value is1.21 arcsec in August and the minimum is 0.66 arcsec in October.The median value of seeing at 12 m is 0.72arcsec in the nighttime and 1.08 arcsec in the daytime.Seeing is a combination of annual and about biannual variations with the same phase as temperature and wind speed indicating that seeing variation with time is influenced by temperature and wind speed.The Richardson number Ri is used to analyze the atmospheric stability and the variations of seeing are consistent with Ri between layers.These quantitative results can provide an important reference for a telescopic observation strategy.

Emerging Trends for Microbiome Analysis： From Single-Cell Functional Imaging to Microbiome Big Data

Method development has always been and will continue to be a core driving force of microbiome science, In this perspective, we argue that in the next decade, method development in microbiome analysis will be driven by three key changes in both ways of thinking and technological platforms： ① a shift from dissecting microbiota structure by sequencing to tracking microbiota state, function, and intercellular interaction via imaging; ② a shift from interrogating a consortium or population of cells to probing individual cells; and ③a shift from microbiome data analysis to microbiome data science. Some of the recent methoddevelopment efforts by Chinese microbiome scientists and their international collaborators that underlie these technological trends are highlighted here. It is our belief that the China Microbiome Initiative has the opportunity to deliver outstanding ＂Made-in-China＂ tools to the international research community, by building an ambitious, competitive, and collaborative program at the forefront of method development for microbiome science.

Preliminary analysis on the noise characteristics of MWISP data

Noise is a significant part within a millimeter-wave molecular line datacube.Analyzing the noise improves our understanding of noise characteristics,and further contributes to scientific discoveries.We measure the noise level of a single datacube from MWISP and perform statistical analyses.We identified major factors which increase the noise level of a single datacube,including bad channels,edge effects,baseline distortion and line contamination.Cleaning algorithms are applied to remove or reduce these noise components.As a result,we obtained the cleaned datacube in which noise follows a positively skewed normal distribution.We further analyzed the noise structure distribution of a 3 D mosaicked datacube in the range l=40°7 to 43°3 and b=-2°3 to 0°3 and found that noise in the final mosaicked datacube is mainly characterized by noise fluctuation among the cells.

Observation and Analysis of VLF Nocturnal Multimode Interference Phenomenon based on Waveguide Mode Theory

Very low frequency(VLF)signals are propagated between the ground-ionosphere.Multimode interference will cause the phase to show oscillatory changes with distance while propagating at night,leading to abnormalities in the received VLF signal.This study uses the VLF signal received in Qingdao City,Shandong Province,from the Russian Alpha navigation system to explore the multimode interference problem of VLF signal propagation.The characteristics of the effect of multimode interference phenomena on the phase are analyzed according to the variation of the phase of the VLF signal.However,the phase of VLF signals will also be affected by the X-ray and energetic particles that are released during the eruption of solar flares,therefore the two phenomena are studied in this work.It is concluded that the X-ray will not affect the phase of VLF signals at night,but the energetic particles will affect the phase change,and the influence of energetic particles should be excluded in the study of multimode interference phenomena.Using VLF signals for navigation positioning in degraded or unavailable GPS conditions is of great practical significance for VLF navigation systems as it can avoid the influence of multimode interference and improve positioning accuracy.

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.

Machine Learning-based Identification of Contaminated Images in Light Curve Data Preprocessing

Attitude is one of the crucial parameters for space objects and plays a vital role in collision prediction and debris removal.Analyzing light curves to determine attitude is the most commonly used method.In photometric observations,outliers may exist in the obtained light curves due to various reasons.Therefore,preprocessing is required to remove these outliers to obtain high quality light curves.Through statistical analysis,the reasons leading to outliers can be categorized into two main types:first,the brightness of the object significantly increases due to the passage of a star nearby,referred to as“stellar contamination,”and second,the brightness markedly decreases due to cloudy cover,referred to as“cloudy contamination.”The traditional approach of manually inspecting images for contamination is time-consuming and labor-intensive.However,we propose the utilization of machine learning methods as a substitute.Convolutional Neural Networks and SVMs are employed to identify cases of stellar contamination and cloudy contamination,achieving F1 scores of 1.00 and 0.98 on a test set,respectively.We also explore other machine learning methods such as ResNet-18 and Light Gradient Boosting Machine,then conduct comparative analyses of the results.

PSRDP:A Parallel Processing Method for Pulsar Baseband Data

To address the problem of real-time processing of ultra-wide bandwidth pulsar baseband data,we designed and implemented a pulsar baseband data processing algorithm(PSRDP)based on GPU parallel computing technology.PSRDP can perform operations such as baseband data unpacking,channel separation,coherent dedispersion,Stokes detection,phase and folding period prediction,and folding integration in GPU clusters.We tested the algorithm using the J0437-4715 pulsar baseband data generated by the CASPSR and Medusa backends of the Parkes,and the J0332+5434 pulsar baseband data generated by the self-developed backend of the Nan Shan Radio Telescope.We obtained the pulse profiles of each baseband data.Through experimental analysis,we have found that the pulse profiles generated by the PSRDP algorithm in this paper are essentially consistent with the processing results of Digital Signal Processing Software for Pulsar Astronomy(DSPSR),which verified the effectiveness of the PSRDP algorithm.Furthermore,using the same baseband data,we compared the processing speed of PSRDP with DSPSR,and the results showed that PSRDP was not slower than DSPSR in terms of speed.The theoretical and technical experience gained from the PSRDP algorithm research in this article lays a technical foundation for the real-time processing of QTT(Qi Tai radio Telescope)ultra-wide bandwidth pulsar baseband data.

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.

Analysis of the wide area differential correction for BeiDou global satellite navigation system

The regional BeiDou Satellite System, or BDS2, broadcasts a differential correction as Equivalent Satellite Clock Correction to correct both orbit and satellite clock errors. For the global BDS, or BDS3, satellite orbit and clock corrections conforming with RTCA standards will be broadcast to authorized users. The hybrid constellation and regional monitoring network pose challenges for the high precision separation of orbit and satellite clock corrections. Three correction models of kinematic,dynamic and Two-way Satellite Time Frequency Transfer(TWSTFT)-based dynamic were studied to estimate the satellite orbit and clock corrections. The correction accuracy of the three models is compared and analyzed based on the BDS observation data. Results show that the accuracies(root mean square, RMS) of dual-frequency real-time positioning for the three models are about 1.76 m, 1.78 m and 2.08 m respectively, which are comparable with the performance of WAAS and EGNOS. With dynamic corrections, the precision of Precise Point Positioning(PPP) experiments may reach about 23 cm after convergence.

Wavelet Analysis of Several Important Periodic Properties in the Relative Sunspot Numbers

We investigate the wavelet transform of yearly mean relative sunspot number series from 1700 to 2002. The curve of the global wavelet power spectrum peaks at 11-yr, 53-yr and 101-yr periods. The evolution of the amplitudes of the three periods is studied. The results show that around 1750 and 1800, the amplitude of the 53-yr period was much higher than that of the the 11-yr period, that the ca. 53-yr period was apparent only for the interval from 1725 to 1850, and was very low after 1850, that around 1750, 1800 and 1900, the amplitude of the 101-yr period was higher than that of the 11-yr period and that, from 1940 to 2000, the 11-yr period greatly dominates over the other two periods.

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.

Data-driven Seeing Prediction for Optics Telescope:from Statistical Modeling,Machine Learning to Deep Learning Techniques

Predicting seeing of astronomical observations can provide hints of the quality of optical imaging in the near future,and facilitate flexible scheduling of observation tasks to maximize the use of astronomical observatories.Traditional approaches to seeing prediction mostly rely on regional weather models to capture the in-dome optical turbulence patterns.Thanks to the developing of data gathering and aggregation facilities of astronomical observatories in recent years,data-driven approaches are becoming increasingly feasible and attractive to predict astronomical seeing.This paper systematically investigates data-driven approaches to seeing prediction by leveraging various big data techniques,from traditional statistical modeling,machine learning to new emerging deep learning methods,on the monitoring data of the Large sky Area Multi-Object fiber Spectroscopic Telescope(LAMOST).The raw monitoring data are preprocessed to allow for big data modeling.Then we formulate the seeing prediction task under each type of modeling framework and develop seeing prediction models through using representative big data techniques,including ARIMA and Prophet for statistical modeling,MLP and XGBoost for machine learning,and LSTM,GRU and Transformer for deep learning.We perform empirical studies on the developed models with a variety of feature configurations,yielding notable insights into the applicability of big data techniques to the seeing prediction task.

The first data release of LAMOST low-resolution single-epoch spectra

LAMOST Data Release 5,covering 17000 deg^(2) from-10°to 80°in declination,contains 9 million co-added low-resolution spectra of celestial objects,each spectrum combined from repeat exposure of two to tens of times during Oct 2011 to Jun 2017.In this paper,we present the spectra of individual exposures for all the objects in LAMOST Data Release 5.For each spectrum,the equivalent width of 60lines from 11 different elements are calculated with a new method combining the actual line core and fitted line wings.For stars earlier than F type,the Balmer lines are fitted with both emission and absorption profiles once two components are detected.Radial velocity of each individual exposure is measured by minimizing χ^(2) between the spectrum and its best template.The database for equivalent widths of spectral lines and radial velocities of individual spectra are available online.Radial velocity uncertainties with different stellar type and signal-to-noise ratio are quantified by comparing different exposure of the same objects.We notice that the radial velocity uncertainty depends on the time lag between observations.For stars observed in the same day and with signal-to-noise ratio higher than 20,the radial velocity uncertainty is below 5 km s^(-1),and increases to 10 km s^(-1) for stars observed in different nights.

Clustering analysis of line indices for LAMOST spectra with AstroStat

The application of data mining in astronomical surveys,such as the Large Sky Area MultiObject Fiber Spectroscopic Telescope(LAMOST)survey,provides an effective approach to automatically analyze a large amount of complex survey data.Unsupervised clustering could help astronomers find the associations and outliers in a big data set.In this paper,we employ the k-means method to perform clustering for the line index of LAMOST spectra with the powerful software Astro Stat.Implementing the line index approach for analyzing astronomical spectra is an effective way to extract spectral features for low resolution spectra,which can represent the main spectral characteristics of stars.A total of 144 340 line indices for A type stars is analyzed through calculating their intra and inter distances between pairs of stars.For intra distance,we use the definition of Mahalanobis distance to explore the degree of clustering for each class,while for outlier detection,we define a local outlier factor for each spectrum.Astro Stat furnishes a set of visualization tools for illustrating the analysis results.Checking the spectra detected as outliers,we find that most of them are problematic data and only a few correspond to rare astronomical objects.We show two examples of these outliers,a spectrum with abnormal continuum and a spectrum with emission lines.Our work demonstrates that line index clustering is a good method for examining data quality and identifying rare objects.

Radio frequency compatibility analysis for BDSBAS and GPS/WAAS

Radio frequency signal compatibility is the basis of interoperability of the Satellite eBased Augmentation System(SBAS).SBAS should abide by relative international radio regulations of International Telecommunication Union(ITU) and meet the compatibility requirements of radio frequency signal between the Global Navigation Satellite System(GNSS)/SBAS,in order to avoid negative mutual interference.According to ITU Proposal and related reference and assumptions,the paper made simulation of signal receiving maximum power in the Bei Dou Satellite-Based Augmentation System(BDSBAS)global signal coverage.And then,interference of BDSBAS to Global Positioning System(GPS)/WAAS(Wide Area Augmentation System) on L1/L5 bands were calculated and analyzed,with equivalent carrierto-noise ratio as the evaluation parameter.The result shows that the carrier-to-noise ratio decrease of GPS/WAAS caused by BDSBAS B1 C and B2 a signals are extremely lower than inter-system interference of GPS/WAAS,and thus can be ignored in practical applications.Therefore,BDSBAS will not cause the service performance degradation of GPS and WAAS.

A GPU based single-pulse search pipeline(GSP)with database and its application to the Commensal Radio Astronomy FAST Survey(CRAFTS)

We developed a GPU based single-pulse search pipeline(GSP)with a candidate-archiving database.Largely based upon the infrastructure of the open source PulsaR Exploration and Search Toolkit(PRESTO),GSP implements GPU acceleration of the de-dispersion and integrates a candidate-archiving database.We applied GSP to the data streams from the Commensal Radio Astronomy FAST Survey(CRAFTS),which resulted in quasi-real-time processing.The integrated candidate database facilitates synergistic usage of multiple machine-learning tools and thus improves efficient identification of radio pulsars such as rotating radio transients(RRATs)and fast radio bursts(FRBs).We first tested GSP on pilot CRAFTS observations with the FAST Ultra-Wide Band(UWB)receiver.GSP detected all pulsars known from the the Parkes multibeam pulsar survey in the corresponding sky area covered by the FAST-UWB.GSP also discovered 13 new pulsars.We measured the computational efficiency of GSP to be~120 times faster than the original PRESTO and~60 times faster than an MPI-parallelized version of PRESTO.

Identifications of RR Lyrae Stars and Quasars from the Simulated Data of Mephisto-W Survey

We have investigated the feasibilities and accuracies of the identifications of RR Lyrae stars and quasars from the simulated data of the Multi-channel Photometric Survey Telescope(Mephisto)W Survey.Based on the variable sources light curve libraries from the Sloan Digital Sky Survey(SDSS)Stripe 82 data and the observation history simulation from the Mephisto-W Survey Scheduler,we have simulated the uvgriz multi-band light curves of RR Lyrae stars,quasars and other variable sources for the first year observation of Mephisto W Survey.We have applied the ensemble machine learning algorithm Random Forest Classifier(RFC)to identify RR Lyrae stars and quasars,respectively.We build training and test samples and extract~150 features from the simulated light curves and train two RFCs respectively for the RR Lyrae star and quasar classification.We find that,our RFCs are able to select the RR Lyrae stars and quasars with remarkably high precision and completeness,with purity=95.4%and completeness=96.9%for the RR Lyrae RFC and purity=91.4%and completeness=90.2%for the quasar RFC.We have also derived relative importances of the extracted features utilized to classify RR Lyrae stars and quasars.

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 search for Mg-poor stars using LAMOST survey data

Most Galactic metal-poor stars exhibit enhanced α-abundances（e.g.[Mg/Fe] ～ ＋0.4） according to previous studies of stellar chemical compositions.However, a handful of metal-poor stars with large deficiencies in Mg（e.g. [Mg/Fe]～-0.2） show severe departures from this α-enhancement trend. The sub-solar[Mg/Fe] ratios of these anomalous stars indicate that they possess different chemical enrichment histories than the majority of Galactic metal-poor stars. In previous work,we presented a method to select Mg-poor metal-poor stars from low-resolution SDSS spectra based on a spectral matching technique. In this paper, a similar method is applied to low-resolution（R ～ 1800） LAMOST spectra. Stellar [Mg/Fe] abundances are determined by using stellar parameters delivered by the LAMOST Data Release2 catalog. From a sample of ～ 60 000 stars with atmospheric parameters in the range Teff = [5500, 6500] K and [Fe/H] = [-2.4, ＋0.5], we select 15 candidate Mg-poor metal-poor stars.

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.