Reminiscences of my life in science

With a great deal of humility I attempt in the following to recall important events in my life and in my scientific career of more than five decades. I am not through yet. I continue to do research though, I admit, not with the energy and fervor I once had. Still, I hope to contribute to science in meaningful ways.     

ALMA observations of the circumstellar envelope around EP Aqr

Atacama Large Millimetre/sub-millimetre Array(ALMA) observations of CO(1–0) and CO(2–1) emissions from the circumstellar envelope of the asymptotic giant branch(AGB) star EP Aqr have been made with four times better spatial resolution than previously available. They are analysed with emphasis on the de-projection in space of the effective emissivity and flux of matter using as input a prescribed configuration of the velocity field, assumed to be radial. The data are found to display an intrinsic axisymmetry with respect to an axis making a small angle with respect to the line of sight. A broad range of wind configurations, from prolate(bipolar) to oblate(equatorial) has been studied and found to be accompanied by significant equatorial emission. Qualitatively, the effective emissivity is enhanced near the equator to produce the central narrow component observed in the Doppler velocity spectra and its dependence on star latitude generally follows that of the wind velocity with the exception of an omni-present depression near the poles. In particular, large equatorial expansion velocities produce a flared disc or a ring of effective emissivity and mass loss. The effect on the determination of the orientation of the star axis of radial velocity gradients, and possibly competing rotation and expansion in the equatorial disc, is discussed. In general,the flux of matter is found to reach a broad maximum at distances of the order of 500 AU from the star.Arguments are given that may be used to favour one wind velocity distribution over another. As a result of the improved quality of the data, a deeper understanding of the constraints imposed on morphology and kinematics has been obtained.     

The symmetry energy and incompressibility constrained by the observations of glitching pulsars

We investigate the masses of glitching pulsars in order to constrain their equation of state(EOS). The observations of glitches(sudden jumps in rotational frequency) may provide information on the interior physics of neutron stars. With the assumption that glitches are triggered by superfluid neutrons, the masses of glitching neutron stars can be estimated using observations of maximum glitches.Together with the observations of thermal emission from glitching pulsars Vela and J1709–4429, the slope of symmetry energy and incompressibility of nuclear matter at saturation density can be constrained.The slope of symmetry energy L should be larger than 67 MeV while the lower limit of incompressibility for symmetric nuclear matter K_0 is 215 MeV. We also obtain a relationship between L and K_0:6.173 MeV + 0.283 K_0≤ L ≤ 7.729 MeV + 0.291 K_0. The restricted EOSs are consistent with the observations of 2-solar-mass neutron stars and gravitational waves from a binary neutron star inspiral.     

CCD astrometric observations of 2017 VR12, Camillo and Midas

We have observed three near-Earth objects(NEOs), 2017 VR12, Camillo and Midas, during 2018.The observations were made with the 1-m telescope, operated by Yunnan Observatories, over two nights.Their precise astrometric positions are derived from 989 CCD observations. The theoretical positions of asteroids are retrieved from the Jet Propulsion Laboratory(JPL) Horizons System and Institut de M′ecanique C′eleste et de Calcul des ′Eph′em′erides(IMCCE). The positions of three asteroids are measured with respect to stars in the Gaia DR2 star catalog. For 2017 VR12, the means(O-C) of right ascension and declination are -0.090′′and-0.623′′respectively based on the published JPL ephemeris, but the corresponding means(O-C) are 3.122′′and-0.636′′based on the published IMCCE ephemeris. The great difference in declination could be explained by several factors.(1) The degraded CCD images caused by the fast apparent motion of the objects lead to a reduction in positioning accuracy.(2) The poor timing system may introduce systematic errors, especially in the high speed direction.(3) The asteroid may be perturbed by Earth when it approaches the Earth too closely. These astrometric results demonstrate that the centroid centering method can reduce the dispersion of non-Gaussian images as compared with the PSF modeling method.For Camillo and Midas, the astrometric results are consistent based on the two ephemerides. Implementing a high-precision timing system, and analyzing some astronomical effects and geometric distortions in CCD images should be carefully considered in future works.     

Transition-region explosive events produced by plasmoid instability

Magnetic reconnection is thought to be a key process in most solar eruptions. Thanks to highresolution observations and simulations, the studied scale of the reconnection process has become smaller and smaller. Spectroscopic observations show that the reconnection site can be very small, which always exhibits a bright core and two extended wings with fast speeds, i.e., transition-region explosive events.In this paper, using the PLUTO code, we perform a 2-D magnetohydrodynamic simulation to investigate small-scale reconnection in double current sheets. Based on our simulation results, such as the line-of-sight velocity, number density and plasma temperature, we can synthesize the line profile of SiIV 1402.77? which is a well known emission line used to study transition-region explosive events on the Sun. The synthetic line profile of Si IV 1402.77? is complex with a bright core and two broad wings which can extend to nearly 200 km s^-1. Our simulation results suggest that the transition-region explosive events on the Sun are produced by plasmoid instability during small-scale magnetic reconnection.     

Numerical simulation and preliminary analysis of typhoon waves during three typhoons in the Yellow Sea and East China Sea

In this study,typhoon waves generated during three typhoons(Damrey(1210),Fung-wong(1416),and Chan-hom(1509))in the Yellow Sea and East China Sea were simulated in a simulating waves nearshore(SWAN)model,and the wind forcing was constructed by combining reanalyzed wind data with a Holland typhoon wind model.Various parameters,such as the Holland fitting parameter(B)and the maximum wind radius?,were investigated in sensitivity experiments in the Holland model that affect the wind field construction.Six different formulations were considered and the parameters determined by comparing the simulated wind results with in-situ wind measurements.The key factors affecting wave growth and dissipation processes from deep to shallow waters were studied,including wind input,whitecapping,and bottom friction.Comparison with in-situ wave measurements suggested that the KOMEN scheme(wind input exponential growth and whitecapping energy dissipation)and the JONSWAP scheme(dissipation of bottom friction)resulted in good reproduction of the significant wave height of typhoon waves.A preliminary analysis of the wave characteristics in terms of wind-sea and swell wave revealed that swell waves dominated with the distance of R to the eye of the typhoon,while wind-sea prevailed in the outer region up to six to eight times the R values despite a clear misalignment between wind and waves.The results support the hypothesis that nonlinear wave-wave interactions may play a key role in the formation of wave characteristics.     

Research on the Method of De-noising the Short-wave Time Signal at Low SNR

A voice enhancement algorithm based on the Empirical Mode Decomposition (EMD) and the improved spectral subtraction is proposed for the low-SNR (Signal Noise Ratio) shortwave time signal. This method is proposed to solve the problem that the shortwave time signal cannot be used for timing in complex noisy environments. The core idea of this method is to use the Hilbert-Huang Transform (HHT) algorithm to make the empirical mode decomposition on the noisy shortwave signal, and to select the intrinsic mode functions containing the shortwave signal information for the signal reconstruction by through the maximum correlation. Then, to make the spectral subtraction on the reconstructed signal to achieve the purpose of noise reduction. The experimental result shows that this method has a better noise reduction than the traditional methods.     

Identification of dissolved and particulate carbonyl compounds produced by marine harmful algal bloom species

Carbonyl compounds, especially polyunsaturated aldehydes (PUAs), are increasingly reported as teratogenic to grazers and deleterious to phytoplankton. While PUAs have been considered to be mainly produced by diatoms after cell wounding, little is known about whether microalgae other than diatoms, particularly harmful algal bloom (HAB)-forming species, produce carbonyl compounds. In this study, we analyzed the carbonyl compounds from eight common HAB-forming species ( Akashiwo sanguinea , Karenia mikimotoi , Karlodinium venefi cum , Margalefi dinium polykrikoides , Prorocentrum donghaiense , P . minimum , Scrippsiella trochoidea , and Heterosigma akashiwo ) using gas chromatography-mass spectrometry in full scan and selected ion monitoring (SIM) modes. Our results show a ubiquitous presence of carbonyl compounds in both dissolved and particulate forms in all species we examined. In the full scan mode, 133 chromatographic peaks were detected from 51 samples altogether. Both the varieties of carbonyl compounds and their quantities were algal species dependent, although most of the carbonyl compounds could not be fully identifi ed according to the mass spectral database only due to the unavailability of enough standards for all analytes. Aided with nine standards in SIM mode, we further identifi ed and quantifi ed all nine aldehydes (2-methyl-2-pentenal, trans-2-nonenal, cis-6-nonenal, 2,6-dimethyl-5-heptenal, trans-2-hexenal, trans-2- decenal, 2,4-heptadienal, trans-trans-2,4-octadienal, and trans-trans-2,4-decadienal). Four of these nine aldehydes were detected in particulate form, which confi rmed that trans-2-nonenal could be produced by K . mikimotoi , Karl . venefi cum , P . donghaiense , P . minimum , S . trochoidea , and H . akashiwo , 2,4-heptadienal and trans-trans-2,4-decadienal by A . sanguinea , M . polykrikoides , and S . trochoidea , and trans, trans-2,4- octadienal by S . trochoidea , respectively. We proved that some dinofl agellate and raphidophyte species could contribute to the pool of carbonyl compounds including PUAs in a marine ecosystem. Some carbonyl compounds, particularly those in high cell quota and/or dissolved concentration, may be related to fi shkilling or allelopathy which needs further identifi cation and quantifi cation.     

Application and comparison of two symptom-based eutrophication-assessment methods in Jiaozhou Bay,China

Increasing numbers of approaches to assess eutrophication, such as estuarine trophic status, the Oslo-Paris Commission Common Procedure, and the Water Framework Directive, focus on the symptoms of eutrophication. In China, however, nutrient index methods dominate the assessment of coastal waters. In this study, an integrated method that includes both water quality and ecological response was compared with the Northwest Pacifi c Action Plan (NOWPAP) Common Procedure. Observation data from Jiaozhou Bay, Shandong, China, were used in a comparison of the two methods in a trophic status study. Overall, both clearly revealed a high level of nutrient enrichment in the bay, indicated by high nutrient concentrations. Though the two methods diff ered in their methodological design in the assessment of the ecological eff ects of nutrient enrichment, they have acquired similar results: the integrated method suggested that the status was good, and the NOWPAP Common Procedure suggested that the status was low (indicating that the bay had no serious eutrophication problem). The introduction of fi lter feeders (shellfi sh aquaculture) into the bay on a reasonable scale may have been eff ective in mitigating eutrophic conditions, and perhaps explains the low ecological impacts there. Our results will be useful to ecosystem-based eutrophication management in the bay and in similar areas.     

机载WIDAS数据的Landsat卫星反照率初步验证

随着精细化监测的需求,中高空间分辨率的地表反照率产品逐渐成为气候模型的主要输入。目前,中高空间分辨率反照率产品的验证主要基于地表站点的通量塔观测数据,区域机载飞行数据的验证依然相对较少。因此,本文基于区域机载数据验证Landsat反照率产品。针对内蒙古自治区根河森林试验区所获取的机载红外广角双模式成像仪(WIDAS)多角度反射率数据,应用BRDF原型反演算法估算其反照率,分析了应用机载数据验证中高空间分辨率反照率产品的潜力。2016年内蒙古根河森林试验区机载WIDAS飞行多角度观测的可用多角度范围为25°,以前的研究表明BRDF原型反照率反演算法表现出对小观测角度的反照率反演结果的鲁棒性。因此,机载WIDAS反照率在一定程度可用于星载反照率的验证。首先,基于核驱动模型和各向异性平整指数(AFX)提取了试验区4种MODIS二向性反射分布函数(BRDF)原型;然后,将其作为先验知识应用到根河森林WIDAS机载数据的反照率反演中;最后,用WIDAS反照率和单个地面通量塔观测的反照率对Landsat卫星数据的反照率进行初步验证。验证结果表明Landsat反照率与WIDAS反照率结果较为一致,但略有低估,总体均方根误差(RMSE)约为0.02,偏差为0.0057。在多角度观测范围较小时,BRDF原型的反照率反演算法可为星载地表反照率的验证提供了一种有效的验证手段。