Local Sensitivity Analysis of a Numerical Model of Volcanic Plinian Columns through Automatic Differentiation

Numerical simulation of geological phenomena are nowadays widely used to understand them. A good example are the numerous codes developed to model mass motion in different processes. Such models are approximations to the real phenomena because the physical conditions are simplified or because the values of the parameters intervening in the modeling are known only approximately. For this reason, an obligatory step in the application of these models is the assessment of their sensitivity to changes in the parameters. Sensitivity analysis is also required to evaluate the influence of the different input variables in the outcome of the model. In general this analysis is carried out by repeated computation of the model to different sets of input variables and settings of the model. A more fruitful procedure to carry out such analysis can be achieved by automatic differentiation. In this paper we illustrate such approach through its application to the model of an eruptive plinian column. The results of the analysis allow to assess, in a systematic and efficient manner, the relevance of the different boundary values and parameters of the model. The analysis reveals the effects on velocity, temperature, gas fraction, radius and ultimately column height, by changes in their boundary value (i.e. at the vent). The effect of changes in the value of the gas constants, density of pyroclasts and specific heat of the erupted material is also evidenced. Alternatively, the results can be regarded as a mean to estimate the values of the variables for different set of boundary values without recurring to new runs of the model. Automatic differentiation is shown to be a convenient scheme to obtain derivatives of model variables. This type of analysis forms in addition the basis for further work in inversion schemes.     

Assessing green gentrification in historically disenfranchised neighborhoods: a longitudinal and spatial analysis of Barcelona

To date, little is known about the extent to which the creation of municipal green spaces over an entire city addresses social or racial inequalities in the distribution of environmental amenities – or whether such an agenda creates contributes to green gentrification. In this study, we evaluate the effects of creating 18 green spaces in socially vulnerable neighborhoods of Barcelona during the 1990s and early 2000s. We examined the evolution over time of six socio-demographic gentrification indicators in the areas close to green spaces in comparison with the entire districts. Our results indicate that new parks in the old town and formerly industrialized neighborhoods seem to have experienced green gentrification. In contrast, most economically depressed areas and working-class neighborhoods with less desirable housing stock and more isolated from the city center gained vulnerable residents as they became greener, indicating a possible redistribution and greater concentration of vulnerable residents through the city.     

Inclusive approaches to urban climate adaptation planning and implementation in the Global South

As cities increasingly engage in climate adaptation planning, many are seeking to promote public participation and facilitate the engagement of different civil society actors. Still, the variations that exist among participatory approaches and the merits and tradeoffs associated with each are not well understood. This article examines the experiences of Quito (Ecuador) and Surat (India) to assess how civil society actors contribute to adaptation planning and implementation. The results showcase two distinct approaches to public engagement. The first emphasizes participation of experts, affected communities, and a wide array of citizens to sustain broadly inclusive programmes that incorporate local needs and concerns into adaptation processes and outcomes. The second approach focuses on building targeted partnerships between key government, private, and civil society actors to institutionalize robust decision-making structures, enhance abilities to raise funds, and increase means to directly engage with local community and international actors. A critical analysis of these approaches suggests more inclusive planning processes correspond to higher climate equity and justice outcomes in the short term, but the results also indicate that an emphasis on building dedicated multi-sector governance institutions may enhance long-term programme stability, while ensuring that diverse civil society actors have an ongoing voice in climate adaptation planning and implementation.

Policy relevance

Many local governments in the Global South experience severe capacity and resource constraints. Cities are often required to devolve large-scale planning and decision-making responsibilities, such as those critical to climate adaptation, to different civil society actors. As a result, there needs to be more rigorous assessments of how civil society participation contributes to the adaptation policy and planning process and what local social, political, and economic factors dictate the way cities select different approaches to public engagement. Also, since social equity and justice are key indicators for determining the effectiveness and sustainability of adaptation interventions, urban adaptation plans and policies must also be designed according to local institutional strengths and civic capacities in order to account for the needs of the poor and most vulnerable. Inclusivity, therefore, is critical for ensuring equitable planning processes and just adaptation outcomes.     

Model reduction and discretization using hybrid finite volumes for flow in porous media containing faults

In this paper, we study two different model reduction strategies for solving problems involving single phase flow in a porous medium containing faults or fractures whose location and properties are known. These faults are represented as interfaces of dimension N ? 1 immersed in an N dimensional domain. Both approaches can handle various configurations of position and permeability of the faults, and one can handle different fracture permeabilities on the two inner sides of the fracture. For the numerical discretization, we use the hybrid finite volume scheme as it is known to be well suited to simulating subsurface flow. Some results, which may be of use in the implementation of the proposed methods in industrial codes, are demonstrated.     

Spatio-temporal evolution of the Choisille River (southern Parisian Basin,France) during the Weichselian and the Holocene as a record of climate trend and human activity in north-western Europe

The morpho-sedimentary evolution of the Choisille floodplain (lowland river, catchment: 288 km²), a tributary of the River Loire in the south-western Parisian Basin, was studied through 61 core drillings along eight transects and a geophysical survey located in four stretches of the river: stretches A and B correspond to two sub-catchments, and stretches C and D are in the main valley. Sixty ¹⁴C and four OSL datings were obtained, and sediments were analysed on seven reference cores. Eight phases of evolution differing markedly from the evolution of more northern areas in the Parisian Basin and north-western Europe were identified from spatio-temporal distribution of nine lithological facies. The deepest incision phase (1) occurred during the first part of the Weichselian, followed by the deposition of a gravelly-sandy unit (phase 2) during the Middle Pleniglacial, which was deeply incised (phase 3), probably during the Bölling. From the Allerød up to the last third of the Boreal (phase 4), sedimentation was continuously dominated by peaty deposits, with no evidence of either increased hydraulic energy during the Younger Dryas, or of incision during the LateGlacial–Holocene transition. This trend seems to reflect the specificity of the south-western Parisian Basin climate from the Late Weichselian up to the end of the Boreal, due to the influence of the Atlantic Ocean, compared to more northern areas where the climate was more continental. The downstream incision trend during the last third of the Boreal up to the Subatlantic (phases 5 and 6) indicates a sharp increase in precipitation and vegetation cover; the lack of peaty sediments, widespread in north-western Europe, and also of precipitated carbonates frequent in the Parisian Basin, seems to be due to local physiographic characteristics. The main part of the sediment filling, which is principally silty and retrograde, began during the Subatlantic (phase 7 and 8) as a result of deforestation of the plateaux for crop farming. High human-induced sediment yield and storage concealed the possible impacts of climate change on fluvial dynamics: lithological facies change from phase 7 to 8 can only indicate the autogenic morphological evolution of the floodplain in accretion. The non-univocal upstream–downstream variation in the start of phase 7 shows that sediment yield varied in space and time in the catchment, particularly in relation to the agricultural potential of the different areas; this observation could be used to testify human-induced sedimentation in other catchments.     

The role of rift‐inherited hyper‐extension in Alpine‐type orogens

Alpine‐type orogens are interpreted to result from the collision of former rifted margins. As many present‐day rifted margins consist of hyper‐extended domains floored by thinned continental crust (<10 km) and/or exhumed mantle, this study explores the influence of rift inheritance on the architecture and final evolution of Alpine‐type orogens. We propose that rift‐related necking zones, separating weakly thinned 25‐ to 30‐km‐thick crust from hyper‐extended domains, may act as buttresses during the transition from subduction to collision. As a result, former necking zones may now be found at the boundary between a highly deformed and overthickened nappe stack, made of relics of hyper‐extended domains, and an external, weakly deformed fold‐and‐thrust belt, which largely escaped significant rift‐related crustal thinning and orogeny‐related thickening. Therefore, the role of rift inheritance is of critical importance and is largely underestimated in controlling the architecture and evolution of Alpine‐type orogens.     

Ten‐year water table recovery after clearcutting and draining boreal forested wetlands of eastern Canada

In boreal forested wetlands, the observed increase in the water table level after clearcutting (watering‐up) is often a threat to sustained ecosystem productivity. Hydrologic recovery refers to the processes by which a water table progressively drops back to its initial level after the cut. In eastern Canada, drainage is used operationally after clearcutting wet sites in order to lower the water table level and accelerate hydrologic recovery. The objective of this study was to evaluate the duration of the watering‐up caused by timber harvesting and the extent to which drainage affected the water table recovery on five peatlands and three hydromorphic mineral sites located in the St. Lawrence Lowlands of Québec (Canada). The mixed wood stands studied are dominated by balsam fir (Abies balsamea (L.) Mill.), eastern white cedar (Thuja occidentalis L.), and red maple (Acer rubrum L). Results indicate that, 10 years after clearcutting, water table levels in undrained plots are still 5 to 7 cm higher than the pre‐cut levels. The slight recovery in water table level plateaued after the third year. Rainfall interception by vegetation was also monitored, and after 10 years had reached nearly 50% of the pre‐cut rate. The immediate water table drawdown following drainage mitigated watering‐up within 40 m of a ditch. The persistent watering‐up observed in this study should encourage using sylvicultural systems adapted to boreal forested wetlands in order to prevent productivity loss and stand conversion. Copyright © 2008 John Wiley & Sons, Ltd.     

Palaeolimnological conditions inferred from fossil diatom assemblages and derivative spectral properties of sediments in thermokarst ponds of subarctic Quebec,Canada

Thermokarst ponds are widespread in arctic and subarctic regions, but little is known about their temporal evolution prior to human observations. This paper presents a pioneer biostratigraphic study conducted at a subarctic site with limnologically contrasted ponds located on the eastern shore of Hudson Bay, Canada. Fossil diatom and visible near infrared (VNIR) derivative spectral analyses were performed on short sediment cores, confirming the occurrence of three distinct stratigraphic facies as already inferred from an anterior sedimentological study: a lacustrine upper facies (UF) and a marine lower facies (LF), separated by an organic‐rich/peat transitional zone (TZ). Diatoms were almost absent from LF, but increased significantly in both TZ and UF. Identified diatom taxa were mainly benthic species (e.g. genera Fragilaria, Pinnularia), and their down‐core distribution appeared to be related to dissolved organic carbon (DOC) and possibly pH conditions. Diatom‐inferred DOC showed a decreasing trend towards the surface (potentially associated with an increase in pH), inverse to the general trend in this region, suggesting the action of other mechanisms on DOC, such as exhaustion of external inputs from limited catchments and the role of discontinuous peat layers (former surfaces of permafrost mounds) during the initial stages of pond formation. These bryophilous substrates in aerophilic habitats probably controlled diatom community composition. The combination of diatom and VNIR data revealed similar trends between (i) opal (amorphous silica) and diatom abundances; (ii) eukaryotic/prokaryotic algae ratio and anoxia or hypoxia in bottom waters; and (iii) limonite (iron oxide) and redox conditions in surface sediments. These findings indicate that diatom community changes and pond limnological evolution in the recent past were controlled mainly by autogenic processes (e.g. local vegetation/soil development, peat accumulation and erosion), rather than by allogenic forcing mechanisms (e.g. precipitation and temperature, geochemical leaching of the surrounding glaciomarine sediments).     

Formation of phyllosilicates in a fault zone affecting deeply buried arkosic sandstones: their influence on petrophysic properties (Annot sandstones, French external Alps)

The Restefond fault, located in the Late Eocene-Early Oligocene Alpine foreland basin, affects the well lithified and low porosity Grès d??Annot. The fault core zone is characterized by the occurrence of highly deformed sandstone lenses. Deformation inside the lenses corresponds to mm to sub-mm-spaced cleavage planes rich in phyllosilicates and up to cm-thick and dm-long quartz-calcite pure extensional veins. The cleavages are mostly composed of newly-formed synkinematic white mica and chlorite. By using thermodynamic thermometers based on the chemical composition of chlorite, a temperature of 200?±?20?°C of fault activity was computed. This temperature shows that the Restefond fault was active at burial conditions comprised between 6.5 and 8?km, assuming a mean geothermal gradient between 25 and 30?°Ckm^?1. The petrophysic properties of sandstones from the core zone and in the hanging and foot wall of the fault were determined on drilled plugs following three spacial directions. The permeability of the highly deformed sandstone from the core zone is about one order of magnitude higher than in the host rock. This increase in permeability occurs in the direction parallel to the S?CC structures and is explained by the occurrence of well-connected micropores localized between platy phyllosilicates. This study shows that the fault petrophysic properties are mostly controlled by the precipitation of synkinematic phyllosilicates under deep burial conditions.