Searches for very high energy gamma rays from blazars with CANGAROO-III telescope in 2005-2009

We have searched for very high energy (VHE) gamma rays from four blazars using the CANGAROO-III imaging atmospheric Cherenkov telescope. We report the results of the observations of H 2356-309, PKS 2155-304, PKS 0537-441, and 3C 279, performed from 2005 to 2009, applying a new analysis to suppress the effects of the position dependence of Cherenkov images in the field of view. No significant VHE gamma ray emission was detected from any of the four blazars. The GeV gamma-ray spectra of these objects were obtained by analyzing Fermi/LAT archival data. Wide range (radio to VHE gamma-ray bands) spectral energy distributions (SEDs) including CANGAROO-III upper limits, GeV gamma-ray spectra, and archival data, even though they are non-simultaneous, are discussed using a one-zone synchrotron self-Compton (SSC) model in combination with a external Compton (EC) radiation. The HBLs (H 2356-309 and PKS 2155-304) can be explained by a simple SSC model, and PKS 0537-441 and 3C 279 are well modeled by a combination of SSC and EC model. We find a consistency with the blazar sequence in terms of strength of magnetic field and component size.     

A Refined Model for Solid Particle Rock Erosion

A procedure for the estimation of distribution parameters of a Weibull distribution model K₁ = f(K_Ic^12/4/σ_C^23/4) for solid particle erosion, as recently suggested in Rock Mech Rock Eng, doi: 10.1007/s00603-014-0658-x, 2014, is derived. The procedure is based on examinations of elastic–plastically responding rocks (rhyolite, granite) and plastically responding rocks (limestone, schist). The types of response are quantified through SEM inspections of eroded surfaces. Quantitative numbers for the distribution parameter K₁ are calculated for 30 rock materials, which cover a wide range of mechanical properties. The ranking according to the parameter K₁ is related to qualitative rock classification schemes. A modified proposal for the erosion of schist due to solid particle impingement at normal incidence is introduced.     

The oxygen isotopic composition of the fraction less than 2 microns (clay fraction) in Recent sediments from western Europe

The oxygen isotopic composition of the clay fraction in Recent argillaceous sediments in western Europe shows regional differences, the δ¹⁸_SMOW values varying between + 15.8 and + 21.7‰. This phenomenon is used for sediment-transport studies. In the Rhine-Meuse estuary a landward transport of marine sediments is observed. The clay fraction of the sediments in the Wadden Sea is only for a minor part derived from the rivers Rhine and Meuse.     

Faunal relationships between the New Zealand plateau and the New Zealand sector of Antarctica based on echinoderm distribution

The problem of New Zealand‐Antarctic faunal relationships is discussed on the basis of echinoderm distribution. The limitations of the data upon which past zoogeographical speculations have been based are pointed out.

Macquarie Island (occupying an intermediate geographical position between the New Zealand Plateau and the Antarctic) shows definite relationships with New Zealand, and the submarine Macquarie Ridge may have provided a connecting migration route. However, only four species of echinoderms are shared between the Ross Sea‐Balleny Islands area (the New Zealand sector of the Antarctic Region) and the New Zealand Pliateau‐Macquarie Island area (the New Zealand Region). Although the as yet unsampled part of the Macquarie‐Balleny Ridge may reveal other faunal similarities, the present systematic sampling has made possible a sounder understanding of the zoogeographical affinities of the two regions.     

Bizarre geochemistry of komatiites from the Crixás greenstone belt,Brazil

Komatiite lava flows in the Crixás greenstone belt, Goiás, Brazil, have textures and volcanic structures typical of Archean komatiites, but are geochemically most unusual. The flows are porphyritic and massive, or layered with spinifex upper parts and olivine cumulate lower parts. MgO contents range from 18 to 40%. In such lavas, only olivine (and minor chromite) can have crystallized, but neither major nor trace elements fall on olivine control lines. In MgO variation diagrams, CaO and Sr fall on lines with slopes steeper than olivine control lines; SiO₂, FeO, Na₂O, K₂O and Y show little systematic variation; Zr shows a large variation that does not correlate with MgO; and Al₂O₃ decreases markedly with decreasing MgO. The aberrant behaviour is highlighted by the REE (rare earth elements) in spinifex and olivine cumulate layers from three flows: in the spinifex layers, chondrite-normalized REE patterns are hump-shaped with maxima at Nd or Sm ((La/Sm)_N=0.6, (Gd/Yb)_N=1.6–2.3), whereas cumulate zones in the same flows have steadily sloping patterns, with LREE enriched relative to HREE ((La/Sm)_N=1.3, (Gd/Yb)_N=1.4). Neither normal magmatic processes acting during emplacement of the komatiites, nor thermal erosion and wall-rock assimilation can explain these effects, and we speculate that elements commonly thought of as “immobile” (e.g. Al, Zr, REE) migrated during hydrothermal alteration or metamorphism. A Pb-Pb whole rock isochron gave an age of 2,728±140 Ma and selected Sm-Nd analyses an apparent isochron age of 2,825±98 Ma (ɛ_Nd≈0). The Pb-Pb age is believed to be the approximate time of emplacement. Interpretation of the Sm-Nd data is complicated by the evidence of mobility of REE.     

Applying essentially non‐oscillatory interpolation to controlled‐source electromagnetic modelling*

Modelling and inversion of controlled‐source electromagnetic (CSEM) fields requires accurate interpolation of modelled results near strong resistivity contrasts. There, simple linear interpolation may produce large errors, whereas higher‐order interpolation may lead to oscillatory behaviour in the interpolated result. We propose to use the essentially non‐oscillatory, piecewise polynomial interpolation scheme designed for piecewise smooth functions that contains discontinuities in the function itself or in its first or higher derivatives. The scheme uses a non‐linear adaptive algorithm to select a set of interpolation points that represent the smoothest part of the function among the sets of neighbouring points. We present numerical examples to demonstrate the usefulness of the scheme. The first example shows that the essentially non‐oscillatory interpolation (ENO) scheme better captures an isolated discontinuity. In the second example, we consider the case of sampling the electric field computed by a finite‐volume CSEM code at a receiver location. In this example, the ENO interpolation performs quite well. However, the overall error is dominated by the discretization error. The other examples consider the comparison between sampling with essentially non‐oscillatory interpolation and existing interpolation schemes. In these examples, essentially non‐oscillatory interpolation provides more accurate results than standard interpolation, especially near discontinuities.     

The biogeochemical sulfur cycle in the marine boundary layer over the Northeast Pacific Ocean

The major components of the marine boundary layer biogeochemical sulfur cycle were measured simultaneously onshore and off the coast of Washington State, U.S.A. during May 1987. Seawater dimethylsulfide (DMS) concentrations on the continental shelf were strongly influenced by coastal upwelling. Concentration further offshore were typical of summer values (2.2 nmol/L) at this latitude. Although seawater DMS concentrations were high on the biologically productive continental shelf (2–12 nmol/L), this region had no measurable effect on atmospheric DMS concentrations. Atmospheric DMS concentrations (0.1–12 nmol/m³), however, were extremely dependent upon wind speed and boundary layer height. Although there appeared to be an appreciable input of non-sea-salt sulfate to the marine boundary layer from the free troposphere, the local flux of DMS from the ocean to the atmosphere was sufficient to balance the remainder of the sulfur budget.     

The gravitational field of topographic-isostatic masses and the hypothesis of mass condensation II-the topographic-isostatic geoid

In Part I we focussed on a convergent representation of the gravitational potential generated bytopographic masses on top of the equipotential surface atMean Sea Level, thegeoid, and by those masses which compensate topography. Topographic masses have also been condensated, namely represented by a single layer. Part II extends the computation of the gravitational field of topographic-isostatic masses by a detailed analysis of itsforce field in terms ofvector-spherical harmonic functions. In addition, the discontinuous mass-condensated topographic gravitational force vector (head force) is given. Once we identify theMoho discontinuity asone interface of isostatically compensated topographical masses, we have computed the topographic potential and the gravitational potential which is generated by isostatically compensated masses atMean Sea Level, the geoid, and illustrated by various figures of geoidal undulations. In comparison to a data oriented global geoid computation ofJ. Engels (1991) the conclusion can be made that the assumption of aconstant crustal mass density, the basic condition for isostatic modeling, does not apply. Insteaddensity variations in the crust, e.g. betweenoceanic and continental crust densities, have to be introduced in order to match the global real geoid and its topographic-isostatic model. The performed analysis documents that thestandard isostatic models based upon aconstant crustal density areunreal.     

Infrasonic detection of a Leonid bolide: 1998 November 17

Abstract— During the early morning hours of the night of the peak of the annual Leonid meteor shower on 1998 November 17, a bright fireball (approximately ?12 to ?14 visual magnitude at 100 km in the zenith) was observed over northern New Mexico with visual sightings as far away from Los Alamos as Albuquerque (～150 km to the south of Los Alamos), including direct persistent trail observations at the U. S. A. F. Starfire Optical Range (SOR), which is also near Albuqerque. This event did not produce any sonic boom reports, presumably because of its high altitude. It was also detected locally by an infrared radiometer at Sandia National Laboratory and by an intensified charge-coupled device (CCD) camera located in Placitas, New Mexico. Subsequent investigations of the data from the six infrasound arrays used by Los Alamos National Laboratory (LANL) and operated for the Department of Energy as a part of the Comprehensive Test Ban Treaty (CTBT) Research and Development program for the International Monitoring System (IMS) showed the presence of an infrasonic signal from the proper direction at the correct time for this bolide from two of our six arrays (both located in Los Alamos). The infrasound recordings (i.e., the wave amplitude and period data) indicated that an explosion occurred in the atmosphere at a source height of ～93.5 km (with respect to sea level) or ～90 km with respect to the altitude of Los Alamos, having its origins slightly to the north and west of Los Alamos. Purely geometric solutions from the ground observers reports combined with direct measurements from the CCD camera at Placitas produced a source height of 91 ± 7 km. The signal characteristics analyzed from 0.5 to 3.0 Hz include a total duration of about 3–4 s for a source directed from Los Alamos toward 353.6 ± 0.4° measured from true north at a maximum elevation arrival angle of ～72.7°. The latter was deduced on the basis of the observed signal trace velocities (for the part of the recording with the highest cross-correlation) and ranged from a constant value of about 920–1150 m/s (depending on the window length used in the analysis) for a ray trajectory along a direct refractive path between the source and the Los Alamos arrays. The dominant signal frequency at maximum amplitude at Los Alamos was ～0.71 Hz. These highly correlated signals had a peak to peak, maximum amplitude of ～2.1 microbars (0.21 Pa). Using several methods that incorporate various observed signal characteristics, total distance traveled, etc., our analysis indicates that the bolide probably had a source energy of ～1.14 t (TNT equivalent) or 4.77 × 10⁹ J. This is ～14.1× smaller than the source energy estimate made using the infrasonic, empirical source energy relationship for low-altitude stationary point sources developed in the 1960s by the Air Force Technical Applications Center (AFTAC), Patrick Air Force Base, Florida. This relation was originally developed, however, for much larger source energies and at much longer ranges.