Source author record

S. Shibata

S. Shibata appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2026arXiv

Observation of the Forbush decrease on 2024 May 10, using the ALPAQUITA air-shower array at the 70-1000 GV rigidity range

The Andes Large area PArticle detector for Cosmic ray and Astronomy (ALPACA) is a new air-shower array experiment under construction in the Bolivian Andes, and its prototype ALPAQUITA surface array has been operating since 2023 April. In addition to the traditional $\ge$3-hit or $\ge$4-hit coincidences to trigger recording air-shower events, ALPAQUITA records the counting rates of the $\ge$1-hit and $\ge$2-hit events (Any1 and Any2, respectively). We report a successful detection of a Forbush decrease occurred on 2024 May 10 caused by a passage of an interplanetary shock formed ahead of the Interplanetary Coronal Mass Ejection. The amplitude detected in the Any1 rate is 4.26$\pm$0.33% at the median primary rigidity of 76GV which is consistent with the observations with the worldwide neutron monitor and muon detector networks. Under the assumption of a power-law rigidity spectrum, we renormalized the errors of the observed amplitude ($A_{obs}$) and fitted them as a function of the median primary rigidity ($R_{m}$) of each detector and observational method. The result $A_{obs} = (10.9\% \pm 0.9\%) \times (R_{m}/10\,GV)^{-0.55 \pm 0.07}$ exhibits a hard nature of this event. Our non-detection in the Any2 rate decrease constrains the amplitude with a 2$σ$ upper limit to be 0.95% at 960GV. This marginally suggests an existence of a spectral softening between 100GV and 1000GV as also suggested by the Misato underground muon detector at 145GV. Although a strong geomagnetic storm was observed during this period, we conclude it does not impact our results. Our novel technique realizes a unique coverage to study the behavior of the Forbush decreases at the highest rigidity.

preprint2022arXiv

Proton Penetration Efficiency over a High Altitude Observatory in Mexico

In association with a large solar flare on November 7, 2004, the solar neutron detectors located at Mt. Chacaltaya (5,250m) in Bolivia and Mt. Sierra Negra (4,600m) in Mexico recorded very interesting events. In order to explain these events, we have performed a calculation solving the equation of motion of anti-protons inside the magnetosphere. Based on these results, the Mt. Chacaltaya event may be explained by the detection of solar neutrons, while the Mt. Sierra Negra event may be explained by the first detection of very high energy solar neutron decay protons (SNDPs) around 6 GeV.

preprint2019arXiv

The angular resolution of GRAPES-3 EAS array after correction for the shower front curvature

The angular resolution of an extensive air shower (EAS) array plays a critical role in determining its sensitivity for the detection of point $γ$-ray sources in the multi-TeV energy range. The GRAPES-3 an EAS array located at Ooty in India (11.4$^{\circ}$N, 76.7$^{\circ}$E, 2200 m altitude) is designed to study $γ$-rays in the TeV-PeV energy range. It comprises of a dense array of 400 plastic scintillators deployed over an area of 25000 m$^2$ and a large area (560 m$^2$) muon telescope. A new statistical method allowed real time determination of the propagation delay of each detector in the GRAPES-3 array. The shape of shower front is known to be curved and here the details of a new method developed for accurate measurement of the shower front curvature is presented. These two developments have led to a sizable improvement in the angular resolution of GRAPES-3 array. It is shown that the curvature depends on the size and age of an EAS. By employing two different techniques, namely, the odd-even and the left-right methods, independent estimates of the angular resolution are obtained. The odd-even method estimates the best achievable resolution of the array. For obtaining the angular resolution, the left-right method is used after implementing the size and age dependent curvature corrections. A comparison of the angular resolution as a function of EAS energy by these two methods shows them be virtually indistinguishable. The angular resolution of GRAPES-3 array is 47$^{\prime}$ for energies E$>$5 TeV and improves to 17$^{\prime}$ at E$>$100 TeV and finally approaching 10$^{\prime}$ at E$>$500 TeV.

preprint2018arXiv

Was the cosmic ray burst detected by the GRAPES-3 on 22 June 2015 caused by transient weakening of geomagnetic field or by an interplanetary anisotropy?

The GRAPES-3 muon telescope in Ooty, India had claimed detection of a 2 hour (h) high-energy ($\sim$20 GeV) burst of galactic cosmic-rays (GCRs) through a $>$50$σ$ surge in GeV muons, was caused by reconnection of the interplanetary magnetic field (IMF) in the magnetosphere that led to transient weakening of Earth's magnetic shield. This burst had occurred during a G4-class geomagnetic storm (storm) with a delay of $\frac{1}{2}$h relative to the coronal mass ejection (CME) of 22 June 2015 (Mohanty et al., 2016). However, recently a group interpreted the occurrence of the same burst in a subset of 31 neutron monitors (NMs) to have been the result of an anisotropy in interplanetary space (Evenson et al., 2017) in contrast to the claim in (Mohanty et al., 2016). A new analysis of the GRAPES-3 data with a fine 10.6$^{\circ}$ angular segmentation shows the speculation of interplanetary anisotropy to be incorrect, and offers a possible explanation of the NM observations. The observed 28 minutes (min) delay of the burst relative to the CME can be explained by the movement of the reconnection front from the bow shock to the surface of Earth at an average speed of 35 km/s, much lower than the CME speed of 700 km/s. This measurement may provide a more accurate estimate of the start of the storm.

preprint2016arXiv

A photon observed by the SEDA-FIB detector on the arrival of the gravitational wave

We reported a photon detection at the arrival time of the gravitational wave on December 26th, 2015. According to the LIGO-Virgo collaboration, the gravitational wave was produced by the merging process of the two black holes. The merged time was determined at 03:38:53.6 UT on December 26th. At 03:38:54.05 GPS time, one of the detectors of SEDA-FIB on board the International Space Station (ISS) detected a photon, arriving from the direction of Corona Australis near the Galactic Bulge (GB). The energy of the photon was about (35+-7) MeV. A 3.0 sigma level detection of the photon is discussed.

preprint2016arXiv

X-ray and Rotational Luminosity Correlation and Magnetic Heating of the Radio Pulsars

Previous works have suggested a correlation between the X-ray luminosity Lx and the rotational luminosity Lrot of the radio pulsars.However, none of the obtained regression lines are statistically acceptable due to large scatters. We construct a statistical model which has an intrinsic Lx-Lrot relation and reproduces the observed Lx distribution about it by using a Monte Carlo simulator, which takes into account the effects obscuring the intrinsic relation,i.e., the anisotropy of radiation, additional heating, uncertainty in distance and detection limit of the instruments. From the ATNF pulsar catalog we collect 57 `ordinary radio pulsars' with significant detection and 42 with upper limits.The sample does not include the high-magnetic field pulsars (>10^{13} G), which are separately analyzed. We obtain a statistically acceptable relation Lx (0.5 - 10 keV)= 10^{31.69} (Lrot / L_0)^{c_1} with c_1 = 1.03 \pm 0.27 and L_0 =10^{35.38}. The distribution about the obtained Lx-Lrot relation is reproduced well by the simulator. Pulsars with abnormally high Lx fall into two types: one is the soft gamma-ray pulsars, and the other is thermally bright pulsars in comparison with the standard cooling curve. On the other hand, pulsars showing low Lx are found to have dim pulsar wind nebulae. We argue that there is an unknown mechanism that governs both the magnetospheric emission and the pulsar wind nebulae, and it might involve the production rate of electron-positron pairs. The high-field pulsars form a distinctive population other than the ordinary pulsars due to their excess luminosities.

preprint2015arXiv

Simultaneous Observation of Solar Neutrons from the ISS and High Mountain Observatories in association with a flare on July 8, 2014

An M6.5-class flare was observed at N12E56 of the solar surface at 16:06 UT on July 8, 2014. In association with this flare, solar neutron detectors located on two high mountains, Mt. Sierra Negra and Chacaltaya and at the space station observed enhancements in the neutral channel. The authors analysed these data and a possible scenario of enhancements produced by high-energy protons and neutrons is proposed, using the data from continuous observation of a solar surface by the ultraviolet telescope onboard the Solar Dynamical Observatory (SDO).

preprint2015arXiv

The Search for >35 MeV Neutrons from the June 3, 2012 Impulsive Flare

We analyzed a highly impulsive solar flare observed on June 3, 2012. In association with this flare, emissions of hard X-rays, high-energy gamma rays, and neutrons were detected by the detectors onboard the FERMI, RHESSI satellites and the International Space Station. We compared those results with the pictures taken by the UV telescope onboard the Solar Dynamics Observatory satellite and found the crossing structure of two magnetic ropes at two positions on the solar surface almost at the same time. High-energy gamma rays were detected by the Fermi Large Area Telescope satellite, implying that the impulsive flare was one of a major source of proton acceleration processes on the solar surface. At the beginning of research, impulsive solar flares were considered to be the main source of particle acceleration processes; our current observations have confirmed this hypothesis.

preprint2013arXiv

A Possible Detection of Solar Gamma-Rays by the Ground Level Detector

On March 7, 2011 from 19:48:00 to 20:03:00 UT, the solar neutron telescope located at Mt. Sierra Negra, Mexico (4,600m) observed a 8.8sigma enhancement. In this paper, we would like to try to explain this enhancement by a hypothesis that a few GeV gamma-rays arrived at the top of the mountain produced by the Sun. We postulate that protons were accelerated at the shock front. They precipitate at the solar surface and produced those gamma-rays. If hypothesis is confirmed, this enhancement is the first sample of GeV gamma-rays observed by a ground level detector.

preprint2013arXiv

Measurement of High-Energy Solar Neutrons with SEDA-FIB onboard the ISS

A new type of solar neutron detector (SEDA-FIB) was launched on board the Space Shuttle Endeavor on July 16 2009, and began collecting data at the International Space Station (ISS) on August 25 2009. This paper summarizes four years of observations with the solar neutron detector SEDA-FIB (Space Environment Data Acquisition using the FIBer detector). The solar neutron detector FIB can determine both the energy and arrival direction of solar neutrons. In this paper, we first present the angular distribution of neutron induced protons obtained in Monte Carlo simulations. The results are compared with the experimental results. Then we provide the angular distribution of background neutrons during one full orbit of the ISS (90 minutes). Next, the angular distribution of neutrons during the flare onset time from 20:02 to 20:10 UT on March 7 2011 is presented. It is compared with the distribution when a solar flare is not occurring. Observed solar neutrons possibly originated from the M-class solar flares that occurred on March 7 (M3.7), June 7 (M2.5), September 24 (M3.0) (weak signal) and November 3 (X1.9) of 2011 and January 23 of 2012 (M8.7). This marked the first time that neutrons have been observed from M-class solar flares. A possible interpretation of the neutron production process will be also provided.

preprint2013arXiv

Solar Activity during Two Millennia as Estimated from Annual Tree Rings

The relationship between solar activity and the global climate is not only an academically interesting issue, but also an important problem for human beings. Lean and Rind have analyzed a considerable amount of climate data from around the world from 1889 to 2006. According to their analysis, the global effect was estimated to be 0.17 plus-minus 0.01K between the solar maximum and minimum. However, they noticed that the effect strongly appeared in the zones between 70N and 30N, and between 25S and 50S. At its peak latitude (near 40), the effect was estimated to be 0.5 K. Therefore, we analyzed a tree that survived at the Murooji temple in Nara prefecture (34d32mN, 136d2mE) Japan, for 391 years. Quite surprisingly, Fourier analysis of the annual growth rate identified two cycles with periodicities of 12 and 25 years during the Maunder minimum. We have continued the analysis, using different samples from the Nagusa shrine in Hyogo Prefecture, Itayanagi City, Aomori Prefecture, Japan, and from Yaku Island in Kyusyu, Japan. An evidence of solar activity was found in all of the samples. In particular, we have discovered a correlation between Swiss glacier fluctuation and the growth rate of the Yaku tree ring.

preprint2012arXiv

Measurement by FIB on the ISS: Two Emissions of Solar Neutrons Detected?

A new type of solar neutron detector (FIB) was launched onboard the Space Shuttle Endeavour on July 16, 2009, and it began collecting data at the International Space Station (ISS) on August 25, 2009. This paper summarizes the three years of observations obtained by the solar neutron detector FIB until the end of July 2012. The solar neutron detector FIB can determine both the energy and arrival direction of neutrons. We measured the energy spectra of background neutrons over the SAA region and elsewhere, and found the typical trigger rates to be 20 counts/sec and 0.22 counts/sec, respectively. It is possible to identify solar neutrons to within a level of 0.028 counts/sec, provided that directional information is applied. Solar neutrons were observed in association with the M-class solar flares that occurred on March 7 (M3.7) and June 7 (M2.5) of 2011. This marked the first time that neutrons were observed in M-class solar flares. A possible interpretaion of the prodcution process is provided.

preprint2006arXiv

Discovery of a Featureless X-Ray Spectrum in the Supernova Remnant Shell of G330.2+1.0

We report here on the first pointed X-ray observation of the supernova remnant (SNR) G330.2+1.0. The X-ray morphology is characterized by an extended shell. Its X-ray spectrum is well represented by a single power-law function with a photon index of $γ\simeq 2.8$ and interstellar absorption of $n_{\rm H}\simeq2.6\times 10^{22}$[cm$^{-2}$]. We interpret this emission as synchrotron radiation from accelerated electrons at the SNR shock, as seen in SN 1006. The surface brightness of the X-ray emission is anti-correlated with the radio emission, and the power-law spectrum is dominated at the western shell where the radio emission is weak. The co-existence of two distinct (radio bright/X-ray faint and radio faint/X-ray bright) shells in a single supernova remnant challenges our understanding of the particle acceleration and radiation mechanisms in different interstellar environments. The object may be a good target for searching TeV gamma-rays and molecular gas surrounding the blast shock. We also report on the nature of a bright point-like source (AX J1601-5143) to the south of the SNR.

preprint2004arXiv

A pulsar outer gap model with trans-field structure

We investigate the electrodynamics of an outer gap in the meridional plane of the aligned-rotator. The charge depletion from the Goldreich-Julian charge density causes a large electric field along the magnetic field line. The electrons or the positrons are accelerated by the field-aligned electric field and radiate the $γ$-rays tangentially to the local magnetic field line. Some of such $γ$-rays collide with $X$-rays to materialize as the electron-positron pairs on different field lines from the field line on which they were emitted. As a result, the electric field structure is expected to change across the field lines. Including these trans-field effects, we solve the formation of the electric field self-consistently with the curvature radiation and the pair creation processes. The $γ$-ray emission and the pair creation are treated by use of Monte Carlo technique. We demonstrate that the distribution of the electric field along the field lines is affected by both the gap geometry and the external currents coming into the gap through the boundaries. In the electrodynamical model, it has been known that the solution disappears if the current density carried by the electron-positron pairs produced in the gap exceeds a critical value. We show that the critical current density is significantly increased when the trans-field structure is taken into account. We also find that the location of the inner boundary of the gap shifts toward the stellar surface from the conventional null surface as the current density increases. The reason for the shift is derived from the stability condition of the inner boundary. We also argue that the ideal-MHD condition holds outside of the gap only when the low energy particles coexist with the high energy particles migrating from the gap.