DOCSLIB.ORG

Physical Oceanography - UNAM, Mexico Lecture 2: the Equations of Ocean Circulation and Ocean Modelling

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Physical Oceanography - UNAM, Mexico Lecture 2: the Equations of Ocean Circulation and Ocean Modelling Physical Oceanography - UNAM, Mexico Lecture 2: The Equations of Ocean Circulation and Ocean Modelling Robin Waldman October 15th 2018 A first taste... ¶uh 1 ¶ ¶uh + (u:r)uh + f k × uh = − rhP + rh:(khurh)uh + (kzu ) ¶t r0 ¶z ¶z Z 0 0 P(z) = r0gh + g rdz z Z z 0 w(z) = − rh:uhdz −H ¶h Z h = −rh: uhdz + P + R − E ¶t −H ¶q ¶ ¶q 1 _ + (u:r)q = rh:(khT rh)q + (kzT ) + Θ ¶t ¶z ¶z rcw ¶S ¶ ¶S + (u:r)S = r :(k r )S + (k ) + S_ ¶t h hS h ¶z zS ¶z r = r(q;S;P0(z)) All the physics of an ocean circulation model is here ! Outline The Equations of Ocean Circulation Ocean modelling Outline The Equations of Ocean Circulation Ocean modelling Conservation of mass : continuity In the following : control volume of zonal, meridional and vertical sizes dx, dy and dz and density r on a fixed Cartesian coordinate system (i;j;k) attached to the ground. d(rdxdydz) = 0 dt dr ddx ddy ddz = dxdydz + r(dydz + dxdz + dxdy ) dt dt dt dt dr = dxdydz + r(dydzdu + dxdzdv + dxdydw) dt Hence dividing by dxdydz : dr ¶u ¶v ¶w + r( + + ) = 0 dt ¶x ¶y ¶z dr () + rr:v = 0 dt ¶ ¶ ¶ with r = ( ¶x ; ¶y ; ¶z ) the space derivative operator. Conservation of mass : continuity Mass conservation : Hence dividing by dxdydz : dr ¶u ¶v ¶w + r( + + ) = 0 dt ¶x ¶y ¶z dr () + rr:v = 0 dt ¶ ¶ ¶ with r = ( ¶x ; ¶y ; ¶z ) the space derivative operator. Conservation of mass : continuity Mass conservation : d(rdxdydz) = 0 dt dr ddx ddy ddz = dxdydz + r(dydz + dxdz + dxdy ) dt dt dt dt dr = dxdydz + r(dydzdu + dxdzdv + dxdydw) dt Conservation of mass : continuity Mass conservation : d(rdxdydz) = 0 dt dr ddx ddy ddz = dxdydz + r(dydz + dxdz + dxdy ) dt dt dt dt dr = dxdydz + r(dydzdu + dxdzdv + dxdydw) dt Hence dividing by dxdydz : dr ¶u ¶v ¶w + r( + + ) = 0 dt ¶x ¶y ¶z dr () + rr:v = 0 dt ¶ ¶ ¶ with r = ( ¶x ; ¶y ; ¶z ) the space derivative operator. dr0 0 dt + (r0 + r )r:v ' r0r:v = 0 () r:v = 0 I To a very good approximation, oceanic currents are non-divergent. Conservation of mass : continuity In the ocean, Boussinesq approximation : relative density variations are small : 0 0 r(x;y;z;t) = r0 + r (x;y;z;t) and r << r0 The continuity equation becomes : Conservation of mass : continuity In the ocean, Boussinesq approximation : relative density variations are small : 0 0 r(x;y;z;t) = r0 + r (x;y;z;t) and r << r0 The continuity equation becomes : dr0 0 dt + (r0 + r )r:v ' r0r:v = 0 () r:v = 0 I To a very good approximation, oceanic currents are non-divergent. Conservation of momentum Newton’s 2nd law : the Lagrangian (material) evolution of momentum is determined by the sum of external (gravity) and body (pressure and friction) forces. ¶u ¶u ¶u ¶u du = dt + dx + dy + dz ¶t ¶x ¶y ¶z du ¶u dt ¶u dx ¶u dy ¶u dz () = + + + dt ¶t dt ¶x dt ¶y dt ¶z dt ¶u ¶u ¶u ¶u = + u + v + w ¶t ¶x ¶y ¶z ¶u = + (u:r)u ¶t (u:r)u is the (nonlinear) advection term. Conservation of momentum ¶u The acceleration in the Eulerian (fixed) framework ¶t is far easier to du observe and model than the Lagrangian one dt . Relation between both : Conservation of momentum ¶u The acceleration in the Eulerian (fixed) framework ¶t is far easier to du observe and model than the Lagrangian one dt . Relation between both : ¶u ¶u ¶u ¶u du = dt + dx + dy + dz ¶t ¶x ¶y ¶z du ¶u dt ¶u dx ¶u dy ¶u dz () = + + + dt ¶t dt ¶x dt ¶y dt ¶z dt ¶u ¶u ¶u ¶u = + u + v + w ¶t ¶x ¶y ¶z ¶u = + (u:r)u ¶t (u:r)u is the (nonlinear) advection term. Conservation of zonal momentum over a control volume ∂u ν (z)δ x δ y - w(z) δxδy u(z) ∂ z ∂ u ν ( y)δ x δ z - v(y) δxδz u(y) ∂ y P(x-δx)/ρ δyδz P(x)/ρ δyδz 0 0 u(x-δx) δyδz u(x-δx) -f v δxδyδz ∂u ∂u - u(x) δyδz u(x) −ν (x−δ x)δ y δ z −ν ( y−δ y)δ x δ z ∂ u ∂ x ∂ y ν (x)δ y δ z δz ∂ x v(y-δy) δxδz u(y-δy) δy w(z-δz) δxδy u(z-δz) ∂ u −ν (z−δ z)δ x δ y k j ∂ z i Figure 1 – Conservation of zonal momentum over a control volume : advection (black), viscous forces (red), pressure forces (blue) and Coriolis acceleration (purple). Conservation of momentum Advection of zonal momentum in all faces of the control volume : Conservation of zonal momentum over a control volume ∂u ν (z)δ x δ y - w(z) δxδy u(z) ∂ z ∂ u ν ( y)δ x δ z - v(y) δxδz u(y) ∂ y P(x-δx)/ρ δyδz P(x)/ρ δyδz 0 0 u(x-δx) δyδz u(x-δx) -f v δxδyδz ∂u ∂u - u(x) δyδz u(x) −ν (x−δ x)δ y δ z −ν ( y−δ y)δ x δ z ∂ u ∂ x ∂ y ν (x)δ y δ z δz ∂ x v(y-δy) δxδz u(y-δy) δy w(z-δz) δxδy u(z-δz) ∂ u −ν (z−δ z)δ x δ y k j ∂ z i Conservation of momentum Advection of zonal momentum in all faces of the control volume : Figure 1 – Conservation of zonal momentum over a control volume : advection (black), viscous forces (red), pressure forces (blue) and Coriolis acceleration (purple). Conservation of zonal momentum over a control volume ∂u ν (z)δ x δ y - w(z) δxδy u(z) ∂ z ∂ u ν ( y)δ x δ z - v(y) δxδz u(y) ∂ y P(x-δx)/ρ δyδz P(x)/ρ δyδz 0 0 u(x-δx) δyδz u(x-δx) -f v δxδyδz ∂u ∂u - u(x) δyδz u(x) −ν (x−δ x)δ y δ z −ν ( y−δ y)δ x δ z ∂ u ∂ x ∂ y ν (x)δ y δ z δz ∂ x v(y-δy) δxδz u(y-δy) δy w(z-δz) δxδy u(z-δz) ∂ u −ν (z−δ z)δ x δ y k j ∂ z i Figure 2 – Conservation of zonal momentum over a control volume : advection (black), viscous forces (red), pressure forces (blue) and Coriolis acceleration (purple). Conservation of momentum Zonal pressure force on the western and eastern faces of the control volume : Conservation of zonal momentum over a control volume ∂u ν (z)δ x δ y - w(z) δxδy u(z) ∂ z ∂ u ν ( y)δ x δ z - v(y) δxδz u(y) ∂ y P(x-δx)/ρ δyδz P(x)/ρ δyδz 0 0 u(x-δx) δyδz u(x-δx) -f v δxδyδz ∂u ∂u - u(x) δyδz u(x) −ν (x−δ x)δ y δ z −ν ( y−δ y)δ x δ z ∂ u ∂ x ∂ y ν (x)δ y δ z δz ∂ x v(y-δy) δxδz u(y-δy) δy w(z-δz) δxδy u(z-δz) ∂ u −ν (z−δ z)δ x δ y k j ∂ z i Conservation of momentum Zonal pressure force on the western and eastern faces of the control volume : Figure 2 – Conservation of zonal momentum over a control volume : advection (black), viscous forces (red), pressure forces (blue) and Coriolis acceleration (purple). Conservation of momentum Zonal pressure force on the western and eastern faces of the control volume : FPx = (P(x − dx) − P(x))dydz hence volumic zonal pressure force : F ¶P Px = − dxdydz ¶x Generalizing to other spatial dimensions : F P = −rP dxdydz Conservation of momentum Friction Ft : I Newton’s law of viscosity for a Boussinesq Newtonian fluid : ¶uj ¶ui tij = +nr( + ) ¶xi ¶xj with tij the viscous stress exerted over the coordinate i on velocity component j, xi ;xj = (x;y;z), ui ;uj = (u;v;w) and n = 8:9 × 10−7m2=s the water kinematic viscosity. I The second term vanishes in a Boussinesq fluid : hence friction behaves just like molecular or heat diffusion ! Conservation of zonal momentum over a control volume ∂u ν (z)δ x δ y - w(z) δxδy u(z) ∂ z ∂ u ν ( y)δ x δ z - v(y) δxδz u(y) ∂ y P(x-δx)/ρ δyδz P(x)/ρ δyδz 0 0 u(x-δx) δyδz u(x-δx) -f v δxδyδz ∂u ∂u - u(x) δyδz u(x) −ν (x−δ x)δ y δ z −ν ( y−δ y)δ x δ z ∂ u ∂ x ∂ y ν (x)δ y δ z δz ∂ x v(y-δy) δxδz u(y-δy) δy w(z-δz) δxδy u(z-δz) ∂ u −ν (z−δ z)δ x δ y k j ∂ z i Figure 3 – Conservation of zonal momentum over a control volume : advection (black), viscous forces (red), pressure forces (blue) and Coriolis acceleration (purple). Conservation of momentum Zonal friction force on all faces of the control volume : Conservation of zonal momentum over a control volume ∂u ν (z)δ x δ y - w(z) δxδy u(z) ∂ z ∂ u ν ( y)δ x δ z - v(y) δxδz u(y) ∂ y P(x-δx)/ρ δyδz P(x)/ρ δyδz 0 0 u(x-δx) δyδz u(x-δx) -f v δxδyδz ∂u ∂u - u(x) δyδz u(x) −ν (x−δ x)δ y δ z −ν ( y−δ y)δ x δ z ∂ u ∂ x ∂ y ν (x)δ y δ z δz ∂ x v(y-δy) δxδz u(y-δy) δy w(z-δz) δxδy u(z-δz) ∂ u −ν (z−δ z)δ x δ y k j ∂ z i Conservation of momentum Zonal friction force on all faces of the control volume : Figure 3 – Conservation of zonal momentum over a control volume : advection (black), viscous forces (red), pressure forces (blue) and Coriolis acceleration (purple).
Load more

Recommended publications
  • Fronts in the World Ocean's Large Marine Ecosystems. ICES CM 2007
    - 1 - This paper can be freely cited without prior reference to the authors International Council ICES CM 2007/D:21 for the Exploration Theme Session D: Comparative Marine Ecosystem of the Sea (ICES) Structure and Function: Descriptors and Characteristics Fronts in the World Ocean’s Large Marine Ecosystems Igor M. Belkin and Peter C. Cornillon Abstract. Oceanic fronts shape marine ecosystems; therefore front mapping and characterization is one of the most important aspects of physical oceanography. Here we report on the first effort to map and describe all major fronts in the World Ocean’s Large Marine Ecosystems (LMEs). Apart from a geographical review, these fronts are classified according to their origin and physical mechanisms that maintain them. This first-ever zero-order pattern of the LME fronts is based on a unique global frontal data base assembled at the University of Rhode Island. Thermal fronts were automatically derived from 12 years (1985-1996) of twice-daily satellite 9-km resolution global AVHRR SST fields with the Cayula-Cornillon front detection algorithm. These frontal maps serve as guidance in using hydrographic data to explore subsurface thermohaline fronts, whose surface thermal signatures have been mapped from space. Our most recent study of chlorophyll fronts in the Northwest Atlantic from high-resolution 1-km data (Belkin and O’Reilly, 2007) revealed a close spatial association between chlorophyll fronts and SST fronts, suggesting causative links between these two types of fronts. Keywords: Fronts; Large Marine Ecosystems; World Ocean; sea surface temperature. Igor M. Belkin: Graduate School of Oceanography, University of Rhode Island, 215 South Ferry Road, Narragansett, Rhode Island 02882, USA [tel.: +1 401 874 6533, fax: +1 874 6728, email: [email protected]].
    [Show full text]
  • A Destabilizing Thermohaline Circulation–Atmosphere–Sea Ice
    642 JOURNAL OF CLIMATE VOLUME 12 NOTES AND CORRESPONDENCE A Destabilizing Thermohaline Circulation±Atmosphere±Sea Ice Feedback STEVEN R. JAYNE MIT±WHOI Joint Program in Oceanography, Woods Hole Oceanographic Institution, Woods Hole, Massachusetts JOCHEM MAROTZKE Center for Global Change Science, Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, Massachusetts 18 November 1996 and 9 March 1998 ABSTRACT Some of the interactions and feedbacks between the atmosphere, thermohaline circulation, and sea ice are illustrated using a simple process model. A simpli®ed version of the annual-mean coupled ocean±atmosphere box model of Nakamura, Stone, and Marotzke is modi®ed to include a parameterization of sea ice. The model includes the thermodynamic effects of sea ice and allows for variable coverage. It is found that the addition of sea ice introduces feedbacks that have a destabilizing in¯uence on the thermohaline circulation: Sea ice insulates the ocean from the atmosphere, creating colder air temperatures at high latitudes, which cause larger atmospheric eddy heat and moisture transports and weaker oceanic heat transports. These in turn lead to thicker ice coverage and hence establish a positive feedback. The results indicate that generally in colder climates, the presence of sea ice may lead to a signi®cant destabilization of the thermohaline circulation. Brine rejection by sea ice plays no important role in this model's dynamics. The net destabilizing effect of sea ice in this model is the result of two positive feedbacks and one negative feedback and is shown to be model dependent. To date, the destabilizing feedback between atmospheric and oceanic heat ¯uxes, mediated by sea ice, has largely been neglected in conceptual studies of thermohaline circulation stability, but it warrants further investigation in more realistic models.
    [Show full text]
  • Glacial Ocean Circulation and Stratification Explained by Reduced
    Glacial ocean circulation and stratification explained by reduced atmospheric temperature Malte F. Jansena,1 aDepartment of the Geophysical Sciences, The University of Chicago, Chicago, IL 60637 Edited by Mark H. Thiemens, University of California, San Diego, La Jolla, CA, and approved November 7, 2016 (received for review June 27, 2016) Earth’s climate has undergone dramatic shifts between glacial and We test the connection between atmospheric temperature interglacial time periods, with high-latitude temperature changes and ocean circulation and stratification changes, using idealized on the order of 5–10 ◦C. These climatic shifts have been asso- numerical simulations, which allow us to isolate the proposed ciated with major rearrangements in the deep ocean circulation mechanism. We use a coupled ocean–sea-ice model, with atmo- and stratification, which have likely played an important role in spheric temperature, winds, and evaporation–precipitation pre- the observed atmospheric carbon dioxide swings by affecting the scribed as boundary conditions (Materials and Methods). The partitioning of carbon between the atmosphere and the ocean. model uses an idealized continental configuration resembling the The mechanisms by which the deep ocean circulation changed, Atlantic and Southern Oceans, where the most elemental circu- however, are still unclear and represent a major challenge to our lation changes have been inferred (3, 5, 6, 12). understanding of glacial climates. This study shows that vari- ous inferred changes in the deep ocean circulation and stratifica- Results tion between glacial and interglacial climates can be interpreted We first focus on the model’s ability to reproduce key features as a direct consequence of atmospheric temperature differences.
    [Show full text]
  • Ice Production in Ross Ice Shelf Polynyas During 2017–2018 from Sentinel–1 SAR Images
    remote sensing Article Ice Production in Ross Ice Shelf Polynyas during 2017–2018 from Sentinel–1 SAR Images Liyun Dai 1,2, Hongjie Xie 2,3,* , Stephen F. Ackley 2,3 and Alberto M. Mestas-Nuñez 2,3 1 Key Laboratory of Remote Sensing of Gansu Province, Heihe Remote Sensing Experimental Research Station, Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences, Lanzhou 730000, China; [email protected] 2 Laboratory for Remote Sensing and Geoinformatics, Department of Geological Sciences, University of Texas at San Antonio, San Antonio, TX 78249, USA; [email protected] (S.F.A.); [email protected] (A.M.M.-N.) 3 Center for Advanced Measurements in Extreme Environments, University of Texas at San Antonio, San Antonio, TX 78249, USA * Correspondence: [email protected]; Tel.: +1-210-4585445 Received: 21 April 2020; Accepted: 5 May 2020; Published: 7 May 2020 Abstract: High sea ice production (SIP) generates high-salinity water, thus, influencing the global thermohaline circulation. Estimation from passive microwave data and heat flux models have indicated that the Ross Ice Shelf polynya (RISP) may be the highest SIP region in the Southern Oceans. However, the coarse spatial resolution of passive microwave data limited the accuracy of these estimates. The Sentinel-1 Synthetic Aperture Radar dataset with high spatial and temporal resolution provides an unprecedented opportunity to more accurately distinguish both polynya area/extent and occurrence. In this study, the SIPs of RISP and McMurdo Sound polynya (MSP) from 1 March–30 November 2017 and 2018 are calculated based on Sentinel-1 SAR data (for area/extent) and AMSR2 data (for ice thickness).
    [Show full text]
  • Antarctic Sea Ice Control on Ocean Circulation in Present and Glacial Climates
    Antarctic sea ice control on ocean circulation in present and glacial climates Raffaele Ferraria,1, Malte F. Jansenb, Jess F. Adkinsc, Andrea Burkec, Andrew L. Stewartc, and Andrew F. Thompsonc aDepartment of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139; bAtmospheric and Oceanic Sciences Program, Geophysical Fluid Dynamics Laboratory, Princeton, NJ 08544; and cDivision of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA 91125 Edited* by Edward A. Boyle, Massachusetts Institute of Technology, Cambridge, MA, and approved April 16, 2014 (received for review December 31, 2013) In the modern climate, the ocean below 2 km is mainly filled by waters possibly associated with an equatorward shift of the Southern sinking into the abyss around Antarctica and in the North Atlantic. Hemisphere westerlies (11–13), (ii) an increase in abyssal stratifi- Paleoproxies indicate that waters of North Atlantic origin were instead cation acting as a lid to deep carbon (14), (iii)anexpansionofseaice absent below 2 km at the Last Glacial Maximum, resulting in an that reduced the CO2 outgassing over the Southern Ocean (15), and expansion of the volume occupied by Antarctic origin waters. In this (iv) a reduction in the mixing between waters of Antarctic and Arctic study we show that this rearrangement of deep water masses is origin, which is a major leak of abyssal carbon in the modern climate dynamically linked to the expansion of summer sea ice around (16). Current understanding is that some combination of all of these Antarctica. A simple theory further suggests that these deep waters feedbacks, together with a reorganization of the biological and only came to the surface under sea ice, which insulated them from carbonate pumps, is required to explain the observed glacial drop in atmospheric forcing, and were weakly mixed with overlying waters, atmospheric CO2 (17).
    [Show full text]
  • DESALINATION: Balancing the Socioeconomic Benefits and Environmental Costs
    DESALINATION: Balancing the Socioeconomic Benefits and Environmental Costs www.research.natixis.com https://gsh.cib.natixis.com executive summary Chapter 1 Making sense of desalination: technological, financial and economic aspects of desalination assets Chapter 2 Sustainability assessment of desalination assets: recognizing the socioeconomic benefits and mitigating environmental costs of desalination Chapter 3 Desalination sustainability performance scorecard acknowledgements appendix biblioghraphy TABLE OF CONTENTS OF TABLE 1. Making sense of desalination: technological, financial and economic aspects of desalination assets 1. DESALINATION TECHNOLOGIES 2. FINANCIAL AND ECONOMIC ASPECTS OF DESALINATION ASSETS 1.1. AN OVERVIEW OF DESALINATION TECHNOLOGIES 2.1.THE DEVELOPMENT AND FINANCING OF DESALINATION ASSETS Thermal desalination: Multistage Flash Distillation and Multieffect Distillation Building and operating desalination assets: complex and evolving value chain Membrane desalination: Reverse Osmosis Project development models: fine-tuning Hybridization of thermal and the appropriate risk-sharing model membrane desalination Bringing capital to desalination assets: A set of parameters to assess the performance an increasingly strategic issue and efficiency of desalination assets Case study of desalination in Israel: innovative 1.2. A BRIEF HISTORY AND financing schemes achieving some of the GEOGRAPHICAL DISTRIBUTION OF lowest desalinated water costs worldwide DESALINATION TECHNOLOGIES Case study of desalination in Singapore: The market
    [Show full text]
  • Strategies for the Simulation of Sea Ice Organic Chemistry: Arctic Tests and Development
    geosciences Article Strategies for the Simulation of Sea Ice Organic Chemistry: Arctic Tests and Development Scott Elliott 1,*, Nicole Jeffery 1, Elizabeth Hunke 1, Clara Deal 2, Meibing Jin 2 ID , Shanlin Wang 1, Emma Elliott Smith 3 and Samantha Oestreicher 4 1 Climate Ocean Sea Ice Modeling (COSIM), Los Alamos National Laboratory, Los Alamos, NM 87545, USA; [email protected] (N.J.); [email protected] (E.H.); [email protected] (S.W.) 2 International Arctic Research Center and University of Alaska, Fairbanks, AK 99775, USA; [email protected] (C.D.), [email protected] (M.J.) 3 Biology Department, University of New Mexico, Albuquerque, NM 87131, USA; [email protected] 4 Applied Mathematics, University of Minnesota, Minneapolis, MN 55455, USA; [email protected] * Correspondence: [email protected]; Tel.: +1-505-606-0118 Received: 11 April 2017; Accepted: 21 June 2017; Published: 14 July 2017 Abstract: A numerical mechanism connecting ice algal ecodynamics with the buildup of organic macromolecules is tested within modeled pan-Arctic brine channels. The simulations take place offline in a reduced representation of sea ice geochemistry. Physical driver quantities derive from the global sea ice code CICE, including snow cover, thickness and internal temperature. The framework is averaged over ten boreal biogeographic zones. Computed nutrient-light-salt limited algal growth supports grazing, mortality and carbon flow. Vertical transport is diffusive but responds to pore structure. Simulated bottom layer chlorophyll maxima are reasonable, though delayed by about a month relative to observations due to uncertainties in snow variability. Upper level biota arise intermittently during flooding events.
    [Show full text]
  • Distribution and Abundance of Select Trace Metals in Chukchi and Beaufort Sea Ice
    Distribution and Abundance of Select Trace Metals in Chukchi and Beaufort Sea Ice Principal Investigators Robert Rember1 Ana M. Aguilar-Islas2 Graduate Student Vincent Domena2 1International Arctic Research Center, University of Alaska Fairbanks 2College of Fisheries and Ocean Sciences, University of Alaska Fairbanks FINAL REPORT December 2016 OCS Study BOEM 2016-079 Contact Information: email: [email protected] phone: 907.474.6782 fax: 907.474.7204 Coastal Marine Institute College of Fisheries and Ocean Sciences University of Alaska Fairbanks P. O. Box 757220 Fairbanks, AK 99775-7220 This study was funded in part by the U.S. Department of the Interior, Bureau of Ocean Energy Management (BOEM) through Cooperative Agreement M13AC00002 between BOEM, Alaska Outer Continental Shelf Region, and the University of Alaska Fairbanks. This report, OCS Study BOEM 2016-079, is available through the Coastal Marine Institute, select federal depository libraries and can be accessed electronically at http://www.boem.gov/Alaska-Scientific-Publications. The views and conclusions contained in this document are those of the authors and should not be interpreted as representing the opinions or policies of the U.S. Government. Mention of trade names or commercial products does not constitute their endorsement by the U.S. Government. TABLE OF CONTENTS LIST OF FIGURES ..................................................................................................................................... iii LIST OF TABLES ......................................................................................................................................
    [Show full text]
  • Numerical Modeling of Brine Disposal for Gaza Central Seawater Desalination Plant
    MEDRC Series of R & D Reports MEDRC Project: 15-DC-003 Numerical Modeling of Brine Disposal for Gaza Central Seawater Desalination Plant MSc. Thesis By Hassan S. Al-Najjar Supervisors Dr. Mazen Abualtayef A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree of Master of Science in Civil Infrastructure Engineering. The Middle East Desalination Research Center Muscat Sultanate of Oman 2015 ISLAMIC UNIVERSITY OF GAZA DEANERY OF HIGH STUDIES FACULTY OF ENGINEERING CIVIL ENGINEERING DEPARTMENT INFRASTRUCTURE ENGINEERING Numerical Modeling of Brine Disposal for Gaza Central Seawater Desalination Plant ا ال ا ا ة ا ه ا Prepared By: Hassan S. Al-Najjar Supervised By: Dr. Mazen Abualtayef A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree of Master of Science in Civil Infrastructure Engineering. 1437 AH, 2015 AD " َوھُ َ اﱠ ِي َ َ َج ْاَ ْ َ ْ ِ ھَ َا َ ْبٌ ُ َ ٌات َوھَ َا ِ ْ ٌ أُ َجٌ َو َ َ+ َ* %َ ْ()َ'ُ َ& %َ ْ َز ً# َو ِ" ْ! ًا َ ْ ُ! ًرا " (ا/.ن:)53 In the name of Allah, Most Gracious, Most Merciful “It is He (Allah) who has let free the two bodies of flowing water: one palatable and sweet, and the other salt and bitter; yet has He made a barrier between them, a partition that is forbidden to be passed.” (Quran 25:53) I Abstract In Gaza, it is planned to construct one of the most important seawater desalination plant in the region of Levantine basin, the plant is named Gaza Central Seawater Desalination Plant (GCDP).
    [Show full text]
  • Vertical Redistribution of Salt and Layered Changes in Global Ocean Salinity
    ARTICLE https://doi.org/10.1038/s41467-019-11436-x OPEN Vertical redistribution of salt and layered changes in global ocean salinity Chao Liu1, Xinfeng Liang 1,6, Rui M. Ponte2, Nadya Vinogradova3,4 & Ou Wang5 Salinity is an essential proxy for estimating the global net freshwater input into the ocean. Due to the limited spatial and temporal coverage of the existing salinity measurements, previous studies of global salinity changes focused mostly on the surface and upper oceans. fl 1234567890():,; Here, we examine global ocean salinity changes and ocean vertical salt uxes over the full depth in a dynamically consistent and data-constrained ocean state estimate. The changes of the horizontally averaged salinity display a vertically layered structure, consistent with the profiles of the ocean vertical salt fluxes. For salinity changes in the relatively well-observed upper ocean, the contribution of vertical exchange of salt can be on the same order of the net surface freshwater input. The vertical redistribution of salt thus should be considered in inferring changes in global ocean salinity and the hydrological cycle from the surface and upper ocean measurements. 1 College of Marine Science, University of South Florida, St Petersburg, FL 33701, USA. 2 Atmospheric and Environmental Research, Lexington, MA 02421, USA. 3 NASA Headquarters, Science Mission Directorate, Washington, DC 20546, USA. 4 Cambridge Climate Institute, Somerville, MA 02145, USA. 5 Jet Propulsion Laboratory, Pasadena, CA 91109, USA. 6Present address: School of Marine Science and Policy, College of Earth, Ocean, and Environment, University of Delaware, 700 Pilottown Road, Lewes, DE 19958, USA. Correspondence and requests for materials should be addressed to X.L.
    [Show full text]
  • Brine Rejection and Hydrate Formation Upon Freezing of Nacl Aqueous
    Brine rejection and hydrate formation upon freezing of NaCl aqueous solutions Ifigeneia Tsironi,a† Daniel Schlesinger,b† Alexander Späha, Lars Erikssonc, Mo Segadc and Fivos Perakis*a Studying the freezing of saltwater on a molecular level is of fundamental importance for improving freeze desalination techniques. Here, we investigate the freezing process of NaCl solutions using a combination of x-ray diffraction and molecular dynamics simulations (MD) for different salt-water concentrations, ranging from seawater conditions to saturation. A linear superposition model reproduces well the brine rejection due to hexagonal ice Ih formation and allows us to quantify the fraction of ice and brine. Furthermore, upon cooling at T = 233 K we observe the formation of NaCl×2H2O hydrates (hydrohalites), which coexist with ice Ih. MD simulations are utilized to model the formation of NaCl crystallites. From the simulations we estimate that the salinity of the newly produced ice is 0.5% mass percent (m/m) due to ion inclusions, which is within the salinity limits of fresh water. In addition, we show the effect of ions on the local ice structure using the tetrahedrality parameter and follow the crystalite formation by using the ion coordination parameter and cluster analysis. 1. Introduction from crystalline ice due to the formation of polycrystalline micro-domains forming upon crystallization22. Experiments Freeze desalination, also known as cryo-desalination, freezing- indicate that frozen NaCl aqueous solutions can form different melting and freeze-thaw desalination, has been suggested as NaCl-water crystal phases23,24. In addition, theoretical an energy effective alternative to distillation processes1. This investigations indicate that on a molecular level, brine is due to the fact that the latent heat of freezing (334 kJ/kg) is rejection progresses through a disordered layer with significantly lower than that of evaporation (2257 kJ/kg)2.
    [Show full text]
  • Intermediate and Deep Water Paleoceanography of the Northern North Atlantic Over the Past 21,000 Years David J
    PALEOCEANOGRAPHY, VOL. 25, PA1211, doi:10.1029/2009PA001833, 2010 Click Here for Full Article Intermediate and deep water paleoceanography of the northern North Atlantic over the past 21,000 years David J. R. Thornalley,1,2 Harry Elderfield,1 and I. Nick McCave1 Received 27 July 2009; revised 22 October 2009; accepted 30 October 2009; published 25 March 2010. [1] Benthic foraminiferal stable isotope records from four high‐resolution sediment cores, forming a depth transect between 1237 m and 2303 m on the South Iceland Rise, have been used to reconstruct intermediate and deep water paleoceanographic changes in the northern North Atlantic during the last 21 ka (spanning Termination I and the Holocene). Typically, a sampling resolution of ∼100 years is attained. Deglacial core chronologies are accurately tied to North Greenland Ice Core Project (NGRIP) ice core records through the correlation of tephra layers and changes in the percent abundance of Neogloboquadrina pachyderma (sinistral) with transitions in NGRIP. The evolution from the glacial mode of circulation to the present regime is punctuated by two periods with low benthic d13Candd18O values, which do not lie on glacial or Holocene water mass mixing lines. These periods correlate with the late Younger Dryas/Early Holocene (11.5–12.2 ka) and Heinrich Stadial 1 (14.7–16.8 ka) during which time freshwater input and sea‐ice formation led to brine rejection both locally and as an overflow exported from the Nordic seas into the northern North Atlantic, as earlier reported by Meland et al. (2008). The export of brine with low d13C values from the Nordic seas complicates traditional interpretations of low d13C values during the deglaciation as incursions of southern sourced water, although the spatial extent of this brine is uncertain.
    [Show full text]