The European Open Science Cloud (EOSC) aims to create a federated environment for hosting and processing research data to support science in all disciplines without geographical boundaries, such that data, software, methods and publications can be shared as part of an Open Science community of practice. This work presents the ongoing activities related to the implementation of visual analytics services, integrated into EOSC, towards addressing the diverse astrophysics user communities needs. These services rely on visualisation to manage the data life cycle process under FAIR principles, integrating data processing for imaging and multidimensional map creation and mosaicing, and applying machine learning techniques for detection of structures in large scale multidimensional maps.

Identification of charged particles in a multilayer detector by the energy loss technique may also be achieved by the use of a neural network. The performance of the network becomes worse when a large fraction of information is missing, for instance due to detector inefficiencies. Algorithms which provide a way to impute missing information have been developed over the past years. Among the various approaches, we focused on normal mixtures models in comparison with standard mean imputation and multiple imputation methods. Further, to account for the intrinsic asymmetry of the energy loss data, we considered skew-normal mixture models and provided a closed form implementation in the Expectation-Maximization (EM) algorithm framework to handle missing patterns. The method has been applied to a test case where the energy losses of pions, kaons and protons in a six-layers Silicon detector are considered as input neurons to a neural network. Results are given in terms of reconstruction efficiency and purity of the various species in different momentum bins.

We present new ATCA multi-wavelength radio measurements (range 2.1-21.2 GHz) of the early-type magnetic star rho Oph A, performed in March 2019 during 3 different observing sessions. These new ATCA observations evidence a clear rotational modulation of the stellar radio emission and the detection of coherent auroral radio emission from rho Oph A at 2.1 GHz. We collected high-resolution optical spectra of rho Oph A acquired by several instruments over a time span of about ten years. We also report new magnetic field measurements of rho Oph A that, together with the radio light curves and the temporal variation of the equivalent width of the HeI line (lambda=5015 Angstrom), were used to constrain the rotation period and the stellar magnetic field geometry. The above results have been used to model the stellar radio emission, modelling that allowed us to constrain the physical condition of rho Oph A's magnetosphere. Past XMM measurements showed periodic X-ray pulses from rho Oph A. We correlate the X-ray light curve with the magnetic field geometry of rho Oph A. The already published XMM data have been re-analyzed showing that the X-ray spectra of rho Oph A are compatible with the presence of a non-thermal X-ray component. We discuss a scenario where the emission phenomena occurring at the extremes of the electromagnetic spectrum, radio and X-ray, are directly induced by the same plasma process. We interpret the observed X-ray and radio features of rho Oph A as having an auroral origin.

We present a catalogue of a large sample of extended radio sources in the SCORPIO field, observed and resolved by the Australia Telescope Compact Array. SCORPIO, a pathfinder project for addressing the early operations of the Australia SKA Pathfinder, is a survey of ~5 square degrees between 1.4 and 3.1 GHz, centered at l=343.5{\deg}, b=0.75{\deg} and with an angular resolution of about 10 arcsec. It is aimed at understanding the scientific and technical challenges to be faced by future Galactic surveys. With a mean sensitivity around 100 $\mu$Jy/beam and the possibility to recover angular scales at least up to 4 arcmin, we extracted 99 extended sources, 35 of them detected for the first time. Among the 64 known sources 55 had at least a tentative classification in literature. Studying the radio morphology and comparing the radio emission with infrared we propose as candidates 6 new H II regions, 2 new planetary nebulae, 2 new luminous blue variable or Wolf--Rayet stars and 3 new supernova remnants. This study provides an overview of the potentiality of future radio surveys in terms of Galactic source extraction and characterization and a discussion on the difficulty to reduce and analyze interferometric data on the Galactic plane.

We present two new radio continuum images from the Australian Square Kilometre Array Pathfinder (ASKAP) survey in the direction of the Small Magellanic Cloud (SMC). These images are part of the Evolutionary Map of the Universe (EMU) Early Science Project (ESP) survey of the Small and Large Magellanic Clouds. The two new source lists produced from these images contain radio continuum sources observed at 960 MHz (4489 sources) and 1320 MHz (5954 sources) with a bandwidth of 192 MHz and beam sizes of 30.0"x30.0" and 16.3"x15.1", respectively. The median Root Mean Squared (RMS) noise values are 186$\mu$Jy beam$^{-1}$ (960 MHz) and 165$\mu$Jy beam$^{-1}$ (1320 MHz). To create point source catalogues, we use these two source lists, together with the previously published Molonglo Observatory Synthesis Telescope (MOST) and the Australia Telescope Compact Array (ATCA) point source catalogues to estimate spectral indices for the whole population of radio point sources found in the survey region. Combining our ASKAP catalogues with these radio continuum surveys, we found 7736 point-like sources in common over an area of 30 deg$^2$. In addition, we report the detection of two new, low surface brightness supernova remnant candidates in the SMC. The high sensitivity of the new ASKAP ESP survey also enabled us to detect the bright end of the SMC planetary nebula sample, with 22 out of 102 optically known planetary nebulae showing point-like radio continuum emission. Lastly, we present several morphologically interesting background radio galaxies.

A new era in radioastronomy will begin with the upcoming large-scale surveys planned at the Australian Square Kilometre Array Pathfinder (ASKAP). ASKAP started its Early Science program in October 2017 and several target fields were observed during the array commissioning phase. The SCORPIO field was the first observed in the Galactic Plane in Band 1 (792-1032 MHz) using 15 commissioned antennas. The achieved sensitivity and large field of view already allow to discover new sources and survey thousands of existing ones with improved precision with respect to previous surveys. Data analysis is currently ongoing to deliver the first source catalogue. Given the increased scale of the data, source extraction and characterization, even in this Early Science phase, have to be carried out in a mostly automated way. This process presents significant challenges due to the presence of extended objects and diffuse emission close to the Galactic Plane. In this context we have extended and optimized a novel source finding tool, named CAESAR , to allow extraction of both compact and extended sources from radio maps. A number of developments have been done driven by the analysis of the SCORPIO map and in view of the future ASKAP Galactic Plane survey. The main goals are the improvement of algorithm performances and scalability as well as of software maintainability and usability within the radio community. In this paper we present the current status of CAESAR and report a first systematic characterization of its performance for both compact and extended sources using simulated maps. Future prospects are discussed in light of the obtained results.

We report the detection of the auroral radio emission from the early-type magnetic star HD142301. New VLA observations of HD142301 detected highly polarized amplified emission occurring at fixed stellar orientations. The coherent emission mechanism responsible for the stellar auroral radio emission amplifies the radiation within a narrow beam, making the star where this phenomenon occurs similar to a radio lighthouse. The elementary emission process responsible for the auroral radiation mainly amplifies one of the two magneto-ionic modes of the electromagnetic wave. This explains why the auroral pulses are highly circularly polarized. The auroral radio emission of HD142301 is characterized by a reversal of the sense of polarization as the star rotates. The effective magnetic field curve of HD142301 is also available making it possible to correlate the transition from the left to the right-hand circular polarization sense (and vice-versa) of the auroral pulses with the known orientation of the stellar magnetic field. The results presented in this letter have implications for the estimation of the dominant magneto-ionic mode amplified within the HD142301 magnetosphere.

At radio wavelengths, solar-type stars emit thermal free-free and gyroresonance, gyrosynchrotron, and impulsive coherent emission. Thermal free-free emission originates at layers where the optical depth is close to unit, while high brightness temperature, variable emission, can be due to flares via gyrosynchrotron emission. We observed the alpha Cen system with the Australian Telescope Compact Array at 2 GHz for three days and 17 GHz for one day. Both stars have been detected at 17 GHz, while only an upper limit has been obtained at low frequency despite the longer integration time. The brightness temperatures are consistent with the temperature of the upper chromosphere of the Sun. Inverting the formulae of the free-free emission, the average electron density of the plasma has been inferred. The same procedure was applied to the data in the millimetre recently acquired with ALMA. A comparison with the atmospheric solar models reveals a higher level of activity in alpha Cen B rather than in alpha Cen A, even if still at quiescent level. The non detection at 2 GHz allows us to put a lower limit in the filling factor of active regions. The claimed detection of an Earth size planet in close orbit to alpha Cen B, although doubtful, opens the opportunity to check the existence of Star-Planet Magnetic Interaction (SPMI). We constructed dynamic spectra in the 1.3 - 2.9 GHz of the 2 stars to search for time-variable coherent emission but obtained a null result.

In this paper, we present the first extended catalogue of far-infrared fluxes of Galactic bubbles. Fluxes were estimated for 1814 bubbles, defined here as the `golden sample', and were selected from the Milky Way Project First Data Release (Simpson et al.) The golden sample was comprised of bubbles identified within the Wide-field Infrared Survey Explorer (WISE) dataset (using 12- and 22-$\mu$m images) and Herschel data (using 70-, 160-, 250-, 350- and 500-$\mu$m wavelength images). Flux estimation was achieved initially via classical aperture photometry and then by an alternative image analysis algorithm that used active contours. The accuracy of the two methods was tested by comparing the estimated fluxes for a sample of bubbles, made up of 126 H II regions and 43 planetary nebulae, which were identified by Anderson et al. The results of this paper demonstrate that a good agreement between the two was found. This is by far the largest and most homogeneous catalogue of infrared fluxes measured for Galactic bubbles and it is a step towards the fully automated analysis of astronomical datasets.

In the next few years the classification of radio sources observed by the large surveys will be a challenging problem, and spectral index is a powerful tool for addressing it. Here we present an algorithm to estimate the spectral index of sources from multiwavelength radio images. We have applied our algorithm to SCORPIO (Umana et al. 2015), a Galactic Plane survey centred around 2.1 GHz carried out with ATCA, and found we can measure reliable spectral indices only for sources stronger than 40 times the rms noise. Above a threshold of 1 mJy, the source density in SCORPIO is 20 percent greater than in a typical extra-galactic field, like ATLAS (Norris et al. 2006), because of the presence of Galactic sources. Among this excess population, 16 sources per square degree have a spectral index of about zero, suggesting optically thin thermal emission such as Hii regions and planetary nebulae, while 12 per square degree present a rising spectrum, suggesting optically thick thermal emission such as stars and UCHii regions.

We present a multi-wavelength study of the Galactic Luminous Blue Variable HR Carinae, based on new high resolution mid-infrared (IR) and radio images obtained with the Very Large Telescope (VLT) and the Australia Telescope Compact Array (ATCA), which have been complemented by far-infrared Herschel-PACS observations and ATCA archive data. The Herschel images reveal the large-scale distribution of the dusty emitting nebula, which extends mainly to the North-East direction, up to 70 arcsec from the central star, and is oriented along the direction of the space motion of the star. In the mid-infrared images, the brightness distribution is characterized by two arc-shaped structures, tracing an inner envelope surrounding the central star more closely. At radio wavelengths, the ionized gas emission lies on the opposite side of the cold dust with respect to the position of the star, as if the ionized front was confined by the surrounding medium in the North-South direction. The comparison with previous data indicates significant changes in the radio nebula morphology and in the mass-loss rate from the central star, which has increased from 6.1$\times10^{-6}$ M$_{\odot}$yr$^{-1}$ in 1994-1995 to $1.17\times10^{-5}$ M$_{\odot}$yr$^{-1}$ in 2014. We investigate possible scenarios that could have generated the complex circumstellar environment revealed by our multi-wavelength data.

We present radio observations of 18 MIPSGAL bubbles performed at 5 GHz (6 cm) with the Karl G. Jansky Very Large Array in configuration B and BnA. The observations were aimed at understanding what kind of information high-resolution and high-sensitivity radio maps can supply on the circumstellar envelopes of different kinds of evolved stars and what their comparison with infrared images with similar resolution can tell us. We found that the 18 bubbles can be grouped into five categories according to their radio morphology. The three bubbles presenting a central point source in the radio images all correspond to luminous blue variable star candidates. Eleven bubbles show an elliptical shape and the total lack of a central object in the radio, and are likely associated with planetary nebulae. Under this assumption we derive their distance, their ionized mass and their distribution on the Galactic plane. We discuss the possibility that the MIPSGAL bubbles catalogue (428 objects) may contain a large fraction of all Galactic planetary nebulae located at a distance between 1.4 kpc and 6.9 kpc and lying in the MIPSGAL field of view. Among the remaining bubbles we identify also a H II region and a proto-planetary nebula candidate.

Automated source extraction and parameterization represents a crucial challenge for the next-generation radio interferometer surveys, such as those performed with the Square Kilometre Array (SKA) and its precursors. In this paper we present a new algorithm, dubbed CAESAR (Compact And Extended Source Automated Recognition), to detect and parametrize extended sources in radio interferometric maps. It is based on a pre-filtering stage, allowing image denoising, compact source suppression and enhancement of diffuse emission, followed by an adaptive superpixel clustering stage for final source segmentation. A parameterization stage provides source flux information and a wide range of morphology estimators for post-processing analysis. We developed CAESAR in a modular software library, including also different methods for local background estimation and image filtering, along with alternative algorithms for both compact and diffuse source extraction. The method was applied to real radio continuum data collected at the Australian Telescope Compact Array (ATCA) within the SCORPIO project, a pathfinder of the ASKAP-EMU survey. The source reconstruction capabilities were studied over different test fields in the presence of compact sources, imaging artefacts and diffuse emission from the Galactic plane and compared with existing algorithms. When compared to a human-driven analysis, the designed algorithm was found capable of detecting known target sources and regions of diffuse emission, outperforming alternative approaches over the considered fields.

The Pierre Auger Observatory has begun a major Upgrade of its already impressive capabilities, with an emphasis on improved mass composition determination using the surface detectors of the Observatory. Known as AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors on top of all 1660 water-Cherenkov detectors, updated and more flexible surface detector electronics, a large array of buried muon detectors, and an extended duty cycle for operations of the fluorescence detectors. This Preliminary Design Report was produced by the Collaboration in April 2015 as an internal document and information for funding agencies. It outlines the scientific and technical case for AugerPrime. We now release it to the public via the arXiv server. We invite you to review the large number of fundamental results already achieved by the Observatory and our plans for the future.

The Cherenkov Telescope Array (CTA) is an initiative to build the next generation, ground-based gamma-ray observatories. We present a prototype workspace developed at INAF that aims at providing innovative solutions for the CTA community. The workspace leverages open source technologies providing web access to a set of tools widely used by the CTA community. Two different user interaction models, connected to an authentication and authorization infrastructure, have been implemented in this workspace. The first one is a workflow management system accessed via a science gateway (based on the Liferay platform) and the second one is an interactive virtual desktop environment. The integrated workflow system allows to run applications used in astronomy and physics researches into distributed computing infrastructures (ranging from clusters to grids and clouds). The interactive desktop environment allows to use many software packages without any installation on local desktops exploiting their native graphical user interfaces. The science gateway and the interactive desktop environment are connected to the authentication and authorization infrastructure composed by a Shibboleth identity provider and a Grouper authorization solution. The Grouper released attributes are consumed by the science gateway to authorize the access to specific web resources and the role management mechanism in Liferay provides the attribute-role mapping.

A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.

Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85{\deg} to +20{\deg}, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

The Pierre Auger Observatory is the world's largest cosmic ray observatory. Our current exposure reaches nearly 40,000 km$^2$ str and provides us with an unprecedented quality data set. The performance and stability of the detectors and their enhancements are described. Data analyses have led to a number of major breakthroughs. Among these we discuss the energy spectrum and the searches for large-scale anisotropies. We present analyses of our X$_{max}$ data and show how it can be interpreted in terms of mass composition. We also describe some new analyses that extract mass sensitive parameters from the 100% duty cycle SD data. A coherent interpretation of all these recent results opens new directions. The consequences regarding the cosmic ray composition and the properties of UHECR sources are briefly discussed.

Joint contributions of the Pierre Auger and Telescope Array Collaborations to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013: cross-calibration of the fluorescence telescopes, large scale anisotropies and mass composition.

Ultra-high energy cosmic rays (UHECRs) interacting with the atmosphere generate extensive air showers (EAS) of secondary particles. The depth corresponding to the maximum development of the shower, $\Xmax$, is a well-known observable for determining the nature of the primary cosmic ray which initiated the cascade process. In this paper, we present an empirical model to describe the distribution of $\Xmax$ for EAS initiated by nuclei, in the energy range from $10^{17}$ eV up to $10^{21}$ eV, and by photons, in the energy range from $10^{17}$ eV up to $10^{19.6}$ eV. Our model adopts the generalized Gumbel distribution motivated by the relationship between the generalized Gumbel statistics and the distribution of the sum of non-identically distributed variables in dissipative stochastic systems. We provide an analytical expression for describing the $\Xmax$ distribution for photons and for nuclei, and for their first two statistical moments, namely $\langle \Xmax\rangle$ and $\sigma^{2}(\Xmax)$. The impact of the hadronic interaction model is investigated in detail, even in the case of the most up-to-date models accounting for LHC observations. We also briefly discuss the differences with a more classical approach and an application to the experimental data based on information theory.

Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013

In the present work we carry out a study of the high energy cosmic rays mass identification capabilities of a hybrid detector employing both fluorescence telescopes and particle detectors at ground using simulated data. It involves the analysis of extensive showers with zenith angles above 60 degrees making use of the joint distribution of the depth of maximum and muon size at ground level as mass discriminating parameters. The correlation and sensitivity to the primary mass are investigated. Two different techniques - clustering algorithms and neural networks - are adopted to classify the mass identity on an event-by-event basis. Typical results for the achieved performance of identification are reported and discussed. The analysis can be extended in a very straightforward way to vertical showers or can be complemented with additional discriminating observables coming from different types of detectors.

We consider multivariate skew-t distributions for modeling composition data of high energy cosmic rays. The model has been validated with simulated data for different primary nuclei and hadronic models focusing on the depth of maximum Xmax and number of muons N{\mu} observables. Further, we consider mixtures of multivariate skew-t distributions for cosmic ray mass composition determination and event-by-event classification. With respect to other approaches in the field, it is based on analytical calculations and allows to incorporate different sets of constraints provided by the present hadronic models. We present some applications to simulated data sets generated with different nuclear abundances assumptions. As it does not fully rely on the hadronic model predictions, the method is particularly suited to the current experimental scenario in which evidences of discrepancies of the measured data with respect to the models have been reported for some shower observables, such as the number of muons at ground level.

The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.

The Pierre Auger Observatory is a facility built to detect air showers produced by cosmic rays above 10^17 eV. During clear nights with a low illuminated moon fraction, the UV fluorescence light produced by air showers is recorded by optical telescopes at the Observatory. To correct the observations for variations in atmospheric conditions, atmospheric monitoring is performed at regular intervals ranging from several minutes (for cloud identification) to several hours (for aerosol conditions) to several days (for vertical profiles of temperature, pressure, and humidity). In 2009, the monitoring program was upgraded to allow for additional targeted measurements of atmospheric conditions shortly after the detection of air showers of special interest, e.g., showers produced by very high-energy cosmic rays or showers with atypical longitudinal profiles. The former events are of particular importance for the determination of the energy scale of the Observatory, and the latter are characteristic of unusual air shower physics or exotic primary particle types. The purpose of targeted (or "rapid") monitoring is to improve the resolution of the atmospheric measurements for such events. In this paper, we report on the implementation of the rapid monitoring program and its current status. The rapid monitoring data have been analyzed and applied to the reconstruction of air showers of high interest, and indicate that the air fluorescence measurements affected by clouds and aerosols are effectively corrected using measurements from the regular atmospheric monitoring program. We find that the rapid monitoring program has potential for supporting dedicated physics analyses beyond the standard event reconstruction.

The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of their shower signals in the water-Cherenkov detectors. In this paper we present for the first time an analysis based on down-going neutrinos. We describe the search procedure, the possible sources of background, the method to compute the exposure and the associated systematic uncertainties. No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. Assuming an E^-2 differential energy spectrum the limit on the single flavour neutrino is (E^2 * dN/dE) < 1.74x10^-7 GeV cm^-2 s^-1 sr^-1 at 90% C.L. in the energy range 1x10^17 eV < E < 1x10^20 eV.

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malarg\"ue and averaged monthly models, the utility of the GDAS data is shown.