Interactions spatiales et auto-organisation des végétations semi-arides

Les recherches récapitulées dans cette thèse de doctorat ont porté sur les causes de l’organisation spatiale des végétations périodiques. Ces structures paysagères aux motifs réguliers, tachetés, tigrés ou labyrinthiques, d’échelle décamétrique à hectométrique, couvrant des étendues considérables sur au moins trois continents, constituent un cas d’école dans l’étude des processus endogènes présidant à l’hétérogénéité du couvert végétal. Ces structures prennent place sur un substrat homogène, mis à part la rétroaction du couvert lui-même, et sont marquées par des écotones abrupts et la persistance d’une proportion considérable de sol nu. Plusieurs modèles ont mis en avant l’existence possible d’un phénomène d’auto-organisation du couvert, qui verrait une structure d’ensemble émerger des interactions locales entre individus. Ces modèles se basent sur le jeu simultané de la consommation de la ressource (compétition) et de l’amélioration de l’un ou l’autre des éléments du bilan de la même ressource par le couvert (facilitation). La condition à l’existence d’une structure d’ensemble spatialement périodique et stable réside dans une différence entre la portée de la compétition (plus grande) et celle de la facilitation. L’apparition de ces structures est modulée par le taux de croissance biologique, qui est le reflet des contraintes extérieures telles que l’aridité, le pâturage ou la coupe de bois. Le modus operandi des interactions spatiales supposées entre individus reste largement à préciser.

Nos recherches ont été menées au sud-ouest de la République du Niger, à l’intérieur et dans les environs du parc Régional du W. Trois axes ont été explorés :(i) Une étude de la dépendance spatiale entre la structure de la végétation (biovolumes cartographiés) et les paramètres du milieu abiotique (relief, sol), sur base d’analyses spectrales et cross-spectrales par transformée de Fourier (1D et 2D). (ii) Une étude diachronique (1956, 1975 et 1996) à large échelle (3000 km²) de l’influence de l’aridité et des pressions d’origine anthropique sur l’auto-organisation des végétations périodiques, basée sur la caractérisation de la structure spatiale des paysages sur photos aériennes via la transformée de Fourier en 2D. (iii) Trois études portant sur les interactions spatiales entre individus :En premier lieu, via l’excavation des systèmes racinaires (air pulsé) ;Ensuite, par un suivi spatio-temporel du bilan hydrique du sol (blocs de gypse) ;Enfin, via le marquage de la ressource par du deutérium.

Nous avons ainsi pu établir que les végétations périodiques constituent bien un mode d’auto-organisation pouvant survenir sur substrat homogène et modulé par les contraintes climatiques et anthropiques. Un ajustement rapide entre l’organisation des végétations périodiques et le climat a pu être montrée en zone protégée. La superficie et l’organisation des végétations périodiques y ont tour à tour progressé et régressé en fonction d’épisodes secs ou humides. Par contre, en dehors de l’aire protégée, la possibilité d’une restauration du couvert semble fortement liée au taux d’exploitation des ressources végétales. Ces résultats ont d’importantes implications quant à la compréhension des interactions entre climat et écosystèmes et à l’évaluation de leurs capacités de charge. La caractérisation de la structure spatiale des végétations arides, notamment par la transformée de Fourier d’images HR, devrait être généralisée comme outil de monitoring de l’état de ces écosystèmes. Nos études portant sur les modes d’interactions spatiales ont permis de confirmer l’existence d’une facilitation à courte portée du couvert végétal sur la ressource. Cependant, cette facilitation ne semble pas s’exercer sur le terme du bilan hydrique traditionnellement avancé, à savoir l’infiltration, mais plutôt sur le taux d’évaporation (deux fois moindre à l’ombre des canopées). Ce mécanisme exclut l’existence de transferts diffusifs souterrains entre sols nu et fourrés. Des transferts inverses semblent d’ailleurs montrés par le marquage isotopique. L’étude du bilan hydrique et la cartographie du micro-relief, ainsi que la profondeur fortement réduite de la zone d’exploitation racinaire, jettent de sérieux doutes quant au rôle communément admis des transferts d’eau par ruissellement/diffusion de surface en tant que processus clé dans la compétition à distance entre les plantes. L’alternative réside dans l’existence d’une compétition racinaire de portée supérieure aux canopées. Cette hypothèse trouve une confirmation tant par les rhizosphères excavées, superficielles et étendues, que dans le marquage isotopique, montrant des contaminations d’arbustes situés à plus de 15 m de la zone d’apport. De même, l’étude du bilan hydrique met en évidence les influences simultanées et contradictoires (facilitation/compétition) des ligneux sur l’évapotranspiration.

/

This PhD thesis gathers results of a research dealing with the causes of the spatial organisation of periodic vegetations. These landscape structures, featuring regular spotted, labyrinthine or banded patterns of decametric to hectometric scale, and extending over considerable areas on at least three continents, constitute a perfect study case to approach endogenous processes leading to vegetation heterogeneities. These patterns occur over homogeneous substratum, except for vegetation’s own feedbacks, and are marked by sharp ecotones and the persistence of a considerable amount of bare soil. A number of models suggested a possible case of self-organized patterning, in which the general structure would emerge from local interactions between individuals. Those models rest on the interplay of competitive and facilitative effects, relating to soil water consumption and to soil water budget enhancement by vegetation. A general necessary condition for pattern formation to occur is that negative interactions (competition) have a larger range than positive interactions (facilitation). Moreover, all models agree with the idea that patterning occurs when vegetation growth decreases, for instance as a result of reduced water availability, domestic grazing or wood cutting, therefore viewing patterns as a self-organised response to environmental constraints. However the modus operandi of the spatial interactions between individual plants remains largely to be specified.

We carried out a field research in South-West Niger, within and around the W Regional Park. Three research lines were explored: (i) The study of the spatial dependency between the vegetation pattern (mapped biovolumes) and the factors of the abiotic environment (soil, relief), on the basis of spectral and cross-spectral analyses with Fourier transform (1D and 2D). (ii) A broad scale diachronic study (1956, 1975, 1996) of the influence of aridity and human induced pressures on the vegetation self-patterning, based on the characterisation of patterns on high resolution remote sensing data via 2D Fourier transform. (iii) Three different approaches of the spatial interactions between individuals: via root systems excavation with pulsed air; via the monitoring in space and time of the soil water budget (gypsum blocks method); and via water resource labelling with deuterated water.

We could establish that periodic vegetations are indeed the result of a self-organisation process, occurring in homogeneous substratum conditions and modulated by climate and human constraints. A rapid adjustment between vegetation patterning and climate could be observed in protected zones. The area and patterning of the periodic vegetations successively progressed and regressed, following drier or wetter climate conditions. On the other hand, outside protected areas, the restoration ability of vegetation appeared to depend on the degree of vegetation resource exploitation. These results have important implications regarding the study of vegetation-climate interactions and the evaluation of ecosystems’ carrying capacities. Spatial pattern characterisation in arid vegetations using Fourier transform of HR remote sensing data should be generalised for the monitoring of those ecosystems. Our studies dealing with spatial interaction mechanisms confirmed the existence of a short range facilitation of the cover on water resource. However, this facilitation does not seem to act through the commonly accepted infiltration component, but rather on the evaporative rate (twice less within thickets). This mechanism excludes underground diffusive transfers between bare ground and vegetation. Inverse transfers were even shown by deuterium labelling. Water budget study and micro-elevation mapping, along with consistent soil shallowness, together cast serious doubts on the traditional mechanism of run-off/diffusion of surface water as a key process of the long range competition between plants. An alternative explanation lies in long range root competition. This hypothesis find support as well in the excavated root systems, shallow and wide, as in isotopic labelling, showing contaminations of shrubs located up to 15 m of the irrigated area. Water budget study also evidenced simultaneous contradictory effects (facilitation/competition) of shrubs on evapotranspiration.

Comparison of the experimental, semi-experimental and ab initio equilibrium structures of acetylene: Influence of relativisitic effects and of the diagonal Born-Oppenheimer corrections

The equilibrium structure of acetylene (also named ethyne) has been reinvestigated to resolve the small discrepancies noted between different determinations. The size of the system as well as the large amount of available experimental data provides the quite unique opportunity to check the magnitude and relevance of various contributions to equilibrium structure as well as to verify the accuracy of experimental results. With respect to pure theoretical investigation, quantum-chemical calculations at the coupled-cluster level have been employed together with extrapolation to the basis set limit, consideration of higher excitations in the cluster operator, inclusion of core correlation effects as well as relativistic and diagonal Born-Oppenheimer corrections. In particular, it is found that the extrapolation to the complete basis set limit, the inclusion of higher excitations in the electronic-correlation treatment and the relativistic corrections are of the same order of magnitude. It also appears that a basis set as large as a core-valence quintuple-zeta set is required for accurately accounting for the inner-shell correlation contribution. From a pure experimental point of view, the equilibrium structure has been determined using very accurate rotational constants recently obtained by a global analysis (that is to say that all non-negligible interactions are explicitely included in the Hamiltonian matrix) of rovibrational spectra. Finally, a semi-experimental equilibrium structure (where the equilibrium rotational constants are obtained from the experimental ground state rotational constants and computed rovibrational corrections) has been obtained from the available experimental ground-state rotational constants for ten isotopic species corrected for computed vibrational corrections. Such a determination led to the revision of the ground-state rotational constants of two isotopologues, thus showing that structural determination is a good method to identify errors in experimental rotational constants. The three structures are found in a very good agreement, and our recommended values are rCC 120.2958(7) pm and rCH 106.164(1) pm. © 2011 American Institute of Physics.

Dinámica y consecuencias del cambio en la cobertura y el uso del suelo en el Chaco semi-árido

La deforestación en Sudamérica afecta principalmente a tres ecosistemas: El Cerrado en Brasil, la selva de Chiquitanos en Bolivia y el Gran Chaco en Bolivia, Paraguay y Argentina, siendo en estos dos últimos países en donde ocurren las mayores transformaciones del paisaje para la producción de commodities para exportación. En la presente tesis, para la porción Noroeste del Chaco Argentino, se analizó la dinámica de cambios de los últimos 30 años; se evaluó la ocurrencia del modelo cambio de uso del suelo denominado "Transición Forestal"; se estudió la influencia de factores locales que controlan la localización de desmontes; y se cuantificó el impacto de dicha transformación sobre el nivel de provisión de servicios ecosistémicos intermedios relativos a la dinámica del C. En el período 1977-2007 ocurrieron cambios que alcanzaron más del 26 por ciento del área de estudio. Los desmontes para actividad agropecuaria alcanzaron un total de 4,5 millones ha, de los cuales el 53 por ciento ocurrieron en el último período (1997-07), siendo los bosques secos y los pastizales las coberturas más afectadas. De continuar con esta tendencia, ocurrirá una inversión del paisaje en un periodo comprendido entre 40 y 100 años, en donde la actividad agropecuaria comenzaría a ser dominante en el paisaje chaqueño. No están ocurriendo ninguno de los modelos de cambios de uso del suelo de compatibilización de producción-conservación discutidos en esta tesis: (a) "Transición Forestal" (Forest Transition), (b) "separación-territorial" (land-sparing); o (c) "integración-territorial" (land-sharing). Se propone un nuevo modelo de cambio de uso de suelo que probablemente se verifica en otras regiones del Gran Chaco Americano, al que se denominó "modelo tipo Tsunami" por sus características de avance en forma de ola, que a su paso homogeniza el paisaje bajo el punto de vista estructural y funcional.

Host species of bees from the tribe Meliponini Nest in semi - arid Chaco forest near Las Lomitas, Formosa (Argentina)

p.85-91

Host species of bees from the tribe Meliponini Nest in semi - arid Chaco forest near Las Lomitas, Formosa (Argentina)

p.85-91

Dinámica y consecuencias del cambio en la cobertura y el uso del suelo en el Chaco semi-árido

La deforestación en Sudamérica afecta principalmente a tres ecosistemas: El Cerrado en Brasil, la selva de Chiquitanos en Bolivia y el Gran Chaco en Bolivia, Paraguay y Argentina, siendo en estos dos últimos países en donde ocurren las mayores transformaciones del paisaje para la producción de commodities para exportación. En la presente tesis, para la porción Noroeste del Chaco Argentino, se analizó la dinámica de cambios de los últimos 30 años; se evaluó la ocurrencia del modelo cambio de uso del suelo denominado "Transición Forestal"; se estudió la influencia de factores locales que controlan la localización de desmontes; y se cuantificó el impacto de dicha transformación sobre el nivel de provisión de servicios ecosistémicos intermedios relativos a la dinámica del C. En el período 1977-2007 ocurrieron cambios que alcanzaron más del 26 por ciento del área de estudio. Los desmontes para actividad agropecuaria alcanzaron un total de 4,5 millones ha, de los cuales el 53 por ciento ocurrieron en el último período (1997-07), siendo los bosques secos y los pastizales las coberturas más afectadas. De continuar con esta tendencia, ocurrirá una inversión del paisaje en un periodo comprendido entre 40 y 100 años, en donde la actividad agropecuaria comenzaría a ser dominante en el paisaje chaqueño. No están ocurriendo ninguno de los modelos de cambios de uso del suelo de compatibilización de producción-conservación discutidos en esta tesis: (a)"Transición Forestal" (Forest Transition), (b)"separación-territorial" (land-sparing); o (c)"integración-territorial" (land-sharing). Se propone un nuevo modelo de cambio de uso de suelo que probablemente se verifica en otras regiones del Gran Chaco Americano, al que se denominó "modelo tipo Tsunami" por sus características de avance en forma de ola, que a su paso homogeniza el paisaje bajo el punto de vista estructural y funcional.

Viscoelastic flow through a planar contraction using a semi-Lagrangian finite volume method

A semi-Lagrangian finite volume scheme for solving viscoelastic flow problems is presented. A staggered grid arrangement is used in which the dependent variables are located at different mesh points in the computational domain. The convection terms in the momentum and constitutive equations are treated using a semi-Lagrangian approach in which particles on a regular grid are traced backwards over a single time-step. The method is applied to the 4 : 1 planar contraction problem for an Oldroyd B fluid for both creeping and inertial flow conditions. The development of vortex behaviour with increasing values of We is analyzed.

Conservative semi-Lagrangian finite volume schemes

Semi-Lagrangian finite volume schemes for the numerical approximation of linear advection equations are presented. These schemes are constructed so that the conservation properties are preserved by the numerical approximation. This is achieved using an interpolation procedure based on area-weighting. Numerical results are presented illustrating some of the features of these schemes.

A semi-Lagrangian finite volume method for Newtonian contraction flows

A new finite volume method for solving the incompressible Navier--Stokes equations is presented. The main features of this method are the location of the velocity components and pressure on different staggered grids and a semi-Lagrangian method for the treatment of convection. An interpolation procedure based on area-weighting is used for the convection part of the computation. The method is applied to flow through a constricted channel, and results are obtained for Reynolds numbers, based on half the flow rate, up to 1000. The behavior of the vortex in the salient corner is investigated qualitatively and quantitatively, and excellent agreement is found with the numerical results of Dennis and Smith [Proc. Roy. Soc. London A, 372 (1980), pp. 393-414] and the asymptotic theory of Smith [J. Fluid Mech., 90 (1979), pp. 725-754].

AC & DC magnetic levitation and semi-levitation modelling

This work presents computation analysis of levitated liquid thermal and flow fields with free surface oscillations in AC and DC magnetic fields. The volume electromagnetic force distribution is continuously updated with the shape and position change. The oscillation frequency spectra are analysed for droplets levitation against gravity in AC and DC magnetic fields at various combinations. For larger volume liquid metal confinement and melting the semi-levitation induction skull melting process is simulated with the same numerical model. Applications are aimed at pure electromagnetic material processing techniques and the material properties measurements in uncontaminated conditions.

The semi-automatic parallelisation of scientific application codes using a computer aided parallelisation toolkit

The shared-memory programming model can be an effective way to achieve parallelism on shared memory parallel computers. Historically however, the lack of a programming standard using directives and the limited scalability have affected its take-up. Recent advances in hardware and software technologies have resulted in improvements to both the performance of parallel programs with compiler directives and the issue of portability with the introduction of OpenMP. In this study, the Computer Aided Parallelisation Toolkit has been extended to automatically generate OpenMP-based parallel programs with nominal user assistance. We categorize the different loop types and show how efficient directives can be placed using the toolkit's in-depth interprocedural analysis. Examples are taken from the NAS parallel benchmarks and a number of real-world application codes. This demonstrates the great potential of using the toolkit to quickly parallelise serial programs as well as the good performance achievable on up to 300 processors for hybrid message passing-directive parallelisations.

Semi-Parametric Analysis of Extreme Forest Fires

Forest fires can cause extensive damage to natural resources and properties. They can also destroy wildlife habitat, affect the forest ecosystem and threaten human lives. In this paper incidences of extreme wildland fires are modelled by a point process model which incorporates time-trend. A model based on a generalised Pareto distribution is used to model data on acres of wildland burnt by extreme fire in the US since 1825. A semi-parametric smoothing approach, which is very useful in exploratory analysis of changes in extremes, is illustrated with the maximum likelihood method to estimate model parameters.