Geant4 simulations of soft proton scattering in X-ray optics

Low energy protons (< 300 keV) can enter the field of view of X-ray telescopes, scatter on their mirror surfaces at small incident angles, and deposit energy on the detector. This phenomenon can cause intense background flares at the focal plane decreasing the mission observing time (e.g. the XMM-Newton mission) or in the most extreme cases, damaging the X-ray detector. A correct modelization of the physics process responsible for the grazing angle scattering processes is mandatory to evaluate the impact of such events on the performance (e.g. observation time, sensitivity) of future X-ray telescopes as the ESA ATHENA mission. The Remizovich model describes particles reflected by solids at glancing angles in terms of the Boltzmann transport equation using the diffuse approximation and the model of continuous slowing down in energy. For the first time this solution, in the approximation of no energy losses, is implemented, verified, and qualitatively validated on top of the Geant4 release 10.2, with the possibility to add a constant energy loss to each interaction. This implementation is verified by comparing the simulated proton distribution to both the theoretical probability distribution and with independent ray-tracing simulations. Both the new scattering physics and the Coulomb scattering already built in the official Geant4 distribution are used to reproduce the latest experimental results on grazing angle proton scattering. At 250 keV multiple scattering delivers large proton angles and it is not consistent with the observation. Among the tested models, the single scattering seems to better reproduce the scattering efficiency at the three energies but energy loss obtained at small scattering angles is significantly lower than the experimental values. In general, the energy losses obtained in the experiment are higher than what obtained by the simulation. The experimental data are not completely representative of the soft proton scattering experienced by current X-ray telescopes because of the lack of measurements at low energies (< 200 keV) and small reflection angles, so we are not able to address any of the tested models as the one that can certainly reproduce the scattering behavior of low energy protons expected for the ATHENA mission. We can, however, discard multiple scattering as the model able to reproduce soft proton funnelling, and affirm that Coulomb single scattering can represent, until further measurements at lower energies are available, the best approximation of the proton scattered angular distribution at the exit of X-ray optics.     

Tracing clinothem geometry and sediment pathways in the prograding Holocene Po Delta system through integrated core stratigraphy

Though clinothem geometry represents a key control on fluid flow in reservoir modelling, tracing clinothem boundaries accurately is commonly limited by the lack of sufficiently precise outcrop or subsurface data. This study shows that in basin systems with strongly heterogeneous compositional signatures, the combination of bulk-sediment geochemistry and benthic foraminiferal distribution can help identify clinothem architecture and generate realistic models of 3D deltaic upbuilding and evolution. Middle-late Holocene deposits in the Po Delta area form an aggradational to progradational parasequence set that reveals the complex interaction of W–E Po Delta progradation, S-directed longshore currents (from Alpine rivers) and Apennines rivers supply. Unique catchment lithologies (ophiolite rocks and dolostones) were used to delineate basin-wide geochemical markers of sediment provenance (Cr and Mg) and to assess distinctive detrital signatures. The geochemical characterization of cored intervals across different components of the sediment routing system enabled a direct linkage between clinothem growth, transport pathways and provenance mixing to be established. On the other hand, abrupt microfaunal variations at clinothem boundaries were observed to reflect the palaeoenvironmental response to sharp changes in sediment flux and fluvial influence. This study documents the ability of an integrated geochemical and palaeoecological approach to delineate three distinct sources (Po, Alps and Apennines) that contributed to coastal progradation and to outline the otherwise lithologically cryptic geometries of clinothems that using conventional sedimentological methods it would be virtually impossible to restore.     

Accuracy assessment of the GPS-based slant total electron content

The main scope of this research is to assess the ultimate accuracy that can be achieved for the slant total electron content (sTEC) estimated from dual-frequency global positioning system (GPS) observations which depends, primarily, on the calibration of the inter-frequency biases (IFB). Two different calibration approaches are analyzed: the so-called satellite-by-satellite one, which involves levelling the carrier-phase to the code-delay GPS observations and then the IFB estimation; and the so-called arc-by-arc one, which avoids the use of code-delay observations but requires the estimation of arc-dependent biases. Two strategies are used for the analysis: the first one compares calibrated sTEC from two co-located GPS receivers that serve to assess the levelling errors; and the second one, assesses the model error using synthetic data free of calibration error, produced with a specially developed technique. The results show that the arc-by-arc calibration technique performs better than the satellite-by-satellite one for mid-latitudes, while the opposite happens for low-latitudes.     

Radio frequency self-interference from a data processing centre at a radio telescope site

Apart from externally generated Radio Frequency Interference (RFI), the occurrence of self-interference is a major concern at any modern radio telescope site. Antenna servo motor controllers, data acquisition processors, and fast computing capabilities operate very close to extremely sensitive and wideband radio astronomical receivers. In this paper, we present a set of measurements of the RFI level generated by a cluster of computers that will be installed at the site of the Sardinia Radio Telescope (SRT). The measured levels are compared to Recommendation ITU-R RA.769-2, which gives the threshold levels for interference detrimental to radio astronomy observations. Our analysis shows that, with proper shielding of the noisiest devices, it will be possible to preserve the present excellent RFI conditions of the SRT site.     

Launching mechanisms of astrophysical jets

The combination of accretion disks and supersonic jets is used to model many active astrophysical objects, viz., young stars, relativistic stars, and active galactic nuclei. However, existing theories on the physical processes by which these structures transfer angular momentum and energy from disks to jets through viscous or magnetic torques are still relatively approximate. Global stationary solutions do not permit understanding the formation and stability of these structures; and global numerical simulations that include both the disk and jet physics are often limited to relatively short time scales and astrophysically out-of-range values of viscosity and resistivity parameters that are instead crucial to defining the coupling of the inflow/outflow dynamics. Along these lines we discuss self-consistent time-dependent simulations of the launching of supersonic jets by magnetized accretion disks, using high resolution numerical techniques. We shall concentrate on the effects of the disk physical parameters, and discuss under which conditions steady state solutions of the type proposed in the self-similar models of Blandford and Payne can be reached and maintained in a self-consistent nonlinear stationary state.     

Environmental and climatic conditions at a potential Glacial refugial site of tree species near the Southern Alpine glaciers. New insights from multiproxy sedimentary studies at Lago della Costa (Euganean Hills, Northeastern Italy)

It has been hypothesized that refugia of thermophilous tree species were located in Northern Italy very close to the Alps, though, this hypothesis has yet to be tested thoroughly. In contrast to Central and Southern Italy with its relative wealth of data, only a few fragmentary records are currently available from Northern Italy for the last Glacial (Würm, Weichselian). Our new study site Lago della Costa lies adjacent to the catchment of the megafans of the Alpine forelands and the braided rivers of the Northeastern Po Plain that have so far inhibited the recovery of continuous Glacial and Late-Glacial records. We analyze pollen, plant macrofossils, charcoal and ostracods to reconstruct the vegetation, fire and lake history for the period 33,000–16,000 cal. BP. We compare our data with Glacial records from Southern Europe to discuss similarities and dissimilarities between these potential refugial areas. A comparison with independent paleoclimatic proxies allows to assess potential linkages between environmental and climatic variability. New macrofossil and pollen data at Lago della Costa unambiguously document the local persistence of boreal tree taxa such as Larix decidua and Betula tree species around the study site during the last Glacial. The regular occurrence of pollen of temperate trees in the organic lake sediments (fine-detritus calcareous gyttja) suggests that temperate taxa such as Corylus avellana, Quercus deciduous, Tilia, Ulmus, Fraxinus excelsior, Carpinus, Abies alba and Fagus sylvatica, most likely survived the Last Glacial Maximum (LGM) at favorable sites in the Euganean Hills. The percentage values of temperate trees are comparable with those from Southern Europe (e.g. Monticchio in Southern Italy). We conclude that the Euganean Hills were one of the northernmost refugial areas of temperate taxa in Europe. However, the relative and absolute abundances of pollen of temperate trees are highly variable. Pollen-inferred declines of temperate tree communities (e.g. Quercetum mixtum) and low ostracod-inferred water levels at Lago della Costa correspond to the cold Heinrich events H-2 (LGM; 23,000–19,000 cal. BP) and H-3 (around 28,000 cal. BP), as recorded in the marine sediments of the North Atlantic. Similar patterns of significant temperate tree population collapses during cold Heinrich events are recorded at southern Mediterranean sites (e.g. Monticchio and the Alboran Sea). These findings suggest close linkages between Northern Atlantic and South-Central European climates during the past Glacial.     

Extension of plasticity theory to debonding,grain dissolution,and chemical damage of calcarenites

The mechanical properties of calcarenites are known to be significantly affected by water saturation: both stiffness and strength decrease for wetting in the short term and for chemical dissolution in the long term. Both processes mainly affect bonds among grains: immediately after inundation depositional bonds fall in suspension, whereas diagenetic bonds dissolve more slowly. In this paper, the authors started from the micro‐structural analysis of the weathering processes to conceive a strain hardening hydro‐chemo‐mechanical coupled elastoplastic constitutive model. The concept of extended hardening rules is here enriched: weathering functions have been determined by employing a micro to macro simplified upscaling procedure. Chemical damage is incorporated into the formulation by means of a scalar damage function. Its evolution is also described by using a multiscale approach. A new term is added to the strain rate tensor in order to incorporate the dissolution induced chemical deformations developing once the soft rock is turned into a granular material. A calibration procedure for the constitutive parameters is suggested, and the model is validated by using both coupled and uncoupled chemo‐mechanical experimental test results. Copyright © 2015 John Wiley & Sons, Ltd.     

Aquifer–peatland connectivity in southern Quebec (Canada)

In areas where peatlands are abundant, they are likely to play a significant role in the hydrological and hydrogeological dynamics of a watershed. Although individual case studies are reported in the literature, there is a large range of aquifer–peatland interactions and there is a need to understand the controls of these interactions. The objectives of this study were (1) to better understand aquifer–peatland connections and how these may be predicted by geology and geomorphic location and (2) to provide a variety of reference sites for glacial geological settings. Slope and depression peatlands were studied in the Abitibi‐Témiscamingue region and in the St. Lawrence Lowlands, two contrasting regions of southern Quebec. A total of 12 transects that span a shallow aquifer–peatland interface were instrumented with piezometers. Field investigations included peatland characterization, monthly water level monitoring, and continuous hydraulic head measurements with pressure transducers. The results indicate that 7 of the 12 transects receive groundwater from the surrounding shallow aquifer. At the peatland margin, four lateral flow patterns were identified and associated with slope peatlands (parallel inflow and divergent flow) and with depression peatlands (convergent flow and parallel outflow). Vertical hydraulic gradients suggest that water flows mainly downwards, i.e. from the peatland to the underlying mineral deposits. Vertical connectivity appears to decrease as the distance from the peatland margin increases. All of these exchanges are important components in the sustainability of peatland hydrogeological functions. The regional comparison of aquifer–peatland flow dynamics performed in this study provides a new set of referenced data for the assessment of aquifer–peatland connectivity. Copyright © 2014 John Wiley & Sons, Ltd.     

Nocturnal Intermittent Coupling Between the Interior of a Pine Forest and the Air Above It

A 1-year set of measurements of CO₂ and energy turbulent fluxes above and within a 25-m pine forest in southern Brazil is analyzed. The study focuses on the coupling state between two levels and its impact on flux determination by the eddy-covariance method. The turbulent series are split in their typical temporal scales using the multiresolution decomposition, a method that allows proper identification of the time scales of the turbulent events. Initially, four case studies are presented: a continually turbulent, a continually calm, a calm then turbulent, and an intermittent night. During transitions from calm to turbulent, large scalar fluxes of opposing signs occur at both levels, suggesting the transference of air accumulated in the canopy during the stagnant period both upwards and downwards. Average fluxes are shown for the entire period as a function of turbulence intensity and a canopy Richardson number, used as an indicator of the canopy coupling state. Above the canopy, CO₂ and sensible heat fluxes decrease in magnitude both at the neutral and at the very stable limit, while below the canopy they increase monotonically with the canopy Richardson number. Latent heat fluxes decrease at both levels as the canopy air becomes more stable. The average temporal scales of the turbulent fluxes at both levels approach each other in neutral conditions, indicating that the levels are coupled in that case. Average CO₂ fluxes during turbulent periods that succeed very calm ones are appreciably larger than the overall average above the canopy and smaller than the average or negative within the canopy, indicating that the transfer of air accumulated during calm portions at later turbulent intervals affects the flux average. The implications of this process for mean flux determination are discussed.