1 / 34

Atmosphere-Sea Hydrodynamic-Ecosystem model study in the sea

International Conference and Young Scientists School on Computational Information Technologies for Environmental Sciences: “CITES-2005” Novosibirsk, Russia, March 13-23, 2005. Atmosphere-Sea Hydrodynamic-Ecosystem model study in the sea. Rein Tamsalu (University of Tartu). Introduction.

dinah
Download Presentation

Atmosphere-Sea Hydrodynamic-Ecosystem model study in the sea

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. International Conference and Young Scientists School on Computational Information Technologies for Environmental Sciences: “CITES-2005”Novosibirsk, Russia, March 13-23, 2005 Atmosphere-Sea Hydrodynamic-Ecosystem model study in the sea Rein Tamsalu (University of Tartu)

  2. Introduction Today the environmental science are very much coupled with everyday life. Management policies need answer to concrete questions concerning the response of nature to both natural and manmade changes in enviromental forcing factors and loding. Numerical Hydro-Ecological modelling is an important tool for a better undrstanding of relations between the processes, and for forecasting these responses.

  3. Atm.-Sea-Hydro-ecological forecasting modelling The goal of our activities is to answer to concrete questions concerning the response of nature to both natural and man-made changes in marine environment. *

  4. Concrete Questions The influence of the Port of Tallinn reconstruction to the Muuga Bay marine environment

  5. The influence of the Port of Tallinn reconstruction to the Muuga Bay marine environment Baltic Sea=3’ =6’ Gulf of Finland=1’ =2’ Gulf of Muuga =0.05’ =0.1’ Talsingi=0.25’ =0.5’

  6. Atm.-Hydro-Ecological forecasting modelling In the Atmosphere-Sea-Hydro-Ecological modelling system are coupled several sub-models: *

  7. Atm.-Sea-Hydro-Ecological Models System

  8. Meso-Scale Atmosph. Model TU-EMHI model ETA model is based on the reference HIRLAM model, ETB is the non-hydrostaticversion. A pre-operational version of nonhydrostaticHIRLAM is used at the Estonian Meteorological Hydrological Institute (EMHI) to test the non-hydrostatic kernel of the model. Two modelling domains, as illustrated in Figure 1, are in use. Grid size of the larger domain ETA is 11km and the smaller domain ETB 3km. As the limited area models require boundary fields from larger models, the ETA model is nested to the FMI operational HIRLAM and ETB to the ETA.

  9. Meso- Scale Atmospheric Model

  10. Marine Circulation Models There are many different models barotropic baroclinic hydrostatic nonhydrostatic .......................

  11. Measured Temperaure Vertical Stucture in the Muuga Bay

  12. Velocity Measurements In Muuga Bay Recording Doppler Current Profiler RDCP 600 ( Aandera Instruments AS, Bergen, Norway.)

  13. Marine Circulation Model It is clear that we need Baroclinic Nonhydrostatic Circulation Model

  14. Marine Circulation Model The governing equations of the circulation model are: Momentum equation for velocity vectorU (u,v,w) Continuity equation for incompressible fluid • Transport-diffusion equations for: • SalinityS TemperatureT State equationfor buoyancyb=f(T,S) Two-equation turbulent model for Kinetic energy kand generic length scale quantity  or  Pressure components sea level fluctuation  , hydrostatic pressure p* and nonhydrostatic pressure p’ are calculated

  15. * Wind wave calculation • Surface wind waves are an integrated effect, in space and time, of driving wind fields. • The wind wave model computes the two-dimensional wave action spectra through integration of the transport equation, where the right hand side consists of several terms describing different evolution mechanisms, such as • energy input from wind ; • the non-linear transfer of energy through the spectrum ; • different kinds of dissipation mechanisms. • Interactive atmospheric input term is used in the Miles’ form. • In this model the so-called narrow-directional approximation is used. This approximation is based on well-known fact that wind waves have a narrow angular spreading function of spectra. The latter circumstance plays a key role in parameterization of nonlinear term.

  16. Size-DependentPankton Community Model • Size- dependent plankton communityfood web is formed by • autotrophs (Ai) i=1,2,…,NP • heterotrophs (Hi) i=1,2,…,NP • bacterioplankton (B) • This plankton communityforms theNPtripletstucture. ESD (m) Picophyto. Bacteriopl. I 0.2 - 2 DIP DIP+DOP Phytoflag. Zooflag. II 2 - 10 DIN DIN+DON Nanophyto. Nanozoopl. 10 - 50 III DIC DOC Netphyto. Microzoopl. 50 - 250 IV Mesozoopl. 250 – 1250

  17. | Growth reactions There are two energy flows to the plankton community The first one is the uptakeof dissolved inorganic nutrientsby autotrophy and it is directed from the autotrophy toward heterotrophy trough grazing. The other is theuptake of dissolved organic and inorganic nutrientsby bacterioplankton and it is directed from bacterioplankton towards heterotrophy trough predation.

  18. | Loss reactions Energy is lost through autotrophyexudation, mortality andrespiration heterotrophyexcretion, mortality andrespiration bacterioplanktonexcretionandrespiration detritusdecay

  19. Different grid resolutions Baltic Sea Gulf of Finland Talsinki area I - 3.0 * 3.0 nm II - 1.0 * 1.0 nm III- ¼ * ¼ nm IV- 1/20 * 1/20 nm Muuga Bay open boundary

  20. Horizontal velocity on the surface layer Muuga Bay Muuga Bay

  21. HORIZONTAL VELOCITY IN THE BOTTOM LAYER Muuga Bay Muuga Bay

  22. Velocity on the cross-section Muuga Bay

  23. Wind Waves In TALSINKI area during SW Storm

  24. Ecological compounds calculation Autotrophs in the beginning of May Heterotrophs in the beginning of May

  25. Suspended Material Calculation Spawning place Reconstruction area

  26. Oil Spill calculation The following oil spill processes are modeled: Transport and deformation of an oil slick due to time and spatially varying winds and currents Oil spreading at the sea surface due to positive buoyancy • Diffusion and dispersion of oil on the sea surface and in the water column • Evaporation of a multi-component mixture of oil Sinking of oil in water, and consequent sedimentation Formation of oil-in-water emulsion Weathering of oil, resulting in changes in density, viscosity, and water content, due to evaporation and emulsification processes Stranding of oil and shoreline interaction

  27. Oil Spill calculation The probability of the oil accident consequence in the NW part of the islandSaaremaa in the summer time during three months.

  28. The influence of the Wind Wavesto theBaroclinic Circulation

  29. Surface Temperature after 30 days calculation No Wind Waves Wind Waves are included

  30. Bottom Temperature after 30 days calculation No Wind Waves Tmax=10 Tmin=6.5 Wind Waves are included Tmax=10 Tmin=6.5

  31. Temperature profile after 30 days calculation No Wind Waves Tmax=10 Tmin=6.5 Wind Waves are included Tmax=10 Tmin=6.5

  32. Eddy Viscosity profile after 30 days calculation No Wind Waves Wind Waves are included Kmax=100cm2/s. Kmin=0.1 cm2/s.

  33. No Wind Waves Surface velocity after 30 days calculation Wind Waves are included

  34. Bottom Velocity after 30 days calculation No Wind Waves Wind Waves are included

More Related