A Parallel 2D Finite Volume Scheme for Solving Systems of Balance Laws

Total Page:16

File Type:pdf, Size:1020Kb

A Parallel 2D Finite Volume Scheme for Solving Systems of Balance Laws Comput. Methods Appl. Mech. Engrg. 195 (2006) 2788–2815 www.elsevier.com/locate/cma A parallel 2d finite volume scheme for solving systems of balance laws with nonconservative products: Application to shallow flows M.J. Castro a,*, J.A. Garcı´a-Rodrı´guez a, J.M. Gonza´lez-Vida b, C. Pare´s a a Dpto. de Ana´lisis Matema´tico, Universidad de Ma´laga, Campus Teatinos s/n, 29071 Ma´laga, Spain b Dpto. de Matema´tica Aplicada, Universidad de Ma´laga, Spain Received 20 January 2005; received in revised form 1 July 2005; accepted 5 July 2005 Abstract The goal of this paper is to construct parallel solvers for 2d hyperbolic systems of conservation laws with source terms and nonconservative products. More precisely, finite volumes solvers on nonstructured grids are considered. The method of lines is applied: at every intercell a projected Riemann problem along the normal direction is considered which is discretized by means of the numerical schemes presented in [M.J. Castro, J. Macı´as, C. Pare´s. A Q-scheme for a class of systems of coupled conservation laws with source term. Application to a two-layer 1-D shallow water system, ESAIM: M2AN 35 (1) (2001) 107–127]. The resulting 2d numerical schemes are explicit and first order accurate. The solver is next parallelized by a domain decomposition technique. The specific application of the scheme to one- and two- layer shallow water systems has been implemented on a PCÕs cluster. An efficient data structure based on OOMPI (C++ object oriented extension of MPI) has been developed to optimize the data exchange among the processors. Some numerical tests are next presented to validate the solver and the performance of its parallel implementation. Finally the two-layer shallow water model is applied to the simulation of the steady exchange through the Strait of Gibraltar. Ó 2005 Elsevier B.V. All rights reserved. Keywords: Parallelization; Domain decomposition; Finite volume schemes; Conservation laws; Source terms; Nonconservative products; Shallow water systems; Two-layer problems; Geophysical flows; Strait of Gibraltar * Corresponding author. E-mail address: [email protected] (M.J. Castro). 0045-7825/$ - see front matter Ó 2005 Elsevier B.V. All rights reserved. doi:10.1016/j.cma.2005.07.007 M.J. Castro et al. / Comput. Methods Appl. Mech. Engrg. 195 (2006) 2788–2815 2789 1. Introduction This article deals with the discretization of 2d hyperbolic systems of conservation laws with source terms and nonconservative products by means of finite volume solvers on nonstructured grids. Problems of this nature arise in computational fluid dynamics. We are concerned in particular with the simulation of free- surface waves in shallow layers of homogeneous fluids or internal waves in stratified fluids composed by two shallow layers of immiscible liquids. The motion of a layer of homogeneous fluid is supposed here to be governed by the shallow water system, formulated under the form of a conservation law with source terms or balance law. In the stratified case, the flow is supposed to be governed by a system composed by two coupled shallow water systems. The global system can be formulated under the form of two coupled balance laws, the coupling terms having the form of nonconservative products. We are mainly interested in the application of these systems to geophysical flows: models based on shal- low water systems are useful for the simulation of rivers, channels, dambreak problems, etc. Stratified fluids are ubiquitous in nature: they appear in atmospheric flows, ocean currents, estuarine systems, etc. This is the situation, for instance, in the Strait of Gibraltar, where surface water from the Atlantic inflows over saltier westward-flowing Mediterranean water. It is well known that numerical methods that are suitable for solving conservation laws can fail in solving systems of balance laws, specially when approaching equilibria or near to equilibria solutions. Bermu´dez and Va´zquez-Cendo´n have shown in [29,3] that methods based on RoeÕs discretization of the flux terms and upwinding the source terms suitably solve these difficulties. For the particular case of the shallow water system, these authors introduced the condition called conservation property or C-property: a first order scheme is said to satisfy this condition if it solves, exactly or up to the second order, the steady state solu- tions corresponding to water at rest. The idea of constructing numerical schemes that preserve some equi- libria, which are called in general well-balanced schemes, has been extended in different ways: see [4] and the references therein, [18,13–15,20]. The presence of nonconservative products add more difficulties: in [12] it was shown that standard meth- ods for conservation laws can fail to convergence when they are applied to coupled systems of conservation or balance laws. In this latter article, some stable numerical schemes for 1d coupled systems of balance laws generalizing those introduced in [29,30,3] were presented. In the particular case of the two-layer shallow water system, the well-balance property required to the numerical schemes was the natural extension of the C-property: equilibria corresponding to water at rest had to be preserved. In [10] this work has been extended to the more complex system corresponding to 1d two-layer shallow water models in symmetric channels with irregular geometries. Besides the difficulties related to source terms and nonconservative products, some specific problems ap- pear in the particular case of shallow water systems like the numerical simulation of wet/dry fronts appear- ing when the thickness of the water layer vanishes. This difficulty also arises in two-layer fluids when the thickness of at least one of the layers vanishes. In real flows these situations appear in coastal simulations, upwelling events, etc. In [21,8,11] some modifications of the numerical scheme presented in [12] were pro- posed to handle with this difficulty. A specific difficulty related to the two-layer shallow water system is the appearance of Kelvin–Helmholtz instabilities due to the destabilizing effect of shear, which may overcome the stabilizing effects of stratifica- tion. When these instabilities develop, the P.D.E. system loses its hyperbolic character and the numerical scheme becomes unstable due to the appearance of complex eigenvalues. In this case, the assumptions of the two-layer model considered here are no longer valid: the interface becomes a rolled-up layer which can- not be represented as the graph of a function. Nevertheless, the model can be useful to simulate flows where Kelvin–Helmholtz instabilities are expected to appear only in a sporadic and local way. In this case a cor- rection has to be added to the numerical scheme to support instability episodes. In [18] a numerical 2790 M.J. Castro et al. / Comput. Methods Appl. Mech. Engrg. 195 (2006) 2788–2815 treatment was proposed consisting in adding locally an extra friction term when complex eigenvalues ap- pear. This extra amount of friction is intended to simulate the loss of momentum induced by the instabil- ities in real flows. The goal of this work is to obtain a first order explicit well-balanced numerical scheme for a general sys- tem of balance laws with nonconservative products generalizing the 1d solver presented in [12] for coupled systems of balance laws. Moreover, this numerical scheme has to be able to handle with wet/dry fronts and Kelvin–Helmholtz instabilities. The construction of such a numerical scheme will be based on the method of lines: following the procedure used in [18] for 2d balance laws, at every edge of the finite volume mesh a projected Riemann problem along the normal direction is considered which is discretized by means of a 1d numerical scheme. Discontinuous Galerkin methods offer an alternative way of building 2d schemes from 1d solvers (see [17] for a recent application to shallow water systems). The 2d solver presented here is original for problems with source terms and nonconservative products, but its application to the particular case of systems of balance laws recovers some known schemes: it coin- cides with the numerical scheme proposed in [29] for the 2d shallow water system and with a particular case of the family of numerical schemes for 2d systems of balance laws introduced in [18]. Due to the first order accuracy of the numerical scheme and to the fact that the space and/or the time scales are frequently very large, the realistic simulation of geophysical flows with the resulting solver can be highly CPU time consuming. Therefore, we study also its parallelization by means of a decomposition do- main technique. In order to minimize the computational time, data exchange among processors has to be minimized. The organization of the article is as follows: in Section 2, we present the general formulation of systems of balance laws with nonconservative products and source terms in 2d domains, and the particular cases corresponding to the one- and two-layer shallow water systems. Next, the construction of the numerical scheme is presented for the general case. Section 4 is devoted to the description of the domain decompo- sition technique used to parallelize the numerical scheme. The application to the one and the two-layer shal- low water systems is presented in Section 5, together with the study of the well-balance properties of the resulting schemes. Some notes about the practical implementation of the numerical scheme are presented in Section 6.We describe briefly the PCÕs cluster used, the finite volume chosen, the mesh partitioner considered and the data structure implemented in OOMPI (a C++ object oriented extension of MPI) in order to optimize the com- munications among processors. Finally, some numerical experiments are presented in order to validate the solver and to assess the per- formance of the one-layer and the two-layer shallow water models.
Recommended publications
  • Interannual Variability of the Mediterranean Outflow Observed in Espartel Sill, Western Strait of Gibraltar J
    JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114, C10018, doi:10.1029/2009JC005496, 2009 Interannual variability of the Mediterranean outflow observed in Espartel sill, western Strait of Gibraltar J. Garcı´a-Lafuente,1 J. Delgado,1 A. Sa´nchez Roma´n,1 J. Soto,1 L. Carracedo,1,2 and G. Dı´az del Rı´o3 Received 6 May 2009; revised 13 July 2009; accepted 24 July 2009; published 22 October 2009. [1] Four-year time series of observations in Espartel sill at the western part of the Strait of Gibraltar have been analyzed in order to investigate the variability of the Mediterranean outflow. It is assumed that the observed variability comes from the changing properties of the dense waters that are located at the maximum depth from where they can be uplifted in the upstream basin (Albora´n Sea, inside the Mediterranean Sea) and evacuated through the strait. From this perspective, the following three mechanisms are investigated: (1) the replenishment of the deep basin by newly formed Western Mediterranean Deep Water that, depending on its density, can either uplift old resident waters or lay above them leaving in any case a cold signature in the temperature series; (2) the presence/absence of the energetic anticyclonic gyres in the Albora´n Sea, particularly the western one, which can transfer momentum to the underlying Mediterranean vein and provide it with additional energy to ascend over the sills of the strait; and (3) the meteorologically enhanced flows that follow the rapid changes of atmospheric pressure over the western Mediterranean basin, which would be able to aspire deeper waters residing in the upstream basin.
    [Show full text]
  • UTMEA, Rome University of Malaga (Spain)
    Symposium on HPC And Data-Intensive Applications In Earth Sciences: Challenges and Opportunities 13 Nov. 2014 Gianmaria Sannino+, A. Carillo+, G. Pisacane+, C. Naranjo° ENEA - UTMEA, Rome University of Malaga (Spain) Energy and Environmental Modelling Unit [email protected] Italian Agency for New Technologies, Energy and Sustainable Economic Development NEMERTE – Numerical Experiment on the Mediterranean model response to Enhanced Resolution and TidE Project leader: Gianmaria Sannino, ENEA, ITALY Collaborators: Patrick Heimbach, Massachusetts Institute of Technology (MIT), UNITED STATES Gabriel Jorda Sanchez, Mediterranean Institute for Advanced Studies (IMEDEA), SPAIN Vincenzo Artale, ENEA, Italy Adriana Carillo, ENEA, Italy Emanuele Lombardi, ENEA, Italy Giovanna Pisacane, ENEA, Italy MariaVittoria Struglia, ENEA, Italy FERMI@CINECA NEMERTE Resource Awarded: 18 Million core hours on FERMI Architecture: 10 BGQ Frames Model: IBM-BG/Q Processor type: IBM PowerA2 @1.6 GHz Computing Cores: 163840 Computing Nodes: 10240 RAM: 1GByte / core (163 PByte total) Internal Network: 5D Torus Disk Space: 2PByte of scratch space Peak Performance: 2PFlop/s N. 7 in Top 500 rank (June 2012) National and PRACE Tier-0 calls The mission of PRACE (Partnership for Advanced Computing in Europe) is to enable high impact scientific discovery and engineering research and development across all disciplines to enhance European competitiveness for the benefit of society. PRACE seeks to realize this mission by offering world class computing and data management resources and services through a peer review process. OUTLINE HOW DID NEMERTE BEGIN? A fascinang story involving ocean models and HPC Mediterranean Thermohaline Circulation (MTHC) The Mediterranean Sea is a semi-enclosed basin displaying an active thermohaline circulation that is sustained by the atmospheric forcing and controlled by the narrow and shallow Strait of Gibraltar The atmospheric forcing drives the Mediterranean basin toward a negative budget of water and heat, and toward a positive budget of salt.
    [Show full text]
  • Oceanographic Processes and Morphosedimentary Products Along the Iberian Margins: a New Multidisciplinary Approach
    Marine Geology 378 (2016) 127–156 Contents lists available at ScienceDirect Marine Geology journal homepage: www.elsevier.com/locate/margo Oceanographic processes and morphosedimentary products along the Iberian margins: A new multidisciplinary approach F. Javier Hernández-Molina a,⁎, Anna Wåhlin b,MiguelBrunoc, Gemma Ercilla d, Estefanía Llave e, Nuno Serra f, Gabriel Rosón g,PerePuigd, Michele Rebesco h, David Van Rooij i,DavidRoquej, César González-Pola k, Francisco Sánchez l,MaríaGómezm,BenedictPreun, Tilmann Schwenk o, Till J.J. Hanebuth p,q, Ricardo F. Sánchez Leal r, Jesús García-Lafuente s, Rachel E. Brackenridge t,u,CarmenJuand,i, Dorrik A.V. Stow t,JoséMaríaSánchez-Gonzálezg a Dept. Earth Sciences, Royal Holloway Univ. London, Egham, Surrey TW20 0EX, UK b Department of Earth Sciences, University of Gothenburg, PO Box 460, SE 405 30 Göteborg, Sweden c CACYTMAR. Univ. Cádiz, Avda República Saharaui S/N, Puerto Real, 11510, Cádiz, Spain d CSIC, ICM, Paseo Marítimo de la Barceloneta, 37-49, 08003 Barcelona, Spain e Instituto Geológico y Minero de España Ríos Rosas, 23, 28003, Madrid, Spain f Institut für Meereskunde, Univ. Hamburg, Bundesstr. 53, 20146, Hamburg, Germany g Facultad de Ciencias do Mar, Univ. Vigo, 36200 Vigo, Spain h OGS, Istituto Nazionale di Oceanografia e di Geofisica Sperimentale, Borgo Grotta Gigante 42/C, 34010 Sgonico, TS, Italy i Renard Centre of Marine Geology, Dept. of Geology and Soil Science. Ghent University, Krijgslaan, 281 S8 B-9000 Gent, Belgium j CSIC, ICMAN, Campus Universitario Rio San Pedro s/n. Puerto Real, Cádiz, CP 11510, Spain k Instituto Español de Oceanografía, C.O. Gijón. c/ Príncipe de Asturias 70 Bis., CP 33212, Gijón, Spain l Instituto Español de Oceanografía, C.O.
    [Show full text]
  • The Mediterranean Outflow Water: Transformations and Pathways Into the Gulf of Cádiz
    The Mediterranean Outflow Water: Transformations and Pathways into the Gulf of Cádiz Marc Gasser i Rubinat Doctoral thesis as a compilation of publications September 2018 Advisor: Josep Lluís Pelegrí Llopart Thesis submitted to the Civil and Environmental Engineering Department, Marine Sciences Program, Universitat Politècnica de Catalunya in partial fulfillment of the requirements for the degree of Doctor of Philosophy Marc Gasser i Rubinat When anxious, uneasy and bad thoughts come, I go to the sea, and the sea drowns them out with its great wide sounds, cleanses me with its noise, and imposes a rhythm upon everything in me that is bewildered and confused". Rainer Maria Rilke III The Mediterranean Outflow Water IV Marc Gasser i Rubinat Abstract The Mediterranean Outflow Water (MOW) is La sortida d'aigua mediterrània (MOW) és un a dense ( r>1028.5 kg/m 3), saline (38.5 g/kg) corrent oceànic dens ( r>1028.5 kg/m3) i salí ocean stream originated in the evaporative (38.5 lg/m3) originat en la conca evaporativa Mediterranean basin flowing westward past de la Mar Mediterrània que flueix passant la Espartel Sill as a fast (>1 m/s) and often baixa d'Espartel en forma d'un corrent de unstable (as indicated by its gradient gravetat molt ràpid (>1 m/s) i sovint inestable Richardson number) gravity current. During (tal i com indica el número de gradient de its descense into the Gulf of Cadiz, the MOW Richardson). Durant el seu descens al Golf entrains the overlying North Atlantic Central de Cadis, la MOW incorpora les aigües Water (NACW), until the density difference atlàntiques (NACW) suprajacents fins que la between both water masses vanishes, and diferència de densitat entre ambdues masses reaches its equilibrium depth.
    [Show full text]
  • The Western Alboran Gyre Helps Ventilate the Western Mediterranean Deep Water Through Gibraltar
    Deep-Sea Research I 63 (2012) 157–163 Contents lists available at SciVerse ScienceDirect Deep-Sea Research I journal homepage: www.elsevier.com/locate/dsri Note The Western Alboran Gyre helps ventilate the Western Mediterranean Deep Water through Gibraltar Cristina Naranjo n, Jesu´ s Garcı´a Lafuente, Jose C. Sa´nchez Garrido, Antonio Sa´nchez Roma´n, Javier Delgado Cabello Physical Oceanography Group, ETSI Telecomunicacio´n, University of Ma´laga, Campus Teatinos, 29071 Ma´laga, Spain article info abstract Article history: Variable properties of the Mediterranean outflow and the variability of the Western Alboran Gyre are Received 27 April 2011 analyzed by means of 5-year long time series of near bottom potential temperature at Espartel sill in Received in revised form the Strait of Gibraltar and altimetry data in the Alboran Sea. Geostrophic velocity at the southern edge 19 October 2011 of the gyre and potential temperature at Espartel sill are significantly correlated (correlation coefficient Accepted 24 October 2011 0.67), suggesting that the intensification of the Alboran Gyre favors the ventilation of Western Available online 29 October 2011 Mediterranean Deep Water. The analysis of historical temperature profiles shows that Western Keywords: Mediterranean Deep Water in the Alboran Sea can be suctioned from a layer between 500 and Western Mediterranean Deep Water 700 m depth for typical changes of the gyre intensification. Deep water ventilation & 2011 Elsevier Ltd. All rights reserved. Strait of Gibraltar Alboran Gyre 1. Introduction approach of the Strait and to bring them west of CS into the Atlantic Ocean. Actually, the ventilation of WMDW would comprise two The question of how the Western Mediterranean Deep Water separated, though linked, set of processes.
    [Show full text]
  • The Mediterranean Outflow in the Strait of Gibraltar and Its Connection
    Ocean Sci., 13, 195–207, 2017 www.ocean-sci.net/13/195/2017/ doi:10.5194/os-13-195-2017 © Author(s) 2017. CC Attribution 3.0 License. The Mediterranean outflow in the Strait of Gibraltar and its connection with upstream conditions in the Alborán Sea Jesús García-Lafuente, Cristina Naranjo, Simone Sammartino, José C. Sánchez-Garrido, and Javier Delgado Physical Oceanography Group, Department of Applied Physics 2, University of Málaga, Málaga, Spain Correspondence to: Jesús García-Lafuente ([email protected]) Received: 23 November 2016 – Discussion started: 13 December 2016 Revised: 24 February 2017 – Accepted: 28 February 2017 – Published: 24 March 2017 Abstract. The present study addresses the hypothesis that the mediate and the Tyrrhenian Dense waters, both of interme- Western Alborán Gyre in the Alborán Sea (the westernmost diate nature. The Winter Intermediate Water is formed along Mediterranean basin adjacent to the Strait of Gibraltar) in- the continental shelf of the Liguro-Provençal sub-basin and fluences the composition of the outflow through the Strait Catalan Sea (Conan and Millot, 1995; Vargas-Yáñez et al., of Gibraltar. The process invoked is that strong and well- 2012), exhibits marked interannual fluctuations that include developed gyres help to evacuate the Western Mediterranean years of no formation (Pinot et al., 2002; Monserrat et al., Deep Water from the Alborán basin, thus increasing its pres- 2008), and is characterized by an absolute minimum of po- ence in the outflow, whereas weak gyres facilitate the outflow tential temperature. Its volume transport is much less than the of Levantine and other intermediate waters. To this aim, in Levantine Intermediate Water and it flows embedded inside situ observations collected at the Camarinal (the main) and this water mass at relatively shallow depths.
    [Show full text]
  • Exchange Flow Through the Strait of Gibraltar As Simulated by a Σ-Coordinate Hydrostatic Model and a Z-Coordinate Nonhydrostatic Model
    3 Exchange Flow through the Strait of Gibraltar as Simulated by a σ-Coordinate Hydrostatic Model and a z-Coordinate Nonhydrostatic Model Gianmaria Sannino1, J. C. Sánchez Garrido2, L. Liberti3, and L. Pratt4 3.1. INTRODUCTION 2007; Sannino et al., 2009a]. The simultaneous presence in the SoG of at least two cross sections in which the The Mediterranean Sea is a semi-enclosed basin dis- exchange is controlled drives the strait dynamics toward playing an active thermohaline circulation (MTHC) that the so-called maximal regime [Bryden and Stommel, 1984; is sustained by the atmospheric forcing and controlled by Armi and Farmer, 1988]. If the exchange is subject to only the narrow and shallow Strait of Gibraltar (hereinafter one hydraulic control, the regime is called submaximal. SoG). The atmospheric forcing drives the Mediterranean The two regimes have different implications for property basin toward a negative budget of water and heat. Over fluxes, response time, and other physical characteristics the basin, evaporation exceeds the sum of precipitation of the coupled circulation in the SoG and Mediterranean and river discharge, while a net heat flux is transferred to Sea. The maximal regime can be expected to have larger the overlying atmosphere through the sea surface. These heat, salt, and mass fluxes and to respond more slowly to fluxes are balanced by the exchange flow that takes place changes in stratification and thermohaline forcing within in Gibraltar. Within the SoG, the MTHC takes the form the Mediterranean Sea and the North Atlantic Ocean of a two-way exchange: an upper layer of fresh and [Sannino et al., 2009a].
    [Show full text]
  • Vertical Structure of Tidal Currents Over Espartel and Camarinal Sills, Strait of Gibraltar ⁎ A
    Available online at www.sciencedirect.com Journal of Marine Systems 74 (2008) 120–133 www.elsevier.com/locate/jmarsys Vertical structure of tidal currents over Espartel and Camarinal sills, Strait of Gibraltar ⁎ A. Sánchez-Román a, , F. Criado-Aldeanueva a, J. García-Lafuente a, J.C. Sánchez a,b a Departamento de Física Aplicada II, Universidad de Málaga, Málaga, Spain b Grupo de Puertos y Costas, Universidad de Granada, Granada, Spain Received 23 July 2007; received in revised form 7 November 2007; accepted 19 November 2007 Available online 5 February 2008 Abstract ADCP velocity data collected in the two main sills of the Strait of Gibraltar (Camarinal and Espartel sills) have been used for analysing the vertical structure of main tidal constituents (M2,S2,O1 and K1) currents in this area. Two different periods (winter and summer) were considered in correspondence to seasonal variations in density profiles. Amplitudes and phases of the various tidal constituents have been compared for both periods and locations. Barotropic and baroclinic parts of the tidal currents have been extracted using the dynamical mode decomposition technique and the relative importance of each mode has been established in terms on the energy associated. In Espartel sill, the barotropic mode is more energetic in wintertime for all constituents except for K1. Baroclinic modes have smaller contribution to total energy. Second and third baroclinic modes tend to be more energetic than the first one, especially during winter, this indicating a relatively complex vertical structure of the tidal currents. Over Camarinal sill, barotropic mode accounts for more than 90% of total energy in all the tidal constituents, the highest value (97%) observed for M2.
    [Show full text]
  • Hydraulic Criticality of the Exchange Flow Through the Strait of Gibraltar
    NOVEMBER 2009 S A N N I N O E T A L . 2779 Hydraulic Criticality of the Exchange Flow through the Strait of Gibraltar GIANMARIA SANNINO Ocean Modelling Unit, ENEA, Rome, Italy LAWRENCE PRATT Woods Hole Oceanographic Institution, Woods Hole, Massachusetts ADRIANA CARILLO Ocean Modelling Unit, ENEA, Rome, Italy (Manuscript received 17 June 2008, in final form 9 January 2009) ABSTRACT The hydraulic state of the exchange circulation through the Strait of Gibraltar is defined using a recently developed critical condition that accounts for cross-channel variations in layer thickness and velocity, applied to the output of a high-resolution three-dimensional numerical model simulating the tidal exchange. The numerical model uses a coastal-following curvilinear orthogonal grid, which includes, in addition to the Strait of Gibraltar, the Gulf of Cadiz and the Alboran Sea. The model is forced at the open boundaries through the specification of the surface tidal elevation that is characterized by the two principal semidiurnal and two diurnal harmonics: M2, S2, O1, and K1. The simulation covers an entire tropical month. The hydraulic analysis is carried out approximating the continuous vertical stratification first as a two-layer system and then as a three-layer system. In the latter, the transition zone, generated by entrainment and mixing between the Atlantic and Mediterranean flows, is considered as an active layer in the hydraulic model. As result of these vertical approximations, two different hydraulic states have been found; however, the simulated behavior of the flow only supports the hydraulic state predicted by the three-layer case. Thus, analyzing the results obtained by means of the three-layer hydraulic model, the authors have found that the flow in the strait reaches maximal exchange about 76% of the tropical monthlong period.
    [Show full text]
  • Science Journals
    CORE Metadata, citation and similar papers at core.ac.uk Provided by Sapientia SCIENCE ADVANCES | RESEARCH ARTICLE OCEANOGRAPHY Copyright © 2017 The Authors, some The Mediterranean Overflow in the Gulf of Cadiz: rights reserved; exclusive licensee A rugged journey American Association for the Advancement 1 1 1 of Science. No claim to Ricardo F. Sánchez-Leal, * María Jesús Bellanco, Luis Miguel Fernández-Salas, original U.S. Government 2 3 4 Jesús García-Lafuente, Marc Gasser-Rubinat, César González-Pola, Works. Distributed 5 3 6 7 8 Francisco J. Hernández-Molina, Josep L. Pelegrí, Alvaro Peliz, Paulo Relvas, David Roque, under a Creative Manuel Ruiz-Villarreal,9 Simone Sammartino,2 José Carlos Sánchez-Garrido2 Commons Attribution NonCommercial The pathways and transformations of dense water overflows, which depend on small-scale interactions between flow License 4.0 (CC BY-NC). dynamics and erosional-depositional processes, are a central piece in the ocean’s large-scale circulation. A novel, high- resolution current and hydrographic data set highlights the intricate pathway travelled by the saline Mediterranean Overflow as it enters the Atlantic. Interaction with the topography constraints its spreading. Over the initial 200 km west of the Gibraltar gateway, distinct channels separate the initial gravity current into several plunging branches depth-sorted by density. Shallow branches follow the upper slope and eventually detach as buoyant plumes. Deeper branches occupy mid slope channels and coalesce upon reaching a diapiric ridge. A still deeper branch, guided by a lower channel wall marked by transverse furrows, experiences small-scale overflows which travel downslope to settle Downloaded from at mid-depths.
    [Show full text]
  • Vertical Density Gradient in the Eastern North Atlantic During the Last 30,000 Years
    Clim Dyn DOI 10.1007/s00382-011-1148-4 Vertical density gradient in the eastern North Atlantic during the last 30,000 years M. Rogerson • G. R. Bigg • E. J. Rohling • J. Ramirez Received: 30 March 2011 / Accepted: 13 July 2011 Ó Springer-Verlag 2011 Abstract Past changes in the density and momentum Regardless, our novel reconstruction of Atlantic density structure of oceanic circulation are an important aspect of structure can be used as the basis for a dynamical measure changes in the Atlantic Meridional Overturning Circulation for validation of model-based AMOC reconstructions. In and consequently climate. However, very little is known addition, our general approach is transferrable to other about past changes in the vertical density structure of the marginal sea outflow plumes, to provide estimates of ocean, even very extensively studied systems such as the oceanic vertical density gradients in other locations. North Atlantic. Here we exploit the physical controls on the settling depth of the dense Mediterranean water plume Keywords Atlantic Ocean Á Meridional overturning Á derived from the Strait of Gibraltar to obtain the first Overflow physics Á Mediterranean outflow Á robust, observations-based, probabilistic reconstruction of Iberian margin Á Paleoceanography the vertical density gradient in the eastern North Atlantic during the last 30,000 years. We find that this gradient was weakened by more than 50%, relative to the present, during 1 Introduction the last Glacial Maximum, and that changes in general are associated with reductions in AMOC intensity. However, Progress in the reconstruction of past Atlantic Meridional we find only a small change during Heinrich Event 1 rela- Overturning Circulation (AMOC) changes (McManus et al.
    [Show full text]
  • Ten Years of Marine Current Measurements in Espartel Sill, 10.1002/2014JC010674 Strait of Gibraltar
    PUBLICATIONS Journal of Geophysical Research: Oceans RESEARCH ARTICLE Ten years of marine current measurements in Espartel Sill, 10.1002/2014JC010674 Strait of Gibraltar Key Points: S. Sammartino1, J. Garcıa Lafuente1, C. Naranjo1,J.C.Sanchez Garrido1,R.Sanchez Leal2, and A 10 year series of ADCP A. Sanchez Roman 3 measurements has been processed The series has been submitted to a 1Physical Oceanography Group, University of Malaga, Malaga, Spain, 2Spanish Institute of Oceanography, Cadiz careful quality control 3 The mean value of the outflow is Oceanography Center, Cadiz, Spain, Mediterranean Institute for Advanced Studies, Esporles, Spain 20.85 6 0.03 Sv Abstract More than 10 year of Acoustic Doppler Current Profiler observations collected at the western- Supporting Information: Supporting Information S1 most sill (Espartel sill) of the Strait of Gibraltar by a monitoring station have been carefully processed to pro- Figure S1 vide the most updated estimation of the Mediterranean outflow. A comprehensive quality control of the Movie S1 factors affecting the uncertainty of the measurements has been carried out and great care has been paid to Movie S2 infer the current at the bottom layer, where direct observations are lacking. The mean outflow in the south- ern channel of the sill section has been estimated as 20.82 Sv (1 Sv 5 1 3 106 m3 s21), with an average con- Correspondence to: S. Sammartino, tribution of the eddy fluxes of 20.04 Sv. This figure is an overestimation, as the mooring measurements, [email protected] assumed valid for the whole section, ignore the lateral friction. On the other hand, it only gives the flow through the southern channel and disregards the fraction flowing through shallower northern part.
    [Show full text]