230 likes | 247 Views
A simple pre-operational model for the Portuguese coast. The INSEA perspective. Summary. Objectives, background Operational model planning Forcing models and external data Portuguese coastal circulation model Validation Automation, access and publishing tools Conclusions, status.
E N D
A simple pre-operational model for the Portuguese coast The INSEA perspective INSEA Athens dec'06
Summary • Objectives, background • Operational model planning • Forcing models and external data • Portuguese coastal circulation model • Validation • Automation, access and publishing tools • Conclusions, status INSEA Athens dec'06
Objectives • INSEA: Task 2.1 ”Implement and improve nesting techniques – Downscalling the Mercator solution for the Portuguese Coast”, Task 4.3 “Develop Remote Data Access” • Me: Thesis “Implement a tridimensional hydrodynamical circulation operational model for the portuguese coast, ready to couple biogeochemical models, as an effective tool to assess eutrophication problems in coastal marine environment such as Harmful Algae Blooms” INSEA Athens dec'06
Background A: Upwelling and large-scale currents near coast Q: What is the link between Tagus estuary’s mouth primary production and Portuguese coastal circulation? INSEA Athens dec'06
Background Q: what does it takes to model hydrodynamic circulation and upwelling on an operational basis? • Accurate winds forecasting • Large scale circulation • 3D mathematical models running on an operational basis INSEA Athens dec'06
Operationalization planning State-process diagram INSEA Athens dec'06
Automation scripts • Perl • Shell scripting • Service/daemon Extract and convert mercator data script flow diagram INSEA Athens dec'06
External data • Open Boundaries: Mercator 14 days forecast • Atmospheric forcing: MM5-IST 7 days forecast • Tide forcing: FES95 and FES2004 • Freshwater land discharges: estimates (so far) • Bathymetry: ETOPO 2’ Calibration/Validation data Remote-sensing imagery, meteorological stations, tidal stations, flow meters, buoys INSEA Athens dec'06
Mercator solution • PSY2V2 solution • resolution: 1/15º ~ 5-7 km • 43 z-coordinate layers • SLA, SST and in-situ data assimilation • Presents mediterranean outflow, meddy permitting (Drillet 2005) • ECMWF atmospheric forcing • 14 days forecast INSEA Athens dec'06
MM5-IST solution • 3 Domains: 81, 27 and 9 km • 25 vertical levels • ECMWF forcings • Winds • Latent heat • Sensible heat • Precipitation • Cloud cover • 7 days forecasts INSEA Athens dec'06
MOHID solution • MOHID Navier-stokes equations solver (Martins et al, 2001), and non-linear seawater equation of state (Jackett et al 1995), under the dimensional analysis considerations: • Beta plane • Hydrostatic approximation • Boussinesq approximation Finite volume discretization, Arakawa C-grid, generic z-level vertical coordinate, partial-step. Time discretization: ADI. Advection: TVD superbee. Turbulent vertical mixing: k-eps model (GOTM using Canuto2001 parameterizations). Wind induced mixing (Craig1994). INSEA Athens dec'06
MOHID solution • W domain: • 2D barotropic model • Lon [-13.7º -5.3º] W; • Lat [46.1º 33.5º] N • Resolution: 0.06º ~ 6.5 km • Sigma vertical coordinate • DT = 180 s • OBC: Radiation using Fes2004 as external solution (Blumberg1985) • Slow start over 5 inertial periods • Tide potential • Biharmonic filter: 1e9 m4/s INSEA Athens dec'06
MOHID solution • P domain: • 3D baroclinic model • Lon [-12.6º -5.5º] W • Lat [45.0º 34.4º] N • Resolution: 0.06º ~ 6.5 km • 42 layers, lagrangian vertical coordinate • DT = 180 s • OBC: Flather1976 radiation over father model + mercator barotropic solution; Flow Relaxation Scheme of S, T, U, V and level with mercator as local solution • Surface boundary forced with MM5 winds ramped over 6 inertial periods and with MM5 radiative forcing. • Tide potential. • Biharmonic filter: 1e10 m4/s • Tagus freshwater discharge estimates taken from 2004 data P domain 3D baroclinic model INSEA Athens dec'06
MOHID solution • E domain: • 3D baroclinic model • Lon [-11.18º -8.78º] W • Lat [40.3º 37.5º] N • Resolution: 0.02º ~ 2.2 km • 42 layers, lagrangian vertical coordinate • DT = 90 s • OBC: Flather1976 radiation over father model barotropic solution; Flow Relaxation Scheme of S, T, U, V and level with mercator as local solution • Surface boundary forced with MM5 winds ramped over 6 inertial periods and with MM5 radiative forcing. • Tide potential. • Biharmonic filter: 1e9 m4/s • Tagus freshwater discharge estimates taken from 2004 data INSEA Athens dec'06
MOHID results INSEA Athens dec'06
MOHID results INSEA Athens dec'06
MOHID results INSEA Athens dec'06
MOHID validation MOHID MODIS INSEA Athens dec'06 Mercator
Publishing process • Opendap server implemented • Live Access Server implemented (experimental) INSEA Athens dec'06
Operational Model Status • Data processing • Mercator fully automated • MM5 50% automated • Running process Works (Task 2.1), 60% automated • Validation process To be developped • Publishing process Implemented (Task 4.3), 40% automated, still experimental INSEA Athens dec'06
Future steps • Implement Primary Production model • Finish automation scripts • Develop validation methods • Improve the used downscalling technique INSEA Athens dec'06
References • Blayo E. and L. Debreu (2005). Revisiting open boundary conditions from the point of view of characteristic variables. Ocean Modelling 9 231–252, 2005. • Coelho Hs, Rjj Neves, M. White, Pc Leitão and Aj Santos. A model for ocean circulation on the Iberian coast Journal of Marine Systems (2002); 32: 153-179. • Adcroft, A.J., Hill, C.N. and Marshall, J. (1997), Representation of Topography by Shaved Cells in a Height Coordinate Ocean Model., Monthly Weather Review,125 (9), 2293, 1997. • Martins, F., R. Neves, P.C. Leitão, and A. Silva, 3D modeling in the Sado estuary using a new generic coordinate approach, oceanologica Acta, 24, S51-S62, 2001. • Flather, R.A., A tidal model of the northwest European continental shelf, Mem. Soc. R. Sci. Liege, 6 (10), 141-164, 1976. • Blumberg, A.F. and L.H. Kantha, 1985. Open boundary condition for circulation models. J. of Hydraulic Engineering, ASCE, 111, 237-2555. • Martinsen, E.A., and H. Engedahl., Implementation and testing of a lateral boundary scheme as an open boundary condition in a barotropic ocean model, Coastal Engineering, 11, 603-627, 1987. • Oey L.Y. (1999) A forcing mechanism for the poleward flow off the southern California coast. Journal of Geophysical Research, vol. 104, nºC6, pages 13529-13539, June 15, 1999. • Burchard, H., and K. Bolding, Comparative analysis of four second-moment turbulence closure models for the oceanic mixed layer, J. Phys. Oceanogr.,31, 1943-1968, 2001. • Canuto, V. M., A. Howard, Y. Cheng, and M. S. Dubovikov. Ocean turbulence Part I - one-point closure model. Momentum and heat vertical diffusivities. J.Phys. Oceanogr., 31,1413-1426, 2001. • Craig, P. D., and M. L. Banner, 1994, Modeling wave-enhanced turbulence in the ocean surface layer. J. Phys. Oceanogr., 24, 2546–2559. • Jackett et al 1995 • Drillet, Y., Bourdallé-Badie, R., Siefridt, L., Le Provost, C., 2005, Meddies in the Mercator North Atlantic and Mediterranean Sea eddy-resolving model, Journal of Geophysical Research, VOL. 110, C03016, doi:10.1029/2003JC002170. • Stevens I., M. Hamann, J. A. Johnson and A.F.G. Fiúza, 2000, Comparisons between a fine resolution model and observations in the Iberian shelf-slope region. Journal of Marine Systems, 26 53–74, INSEA Athens dec'06