Modelling the Channel and the Bay of Biscay using the MARS model. Contributions: Dyneco Physed / Pelagos / Benthos EMH

Transcription

1 Modelling the Channel and the Bay of Biscay using the MARS model Contributions: Dyneco Physed / Pelagos / Benthos EMH

2 Overview 1. The MANGA Configuration 2. Model validation: Sea Surface Temperature 3. Sensitivity studies 4. Numerical Methods: Advection schemes and diapycnal diffusion 5. Future developments 6. Model applications for the biology

3 MANGA configuration A Bay of Biscay / Channel regional modelling

4 25/12/ /07/ /07/2006 Physical spin-up Analyzed period Heat fluxes. From Luyten and Mulder (1992) to Large and Yeager (2004) and Fairall et al. (2003) Turbulence. From Gaspard et al (1990) to the Generic length scale turbulence closure model (Umlauf and Burchard, 2003)

5 MANGA Model MARS v. 8.18* Resolution Grid Coordinates 4km, 2.5km / 30 levels Spherical coordinates, Arakawa C Generalized sigma Bathymetry Composite (filtered bathymetry, Δh/h < 0.2) Vertical mixing TKE vertical mixing model (Gaspard, 1990) Rivers Atmospheric forcings Oceanic forcings Bottom friction Solar penetration 72 rivers Meteorological forcing from Arpege (Météo-France) OGCM : ORCA025 B83 in hindcast mode Tide : Sea level calculated by a barotropic 2D model (Mars2D) of the Western European Shelf forced by tides, atmospheric pressure and surface wind stress. Quadratic bottom friction (Z 0 formulation) Jerlov classification + monthly climatology (Gohin, OSD, 2011) * In this MARS version: - New numeric temporal scheme possibility to take = 0 : no explicit viscosity - Better numerical dispersion (better propagation of physical wave on a wider part of the spectrum) - Improvement of the standard ADI scheme on the second order error

6 Characteristics of MANGA configuration 4 km et 2.5 km horizontal resolution Two spatial extensions Bathy manga étendu

7 Characteristics of MANGA configuration 30 vertical generalized sigma layers z (1 ) h c (h h c )C( ) C( ) (1 b) sinh( ) sinh( ) b tanh[ ( 0.5)] tanh(0.5 ) 2tanh(0.5 ) b= 0.0; theta = 6.0 ; hc=9.0

8 Open boundary conditions Dynamics (dirichlet) Tracers

9 Model validation From Sea Surface Temperature

10 Observation coverage From Ocean and Sea Ice Satellite Application Facility (OSI SAF): SEVIRI SST Hourly data, midnight considered 0.05 grid

11 Mean bias Weak bias (<0.4 C) over 2 years except: above the continental slope Model along the western northern boundaries close to rivers Satellite

12 Seasonal mean bias DJF MAM Larger biases in summer Weak stratification too much vertical mixing JJA SON

13 Vertical section: 6 W

14 Sensitivity studies A set of parameters to tune

15 Comparison TKE and k-ε TKE k-ε Using k-ε: Improvement along the western boundary and in the Ushant front larger bias over the continental slope Further investigations needed for a final conclusion In the meantime, other experiments in progress (for example, using different bulk formulations Luyten & Mulder or Fairall)

16 Numerical methods Advection schemes and diapycnal diffusion

17 Description of the experiment (1/2)

18 Description of the experiment (2/2)

19 Results: sensitivity to the temporal integration method Methods efficient to reduce errors from implicit diffusion (x,y)/z t Diffusion increases with the CFL value

20 Results: sensitivity to the temporal integration method Methods efficient to reduce errors from implicit diffusion (x,y)/z t Diffusion increases with the CFL value

21 Future developments

22 MARS evolution in progress Few planed actions Assessment of turbulence closure schemes with 2 equations (k-ε, k-w, k-kl, generic) Parametrization of air/sea exchanges Open boundary conditions Influence of waves on the littoral circulation Generalised coordinates

23 Model applications for the biology From the biogeochemistry to the fish ecology

24 The derived indices from the Hindcast Stratification Fronts Plumes Eddies Primary Prod. A set of ~20 mesoscale and biological indices to characterise the oceanography in the BoB. and relate to fish ecology

25 Use in fisheries oceanography Survival maps of anchovy larvae (from homogeneous spawning (l), from observed spawning (c), at the end of drift) Seasonal pattern of drift for anchovy larvae from dark (beginning) to light (end of season) Drift and survival of anchovy larvae from an IBM forced by the hindcast

26 Hindcast use in fisheries oceanography Maps of indices Time-series River plumes Temperature, Currents Primary Prod. Environmental Monitoring Indicators Dynamic models of fish (IBM larvae / DEB adults) Habitat mapping

27 Development of the BGC model for the openwaters in the South of the BoB Impact of the offshore circulation on the biological production over the shelf WATER Silicate Nitrate Diatoms Zoo. P ad. Ammo. Phosphate Dino. µzoo. Equilibrium N 30% 100% P, Si Nano & Pico P Diss. N Diss. Si Diss. SEDIMENT Variable 2 (OW / C) P part. N part. Si part.

28 Mean August Chl.a 1- Fast remineralisation of the organic matter into the sediment (importance of the benthic activities and our limitation in the evaluation of the rivers impact). Coastal N Open Waters N

29 Mean August Chl.a 2 Improvement du to the hydrodynamic model evolution oN oN 9 52oN oN 48oN oN 5 48oN oN 3 46oN 8oW 6oW 4oW 2oW 0o 2oE 4oE oN 8oW 6oW 4oW 2oW 0o 2oE 4oE 0

30 ECO-MARS 3D in the Channel and the Bay of Biscay Different configurations

31 The biogeochemical model

32 Results example: Pseudo-Nitzchia Mean Pseudo-Nitzschia from (REPHY) Simulated biomass in 2003 Pseudo-Nitzschia Comptages REPHY Biomasse azotée simulée (strate 0-50m)