We report on the highest energy photons from the Crab Nebula observed by the Tibet air shower array with the underground water-Cherenkov-type muon detector array. Based on the criterion of muon number measured in an air shower, we successfully suppress 99.92% of the cosmic-ray background events with energies $E>100$ TeV. As a result, we observed 24 photon-like events with $E>100$ TeV against 5.5 background events, which corresponds to 5.6$\sigma$ statistical significance. This is the first detection of photons with $E>100$ TeV from an astrophysical source.

The progenitor stars of core-collapse supernovae (CCSNe) are asymmetrically fluctuating due to turbulent convections in the late stages of their lives. The progenitor asymmetry at the pre-supernova stage has recently caught the attention as a new ingredient to facilitate shock revival in the delayed neutrino-heating mechanism. In this paper, we investigate the importance of the progenitor asymmetries to shock revival with a semi-analytical approach. Free parameters were chosen such that the time evolution of shock radii and mass accretion rates are compatible with the results of detailed numerical simulations of CCSNe in spherical symmetry. We first estimate the amplitude of asymmetries required for the shock revival by the impulsive change of pre-shock flows in the context of neutrino heating mechanism, and then convert the amplitude to the corresponding amplitude in the pre-supernova phase by taking into account the growth of asymmetries during infall. We apply our model to various types of progenitors and find that the requisite amplitude of pre-supernova asymmetry is roughly three times larger than the prediction by current stellar evolution models unless other additional physical ingredients such as multi-dimensional fluid instabilities and turbulent convections in post-shock flows aid shock revival. We thus conclude that progenitor asymmetries can not trigger the shock revival by the impulsive way but rather play a supplementary role in reality.

We investigate axisymmetric steady solutions of (magneto)hydrodynamics equations that describe approximately accretion flows through a standing shock wave and discuss the effects of rotation and magnetic field on the revival of the stalled shock wave in supernova explosions. We develop a new powerful numerical method to calculate the 2-dimensional (2D) steady accretion flows self-consistently. We first confirm the results of preceding papers that there is a critical luminosity of irradiating neutrinos, above which there exists no steady solution in spherical models. If a collapsing star has rotation and/or magnetic field, the accretion flows are no longer spherical owing to the centrifugal force and/or Lorentz force and the critical luminosity is modified.In fact we find that the critical luminosity is reduced by about 50% - 70% for rapid rotations and about 20% - 50% for strong toroidal magnetic fields, depending on the mass accretion rate. These results may be also interpreted as an existence of the critical specific angular momentum or critical magnetic field, above which there exists no steady solution and the standing shock wave will revive for a given combination of mass accretion rate and neutrino luminosity.

We propose a new class of method for solving nonlinear systems of equations, which, among other things,has four nice features: (i) it is inspired by the mathematical property of damped oscillators, (ii) it can be regarded as a simple extention to the Newton-Raphson(NR) method, (iii) it has the same local convergence as the NR method does, (iv) it has a significantly wider convergence region or the global convergence than that of the NR method. In this article, we present the evidence of these properties, applying our new method to some examples and comparing it with the NR method.

We examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

The existence of neutron stars with $2M_\odot$ requires the strong stiffness of the equation of state (EoS) of neutron-star matter. We introduce a multi-pomeron exchange potential (MPP) working universally among 3- and 4-baryons to stiffen the EoS. Its strength is restricted by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The EoSs are derived using the Brueckner-Hartree-Fock (BHF) and the cluster variational method (CVM) with the nuclear interactions ESC and AV18. The mass-radius relations are derived by solving the Tolmann-Oppenheimer-Volkoff (TOV) equation, where the maximum masses over $2M_\odot$ are obtained on the basis of the terrestrial data. Neutron-star radii $R$ at a typical mass $1.5M_\odot$ are predicted to be $12.3\!\sim\!13.0$ km. The uncertainty of calculated radii is mainly from the ratio of 3- and 4-pomeron coupling constants, which cannot be fixed by any terrestrial experiment. Though values of $R(1.5M_\odot)$ are not influenced by hyperon-mixing effects, finely-observed values for them indicate degrees of EoS softening by hyperon mixing in the region of $M\!\sim\!2M_\odot$. If $R(1.5M_\odot)$ is less than about 12.4 km, the softening of EoS by hyperon mixing has to be weak. Useful information can be expected by the space mission NICER offering precise measurements for neutron-star radii within $\pm 5\%$.

The explosion mechanism of core-collapse supernovae has not been fully understood yet but multi-dimensional fluid instabilities such as standing accretion shock instability (SASI) and convection are now believed to be crucial for shock revival. Another multi-dimensional effect that has been recently argued is the asymmetric structures in progenitors, which are induced by violent convections in silicon/oxygen layers that occur before the onset of collapse, as revealed by recent numerical simulations of the last stage of massive star evolutions. Furthermore, it has been also demonstrated numerically that accretions of such non-spherical envelopes could facilitate shock revival. These two multi-dimensional may hence hold a key to successful explosions. In this paper, we performed a linear stability analysis of the standing accretion shock in core-collapse supernovae, taking into account non-spherical, unsteady accretion flows onto the shock to clarify the possible links between the two effects. We found that such pre-shock perturbations can excite the fluid instabilities efficiently and hence help the shock revive in core-collapse supernovae.

We investigate the effects of the odd-state part of bare $\Lambda \Lambda$ interactions on the structure of neutron stars (NSs) by constructing equations of state (EOSs) for uniform nuclear matter containing $\Lambda$ and $\Sigma^-$ hyperons with use of the cluster variational method. The isoscalar part of the Argonne v18 two-nucleon potential and the Urbana IX three-nucleon potential are employed as the interactions between nucleons, whereas, as the bare $\Lambda N$ and even-state $\Lambda \Lambda$ interactions, two-body central potentials that are determined so as to reproduce the experimental data on single- and double-$\Lambda$ hypernuclei are adopted. In addition, the $\Sigma^- N$ interaction is constructed so as to reproduce the empirical single-particle potential of $\Sigma^-$ in symmetric nuclear matter. Since the odd-state part of the $\Lambda \Lambda$ interaction is not known owing to lack of experimental data, we construct four EOSs of hyperonic nuclear matter, each with a different odd-state part of the $\Lambda \Lambda$ interaction. The EOS obtained for NS matter becomes stiffer as the odd-state $\Lambda \Lambda$ interaction becomes more repulsive, and correspondingly the maximum mass of NSs increases. It is interesting that the onset density of $\Sigma^-$ depends strongly on the repulsion of the odd-state $\Lambda \Lambda$ interaction. Furthermore, we take into account the three-baryon repulsive force to obtain results that are consistent with observational data on heavy NSs.

Neutron stars with large masses $\sim 2M_{\odot}$ require the hard stiffness of equation of state (EoS) of neutron-star matter. On the other hand, hyperon mixing brings about remarkable softening of EoS. In order to solve this problem, a multi-pomeron exchange potential (MPP) is introduced as a model for the universal many-body repulsion in baryonic systems on the basis of the Extended Soft Core (ESC) baryon-baryon interaction. The strength of MPP is determined by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The interactions in $\Lambda\!N$, $\Sigma\!N$ and $\Xi\!N$ channels are shown to be consistent with experimental indications. The EoS in neutron-star matter with hyperon mixing is obtained from ESC in addition of MPP, and mass-radius relations of neutron stars are derived. The maximum mass is shown to reach $2M_{\odot}$ even in the case of including hyperon mixing on the basis of model-parameters determined by terrestrial experiments.

A multi-pomeron exchange potential (MPP) is proposed as a model for the universal many-body repulsion in baryonic systems on the basis of the Extended Soft Core (ESC) bryon-baryon interaction. The strength of MPP is determined by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The interaction in $\Lambda N$ channels is shown to reproduce well the experimental $\Lambda$ binding energies. The equation of state (EoS) in neutron matter with hyperon mixing is obtained including the MPP contribution, and mass-radius relations of neutron stars are derived. It is shown that the maximum mass can be larger than the observed one $2M_{\odot}$ even in the case of including hyperon mixing on the basis of model-parameters determined by terrestrial experiments.

In light of positive signals reported by the CDMS-II Si experiment and the recent results of the LUX and SuperCDMS experiments, we study isospin-violating dark matter scenarios assuming that the interaction of the dark matter is mediated by colored particles. We investigate the phenomenology of the model, including collider searches, flavor and CP phenomenology. A minimal possible scenario includes scalar dark matter and new vector-like colored fermions with masses of O(1) TeV as mediators. Such a scenario may be probed at the 14 TeV LHC, while flavor and CP constraints are stringent and severe tuning in the couplings is unavoidable. We also found that, as an explanation of the CDMS-II Si signal, isospin-violating fermionic dark matter models with colored scalar mediators are disfavored by the LHC constraints.

A multi-pomeron exchange potential (MPP) is proposed as a model for the three-body repulsion indicated in neutron-star matter, which works universally among three- and four-baryons. Its strength is determined by analyzing the nucleus-nucleus scattering with the G-matrix folding model. The EoS in neutron matter is obtained including the MPP contribution. The neutron-star mass is calculated by solving the TOV equation. The maximum mass is obtained to be larger than the observed one $1.97 M_{solar}$ on the basis of the experimental data.

We report on a clear solar-cycle variation of the Sun's shadow in the 10 TeV cosmic-ray flux observed by the Tibet air shower array during a full solar cycle from 1996 to 2009. In order to clarify the physical implications of the observed solar cycle variation, we develop numerical simulations of the Sun's shadow, using the Potential Field Source Surface (PFSS) model and the Current Sheet Source Surface (CSSS) model for the coronal magnetic field. We find that the intensity deficit in the simulated Sun's shadow is very sensitive to the coronal magnetic field structure, and the observed variation of the Sun's shadow is better reproduced by the CSSS model. This is the first successful attempt to evaluate the coronal magnetic field models by using the Sun's shadow observed in the TeV cosmic-ray flux.

We evaluate the exposure during nadir observations with JEM-EUSO, the Extreme Universe Space Observatory, on-board the Japanese Experiment Module of the International Space Station. Designed as a mission to explore the extreme energy Universe from space, JEM-EUSO will monitor the Earth's nighttime atmosphere to record the ultraviolet light from tracks generated by extensive air showers initiated by ultra-high energy cosmic rays. In the present work, we discuss the particularities of space-based observation and we compute the annual exposure in nadir observation. The results are based on studies of the expected trigger aperture and observational duty cycle, as well as, on the investigations of the effects of clouds and different types of background light. We show that the annual exposure is about one order of magnitude higher than those of the presently operating ground-based observatories.

We perform some experimental simulations in spherical symmetry and axisymmetry to understand the post-shock-revival evolution of core-collapse supernovae. Assuming that the stalled shock wave is relaunched by neutrino heating and employing the so-called light bulb approximation, we induce shock revival by raising the neutrino luminosity by hand up to the critical value, which is also de- termined by dynamical simulations. A 15M_{sun} progenitor model is employed. We incorporate nuclear network calculations with a consistent equation of state in the simulations to account for the energy release by nuclear reactions and their feedback to hydrodynamics. Varying the shock-relaunch time rather arbitrarily, we investigate the ensuing long-term evolutions systematically, paying particular attention to the explosion energy and nucleosynthetic yields as a function of this relaunch time, or equivalently the accretion rate at shock revival. We study in detail how the diagnostic explosion energy approaches the asymptotic value and which physical processes contribute to the explosion energy in what proportions as well as their dependence on the relaunch time and the dimension of dynamics. We find that the contribution of nuclear reactions to the explosion energy is comparable to or greater than that of neutrino heating. In particular, recombinations are dominant over burnings in the contributions of nuclear reactions. Interestingly 1D models studied in this paper cannot produce the appropriate explosion energy and nickel mass simultaneously, overproducing nickels, whereas this problem is resolved in 2D models if the shock is relaunched at 300-400ms after bounce.

A new hybrid experiment has been started by AS{\gamma} experiment at Tibet, China, since August 2011, which consists of a low threshold burst-detector-grid (YAC-II, Yangbajing Air shower Core array), the Tibet air-shower array (Tibet-III) and a large underground water Cherenkov muon detector (MD). In this paper, the capability of the measurement of the chemical components (proton, helium and iron) with use of the (Tibet-III+YAC-II) is investigated by means of an extensive Monte Carlo simulation in which the secondary particles are propagated through the (Tibet-III+YAC-II) array and an artificial neural network (ANN) method is applied for the primary mass separation. Our simulation shows that the new installation is powerful to study the chemical compositions, in particular, to obtain the primary energy spectrum of the major component at the knee.

We develop a new semi-dynamical method to study shock revival by neutrino heating in core- collapse supernovae. Our new approach is an extension of the previous studies that employ spherically symmetric, steady, shocked accretion flows together with the light bulb approximation. The latter has been widely used in the supernova community for the phenomenological investigation of the criteria for successful supernova explosions. In the present approach, on the other hand, we get rid of the steady-state condition and take into account shock wave motions instead. We have in mind the scenario that not the critical luminosity but the critical fluctuation generated by hydrodynamical instabilities such as SASI and neutrino-driven convection in the post-shock region determines the onset of shock revival. After confirming that the new approach indeed captures the dynamics of revived shock wave qualitatively, we then apply the method to various initial conditions and find that there is a critical fluctuation for shock revival, which can be well fit by the following formula: f_crit ~ 0.8 * (M_in/1.4M_sun) * {1- (rsh/10^8cm)}, in which fcrit denotes the critical pressure fluctuation normalized by the unperturbed post-shock value. Min and rsh stand for the mass of the central compact object and the shock radius, respectively. The critical fluctuation decreases with the shock radius, whereas it increases with the mass of the central object. We discuss the possible implications of our results for 3D effects on shock revival, which is currently controversial in the supernova community.

Contributions of the JEM-EUSO Collaboration to the 32nd International Cosmic Ray Conference, Beijing, August, 2011.

The AKARI All-Sky Catalogues are an important infrared astronomical database for next-generation astronomy that take over the IRAS catalog. We have developed an online service, AKARI Catalogue Archive Server (AKARI-CAS), for astronomers. The service includes useful and attractive search tools and visual tools. One of the new features of AKARI-CAS is cached SIMBAD/NED entries, which can match AKARI catalogs with other catalogs stored in SIMBAD or NED. To allow advanced queries to the databases, direct input of SQL is also supported. In those queries, fast dynamic cross-identification between registered catalogs is a remarkable feature. In addition, multiwavelength quick-look images are displayed in the visualization tools, which will increase the value of the service. In the construction of our service, we considered a wide variety of astronomers' requirements. As a result of our discussion, we concluded that supporting users' SQL submissions is the best solution for the requirements. Therefore, we implemented an RDBMS layer so that it covered important facilities including the whole processing of tables. We found that PostgreSQL is the best open-source RDBMS products for such purpose, and we wrote codes for both simple and advanced searches into the SQL stored functions. To implement such stored functions for fast radial search and cross-identification with minimum cost, we applied a simple technique that is not based on dividing celestial sphere such as HTM or HEALPix. In contrast, the Web application layer became compact, and was written in simple procedural PHP codes. In total, our system realizes cost-effective maintenance and enhancements.

We analyze the large-scale two-dimensional sidereal anisotropy of multi-TeV cosmic rays by Tibet Air Shower Array, with the data taken from 1999 November to 2008 December. To explore temporal variations of the anisotropy, the data set is divided into nine intervals, each in a time span of about one year. The sidereal anisotropy of magnitude about 0.1% appears fairly stable from year to year over the entire observation period of nine years. This indicates that the anisotropy of TeV Galactic cosmic rays remains insensitive to solar activities since the observation period covers more than a half of the 23rd solar cycle.

Using the Tibet-III air shower array, we search for TeV gamma-rays from 27 potential Galactic sources in the early list of bright sources obtained by the Fermi Large Area Telescope at energies above 100 MeV. Among them, we observe 7 sources instead of the expected 0.61 sources at a significance of 2 sigma or more excess. The chance probability from Poisson statistics would be estimated to be 3.8 x 10^-6. If the excess distribution observed by the Tibet-III array has a density gradient toward the Galactic plane, the expected number of sources may be enhanced in chance association. Then, the chance probability rises slightly, to 1.2 x 10^-5, based on a simple Monte Carlo simulation. These low chance probabilities clearly show that the Fermi bright Galactic sources have statistically significant correlations with TeV gamma-ray excesses. We also find that all 7 sources are associated with pulsars, and 6 of them are coincident with sources detected by the Milagro experiment at a significance of 3 sigma or more at the representative energy of 35 TeV. The significance maps observed by the Tibet-III air shower array around the Fermi sources, which are coincident with the Milagro >=3sigma sources, are consistent with the Milagro observations. This is the first result of the northern sky survey of the Fermi bright Galactic sources in the TeV region.

We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the "skewed" feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.

We have designed and demonstrated a Superconducting Quantum Interference Device (SQUID) array linearized with cryogenic feedback. To achieve the necessary loop gain a 300 element series array SQUID is constructed from three monolithic 100-element series arrays. A feedback resistor completes the loop from the SQUID output to the input coil. The short feedback path of this Linearized SQUID Array (LISA) allows for a substantially larger flux-locked loop bandwidth as compared to a SQUID flux-locked loop that includes a room temperature amplifier. The bandwidth, linearity, noise performance, and dynamic range of the LISA are sufficient for its use in our target application: the multiplexed readout of transition-edge sensor bolometers.

The amplitude of the Compton-Getting (CG) anisotropy contains the power-law index of the cosmic-ray energy spectrum. Based on this relation and using the Tibet air-shower array data, we measure the cosmic-ray spectral index to be $-3.03 \pm 0.55_{stat} \pm < 0.62_{syst}$ between 6 TeV and 40 TeV, consistent with $-$2.7 from direct energy spectrum measurements. Potentially, this CG anisotropy analysis can be utilized to confirm the astrophysical origin of the ``knee'' against models for non-standard hadronic interactions in the atmosphere.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We report here on the detailed procedure of our Monte Carlo simulation.

The Tibet air shower array, which has an effective area of 37,000 square meters and is located at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few TeV. We are planning to add a large muon detector array to it for the purpose of increasing its sensitivity to cosmic gamma rays in the 100 TeV (10 - 1000 TeV) energy region by discriminating them from cosmic-ray hadrons. We report on the possibility of detection of gamma rays in the 100 TeV energy region in our field of view, based on the improved sensitivity of our air shower array deduced from the full Monte Carlo simulation.

We performed a search for cosmic-ray antideuterons using data collected during four BESS balloon flights from 1997 to 2000. No candidate was found. We derived, for the first time, an upper limit of 1.9E-4 (m^2 s sr GeV/nucleon)^(-1) for the differential flux of cosmic-ray antideuterons, at the 95% confidence level, between 0.17 and 1.15 GeV/nucleon at the top of the atmosphere.

Primary and atmospheric cosmic-ray spectra were precisely measured with the BESS-TeV spectrometer. The spectrometer was upgraded from BESS-98 to achieve seven times higher resolution in momentum measurement. We report absolute fluxes of primary protons and helium nuclei in the energy ranges, 1-540 GeV and 1-250 GeV/n, respectively, and absolute flux of atmospheric muons in the momentum range 0.6-400 GeV/c.

We report new measurements of the atmospheric muons at mountain altitude. The measurement was carried out with the BESS detector at the top of Mt. Norikura, Japan. The altitude is 2,770 m above sea level. Comparing our results and predictions given by some interaction models, a further appropriate model has been investigated. These studies would improve accuracy of atmospheric neutrino calculations.

Cosmic-ray proton and antiproton spectra were measured at mountain altitude, 2770 m above sea level. We observed more than 2 x 10^5 protons and 10^2 antiprotons in a kinetic energy range between 0.25 and 3.3 GeV. The zenith-angle dependence of proton flux was obtained. The observed spectra were compared with theoretical predictions.

Measurement of cosmic-ray proton, antiproton and muon spectra was carried out at mountain altitude. We observed 2 x 10^5 protons and 10^2 antiprotons in a kinetic energy region of 0.25 -- 3.3 GeV. Zenith-angle dependence of proton fluxes was obtained. Atmospheric muon spectra were measured simultaneously. The observed antiproton spectrum showed some deviation from theoretical predictions particularly in a low energy region.

The cosmic-ray proton, helium, and muon spectra at small atmospheric depths of 4.5 -- 28 g/cm^2 were precisely measured during the slow descending period of the BESS-2001 balloon flight. The variation of atmospheric secondary particle fluxes as a function of atmospheric depth provides fundamental information to study hadronic interactions of the primary cosmic rays with the atmosphere.

The vertical absolute fluxes of atmospheric muons and muon charge ratio have been measured precisely at different geomagnetic locations by using the BESS spectrometer. The observations had been performed at sea level (30 m above sea level) in Tsukuba, Japan, and at 360 m above sea level in Lynn Lake, Canada. The vertical cutoff rigidities in Tsukuba (36.2 N, 140.1 E) and in Lynn Lake (56.5 N, 101.0 W) are 11.4 GV and 0.4 GV, respectively. We have obtained vertical fluxes of positive and negative muons in a momentum range from 0.6 to 20 GeV/c with systematic errors less than 3 % in both measurements. By comparing the data collected at two different geomagnetic latitudes, we have seen an effect of cutoff rigidity. The dependence on the atmospheric pressure and temperature, and the solar modulation effect have been also clearly observed. We also clearly observed the decrease of charge ratio of muons at low momentum side with at higher cutoff rigidity region.

The energy spectra of cosmic-ray low-energy antiprotons and protons have been measured by BESS in 1999 and 2000, during a period covering the solar magnetic field reversal. Based on these measurements, a sudden increase of the antiproton to proton flux ratio following the solar magnetic field reversal was observed, and it generally agrees with a drift model of the solar modulation.

A search for ultra-high energy gamma-rays emitted by the young, nearby supernova remnant that was discovered recently by the COMPTEL and ROSAT satellites was made using the JANZOS database for the period 1987-1993. A 95% confidence upper limit on the flux above 100 TeV of 3x10^(-13)cm^(-2)sec^(-1) was obtained. This is an order of magnitude below the expected flux based on the standard model of cosmic ray acceleration in supernova shocks. An optical survey of the region that has been commenced is also reported. This uses UK and ESO Schmidt plates, and CCD images by a NZ/Japan microlensing group.

The Telescope Array prototype detectors were installed at Akeno Observatory and at the Utah Fly's Eye site. Using these detectors, we have observed the Crab Nebula and AGN's since the end of 1995. The successful detections of TeV gamma rays from Crab Nebula and Mkn501 are reported.