Quantifying bioirrigation using ecological parameters: a stochastic approach†

Irrigation by benthic macrofauna has a major influence on the biogeochemistry and microbial community structure of sediments. Existing quantitative models of bioirrigation rely primarily on chemical, rather than ecological, information and the depth-dependence of bioirrigation intensity is either imposed or constrained through a data fitting procedure. In this study, stochastic simulations of 3D burrow networks are used to calculate mean densities, volumes and wall surface areas of burrows, as well as their variabilities, as a function of sediment depth. Burrow networks of the following model organisms are considered: the polychaete worms Nereis diversicolor and Schizocardium sp., the shrimp Callianassa subterranea, the echiuran worm Maxmuelleria lankesteri, the fiddler crabs Uca minax, U. pugnax and U. pugilator, and the mud crabs Sesarma reticulatum and Eurytium limosum. Consortia of these model organisms are then used to predict burrow networks in a shallow water carbonate sediment at Dry Tortugas, FL, and in two intertidal saltmarsh sites at Sapelo Island, GA. Solute-specific nonlocal bioirrigation coefficients are calculated from the depth-dependent burrow surface areas and the radial diffusive length scale around the burrows. Bioirrigation coefficients for sulfate obtained from network simulations, with the diffusive length scales constrained by sulfate reduction rate profiles, agree with independent estimates of bioirrigation coefficients based on pore water chemistry. Bioirrigation coefficients for O(2 )derived from the stochastic model, with the diffusion length scales constrained by O(2 )microprofiles measured at the sediment/water interface, are larger than irrigation coefficients based on vertical pore water chemical profiles. This reflects, in part, the rapid attenuation with depth of the O(2 )concentration within the burrows, which reduces the driving force for chemical transfer across the burrow walls. Correction for the depletion of O(2 )in the burrows results in closer agreement between stochastically-derived and chemically-derived irrigation coefficient profiles

Observations of Sagittarius A* during the pericenter passage of the G2 object with MAGIC

Context. We present the results of a multi-year monitoring campaign of the Galactic center (GC) with the MAGIC telescopes. These observations were primarily motivated by reports that a putative gas cloud (G2) would be passing in close proximity to the super-massive black hole (SMBH), associated with Sagittarius A*, located at the center of our galaxy. This event was expected to give astronomers a unique chance to study the effect of in-falling matter on the broad-band emission of a SMBH.Aims. We search for potential flaring emission of very-high-energy (VHE; >= 100 GeV) gamma rays from the direction of the SMBH at the GC due to the passage of the G2 object. Using these data we also study the morphology of this complex region.Methods. We observed the GC region with the MAGIC Imaging Atmospheric Cherenkov Telescopes during the period 2012-2015, collecting 67 h of good-quality data. In addition to a search for variability in the flux and spectral shape of the GC gamma-ray source, we use a point-source subtraction technique to remove the known gamma-ray emitters located around the GC in order to reveal the TeV morphology of the extended emission inside that region.Results. No effect of the G2 object on the VHE gamma-ray emission from the GC was detected during the 4 yr observation campaign. We confirm previous measurements of the VHE spectrum of Sagittarius A*, and do not detect any significant variability of the emission from the source. Furthermore, the known VHE gamma-ray emitter at the location of the supernova remnant G0.9+0.1 was detected, as well as the recently discovered VHE source close to the GG radio arc

Investigating the peculiar emission from the new VHE gamma-ray source H1722+119

The Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes observed the BL Lac object H1722+119 (redshift unknown) for six consecutive nights between 2013 May 17 and 22, for a total of 12.5 h. The observations were triggered by high activity in the optical band measured by the KVA (Kungliga Vetenskapsakademien) telescope. The source was for the first time detected in the very high energy (VHE, E > 100 GeV) gamma-ray band with a statistical significance of 5.9s. The integral flux above 150 GeV is estimated to be (2.0 +/- 0.5) per cent of the Crab nebula flux. We used contemporaneous high energy (HE, 100 MeV < E < 100 GeV) gamma-ray observations from Fermi-Large Area Telescope to estimate the redshift of the source. Within the framework of the current extragalactic background light models, we estimate the redshift to be z = 0.34 +/- 0.15. Additionally, we used contemporaneous X-ray to radio data collected by the instruments on board the Swift satellite, the KVA, and the Owens Valley Radio Observatory telescope to study multifrequency characteristics of the source. We found no significant temporal variability of the flux in the HE and VHE bands. The flux in the optical and radio wavebands, on the other hand, did vary with different patterns. The spectral energy distribution of H1722+119 shows surprising behaviour in the similar to 3 x 10(14)-10(18) Hz frequency range. It can be modelled using an inhomogeneous helical jet synchrotron self-Compton model

Testing emission models on the extreme blazar 2WHSP J073326.7+515354 detected at very high energies with the MAGIC telescopes

Extreme high-energy-peaked BL Lac objects (EHBLs) are an emerging class of blazars. Their typical two-hump-structured spectral energy distribution (SED) peaks at higher energies with respect to conventional blazars. Multiwavelength (MWL) observations constrain their synchrotron peak in the medium to hard X-ray band. Their gamma-ray SED peaks above the GeV band, and in some objects it extends up to several TeV. Up to now, only a few EHBLs have been detected in the TeV gamma-ray range. In this paper, we report the detection of the EHBL 2WHSP J073320,7+515354, observed and detected during 2018 in TeV gamma rays with the MAGIC telescopes. The broad-band SED is studied within an MWL context, including an analysis of the Fermi-LAT data over 10 yr of observation and with simultaneous Swift-XRT, Swift-UVOT, and KVA data. Our analysis results in a set of spectral parameters that confirms the classification of the source as an EIME. In order to investigate the physical nature of this extreme emission, different theoretical frameworks were tested to model the broadband SED. The hard TeV spectrum of 2WHSP J073326.7+515354 sets the SED far from the energy equipartition regime in the standard one-zone leptonic scenario of blazar emission. Conversely, more complex models of the jet, represented by either a two-zone spine-layer model or a hadronic emission model, better represent the broad-hand SED

Multiband variability studies and novel broadband SED modeling of Mrk 501 in 2009

Aims. We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1, which includes, among other instruments, MAGIC, VERITAS, Whipple 10 m, and Fermi-LAT to cover the gamma-ray range from 0.1 GeV to 20 TeV; RXTE and Swift to cover wavelengths from UV to hard X-rays; and GASP-WEBT, which provides coverage of radio and optical wavelengths. Optical polarization measurements were provided for a fraction of the campaign by the Steward and St. Petersburg observatories. We evaluate the variability of the source and interband correlations, the gamma-ray flaring activity occurring in May 2009, and interpret the results within two synchrotron self-Compton (SSC) scenarios.Methods. The multiband variability observed during the full campaign is addressed in terms of the fractional variability, and the possible correlations are studied by calculating the discrete correlation function for each pair of energy bands where the significance was evaluated with dedicated Monte Carlo simulations. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters, hence providing a less biased interpretation than the "single-curve SSC model adjustment" typically reported in the literature.Results. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found on the basis of the acquired data set. The SSC model grid-scan shows that the flaring activity around May 22 cannot be modeled adequately with a one-zone SSC scenario (using an electron energy distribution with two breaks), while it can be suitably described within a two (independent) zone SSC scenario. Here, one zone is responsible for the quiescent emission from the averaged 4.5-month observing period, while the other one, which is spatially separated from the first, dominates the flaring emission occurring at X-rays and very-high-energy (> 100 GeV, VHE) gamma-rays. The flaring activity from May 1, which coincides with a rotation of the electric vector polarization angle (EVPA), cannot be satisfactorily reproduced by either a one-zone or a two-independent-zone SSC model, yet this is partially affected by the lack of strictly simultaneous observations and the presence of large flux changes on sub-hour timescales (detected at VHE gamma rays).Conclusions. The higher variability in the VHE emission and lack of correlation with the X-ray emission indicate that, at least during the 4.5-month observing campaign in 2009, the highest energy (and most variable) electrons that are responsible for the VHE gamma rays do not make a dominant contribution to the similar to 1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the SSC scenario. The studies with our dedicated SSC grid-scan show that there is some degeneracy in both the one-zone and the two-zone SSC scenarios probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The observed gamma-ray flaring activity, with the EVPA rotation coincident with the first gamma-ray flare, resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties

MAGIC detection of very high energy gamma-ray emission from the low-luminosity blazar 1ES 1741+196

We present the first detection of the nearby (z = 0.084) low-luminosity BL Lac object 1ES 1741+196 in the very high energy (E > 100 GeV) band. This object lies in a triplet of interacting galaxies. Early predictions had suggested 1ES 1741+196 to be, along with several other high-frequency BL Lac sources, within the reach of MAGIC detectability. Its detection by MAGIC, later confirmed by VERITAS, helps to expand the small population of known TeV BL Lacs. The source was observed with the MAGIC telescopes between 2010 April and 2011 May, collecting 46 h of good quality data. These observations led to the detection of the source at 6.0 s confidence level, with a steady flux F(> 100 GeV) = (6.4 +/- 1.7(stat) +/- 2.6(syst)) x 10-12 ph cm(-2)s(-1) and a differential spectral photon index Gamma = 2.4 +/- 0.2(stat) +/- 0.2(syst) in the range of similar to 80 GeV-3 TeV. To study the broad-band spectral energy distribution (SED) simultaneous with MAGIC observations, we use KVA, Swift/UVOT and XRT and Fermi/LAT data. One-zone synchrotron-self-Compton (SSC) modelling of the SED of 1ES 1741+196 suggests values for the SSC parameters that are quite common among known TeV BL Lacs except for a relatively low Doppler factor and slope of electron energy distribution. A thermal feature seen in the SED is well matched by a giant elliptical's template. This appears to be the signature of thermal emission from the host galaxy, which is clearly resolved in optical observations

Investigating the peculiar emission from the new VHE gamma-ray source H1722+119

The Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes observed the BL Lac object H1722+119 (redshift unknown) for six consecutive nights between 2013 May 17 and 22, for a total of 12.5 h. The observations were triggered by high activity in the optical band measured by the KVA (Kungliga Vetenskapsakademien) telescope. The source was for the first time detected in the very high energy (VHE, E > 100 GeV) gamma-ray band with a statistical significance of 5.9s. The integral flux above 150 GeV is estimated to be (2.0 +/- 0.5) per cent of the Crab nebula flux. We used contemporaneous high energy (HE, 100 MeV < E < 100 GeV) gamma-ray observations from Fermi-Large Area Telescope to estimate the redshift of the source. Within the framework of the current extragalactic background light models, we estimate the redshift to be z = 0.34 +/- 0.15. Additionally, we used contemporaneous X-ray to radio data collected by the instruments on board the Swift satellite, the KVA, and the Owens Valley Radio Observatory telescope to study multifrequency characteristics of the source. We found no significant temporal variability of the flux in the HE and VHE bands. The flux in the optical and radio wavebands, on the other hand, did vary with different patterns. The spectral energy distribution of H1722+119 shows surprising behaviour in the similar to 3 x 10(14)-10(18) Hz frequency range. It can be modelled using an inhomogeneous helical jet synchrotron self-Compton model

VERY HIGH ENERGY gamma-RAYS FROM THE UNIVERSE'S MIDDLE AGE: DETECTION OF THE z=0.940 BLAZAR PKS 1441+25 WITH MAGIC

The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 sigma using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO. B0218 +357 (z = 0.944), PKS. 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies with the Large Area Telescope on board Fermi. Multi-wavelength observations suggest a subdivision of the high state into two distinct flux states. In the band covered by MAGIC, the variability timescale is estimated to be 6.4 +/- 1.9 days. Modeling the broadband spectral energy distribution with an external Compton model, the location of the emitting region is understood as originating in the jet outside the broad-line region (BLR) during the period of high activity, while being partially within the BLR during the period of low (typical) activity. The observed VHE spectrum during the highest activity is used to probe the extragalactic background light at an unprecedented distance scale for ground-based gamma-ray astronomy