Gap Filling and Quality Assessment of CO2 and Water Vapour Fluxes above an Urban Area with Radial Basis Function Neural Networks

Vertical turbulent fluxes of water vapour, carbon dioxide, and sensible heat were measured from 16 August to the 28 September 2006 near the city centre of Münster in north-west Germany. In comparison to results of measurements above homogeneous ecosystem sites, the CO₂ fluxes above the urban investigation area showed more peaks and higher variances during the course of a day, probably caused by traffic and other varying, anthropogenic sources. The main goal of this study is the introduction and establishment of a new gap filling procedure using radial basis function (RBF) neural networks, which is also applicable under complex environmental conditions. We applied adapted RBF neural networks within a combined modular expert system of neural networks as an innovative approach to fill data gaps in micrometeorological flux time series. We found that RBF networks are superior to multi-layer perceptron (MLP) neural networks in the reproduction of the highly variable turbulent fluxes. In addition, we enhanced the methodology in the field of quality assessment for eddy covariance data. An RBF neural network mapping system was used to identify conditions of a turbulence regime that allows reliable quantification of turbulent fluxes through finding an acceptable minimum of the friction velocity. For the data analysed in this study, the minimum acceptable friction velocity was found to be 0.15 m s⁻¹. The obtained CO₂ fluxes, measured on a tower at 65 m a.g.l., reached average values of 12 μmol m⁻² s⁻¹ and fell to nighttime minimum values of 3 μmol m ⁻² s⁻¹. Mean daily CO₂ emissions of 21 g CO₂ m⁻²d ⁻¹ were obtained during our 6-week experiment. Hence, the city centre of Münster appeared to be a significant source of CO₂. The half-hourly average values of water vapour fluxes ranged between 0.062 and 0.989 mmol m⁻² s⁻¹and showed lower variances than the simultaneously measured fluxes of CO₂.     

Tidal effects and the environment dependence of halo assembly

We explore a possible origin for the puzzling anti-correlation between the formation epoch of galactic dark-matter haloes and their environment density. This correlation has been revealed from cosmological N -body simulations and is in conflict with the extended Press–Schechter model of halo clustering. Using similar simulations, we first quantify the straightforward association of an early formation epoch with a reduced mass-growth rate at late times. We then find that a primary driver of suppressed growth, by accretion and mergers, is tidal effects dominated by a neighbouring massive halo. The tidal effects range from a slowdown of the assembly of haloes due to the shear along the large-scale filaments that feed the massive halo to actual mass loss in haloes that pass through the massive halo. Using the restricted three-body problem, we show that haloes are prone to tidal mass loss within 1.5 virial radii of a larger halo. Our results suggest that the dependence of the formation epoch on environment density is a secondary effect induced by the enhanced density of haloes in filaments near massive haloes where the tides are strong. Our measures of assembly rate are particularly correlated with the tidal field at high redshifts   z ∼ 1  .     

Origin of fault domains and fault-domain boundaries (transfer zones and accommodation zones) in extensional provinces: Result of random nucleation and self-organized fault growth

In many extensional provinces, large normal faults dip in the same direction forming fault domains. Features variously named transfer faults, transfer zones, and accommodation zones (hereafter non-genetically referred to as fault-domain boundaries) separate adjacent fault domains. Experimental modeling of distributed extension provides insights on the origin, geometry, and evolution of these fault domains and fault-domain boundaries. In our scaled models, a homogeneous layer of wet clay or dry sand overlies a latex sheet that is stretched orthogonally or obliquely between two rigid sheets. Fault domains and fault-domain boundaries develop in all models in both map view and cross-section. The number, size, and arrangement of fault domains as well as the number and orientation of fault-domain boundaries are variable, even for models with identical boundary conditions. The fault-domain boundaries in our models differ profoundly from those in many published conceptual models of transfer/accommodation zones. In our models, fault-domain boundaries are broad zones of deformation (not discrete strike-slip or oblique-slip faults), their orientations are not systematically related to the extension direction, and they can form spontaneously without any prescribed pre-existing zones of weakness. We propose that fault domains develop because early-formed faults perturb the stress field, causing new nearby faults to dip in the same direction (self-organized growth). As extension continues, faults from adjacent fault domains propagate toward each another. Because opposite-dipping faults interfere with one another in the zone of overlap, the faults stop propagating. In this case, the geometry of the domain boundaries depends on the spatial arrangement of the earliest formed faults, a result of the random distribution of the largest flaws at which the faults nucleate.     

Properties of dusty tori in active galactic nuclei – I. The case of SWIRE/SDSS quasars

We derive the properties of dusty tori in active galactic nuclei from the comparison of observed spectral energy distributions (SEDs) of SDSS quasars and a precomputed grid of torus models. The observed SEDs comprise SDSS photometry, Two-Micron All-Sky Survey J , H and K data, whenever available, and mid-infrared (mid-IR) data from the Spitzer Wide-area InfraRed Extragalactic Survey. The adopted model is that of Fritz, Franceschini & Hatziminaoglou. The fit is performed by standard  χ²  -minimization; the model, however, can be a multicomponent comprising a stellar and a starburst component, whenever necessary. Models with low equatorial optical depth, τ_9.7, were allowed as well as 'traditional' models with  τ_9.7≥ 1.0  , corresponding to   A _V≥ 22  and the results were compared. Fits using high optical depth tori models only produced dust more compactly distributed than in the configuration where all τ_9.7 models were permitted. Tori with decreasing dust density with the distance from the centre were favoured while there was no clear preference for models with or without angular variation of the dust density. The computed outer radii of the tori are of some tens of parsecs large but can reach, in a few cases, a few hundreds of parsecs. The mass of dust, M _Dust, and IR luminosity, L _IR, integrated in the wavelength range between 1 and 1000 μm, do not show significant variations with redshift, once the observational biases are taken into account. Objects with 70-μm detections, representing 25 per cent of the sample, are studied separately and the starburst contribution (whenever present) to the IR luminosity can reach, in the most extreme but very few cases, 80 per cent.     

Contrasting degrees of recrystallization of carbonaceous material in the Nelson aureole,British Columbia and Ballachulish aureole,Scotland, with implications for thermometry based on Raman spectroscopy of carbonaceous material

The degree of recrystallization of carbonaceous material (CM), as monitored by Raman microspectroscopy, was examined as a function of metamorphic grade in two well‐studied contact aureoles containing carbonaceous pelites: the Nelson aureole, British Columbia and the Ballachulish aureole, Scotland. Here, we use (a) the R2 ratio extracted from the Raman spectrum of CM as a proxy for the degree of graphitization (0.0 in perfect graphite then increasing with structural defects) and (b) the second‐order S1 band (~2,700 cm^?1) as a marker for the tridimensional ordering of CM. The Nelson aureole (garnet–staurolite–andalusite–sillimanite–K‐feldspar sequence, ~550–650°C, 3.5–4.0 kbar) was developed in rocks that were unmetamorphosed prior to contact metamorphism, whereas the Ballachulish aureole (cordierite–andalusite–K‐feldspar–sillimanite sequence, ~550–700°C, ~3.0 kbar) was developed in rocks that had been metamorphosed to garnet grade conditions (~7 kbar, ~500°C) c. 45 Ma before contact metamorphism. Thirty‐one samples were examined from Nelson and 29 samples from Ballachulish. At Nelson, the R2 ratio steadily decreases from ~0.25 to 0.0 as the igneous contact is approached, whereas at Ballachulish, the R2 ratio remains largely unchanged from regional values (~0.20–0.25) until less than 100 m from the igneous contact. The second‐order S1 band reveals that carbonaceous material (CM) was transformed to highly “ordered” locally tridimensional graphitic carbon at Ballachulish by regional metamorphism prior to contact metamorphism, whereas CM was still a disordered turbostratic (bidimensional) material before contact metamorphism in the case of Nelson. Pretexturation of CM likely induced sluggish recrystallization of CM and delayed graphitization in the Ballachulish aureole. Temperatures of recrystallization of the CM in the two aureoles were estimated using different published calibrations of the thermometry based on Raman Spectroscopy of Carbonaceous Material (RSCM), with differences among the calibrations being minor. In the Nelson aureole, temperatures are in reasonable agreement with those indicated by the metapelitic phase equilibria (all within 50°C, most within 25°C). In the Ballachulish aureole, the retarded crystallization noted above results in increasing underestimates of temperatures compared to the metapelitic phase equilibria (up to ~75°C too low within 200 m of the igneous contact). Our study calls for careful attention when using RSCM thermometry in complexly polymetamorphosed rocks to assess properly the meaning of the calculated temperature.     

Testing the fidelity of thermometers at ultrahigh temperatures

A highly residual granulite facies rock (sample RG07‐21) from Lunnyj Island in the Rauer Group, East Antarctica, presents an opportunity to compare different approaches to constraining peak temperature in high‐grade metamorphic rocks. Sample RG07‐21 is a coarse‐grained pelitic migmatite composed of abundant garnet and orthopyroxene along with quartz, biotite, cordierite, and plagioclase with accessory rutile, ilmenite, zircon, and monazite. The inferred sequence of mineral growth is consistent with a clockwise pressure–temperature (P–T) evolution when compared with a forward model (P–T pseudosection) for the whole‐rock chemical composition. Peak metamorphic conditions are estimated at 9 ± 0.5 kbar and 910 ± 50°C based on conventional Al‐in‐orthopyroxene thermobarometry, Zr‐in‐rutile thermometry, and calculated compositional isopleths. U–Pb ages from zircon rims and neocrystallized monazite grains yield ages of c. 514 Ma, suggesting that crystallization of both minerals occurred towards the end of the youngest pervasive metamorphic episode in the region known as the Prydz Tectonic Event. The rare earth element compositions of zircon and garnet are consistent with equilibrium growth of these minerals in the presence of melt. When comparing the thermometry methods used in this study, it is apparent that the Al‐in‐orthopyroxene thermobarometer provides the most reliable estimate of peak conditions. There is a strong textural correlation between the temperatures obtained using the Zr‐in‐rutile thermometer––maximum temperatures are recorded by a single rutile grain included within orthopyroxene, whereas other grains included in garnet, orthopyroxene, quartz, and biotite yield a range of temperatures down to 820°C. Ti‐in‐zircon thermometry returns significantly lower temperature estimates of 678–841°C. Estimates at the upper end of this range are consistent with growth of zircon from crystallizing melt at temperatures close to the elevated (H₂O undersaturated) solidus. Those estimates, significantly lower than the calculated temperature of this residual solidus, may reflect isolation of rutile from the effective equilibration volume leading to an activity of TiO₂ that is lower than the assumed value of unity.