Fracture sealing of quartzite under a temperature gradient: experimental results

Veined lithologies are formed by fracturing and sealing processes, with the veins representing former fluid conduits through the rock. Although detailed fieldwork and numerical simulations have provided a better understanding of vein growth, few studies have attempted to seal fractures and generate veins experimentally. In this pilot study, we subjected fractured quartzite to temperature gradients of 45–125 °C under hydrothermal conditions in a static fluid, with the aim of precipitating secondary quartz in the cooler portions of the fracture. Results show that secondary quartz precipitates due to the imposed temperature gradient, causing the initial fracture to seal locally. Although no systematic sealing pattern was observed along the fracture, samples subjected to higher temperatures exhibit a smaller fracture width and appear to have reacted more extensively. Electron microprobe mapping visualizes the spatial distribution of secondary quartz, which contains elevated concentrations of aluminium.     

Variability associated with α accretion disc theory for standard and advection-dominated discs

The α turbulent viscosity formalism for accretion discs must be interpreted as a mean field theory, modelling a steady state only on spatial or time-scales greater than those of the turbulence. The extent of the scale separation determines the relative precision error (RPE) of the predicted luminosity L _ν. Turbulence and the use of α implies that (1) field line stretching gives a magnetic pressure  α²/6 of the total pressure generally, and a one-to-one relation between α and the pressure ratio for thin discs, and (2) large turbulent scales in advection-dominated accretion flows (ADAFs) predict a lower L _ν precision than thin discs for a given observation duration and central mass. The allowed variability (or RPE) at frequency ν increases with the size of the contributing region. For X-ray binary ADAFs, the RPE ∼ 5 per cent at R  ≤ 1000 Schwarzchild radii ( R _s) for averages over  1000 s. However, current data for galaxies like NGC 4258 and M87 give RPEs in L _ν of 50–100 per cent even at R  ≤ 100  R _S. More data are required, but systematic deviations from ADAF predictions are more significant than random deviations, and may constrain properties of the turbulence, the accretion mode, the assumption of a steady state or the accretion rate.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

Understanding the implications of the EU-LULUCF regulation for the wood supply from EU forests to the EU

Background

In June 2018, the European Parliament and Council of the European Union adopted a legislative regulation for incorporating greenhouse gas emissions and removals from Land Use, Land Use Change and Forestry (EU-LULUCF) under its 2030 Climate and Energy Framework. The LULUCF regulation aim to incentivise EU Member States to decrease greenhouse gas emissions and increase removals in the LULUCF sector. The regulation, however, does not set a target for increasing the LULUCF carbon sink, but rather includes a ‘no net debit’ target for LULUCF (Forests and Agricultural soils). For Managed Forest Land (MFL) an accounting framework with capped credits for additional mitigation against a set forest reference level (FRL) was agreed for 2021–2030. The FRL gives the projected future carbon sink in the two compliance periods 2021–2025 and 2026–2030 under “continuation of forest management practices as they were in the reference period 2000–2009”. This FRL was disputed by some Member States as it was perceived to put a limit on their future wood harvesting from MFL. Here we simulated with the EFISCEN European forest model the “continuation of forest management practices” and determined the corresponding wood harvest for 26 EU countries under progressing age classes.

Results

The simulations showed that under “continuation of forest management practices” the harvest (wood removals) in the 26 EU countries as a whole can increase from 420 million m³/year in 2000–2009 to 560 million m³/year in 2050 due to progressing age classes. This implies there is a possibility to increase absolute wood harvests without creating debits compared to the forest reference level. However, the manner in which ‘continuation of forest management’ developed with a progressing age class development over time, meant that in some countries the future harvesting exceeded 90% of the increment. Since this generally is considered to be unsustainable we additionally set a harvesting cut-off as max 90% of increment to be harvested for each individual country as a possible interpretation of sustainability criteria that are included in the regulation. Using this additional limit the projected harvest will only increase to 493 million m³/year.

Conclusions

The worry from Member States (MS) that the FRL will prevent any additional harvesting seems unwarranted. Due to differences between Member States concerning the state of their forest resources, the FRL as a baseline for harvesting works out very differently for the different Member States. The FRL may have other unforeseen consequences which we discuss. Under all scenarios the living forest biomass sink shows a decline. This can be counteracted through incentivising measures under Climate Smart Forestry.

     

Stochastic wobble of accretion discs and jets from turbulent rocket torques

Models of accretion discs and their associated outflows often incorporate assumptions of axisymmetry and symmetry across the disc plane. However, for turbulent discs these symmetries only apply to averaged quantities and do not apply locally. The local asymmetries can induce local imbalances in outflow power across the disc mid-plane, which can in turn induce local tilting torques. Here we calculate the effect of the resulting stochastic torques on disc annuli that are a consequence of standard mean field accretion disc models. The torques induce a random walk of the vector perpendicular to the plane of each averaged annulus. This random walk is characterized by a radially dependent diffusion coefficient which we calculate for small angle tilt. We use the coefficient to calculate a radially dependent time-scale for annular tilt and associated jet wobble. The wobble time depends on the square of the wander angle so the age of a given system determines the maximum wobble angle. We apply this to examples of blazars, young stellar objects and binary engines of pre-planetary nebulae and microquasars. It is noteworthy that for an averaging time   t _w∼ 3 d  , we estimate a wobble angle for jets in SS 433 of  θ∼ 0.8°  , not inconsistent with observational data. In general the non-periodic nature of the stochastic wobble could distinguish it from faster periodic jet precession.     

A GIS-based method for defining snow zones: application to the western United States

This study maps the geographic extent of intermittent and seasonal snow cover in the western United States using thresholds of 2000–2010 average snow persistence derived from moderate resolution imaging spectroradiometer snow cover area data from 1 January to 3 July. Results show seasonal snow covers 13% of the region, and intermittent snow covers 25%. The lower elevation boundaries of intermittent and seasonal snow zones increase from north-west to south-east. Intermittent snow is primarily found where average winter land surface temperatures are above freezing, whereas seasonal snow is primarily where winter temperatures are below freezing. However, temperatures at the boundary between intermittent and seasonal snow exhibit high regional variability, with average winter seasonal snow zone temperatures above freezing in west coast mountain ranges. Snow cover extent at peak accumulation is most variable at the upper elevations of the intermittent snow zone, highlighting the sensitivity of this snow zone boundary to climate conditions.     

Mapping collective human activity in an urban environment based on mobile phone data

Identifying and characterizing variations of human activity – specifically changes in intensity and similarity – in urban environments provide insights into the social component of those eminently complex systems. Using large volumes of user-generated mobile phone data, we derive mobile communication profiles that we use as a proxy for the collective human activity. In this article, geocomputational methods and geovisual analytics such as self-organizing maps (SOM) are used to explore the variations of these profiles, and its implications for collective human activity. We evaluate the merits of SOM as a cross-dimensional clustering technique and derived temporal trajectories of variations within the mobile communication profiles. The trajectories’ characteristics such as length are discussed, suggesting spatial variations in intensity and similarity in collective human activity. Trajectories are linked back to the geographic space to map the spatial and temporal variation of trajectory characteristics. Different trajectory lengths suggest that mobile phone activity is correlated with the spatial configuration of the city, and so at different times of the day. Our approach contributes to the understanding of the space-time social dynamics within urban environments.     

Conjoint modeling of residential group preferences: A comparison of the internal validity of hierarchical information integration approaches

 In this paper, two approaches for measuring residential group preferences, based on the method of Hierarchical Information Integration (HII), are compared. In particular, the hypothesis that group-based preference models estimated from integrated HII experiments better predict group preferences than part individual-based group models estimated from classical HII experiments is tested. To that effect, the models' ability to predict group preferences for new residential alternatives is compared in a study of residential preferences of co-ops. Results indicate that integrated HII group experiments indeed result in better predictions of residential preferences.