On September 10 2017, a fast coronal mass ejection (CME) erupted from the active region (AR) 12673, leading to a ground level enhancement (GLE) event at Earth. Using the 2D improved Particle Acceleration and Transport in the Heliosphere (iPATH) model, we model the large solar energetic particle (SEP) event of 10 September 2017 observed at Earth, Mars and STEREO-A. Based on observational evidence,we assume that the CME-driven shock experienced a large lateral expansion shortly after the eruption, which is modelled by a double Gaussian velocity profile in this simulation. We use the in-situ shock arrival times and the observed CME speeds at multiple spacecraft near Earth and Mars as constraints to adjust the input model parameters. The modelled time intensity profiles and fluence for energetic protons are then compared with observations. Reasonable agreements with observations at Mars and STEREO-A are found. The simulated results at Earth differ from observations of GOES-15. Instead, the simulated results at a heliocentric longitude 20 degree west to Earth fit reasonably well with the GOES observation. This can be explained if the pre-event solar wind magnetic field at Earth is not described by a nominal Parker field. Our results suggest that a large lateral expansion of the CME-driven shock and a distorted interplanetary magnetic field due to previous events can be important in understanding this GLE event.

An observation-based Galactic Cosmic Ray (GCR) spectral model for heavy nuclei is developed. Zhao and Qin (J. Geophys. Res. Space Phys.118, 1837(2013)) proposed an empirical elemental GCR spectra model for nuclear charge 5-28 over the energy range from 30 to 500 MeV/nuc, which is proved to be successful in predicting yearly averaged GCR heavy nuclei spectra.Based on the latest highly statistically precise measurements from ACE/CRIS,a further elemental GCR model with monthly averaged spectra is presented. The model can reproduce the past and predict the futureGCR intensity monthly by correlating model parameters with thecontinuous sunspot number (SSN) record. The effects of solar activity on GCR modulation are considered separately for odd and even solar cycles. Compared with other comprehensive GCR models, our modeling results are satisfyingly consistent with the GCR spectral measurements from ACE/SIS and IMP-8, and have comparable prediction accuracy as the Badhwar & O'Neill 2014 model.A detailed error analysis is also provided.Finally, the GCR carbon and iron nuclei fluxes for the subsequent two solar cycles (SC 25 and 26) are predicted and they show a potential trend in reduced flux amplitude, which is suspected to be relevant to possible weak solar cycles.

Recent detection of superflares on solar-type stars by Kepler mission raised a possibility that they can be associated with energetic coronal mass ejections (CMEs) and energetic particle events (SEPs). These space weather events can impact habitability of exoplanets around these stars. Here we use the improved Particle Acceleration and Transport in the Heliosphere (iPATH) model, to model the time intensity profile and spectrum of SEPs accelerated at CME-driven shocks from stars of different ages traced by their rotation rates. We consider a solar-like (G-type) star with 6 different rotation rates varying from 0.5 to 3.0 times rotation speed of the sun. In all 6 cases, a fast CME is launched with the same speed of 1500 km/sec and the resulting time intensity profiles at 3 locations and and energy spectra at 5 locations at 1 AU are obtained. The maximum particle energy at the shock front as a function of r is also shown. Our results suggest that within 0.8 AU the maximum particle energy at the shock front increases with the rotation rate of the star. However, event integrated spectra for the five selected locations along the CME path show complicated patterns. This is because the Parker magnetic field for rapidly rotating stars is more tightly winded. Our results can be used in estimating the radiation environments of terrestrial-type exoplanets around solar-type stars.

We explore methods for robust estimation of the 21 cm signal from the Epoch of Reionisation (EoR). A Kernel Density Estimator (KDE) is introduced for measuring the spatial temperature fluctuation power spectrum from the EoR. The KDE estimates the underlying probability distribution function of fluctuations as a function of spatial scale, and contains different systematic biases and errors to the typical approach to estimating the fluctuation power spectrum. Extraction of histograms of visibilities allows moments analysis to be used to discriminate foregrounds from 21 cm signal and thermal noise. We use the information available in the histograms, along with the statistical dis-similarity of foregrounds from two independent observing fields, to robustly separate foregrounds from cosmological signal, while making no assumptions about the Gaussianity of the signal. Using two independent observing fields to robustly discriminate signal from foregrounds is crucial for the analysis presented in this paper. We apply the techniques to 13 hours of Murchison Widefield Array (MWA) EoR data over two observing fields. We compare the output to that obtained with a comparative power spectrum estimation method, and demonstrate the reduced foreground contamination using this approach. Using the second moment obtained directly from the KDE distribution functions yields a factor of 2-3 improvement in power for k < 0.3hMpc^{-1} compared with a matched delay space power estimator, while weighting data by additional statistics does not offer significant improvement beyond that available for thermal noise-only weights.

Understanding the sources of lunar water is crucial for studying the history of lunar evolution, and also the solar wind interaction with the Moon and other airless bodies. Recent observations revealed lunar hydration is very likely a surficial dynamic process driven by solar wind. Solar wind is shielded over a period of 3-5 days as the Moon passes through the Earth's magnetosphere, during which a significant loss of hydration is expected from previous works.Here we study lunar hydration inside the magnetosphere using orbital spectral data, which unexpectedly found that the polar surficial OH/H2O abundance remains at the same level when in the solar wind and in the magnetosphere. We suggest that particles from the magnetosphere (Earth wind, naturally different from solar wind) contribute to lunar hydration. From lunar orbital plasma observations, we find the existence of optimal energy ranges, other than 1 keV as previously thought, for surface hydration formation. These optimal energy ranges deduced from space observations may provide strong implications for laboratory experiments simulating lunar hydration processes.

We investigate the possibility that the dwarf galaxies Crater II and Hercules have previously been tidally stripped by the Milky Way. We present Magellan/IMACS spectra of candidate member stars in both objects. We identify 37 members of Crater II, 25 of which have velocity measurements in the literature, and we classify 3 stars within that subset as possible binaries. We find that including or removing these binary candidates does not change the derived velocity dispersion of Crater II. Excluding the binary candidates, we measure a velocity dispersion of $\sigma_{V_{los}} = 2.7^{+0.5}_{-0.4}$ \kms, corresponding to $M/L = 47^{+17}_{-13}$ M$_{\odot}/$L$_{\odot}$. We measure a mean metallicity of [Fe/H] = $-1.95^{+0.06}_{-0.05}$, with a dispersion of $\sigma_{\mbox{[Fe/H]}} = 0.18^{+0.06}_{-0.08}$. Our velocity dispersion and metallicity measurements agree with previous measurements for Crater II, and confirm that the galaxy resides in a kinematically cold dark matter halo. We also search for spectroscopic members stripped from Hercules in the possible extratidal stellar overdensities surrounding the dwarf. For both galaxies, we calculate proper motions using \textit{Gaia} DR2 astrometry, and use their full 6D phase space information to evaluate the probability that their orbits approach sufficiently close to the Milky Way to experience tidal stripping. Given the available kinematic data, we find a probability of $\sim$40% that Hercules has suffered tidal stripping. The proper motion of Crater II makes it almost certain to be stripped.

We present a search for new members of the 300 km s$^{-1}$ stream (300S) near the dwarf galaxy Segue 1 using wide-field survey data. We identify 11 previously unknown bright stream members in the APOGEE-2 and SEGUE-1 and 2 spectroscopic surveys. Based on the spatial distribution of the high-velocity stars, we confirm for the first time that this kinematic structure is associated with a 24$^{\circ}$-long stream seen in SDSS and Pan-STARRS imaging data. The 300S stars display a metallicity range of $-2.17 < {\rm [Fe/H]} < -1.24$, with an intrinsic dispersion of 0.21$_{-0.09}^{+0.12}$ dex. They also have chemical abundance patterns similar to those of Local Group dwarf galaxies, as well as that of the Milky Way halo. Using the open-source code galpy to model the orbit of the stream, we find that the progenitor of the stream passed perigalacticon about 70 Myr ago, with a closest approach to the Galactic Center of about 4.1 kpc. Using Pan-STARRS DR1 data, we obtain an integrated stream luminosity of $4 \times 10^3$ L$_{\odot}$. We conclude that the progenitor of the stream was a dwarf galaxy that is probably similar to the satellites that were accreted to build the present-day Milky Way halo.

In a two-component jet model, the emissions are the sum of the core and extended emissions: $S^{\rm ob}=S_{\rm core}^{\rm ob}+S_{\rm ext}^{\rm ob}$, with the core emissions, $S_{\rm core}^{\rm ob}= f S_{\rm ext}^{\rm ob}\delta^{q}$, being a function of the Doppler factor, $\delta$, the extended emission, $S_{\rm ext}^{\rm ob}$, jet type dependent factor, $q$, and the ratio of the core to the extended emissions in the comoving frame, $f$. The $f$ is an unobservable but important parameter. Following our previous work, we collect 65 blazars with available Doppler factor, $\delta$, superluminal velocity, $\beta_{app}$, and core-dominance parameter, $R$, calculate the ratio, $f$, and peform statistical analyses. We find that the ratio, $f$, in BL Lacs is on average larger than that in FSRQs. We suggest that the difference of the ratio $f$ between FSRQs and BL Lacs is one of the possible reasons that cause the difference of other observed properties between them. We also find some significant correlations between $\log f$ and other parameters, including intrinsic (de-beamed) peak frequency, $\log \nu _{\rm p}^{\rm in}$, intrinsic polarization, $\log P^{\rm in}$, and core-dominance parameter, $\log R$, for the whole sample. In addition, we show that the ratio, $f$, can be estimated by $R$.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is a co-planar interferometer array operating at a wavelength of 3mm to measure the Sunyaev-Zeldovich effect (SZE) of galaxy clusters. In the first phase of operation -- with a compact 7-element array with 0.6m antennas (AMiBA-7) -- we observed six clusters at angular scales from 5\arcmin to 23\arcmin. Here, we describe the expansion of AMiBA to a 13-element array with 1.2m antennas (AMiBA-13), its subsequent commissioning, and our cluster SZE observing program. The most important changes compared to AMiBA-7 are (1) array re-configuration with baselines ranging from 1.4m to 4.8m covering angular scales from 2\arcmin to 11.5\arcmin, (2) thirteen new lightweight carbon-fiber-reinforced plastic (CFRP) 1.2m reflectors, and (3) additional correlators and six new receivers. From the AMiBA-13 SZE observing program, we present here maps of a subset of twelve clusters. In highlights, we combine AMiBA-7 and AMiBA-13 observations of Abell 1689 and perform a joint fitting assuming a generalized NFW pressure profile. Our cylindrically integrated Compton-y values for this cluster are consistent with the BIMA/OVRA, SZA, and Planck results. We report the first targeted SZE detection towards the optically selected galaxy cluster RCS J1447+0828, and we demonstrate the ability of AMiBA SZE data to serve as a proxy for the total cluster mass. Finally, we show that our AMiBA-SZE derived cluster masses are consistent with recent lensing mass measurements in the literature.

We present a new way to solve the platform deformation problem of co-planar interferometers. The platform of a co-planar interferometer can be deformed due to driving forces and gravity. A deformed platform will induce extra components into the geometric delay of each baseline, and change the phases of observed visibilities. The reconstructed images will also be diluted due to the errors of the phases. The platform deformations of The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) were modelled based on photogrammetry data with about 20 mount pointing positions. We then used the differential optical pointing error between two optical telescopes to fit the model parameters in the entire horizontal coordinate space. With the platform deformation model, we can predict the errors of the geometric phase delays due to platform deformation with given azimuth and elevation of the targets and calibrators. After correcting the phases of the radio point sources in the AMiBA interferometric data, we recover 50% - 70% flux loss due to phase errors. This allows us to restore more than 90% of a source flux. The method outlined in this work is not only applicable to the correction of deformation for other co-planar telescopes but also to single dish telescopes with deformation problems. This work also forms the basis of the upcoming science results of AMiBA-13.