Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
In this paper we present a new local remeshing algorithm that is dedicated to the problem of erosion in finite element models whose grid follows the movement of the free surface. The method, which we name Surface Lagrangian Remeshing (SLR), is adapted to 2D Lagrangian models which couple surface erosion with deformation of Earth materials. The remeshing procedure preserves nodes defining the surface submitted to erosion and removes nodes belonging to surface elements whose internal angles or area is critically low. This algorithm is ideally suited to track long term surface evolution. To validate the method we perform a set of numerical tests, using triangular finite elements, which compare the results obtained with the SLR algorithm with global remeshing and with analytical results. The results show good agreements with analytical solutions. Interpolation errors associated with remeshing are generated locally and numerical diffusion is restricted to the remeshed domain itself. In addition this method is computationally costless compared to classical global remeshing algorithm. We propose to couple the SLR method with the Dynamical Lagrangian Remeshing (DLR) algorithm to enable local remeshing only of Lagrangian models coupling large deformation of Earth materials with large erosion.
ISSN: 0098-3004

hal-00639472
https://hal.archives-ouvertes.fr/hal-00639472
DOI : 10.1016/j.cageo.2010.05.023

Éditeur(s) : HAL CCSD Inter-Research
Résumé : International audience
Ecological connectivity is important for effective marine planning and biodiversity conservation. Our aim was to identify factors important in influencing variation in benthic community structure on shallow rocky reefs in 2 regions of the Mediterranean Sea with contrasting oceanographic regimes. We assessed beta (β) diversity at 146 sites in the littoral and shallow sub-littoral from the Adriatic/Ionian Seas (eastern region) and Ligurian/Tyrrhenian Seas (western region) using a null modelling approach to account for variation in species richness. The distance decay relationship between species turnover within each region and geographic distance by sea was determined using generalised linear models. Mantel tests were used to examine correlations between β diversity and connectivity by ocean currents, estimated from Lagrangian dispersal simulations. Variation in β diversity between sites was partitioned according to environmental and spatial components using a distance-based redundancy approach. Species turnover along a gradient of geographic distance was greater by a factor of 3 to 5 in the western region than the eastern region, suggesting lower con-nectivity between sites. β diversity was correlated with connectivity by ocean currents at both depths in the eastern region but not in the western region. The influence of spatial and environmental predictors of β diversity varied considerably between regions, but was similar between depths. Our results highlight the interaction of oceanographic, spatial and environmental processes influencing benthic marine β diversity. Persistent currents in the eastern region may be responsible for lower observed β diversity compared to the western region, where patterns of water circulation are more variable.
Marine Ecology Progress Series

Droits : http://creativecommons.org/licenses/by/
hal-01357935
https://hal-univ-perp.archives-ouvertes.fr/hal-01357935
https://hal-univ-perp.archives-ouvertes.fr/hal-01357935/document
https://hal-univ-perp.archives-ouvertes.fr/hal-01357935/file/m553p001.pdf
DOI : 10.3354/meps11783

Éditeur(s) : HAL CCSD
Résumé : Among erosion processes, river incision is classicaly described as a key process controlling erosion of landscapes. Incision efficiency is mainly influenced by climate and erodibility. This latter is not only dependent on the nature of the bedrock, but also on its past deformation, which affects its rheological effective properties, such as fracture density. The main objectives of this thesis are: (1) to better constrain the relationship between effective properties and erodibility, and (2) to quantify the influence of erodibility and erosion on both the temporal and spatial building or decay of the topography. Several numerical tools are developed. A 1D formalism of landscape evolution is introduced, including river incision with stochastic distribution of water discharge and hillslope landsliding. A new remeshing algorithm called Surface Lagrangian Remeshing (SLR) is developed as a complement to remeshing algorithms dealing with internal elements. It allows one to take into account long-term erosion into 2D Lagrangian numerical codes based on triangular finite elements. Then the potentiality of measuring erodibility in-situ using a Schmidt hammer (R) is assessed for the active orogen of Taiwan, the diagenetic Annot sandstones and St Clement fault zone. Results suggest a strong control of R by effective properties. A linear model based on effective medium theory is applied to a fault zone with an unmacthed resolution (750 measures, 25 measures per square meter). The model successfully correlates R to fracture density. These results demonstrate that effective elasticity as well as erodibility are sensitive to the density and type of fractures. Finally I focus on the erosional and rheological conditions that allows reproducing post-orogenic evolution of mountain belts. A model coupling surface erosion and regional isostatic uplift is consistent with observations. The topographic decay and decrease of the ratio of surface elevation over crustal root thickness is at first order controlled by the initial geometry of the mountain belt and erosion efficiency. This new model highlights the control ofclimate and erodibility on the topographic decay and of lithospheric rheology on the perseverance of crustal roots.
Parmi les processus d'érosion, l'incision des rivières est classiquement décrit comme un processus clé contrôlant l'érosion des paysages. L'efficacité de l'incision est principalement influencée par le climat et par l'érodabilité. Ce dernier paramètre ne dépend pas seulement de la nature du substratum rocheux, mais aussi de sa déformation passée qui affecte ses propriétés rhéologiques équivalentes. Les principaux objectifs de cette thèse sont: (1) de mieux contraindre les relations entre propriétés équivalentes et érodabilité, et (2) de quantifier l'influence de l'érosion et de l'érodabilité sur la formation ou la décadence, spatiale et temporelle, de la topographie. Plusieurs outils numériques sont développés. Un formalisme 1D d'évolution des paysages est proposé, prenant en compte simultanément l'incision des rivières avec une distribution stochastique des débits en eau et l'érosion des versants par glissements de terrain. Un nouvel algorithme de remaillage appelé Surface Lagrangian Remeshing (SLR) est développé. Il permet de prendre en compte l'érosion à long terme dans les codes numériques 2D Lagrangiens basés sur des éléments finis triangulaires. Ensuite, la possibilité de mesurer in-situ l'érodabilité avec un marteau de Schmidt (R) est évaluée pour: l'orogène actif de Taiwan, les grès diagénétiques d'Annot et la zone de faille de St Clément. Les résultats suggèrent un fort contrôle de R par les propriétés équivalentes des roches testées. Un modèle linéaire, basé sur la théorie des milieux équivalents, est appliqué à une zone de faille avec une résolution inégalée (750 mesures, 25 mesures par mètre carré). Le modèle permet de corréler avec succès la densité de fracture et R. Ces résultats démontrent empiriquement que l'élasticité équivalente ainsi que l'érodabilité sont sensibles à la densité et à la nature des fractures. Enfin j'étudie les conditions d'érosion et de rhéologie qui permettent de reproduire l'évolution des chaînes de montagnes post-orogéniques. Un modèle couplant érosion en surface et soulèvement isostatique régional est compatible avec les observations. Les taux de décroissance topographique et de diminution du rapport de l'élévation de surface sur l'épaisseur de racine crustale sont contrôlés au premier ordre par la géométrie initiale de la chaîne de montagne et par l'efficacité de l'érosion. Ce nouveau modèle met en évidence le contrôle du climat et de l'érodabilité sur la décroissance topographique et de la rhéologie lithosphérique sur la persévérance des racines crustales.
https://tel.archives-ouvertes.fr/tel-00602341

tel-00602341
https://tel.archives-ouvertes.fr/tel-00602341
https://tel.archives-ouvertes.fr/tel-00602341/document
https://tel.archives-ouvertes.fr/tel-00602341/file/thesephilippesteer.pdf

Éditeur(s) : HAL CCSD European Geosciences Union
Résumé : We designed a thermo-mechanical numerical model for fast-spreading mid-ocean ridge with variable viscosity, hydrothermal cooling, latent heat release, sheeted dyke layer, and variable melt intrusion possibilities. The model allows for modulating several accretion possibilities such as the "gabbro glacier" (G), the "sheeted sills" (S) or the "mixed shallow and MTZ lenses" (M). These three crustal accretion modes have been explored assuming viscosity contrasts of 2 to 3 orders of magnitude between strong and weak phases and various hydrothermal cooling conditions depending on the cracking temperatures value. Mass conservation (stream-function), momentum (vorticity) and temperature equations are solved in 2-D cartesian geometry using 2-D, alternate direction, implicit and semi-implicit finite-difference scheme. In a first step, an Eulerian approach is used solving iteratively the motion and temperature equations until reaching steady states. With this procedure, the temperature patterns and motions that are obtained for the various crustal intrusion modes and hydrothermal cooling hypotheses display significant differences near the mid-ocean ridge axis. In a second step, a Lagrangian approach is used, recording the thermal histories and cooling rates of tracers travelling from the ridge axis to their final emplacements in the crust far from the mid-ocean ridge axis. The results show that the tracer's thermal histories are depending on the temperature patterns and the crustal accretion modes near the mid-ocean ridge axis. The instantaneous cooling rates obtained from these thermal histories betray these discrepancies and might therefore be used to characterize the crustal accretion mode at the ridge axis. These deciphering effects are even more pronounced if we consider the average cooling rates occurring over a prescribed temperature range. Two situations were tested at 1275-1125 °C and 1050-850 °C. The first temperature range covers mainly the crystallization range that is characteristic of the high temperature areas in the model (i.e. the near-mid-oceanic-ridge axis). The second temperature range corresponds to areas in the model where the motion is mainly laminar and the vertical temperature profiles are closer to conductive. Thus, this situation results in less discriminating efficiency among the crustal accretion modes since the thermal and dynamic properties that are described are common to all the crustal accretion modes far from the ridge axis. The results show that numerical modeling of thermo-mechanical properties of the lower crusts may bring useful information to characterize the ridge accretion structure, hydrothermal cooling and thermal state at the fast-spreading ridges and may open discussions with petrological cooling rate results.
ISSN: 1991-959X

hal-00907012
https://hal.archives-ouvertes.fr/hal-00907012
DOI : 10.5194/gmd-6-1659-2013

Éditeur(s) : American Geophysical Union
Résumé : Lagrangian analyses of a global ocean circulation model quantify the salinity changes experienced by the warm limb of the thermohaline circulation during the northward flow to the Atlantic deep convection regions. 6 Sv out of the estimated 10-Sv transfer from 45 degrees S to 47 degrees N flow through regions of prevailing surface evaporation: the southern and northern formation regions of Salinity Maximum Water and the Gulf of Cadiz/Mediterranean Sea domain. The remaining transport gains salinity through mixing with adjacent waters. As much as 6 Sv flow through the low-salinity surface mixed layer at the latitudes of the ITCZ whose effect annihilates that of the southern region of Salinity Maximum Water. Most of the salinity increase corresponds to the transformation of South to North Atlantic Central Water, with strong diapycnal transfers for the water that intersects the high and low salinity regions, and nearly isopycnal modifications for the water that avoids these regions.
Geophysical Research Letters ( GRL ) (0094-8276) (American Geophysical Union), 2006-03 , Vol. 33 , N. 6 , P. NIL_44-NIL_47

Droits : 2006 by the American Geophysical Union
http://archimer.ifremer.fr/doc/2006/publication-1213.pdf
DOI:10.1029/2005GL024938
http://archimer.ifremer.fr/doc/00000/1213/

Éditeur(s) : Elsevier Science Bv
Résumé : A Lagrangian approach based on a physical-biogeochemical modeling was used to compare the potential transfer of cadmium (Cd) from natural and anthropogenic sources to plankton communities (Cd-uptake) in the NorthWest African upwelling. In this region, coastal upwelling was estimated to be the main natural source of Cd while the most significant anthropogenic source for marine ecosystem is provided by phosphate industry. In our model experiment, Cd-uptake (natural or anthropogenic) in the North-West African upwelling is the result of an interplay between the Cd dispersion (by advection processes) and the simulated biological productivity. In the Moroccan waters, advection processes limit the residence time of water masses resulting in a low natural Cd-uptake by plankton communities while anthropogenic Cd-uptake is high. As expected, the situation is reversed in the Senegalo-Mauritanian upwelling where natural Cd-uptake is higher than anthropogenic Cd-uptake. Based upon an estimate of Cd sources, our modeling study shows, unexpectedly, that the anthropogenic signal of potential Cd-bioaccumulation in the Moroccan upwelling is of the same order of magnitude as the natural signal mainly present in the Senegalo-Mauritanian upwelling region. A comparison with observed Cd levels in mollusk and fishes, which shows overall agreement with our simulations, is confirming our estimates.
Science Of The Total Environment (0048-9697) (Elsevier Science Bv), 2015-02 , Vol. 505 , P. 870-888

Droits : 2014 Elsevier B.V. All rights reserved.
http://archimer.ifremer.fr/doc/00252/36310/35200.pdf
DOI:10.1016/j.scitotenv.2014.10.045
http://archimer.ifremer.fr/doc/00252/36310/

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Tropical islands are characterized by thermal and orographical forcings which may generate microscale air mass circulations. The Lesser Antilles Arc includes small tropical islands (width lower than 50 km) where a total of one-and-a-half million people live. Air quality over this region is affected by anthropogenic and volcanic emissions, or saharan dust. To reduce risks for the population health, the atmospheric dispersion of emitted pollutants must be predicted. In this study, the dispersion of anthropogenic nitrogen oxides (NOx) is numerically modelled over the densely populated area of the Guadeloupe archipelago under weak trade winds, during a typical case of severe pollution. The main goal is to analyze how microscale resolutions affect air pollution in a small tropical island. Three resolutions of domain grid are selected: 1 km, 333 m and 111 m. The Weather Research and Forecasting model (WRF) is used to produce real nested microscale meteorological fields. Then the weather outputs initialize the Lagrangian Particle Dispersion Model (FLEXPART). The forward simulations of a power plant plume showed good ability to reproduce nocturnal peaks recorded by an urban air quality station. The increase in resolution resulted in an improvement of model sensitivity. The nesting to subkilometer grids helped to reduce an overestimation bias mainly because the LES domains better simulate the turbulent motions governing nocturnal flows. For peaks observed at two air quality stations, the backward sensitivity outputs identified realistic sources of NOx in the area. The increase in resolution produced a sharper inverse plume with a more accurate source area. This study showed the first application of the FLEXPART-WRF model to microscale resolutions. Overall, the coupling model WRF-LES-FLEXPART is useful to simulate the pollutant dispersion during a real case of calm wind regime over a complex terrain area. The forward and backward simulation results showed clearly that the subkilometer resolution of 333 m is necessary to reproduce realistic air pollution patterns in this case of short-range transport over a complex terrain area. Globally, this work contributes to enrich the sparsely documented domain of real nested microscale air pollution modelling. This study dealing with the determination of the proper resolution grid and proper turbulence scheme, is of significant interest to the near-source and complex terrain air quality research community.
ISSN: 1352-2310

hal-01349208
https://hal.archives-ouvertes.fr/hal-01349208
DOI : 10.1016/j.atmosenv.2016.05.028

Résumé : Ce travail est l'expression d'une volonté de chercheurs, de la Carabe d'améliorer les connaissances scientifiques méso- et micro-météorologiques appliquées aux milieux insulaires sous influence des alizés et de développer la recherche dans ces domaines.On sait que le phénomène météorologique le plus remarquable impactant les îles de la Caraïbe reste le cyclone tropical. Mais d'autres phénomènes, a des échelles inferieures, telles que les pluies intenses, les houles, la dégradation de la qualité de l'air ont une importance aigüe en termes de risques naturels ou de risques sur la sante. Ces exemples attestent la nécessité d'utiliser des méthodes de descente d'échelle pour exploiter l'information météorologique et/ou climatique de grande échelle et dériver des scenarios locaux et régionaux appliques aux territoires insulaires. Ce défi est important car l'attente d'analyses scientifiques pertinentes dans ces domaines est grande.Les travaux présentés dans ce mémoire ont pour principaux objectifs la simulation numérique puis l'analyse des mécanismes méso- et micro-échelles qui induisent des circulations locales diurnes et nocturnes sur l'archipel de la Guadeloupe à l'aide de codes numériques météorologiques éprouves car largement utilises en recherche et en prévision opérationnelle.Ils constituent donc la première étude de modélisation numérique à haute résolution en basse atmosphère, par descente d'échelle dynamique, pour des intervalles d'espace compris entre 1 km et 111 m sur cet archipel.Le modelé atmosphérique Weather Research and Forecasting ARW V3.4 (WRF ARW), non-hydrostatique, a été utilise pour l'ensemble des simulations pour modéliser la troposphère depuis l'échelle globale à l'échelle de la turbulence.Trois situations météorologiques classiques réelles d'une durée de 48 heures, correspondant à 80 % des situations météorologiques observées dans la zone, alizes soutenus (AS), alizes moyens (AM), alizes faibles (AF) ont été examinées. Ces situations sont caractérisées par les valeurs du nombre de Froude local suivantes : 0,82 (AS), 0,41 (AM) et 0,21 (AF). Six domaines de modélisation ont été sélectionnés pour effectuer les descentes d'échelle dynamiques : D01 (maille de 27 km), D02 (maille de 9 km), D03 (maille de 3 km), D04 (maille de 1 km), D05 (maille de 333 m) et D06 (maille de 111 m) avec soixante-dix niveaux verticaux. Les quatre premiers domaines (D01 à D04) couvrent l'archipel de la Guadeloupe et sont utilisés en mode méso-échelle à l'aide d'un schéma d'ensemble de couche limite planétaire YSU. Les domaines D05 (couvrant l'île de la Basse-Terre et le centre de l'archipel) et D06 (couvrant la zone littorale et rurale du Petit Cul-de-Sac Marin et la zone urbaine de l'agglomération pointoise) sont utilisés en mode Large Eddy Simulation avec une fermeture de la turbulence 1,5 TKE 3D. Le modèle WRF a été forcé toutes les six heures par l'assimilation des champs d'analyses globales du modèle NCEP FNL (1° de résolution). Les simulations effectuées ont permis d'obtenir des champs de variables météorologiques 10-minutes à très hautes résolutions spatiales.Les résultats des simulations méso- et micro-échelles ont été confrontés aux valeurs expérimentales obtenues à l'aide de capteurs places sur des mâts météorologiques (campagne Atmo-Mangrov et réseau de mesure Météo-France). Il s'agissait d'optimiser l'utilisation des couplages de codes numériques tout en conservant la possibilité de les confronter aux observations expérimentales.Les résultats des simulations numériques micro-échelles, des différents cas (AS, AM, AF) sont utilisées pour forcer (c'est-à-dire définir les conditions limites) un modèle lagrangien de dispersion de particules : FLEXPART. Le système couple FLEXPART-WRF a été employé dans le domaine D06 pour étudier la dispersion du panache d'oxydes d'azote émis par la principale usine de production d'électricité de l'île.
The present work expresses the will of Caribbean researchers to improve the meso- and micro-meteorological scientific knowledge of the trade winds influenced island areas, and to develop research in these domains.It is well known that tropical hurricanes remain the most remarkable meteorological phenomena that affect the Caribbean islands. But some other phenomena, of smaller scale, such as intense rainfall events, swells, or air quality degradation, are of extreme importance for natural or health hazards. These examples show the need to use downscaling methods to exploit large scale meteorological or climatic information, and to obtain local and regional scenarios for the island areas. This is an important challenge, as sound scientific studies in these matters are eagerly expected. The aim of the research works exposed in the present dissertation is numerical simulation and analysis of the meso- and microscale mechanisms that induce diurnal and nocturnal local circulations in the Guadeloupe archipelago, using numerical meteorological models that are widely used in research and in operational forecasting.These works represent the first high-resolution (1 km to 111 m) numerical study of the lower atmosphere over the Guadeloupe archipelago. The Weather Research and Forecasting ARW 3.4 (WRF-ARW) model is used to simulate the troposphere from global scale to turbulence scale. Real cases of three typical weather types (80% of cases during a year) are examined during 48 hours : strong trade winds (STW), medium trade winds (MTW) and weak trade winds (WTW). These weather types are characterized by values of the local Froude number : 0.82 (STW), 0.41 (MTW) and 0.21 (WTW). Six domains have been selected for the dynamical downscaling : D01 (grid spacing of 27 km), D02 (grid spacing of 9 km), D03 (grid spacing of 3 km), D04 (grid spacing of 1 km), D05 (grid spacing of 333 m) and D06 (grid spacing of 111 m), including 70 vertical levels. The first four domains (D01 to D04) cover the Guadeloupe archipelago and are used in the meso-scale simulations with the planetary boundary layer scheme YSU (ensemble mean). Domain D05 (covering the Basse-Terre island and the middle of the archipelago) and domain D06 (covering the coastal and rural area of Le Petit Cul-de-Sac Marin and the urban area of Pointe-à-Pitre), are employed in the micro-scale simulation (LES) with the 3D TKE 1.5 order closure scheme. WRF has been 6 hourly reinitialized with the NCEP FNL global analyses (resolution of 1°). These simulations permitted to obtain 10-minutes meteorological variable fields with a very high resolution (111 m).Meso-scale and micro-scale model results have been evaluated with observational data from meteorological stations (field campaign Atmo-Mangrov, French Met Office).Once validated, the micro-scale model outputs have been used for the assimilation of the lagrangian particle dispersion model : FLEXPART. The coupling FLEXPART-WRF has been employed in domain D06, in order to analyze the dispersion of the nitrogen oxide plume emitted by the main power plant of the archipelago.

http://www.theses.fr/2014AGUY0735/document

Éditeur(s) : Wiley-blackwell Publishing, Inc
Résumé : Oceanic larvae of the European (Anguilla anguilla) and American (A. rostrata) eels have to cross the Atlantic Ocean from the Sargasso Sea to European or North American coasts before entering continental habitats. In some European rivers, eel recruitment is now < 1% of levels in the 1980s. A better understanding of the effects of anthropogenic pressures and environmental fluctuations on eel larvae and subsequent recruitment is a prerequisite to build efficient management plans. The present paper provides insight into the critical oceanic phase of the eel life cycle with a focus on the duration of the larval migration whose estimates varies between 7 months and more than 2 years in both species. Does this range correspond to a natural variability in larval duration or does it stem from methodological artefacts? We first review the different methods used to estimate the duration of larval migration and critically describe their possible sources of misinterpretation. We then evaluate the consistency of these methods with the current knowledge about the ecology and physiology of eel larvae and the physical oceanography. While a moderate discrepancy in migration duration was found between methods for the American eel, the discrepancy was large in the European eel. In this species, otolith microstructure studies indicated migration durations between 7 and 9 months, while other methods pointed to durations of about 2 years. We show that estimates in favour of a long migration duration seem more robust to methodological caveats than methods estimating short durations of migration.
Fish And Fisheries (1467-2960) (Wiley-blackwell Publishing, Inc), 2010-09 , Vol. 11 , N. 3 , P. 289-306

Droits : 2010 Blackwell Publishing Ltd.
http://archimer.ifremer.fr/doc/00011/12208/9033.pdf
DOI:10.1111/j.1467-2979.2010.00362.x
http://archimer.ifremer.fr/doc/00011/12208/

Éditeur(s) : HAL CCSD
Résumé : Meteorological phenomena with scales lower than hurricanes, from extreme rainfalls to pollution peaks, represent important environmental and health risks in the Caribbean. These examples show the need to downscale large-scale meteorological data to island-scale fields.The main goal of the present numerical study is the analysis of the meso- and micro-scale mechanisms inducing diurnal and nocturnal circulations in the Guadeloupe archipelago. The Weather Research and Forecasting model is employed. The micro-scale results are obtained by the coupling of the 1-km simulated fields and the Large Eddy Simulation method. Then, in order to analyze the dispersion of a pollutants plume, the micro-scale outputs with 111 m of resolution are used to force the lagrangian particles dispersion model FLEXPART. The study has been conducted for three real cases of weather types linked with strong, medium and weak trade winds.The results of simulations, validated with available observed data, permitted to simulate all local circulations, including a nocturnal counter-current flow which may induce a pollution peak in the urban area of Pointe-à-Pitre city.These works, which may be applied to other Lesser Antilles islands, suggest many environmental and energy applications.
Les phénomènes météorologiques d'échelles inférieures aux cyclones, depuis les pluies intenses, jusqu'à la dégradation de la qualité de l'air, ont pour les Caraïbes une importance aiguë en termes de risques naturels ou de risques sur la santé. Ces exemples attestent de la nécessité d'utiliser des méthodes de descentes d'échelle pour exploiter l'information météorologique de grande échelle et dériver des scénarios locaux appliqués aux territoires insulaires. Les travaux présentés ont pour principaux objectifs la simulation numérique et l'analyse des mécanismes méso- et micro-échelles qui induisent des circulations locales diurnes et nocturnes sur l'archipel de la Guadeloupe, à l'aide du modèle numérique météorologique Weather Research and Forecasting. Les résultats des simulations micro-échelles sont obtenus par le couplage des champs modélisés à 1 km avec la méthode Large Eddy Simulation. Les sorties micro-échelles, à 111 m, sont ensuite utilisées pour forcer le modèle lagrangien de dispersion de particules FLEXPART afin d'étudier la dispersion d'un panache de polluants. L'étude a été menée pour trois cas réels de situations météorologiques correspondant à des alizés soutenus, moyens et faibles. Les résultats des simulations, après validation, et comparaison aux données d'observation disponibles, permettent de restituer l'ensemble des circulations locales et soulignent l'existence d'un écoulement nocturne à contre courant susceptible d'induire une forte pollution sur la zone urbaine de Pointe-à-Pitre. Cette étude, transposable aux autres îles des petites Antilles, laisse envisager de nombreuses applications environnementales et énergétiques.
https://hal.archives-ouvertes.fr/tel-01074600

Droits : info:eu-repo/semantics/OpenAccess
tel-01074600
https://hal.archives-ouvertes.fr/tel-01074600
https://hal.archives-ouvertes.fr/tel-01074600/document
https://hal.archives-ouvertes.fr/tel-01074600/file/rcece_thesis.pdf

Éditeur(s) : American Geophysical Union
Résumé : The origins and fates of the Atlantic salinity maximum water (SMW), formed through excess evaporation in the tropics and subtropics of both hemispheres, are studied using monthly mean outputs of a numerical simulation of the world ocean climatological circulation. After defining formation domains from the surface salinity field and the vertical stratification, a Lagrangian technique is used to estimate the formation rates and main pathways in each hemisphere and the role of this water in the framework of the warm water return flow of the meridional overturning cell. Formation rates around 9 and 11 Sv are found in the Southern and Northern Hemispheres, respectively. While the export of the southern SMW from its formation area is realized by the western boundary currents, that of the northern SMW mainly results from interior subduction. Equatorward of the formation regions, a fraction of each SMW variety is entrained in the subtropical cells that connect the subtropics to the equatorial region. Poleward of them, both varieties are seen to feed the regions of subtropical mode water (STMW) formation around 35degrees of latitude in both hemispheres. The bulk of the transport associated with each variety eventually turns northward: This amounts to similar to6 Sv of southern SMW gathered in the North Brazil Undercurrent (NBUC), and similar to10 Sv of northern SMW found in the Gulf Stream at 35degreesN, of which 8 Sv have gone through the Caribbean Sea. Of the 13.4-Sv northward transport of the meridional overturning cell estimated by the model at 47degreesN, more than 50% (6.9 Sv) is found to have transited through at least one of the SMW regions. This gives an indication of the likely important role of SMW formation in the observed northward salinity increase of the upper Atlantic Ocean.
Journal of Geophysical Union - Research C - Oceans (0148-0227) (American Geophysical Union), 2002-09 , Vol. 107 , N. C10 , P. NIL_610-NIL_624

Droits : 2002 by the American Geophysical Union
http://archimer.ifremer.fr/doc/2002/publication-769.pdf
DOI:10.1029/2002JC001318
http://archimer.ifremer.fr/doc/00000/769/

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : Some major challenges for geophysicists and structural geologists using three-dimensional boundary element method codes (3D-BEM) are: (1) reducing the amount of memory required to solve large and dense systems and (2) incorporation of inequality constraints such as traction inequality constraints (TIC) and displacement inequality constraints (DIC). The latter serves two purposes. First, for example, inequality constraints can be used to simulate frictional slip (using TIC). Second, these constraints can prevent element interpenetration while allowing opening mode (using DIC). We have developed a method that simultaneously incorporates both types of functionality of the inequality constraints. We show that the use of an appropriate iterative solver not only avoids the allocation of significant memory for solving the system (allowing very large model computation and simplifying parallelization on multi-core processors), but also admits interesting features such as natural incorporation of TICs and DICs. Compared to other techniques of contact management (e.g., Lagrange multipliers, penalty method, or complementarity problem), this new simple methodology, which does not use any incremental trial-and-error procedures, brings more flexibility, while making the system more stable and less subject to round-off errors without any computational overhead. We provide validations and comparisons of the inequality constraints implementation using 2D analytical and numerical solutions.
ISSN: 1420-0597

hal-00561035
https://hal.archives-ouvertes.fr/hal-00561035
DOI : 10.1007/s10596-009-9170-x