Reionization and the Cosmic Dawn with the Square Kilometre Array

The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21 cm line from the earliest phases of star and galaxy formation in the Universe. This 21 cm signal provides a new and unique window both on the time of the formation of the first stars and accreting black holes and the subsequent period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.     

A first deep seismic survey in the Sea of Marmara: Deep basins and whole crust architecture and evolution

Increased source strength, streamer length and dense spatial coverage of seismic reflection profiles of the SEISMARMARA Leg 1 allow to image the deep structure of the marine North Marmara Trough (NMT) on the strike-slip North Anatolian Fault (NAF) west of the destructive Izmit 1999 earthquake. A reflective lower crust and the Moho boundary are detected. They appear upwarped on an E-W profile from the southern Central Basin eastwards, towards more internal parts of the deformed region. Thinning of the upper crust could use a detachment suggested from an imaged dipping intracrustal reflector that would allow upper crustal material to be dragged from beneath it and above the lower crust, accounting for the extensional component but also southwest motion of the southern margin of the NMT. Sections across the eastern half of the NMT, crossing the Cinarcik and Imrali basins, reveal several faults that are active reaching into the basement and have varying strike and proportions of normal and strike-slip displacement. They might be viewed as petals of a large scale negative flower-structure that spreads over a width of 30 km at surface and is rooted deeper in the lithosphere. Under the Central Basin a very thick sediment infill is revealed and its extensional bounding faults are active and imaged as much as 8 km apart down to 6 km depth. We interpret them as two deep-rooted faults encompassing a foundering basement block, rather than being merely pulled-apart from a jog in a strike-slip above a décollement. The deep-basin lengthening would account for only a modest part of the proposed 60 km finite motion since 4 Myr along the same direction oblique to the NMT that sidesteps the shear motion from its two ends. Thus differential motion occurred much beyond the deep basins, like subsidence involving the NMT bounding faults and the intracrustal detachments. The complex partitioned motion localized on active faults with diverse natures and orientations is suggested to represent the overburden deformation induced from horizontal plane simple shear occurring in depth at lithospheric scale, and in front of the North Anatolian Fault when it propagated through the region.     

Multi-stage evolution of a sub-aerial volcanic ridge over the last 1.3 Myr: S. Jorge Island, Azores Triple Junction

New K/Ar dating and geochemical analyses have been carried out on the WNW–ESE elongated oceanic island of S. Jorge to reconstruct the volcanic evolution of a linear ridge developed close to the Azores triple junction. We show that S. Jorge sub-aerial construction encompasses the last 1.3 Myr, a time interval far much longer than previously reported. The early development of the ridge involved a sub-aerial building phase exposed in the southeast end of the island and now constrained between 1.32 ± 0.02 and 1.21 ± 0.02 Ma. Basic lavas from this older stage are alkaline and enriched in incompatible elements, reflecting partial melting of an enriched mantle source. At least three differentiation cycles from alkaline basalts to mugearites are documented within this stage. The successive episodes of magma rising, storage and evolution suggest an intermittent re-opening of the magma feeding system, possibly due to recurrent tensional or trans-tensional tectonic events. Present data show a gap in sub-aerial volcanism before a second main ongoing building phase starting at about 750 ka. Sub-aerial construction of the S. Jorge ridge migrated progressively towards the west, but involved several overlapping volcanic episodes constrained along the main WNW–ESE structural axis of the island. Mafic magmas erupted during the second phase have been also generated by partial melting of an enriched mantle source. Trace element data suggest, however, variable and lower degrees of partial melting of a shallower mantle domain, which is interpreted as an increasing control of lithospheric deformation on the genesis and extraction of primitive melts during the last 750 kyr. The multi-stage development of the S. Jorge volcanic ridge over the last 1.3 Myr has most likely been greatly influenced by regional tectonics, controlled by deformation along the diffuse boundary between the Nubian and the Eurasian plates, and the increasing effect of sea-floor spreading at the Mid-Atlantic Ridge.     

Pb,Hf and Nd isotope compositions of the two Réunion volcanoes (Indian Ocean): A tale of two small-scale mantle “blobs”?

Pb, Hf, Nd and Sr isotopes of basaltic lavas from the two Réunion Island volcanoes are reported in order to examine the origin of the sources feeding these volcanoes and to detect possible changes through time. Samples, chosen to cover the whole lifetime of the two volcanoes (from 2 Ma to present), yield a chemically restricted (compared to OIB lavas) but complex distribution. Réunion plume isotopic characteristics have been defined on the basis of the composition of uncontaminated shield-building lavas from the Piton de la Fournaise volcano. The average ?_Nd, ?_Hf, ⁸⁷Sr/⁸⁶Sr and ²⁰⁶Pb/²⁰⁴Pb, ²⁰⁷Pb/²⁰⁴Pb, and ²⁰⁸Pb/²⁰⁴Pb isotope ratios calculated for this component are + 4.4, + 9.1, 0.70411, 18.97, 15.59 and 39.03, respectively. In Pb–Pb isotope space, each volcano defines a distinct linear trend but slight variations are also detected within the various volcanic sequences. The Piton des Neiges volcano yields a distinct and significantly more scattered isotopic distribution than Piton de la Fournaise for both Pb, Hf and Nd isotope tracers. A principal component analysis of the Pb isotope data from Piton de la Fournaise reveals a major contribution of the C and EM-1 components (with a clear Dupal flavor) as main components for the modern Réunion plume. The same components have been identified for Piton des Neiges but with a stronger participation of a depleted mantle component and a weaker EM-1 contribution. The compositional change of the lavas erupted by the Piton des Neiges and Piton de la Fournaise volcanoes is attributed to the impingement of two small-scale blobs of plume material at the base of the Réunion lithosphere. Compared to other hot-spots worldwide, in particular Hawaii and Kerguelen, magmas beneath Réunion are generated from a considerably more homogeneous, compositionally more primitive plume higher in ²⁰⁶Pb. Although shallow-level contamination processes have been locally detected they did not alter significantly the composition of the plume magmas. This is tentatively attributed to mantle dynamics producing small, high-velocity blobs that ascend rapidly through the lithosphere, and to the lack of a well-developed magma chamber at depth in the lithosphere.     

Inversion of time-dependent nuclear well-logging data using neural networks

The purpose of this work was to investigate a new and fast inversion methodology for the prediction of subsurface formation properties such as porosity, salinity and oil saturation, using time‐dependent nuclear well logging data. Although the ultimate aim is to apply the technique to real‐field data, an initial investigation as described in this paper, was first required; this has been carried out using simulation results from the time‐dependent radiation transport problem within a borehole. Simulated neutron and γ‐ray fluxes at two sodium iodide (NaI) detectors, one near and one far from a pulsed neutron source emitting at ～14 MeV, were used for the investigation. A total of 67 energy groups from the BUGLE96 cross section library together with 567 property combinations were employed for the original flux response generation, achieved by solving numerically the time‐dependent Boltzmann radiation transport equation in its even parity form. Material property combinations (scenarios) and their correspondent teaching outputs (flux response at detectors) are used to train the Artificial Neural Networks (ANNs) and test data is used to assess the accuracy of the ANNs. The trained networks are then used to produce a surrogate model of the expensive, in terms of computational time and resources, forward model with which a simple inversion method is applied to calculate material properties from the time evolution of flux responses at the two detectors. The inversion technique uses a fast surrogate model comprising 8026 artificial neural networks, which consist of an input layer with three input units (neurons) for porosity, salinity and oil saturation; and two hidden layers and one output neuron representing the scalar photon or neutron flux prediction at the detector. This is the first time this technique has been applied to invert pulsed neutron logging tool information and the results produced are very promising. The next step in the procedure is to apply the methodology to real data.     

Using a groundwater quality negotiation support system to change land-use management near a drinking-water abstraction in the Netherlands

A negotiation support system (NSS) was developed to solve groundwater conflicts that arose during land-use management. It was set up in cooperation with the stakeholders involved to provide information on the impact of land use, e.g., agriculture, nature (forested areas), recreation, and urban areas, on the quality of both infiltrating and abstracted groundwater. This NSS combined simulation programs that calculate (1) the concentrations of nitrate in shallow groundwater for each land-use area and (2) the transport of nitrate in the groundwater-saturated zone. The user interface of the NSS enabled scenario analyses. The NSS was validated at a drinking-water abstraction near Holten (the Netherlands) using a spatial planning process aimed at sustainable land-use and groundwater-resource management. Two land-use scenarios were considered: a base scenario reflecting the autonomous development and an adapted land-use scenario. The calculated results for shallow groundwater provided an explicit spatial overview of the impact of historical land use and N application on the quality of abstracted groundwater as well as insight into the impact of changes in land use and N application. Visualization of the conflicting interests of agriculture and the drinking-water abstraction helped all stakeholders accept the necessary changes in land use identified by the adapted land-use scenario of the NSS. These changes were included in the preferred land-use management option in the regional planning process, which has since been formalized. The NSS provided system insight, scoping analyses, and education, in addition to generating quantitative information on the impact of land-use functions on groundwater quality.     

Anthropogenic debris on beaches in the SE Pacific (Chile): Results from a national survey supported by volunteers

Anthropogenic marine debris (AMD) is an ubiquitous problem, which has motivated public participation in activities such as beach surveys and clean-up campaigns. While it is known that beaches in the SE Pacific are also affected by this problem, the quantities and types of AMD remain largely unknown. In the context of an outreach project, volunteers (1500 high-school students) participated in a nation-wide survey of AMD on 43 beaches distributed randomly along the entire Chilean coast (18°S to 53°S). The mean density of AMD was 1.8 items m⁻² and the major types were plastics, cigarette butts and glass. Densities in central Chile were lower than in northern and southern Chile, which could be due to different attitudes of beach users or to intense beach cleaning in central regions. We suggest that public participation in surveys and cleaning activities will raise awareness and thereby contribute to an improvement of the situation.     

Influence of Degree of Saturation in the Electric Resistivity–Hydraulic Conductivity Relationship

The relationship between aquifer hydraulic conductivity and aquifer resistivity, either measured on the ground surface by vertical electrical sounding (VES) or from resistivity logs, or measured in core samples have been published for different types of aquifers in different locations. Generally, these relationships are empirical and semi-empirical, and confined in few locations. This relation has a positive correlation in some studies and negative in others. So far, there is no potentially physical law controlling this relation, which is not completely understood. Electric current follows the path of least resistance, as does water. Within and around pores, the model of conduction of electricity is ionic and thus the resistivity of the medium is controlled more by porosity and water conductivity than by the resistivity of the rock matrix. Thus, at the pore level, the electrical path is similar to the hydraulic path and the resistivity should reflect hydraulic conductivity. We tried in this paper to study the effect of degree of groundwater saturation in the relation between hydraulic conductivity and bulk resistivity via a simple numerical analysis of Archie’s second law and a simplified Kozeny-Carmen equation. The study reached three characteristic non-linear relations between hydraulic conductivity and resistivity depending on the degree of saturation. These relations are: (1) An inverse power relation in fully saturated aquifers and when porosity equals water saturation, (2) An inverse polynomial relation in unsaturated aquifers, when water saturation is higher than 50%, higher than porosity, and (3) A direct polynomial relation in poorly saturated aquifers, when water saturation is lower than 50%, lower than porosity. These results are supported by some field scale relationships.     

Contrasting regional discharge evolutions in the Amazon basin (1974–2004)

Former hydrological studies in the Amazon Basin generally describe annual discharge variability on the main stem. However, the downstream Amazon River only represents the mean state of the Amazonian hydrological system. This study therefore uses a new data set including daily discharge in 18 sub-basins to analyze the variability of regional extremes in the Amazon basin, after recalling the diversity of the hydrological annual cycles within the Amazon basin. Several statistical tests are applied in order to detect trends and breaks in the time series. We show that during the 1974–2004 period, the stability of the mean discharge on the main stem in Óbidos is explained by opposite regional features that principally involve Andean rivers: a decrease in the low stage runoff, particularly important in the southern regions, and an increase in the high stage runoff in the northwestern region. Both features are observed from the beginning of the nineties. These features are also observed in smaller meridian sub-basins in Peru and Bolivia. Moreover we show that the changes in discharge extremes are related to the regional pluriannual rainfall variability and the associated atmospheric circulation as well as to tropical large-scale climatic indicators.