Refine
Year of publication
Language
- English (76)
Is part of the Bibliography
- yes (76)
Keywords
- gamma rays: general (24)
- radiation mechanisms: non-thermal (20)
- ISM: supernova remnants (14)
- gamma rays: galaxies (13)
- galaxies: active (12)
- cosmic rays (11)
- acceleration of particles (8)
- ISM: clouds (6)
- astroparticle physics (6)
- gamma rays: ISM (6)
Transverse dispersion, or tracer spreading orthogonal to the mean flow direction, which is relevant e.g, for quantifying bio-degradation of contaminant plumes or mixing of reactive solutes, has been studied in the literature less than the longitudinal one. Inferring transverse dispersion coefficients from field experiments is a difficult and error-prone task, requiring a spatial resolution of solute plumes which is not easily achievable in applications. In absence of field data, it is a questionable common practice to set transverse dispersivities as a fraction of the longitudinal one, with the ratio 1/10 being the most prevalent. We collected estimates of field-scale transverse dispersivities from existing publications and explored possible scale relationships as guidance criteria for applications. Our investigation showed that a large number of estimates available in the literature are of low reliability and should be discarded from further analysis. The remaining reliable estimates are formation-specific, span three orders of magnitude and do not show any clear scale-dependence on the plume traveled distance. The ratios with the longitudinal dispersivity are also site specific and vary widely. The reliability of transverse dispersivities depends significantly on the type of field experiment and method of data analysis. In applications where transverse dispersion plays a significant role, inference of transverse dispersivities should be part of site characterization with the transverse dispersivity estimated as an independent parameter rather than related heuristically to longitudinal dispersivity.
High-performance numerical codes are an indispensable tool for hydrogeologists when modeling subsurface flow and transport systems. But as they are written in compiled languages, like C/C++ or Fortran, established software packages are rarely user-friendly, limiting a wider adoption of such tools. OpenGeoSys (OGS), an open-source, finite-element solver for thermo-hydro-mechanical-chemical processes in porous and fractured media, is no exception. Graphical user interfaces may increase usability, but do so at a dramatic reduction of flexibility and are difficult or impossible to integrate into a larger workflow. Python offers an optimal trade-off between these goals by providing a highly flexible, yet comparatively user-friendly environment for software applications. Hence, we introduceogs5py, a Python-API for the OpenGeoSys 5 scientific modeling package. It provides a fully Python-based representation of an OGS project, a large array of convenience functions for users to interact with OGS and connects OGS to the scientific and computational environment of Python.
Most large-scale hydrologic models fall short in reproducing groundwater head dynamics and simulating transport process due to their oversimplified representation of groundwater flow. In this study, we aim to extend the applicability of the mesoscale Hydrologic Model (mHM v5.7) to subsurface hydrology by coupling it with the porous media simulator OpenGeoSys (OGS). The two models are one-way coupled through model interfaces GIS2FEM and RIV2FEM, by which the grid-based fluxes of groundwater recharge and the river-groundwater exchange generated by mHM are converted to fixed-flux boundary conditions of the groundwater model OGS. Specifically, the grid-based vertical reservoirs in mHM are completely preserved for the estimation of land-surface fluxes, while OGS acts as a plug-in to the original mHM modeling framework for groundwater flow and transport modeling. The applicability of the coupled model (mHM-OGS v1.0) is evaluated by a case study in the central European mesoscale river basin - Nagelstedt. Different time steps, i.e., daily in mHM and monthly in OGS, are used to account for fast surface flow and slow groundwater flow. Model calibration is conducted following a two-step procedure using discharge for mHM and long-term mean of groundwater head measurements for OGS. Based on the model summary statistics, namely the Nash-Sutcliffe model efficiency (NSE), the mean absolute error (MAE), and the interquartile range error (QRE), the coupled model is able to satisfactorily represent the dynamics of discharge and groundwater heads at several locations across the study basin. Our exemplary calculations show that the one-way coupled model can take advantage of the spatially explicit modeling capabilities of surface and groundwater hydrologic models and provide an adequate representation of the spatiotemporal behaviors of groundwater storage and heads, thus making it a valuable tool for addressing water resources and management problems.
Most large-scale hydrologic models fall short in reproducing groundwater head dynamics and simulating transport process due to their oversimplified representation of groundwater flow. In this study, we aim to extend the applicability of the mesoscale Hydrologic Model (mHM v5.7) to subsurface hydrology by coupling it with the porous media simulator OpenGeoSys (OGS). The two models are one-way coupled through model interfaces GIS2FEM and RIV2FEM, by which the grid-based fluxes of groundwater recharge and the river-groundwater exchange generated by mHM are converted to fixed-flux boundary conditions of the groundwater model OGS. Specifically, the grid-based vertical reservoirs in mHM are completely preserved for the estimation of land-surface fluxes, while OGS acts as a plug-in to the original mHM modeling framework for groundwater flow and transport modeling. The applicability of the coupled model (mHM-OGS v1.0) is evaluated by a case study in the central European mesoscale river basin - Nagelstedt. Different time steps, i.e., daily in mHM and monthly in OGS, are used to account for fast surface flow and slow groundwater flow. Model calibration is conducted following a two-step procedure using discharge for mHM and long-term mean of groundwater head measurements for OGS. Based on the model summary statistics, namely the Nash-Sutcliffe model efficiency (NSE), the mean absolute error (MAE), and the interquartile range error (QRE), the coupled model is able to satisfactorily represent the dynamics of discharge and groundwater heads at several locations across the study basin. Our exemplary calculations show that the one-way coupled model can take advantage of the spatially explicit modeling capabilities of surface and groundwater hydrologic models and provide an adequate representation of the spatiotemporal behaviors of groundwater storage and heads, thus making it a valuable tool for addressing water resources and management problems.
Geostatistics as a subfield of statistics accounts for the spatial correlations encountered in many applications of, for example, earth sciences. Valuable information can be extracted from these correlations, also helping to address the often encountered burden of data scarcity. Despite the value of additional data, the use of geostatistics still falls short of its potential. This problem is often connected to the lack of user-friendly software hampering the use and application of geostatistics. We therefore present GSTools, a Python-based software suite for solving a wide range of geostatistical problems. We chose Python due to its unique balance between usability, flexibility, and efficiency and due to its adoption in the scientific community. GSTools provides methods for generating random fields; it can perform kriging, variogram estimation and much more. We demonstrate its abilities by virtue of a series of example applications detailing their use.
A hotspot at a position compatible with the BL. Lac object 1ES 2322-409 was serendipitously detected with H.E.S.S. during observations performed in 2004 and 2006 on the blazar PKS 2316-423. Additional data on 1ES 2322-409 were taken in 2011 and 2012, leading to a total live-time of 22.3 h. Point-like very-high-energy (VHE; E > 100 GeV) gamma-ray emission is detected from a source centred on the IFS 2322-409 position, with an excess of 116.7 events at a significance of 6.0 sigma. The average VHE gamma-ray spectrum is well described with a power law with a photon index Gamma = 3.40 +/- 0.66(stat) +/- 0.20(sys) and an integral flux Phi(E > 200 GeV) = (3.11 +/- 0.71(stat) 0.62(sys)) x 10(-2)cm(-2)s(-1), which corresponds to 1.1 per cent of the Crab nebula flux above 200 GeV. Multiwavelength data obtained with Fermi LAT, Swift XRT and UVOT, RXTE PCA, ATOM, and additional data from WISE, GROND, and Catalina are also used to characterize the broad-band non-thermal emission of lES 2322-409. The multiwavelength behaviour indicates day-scale variability. Swift UVOT and XRT data show strong variability at longer scales. A spectral energy distribution (SED) is built from contemporaneous observations obtained around a high state identified in Swift data. A modelling of the SED is performed with a stationary homogeneous one-zone synchrotronself-Compton leptonic model. The redshift of the source being unknown, two plausible values were tested for the modelling. A systematic scan of the model parameters space is performed, resulting in a well-constrained combination of values providing a good description of the broad-band behaviour of 1ES 2322-409.
The blazar Mrk 501 (z = 0.034) was observed at very-high-energy (VHE, E greater than or similar to 100 GeV) gamma-ray wavelengths during a bright flare on the night of 2014 June 23-24 (MJD 56832) with the H.E.S.S. phase-II array of Cherenkov telescopes. Data taken that night by H.E.S.S. at large zenith angle reveal an exceptional number of gamma-ray photons at multi-TeV energies, with rapid flux variability and an energy coverage extending significantly up to 20 TeV. This data set is used to constrain Lorentz invariance violation (LIV) using two independent channels: a temporal approach considers the possibility of an energy dependence in the arrival time of gamma-rays, whereas a spectral approach considers the possibility of modifications to the interaction of VHE gamma-rays with extragalactic background light (EBL) photons. The non-detection of energy-dependent time delays and the non-observation of deviations between the measured spectrum and that of a supposed power-law intrinsic spectrum with standard EBL attenuation are used independently to derive strong constraints on the energy scale of LIV (E-QG) in the subluminal scenario for linear and quadratic perturbations in the dispersion relation of photons. For the case of linear perturbations, the 95% confidence level limits obtained are E-QG,E-1 > 3.6 x 10(17) GeV using the temporal approach and E-QG,E-1 > 2.6 x 10(19) GeV using the spectral approach. For the case of quadratic perturbations, the limits obtained are E-QG,E-2 > 8.5 x 10(10) GeV using the temporal approach and E-QG,E-2 > 7.8 x 10(11) GeV using the spectral approach.
First ground-based measurement of sub-20 GeV to 100 GeV gamma-Rays from the Vela pulsar with HESS II
(2018)
Aims. We report on the measurement and investigation of pulsed high-energy y-ray emission from the Vela pulsar, PSR B0833-45, based on observations with the largest telescope of H.E.S.S., CT5, in monoscopic mode, and on data obtained with the Fermi-LAT. Methods. Data from 40.3 h of observations carried out with the H.E.S.S. II array from 2013 to 2015 have been used. A dedicated very low-threshold event reconstruction and analysis pipeline was developed to achieve the lowest possible energy threshold. Eight years of Fermi-LAT data were analysed and also used as reference to validate the CT5 telescope response model and analysis methods. Results. A pulsed gamma-ray signal at a significance level of more than 15 sigma is detected from the P2 peak of the Vela pulsar light curve. Of a total of 15 835 events, more than 6000 lie at an energy below 20 GeV, implying a significant overlap between H.E.S.S. II-CT5 and the Fermi-LAT. While the investigation of the pulsar light curve with the LAT confirms characteristics previously known up to 20 GeV in the tens of GeV energy range, CT5 data show a change in the pulse morphology of P2, i.e. an extreme sharpening of its trailing edge, together with the possible onset of a new component at 3.4 sigma significance level. Assuming a power-law model for the P2 spectrum, an excellent agreement is found for the photon indices (Gamma similar or equal to 4.1) obtained with the two telescopes above 10 GeV and an upper bound of 8% is derived on the relative offset between their energy scales. Using data from both instruments, it is shown however that the spectrum of P2 in the 10-100 GeV has a pronounced curvature; this is a confirmation of the sub-exponential cut-off form found at lower energies with the LAT. This is further supported by weak evidence of an emission above 100 GeV obtained with CT5. In contrast, converging indications are found from both CT5 and LAT data for the emergence of a hard component above 50 GeV in the leading wing (LW2) of P2, which possibly extends beyond 100 GeV. Conclusions. The detection demonstrates the performance and understanding of CT5 from 100 GeV down to the sub-20 GeV domain, i.e. unprecedented low energy for ground-based gamma-ray astronomy. The extreme sharpening of the trailing edge of the P2 peak found in the H.E.S.S. II light curve of the Vela pulsar and the possible extension beyond 100 GeV of at least one of its features, LW2, provide further constraints to models of gamma-Ray emission from pulsars.
Context. NGC 253 is one of only two starburst galaxies found to emit gamma-rays from hundreds of MeV to multi-TeV energies. Accurate measurements of the very-high-energy (VHE; E> 100 GeV) and high-energy (HE; E > 60 MeV) spectra are crucial to study the underlying particle accelerators, probe the dominant emission mechanism(s) and to study cosmic-ray interaction and transport. Aims. The measurement of the VHE gamma-ray emission of NGC 253 published in 2012 by H.E.S.S. was limited by large systematic uncertainties. Here, the most up to date measurement of the gamma-ray spectrum of NGC 253 is investigated in both HE and VHE gamma-rays. Assuming a hadronic origin of the gamma-ray emission, the measurement uncertainties are propagated into the interpretation of the accelerated particle population. Methods. The data of H.E.S.S. observations are reanalysed using an updated calibration and analysis chain. The improved Fermi-LAT analysis employs more than 8 yr of data processed using pass 8. The cosmic-ray particle population is evaluated from the combined HE-VHE gamma-ray spectrum using NAIMA in the optically thin case. Results. The VHE gamma-ray energy spectrum is best fit by a power-law distribution with a flux normalisation of (1.34 +/- 0.14(stat) +/- 0.27(sys)) x 10(-13) cm(-2) s(-1) TeV-1 at 1 TeV - about 40% above, but compatible with the value obtained in Abramowski et al. (2012). The spectral index Gamma = 2.39 +/- 0.14(stat) +/- 0.25(sys) is slightly softer than but consistent with the previous measurement within systematic errors. In the Fermi energy range an integral flux of F(E > 60 MeV) = (1.56 +/- 0.28(stat) +/- 0.15(sys)) x 10(-8) cm(-2) s(-1) is obtained. At energies above similar to 3 GeV the HE spectrum is consistent with a power-law ranging into the VHE part of the spectrum measured by H.E.S.S. with an overall spectral index Gamma = 2.22 +/- 0.06(stat). Conclusions. Two scenarios for the starburst nucleus are tested, in which the gas in the starburst nucleus acts as either a thin or a thick target for hadronic cosmic rays accelerated by the individual sources in the nucleus. In these two models, the level to which NGC 253 acts as a calorimeter is estimated to a range of f(cal) = 0.1 to 1 while accounting for the measurement uncertainties. The presented spectrum is likely to remain the most accurate measurements until the Cherenkov Telescope Array (CTA) has collected a substantial set of data towards NGC 253.
Dwarf spheroidal galaxies are among the most promising targets for detecting signals of Dark Matter (DM) annihilations. The H.E.S.S. experiment has observed five of these systems for a total of about 130 hours. The data are re-analyzed here, and, in the absence of any detected signals, are interpreted in terms of limits on the DM annihilation cross section. Two scenarios are considered: i) DM annihilation into mono-energetic gamma-rays and ii) DM in the form of pure WIMP multiplets that, annihilating into all electroweak bosons, produce a distinctive gamma-ray spectral shape with a high-energy peak at the DM mass and a lower-energy continuum. For case i), upper limits at 95% confidence level of about <sigma upsilon > less than or similar to 3 x 10(-25) cm(3) s(-1) are obtained in the mass range of 400 GeV to 1TeV. For case ii), the full spectral shape of the models is used and several excluded regions are identified, but the thermal masses of the candidates are not robustly ruled out.