Reduced-order optimal interpolation for biomass assimilation

G. Crispi, M. Pacciaroni, D. Viezzoli
2006 Ocean Science Discussions  
The analysis of the evolution of the chemical and biological characteristics of the Mediterranean marine ecosystem requires an integrated approach. Consistently, an ecosystem description, embedded in the MFSTEP one eighth degree threedimensional general circulation model, is tested and used following a twin experiment 5 approach. The ecosystem model is based on the NPZD trophic chain: inorganic nitrogen, N, phytoplankton, P, zooplankton, Z and detritus, D. Assimilation of synthetic biomass data
more » ... thetic biomass data is performed by means of the reduced-order optimal interpolation system SOFA. The synthetic "sea-truth" data are daily averages obtained from a sixty-nine days 10 reference run (RR). The twin experiments consist in performing two runs: a free run (FR) with wrong summer-depleted phytoplankton initial conditions and an assimilated run (AR), in which, starting from the same FR wrong phytoplankton concentrations, weekly averaged surface biomasses extracted from the RR results are assimilated. The comparison of the FR results with the AR ones shows a good convergence, on a 15 basin Mediterranean scale, confirming improvements of the forecasting in each of the four ecological compartments. Regional trophic regimes are analysed and interpreted in the western and eastern Mediterranean subbasins, for explaining the deteriorating behaviour of the total nitrogen. 20 ers of the ocean. Since polar-orbiting satellites' swaths cover the globe at high spatial 504 OSD 3, 503-539, 2006 Abstract 25 (Demirov and Pinardi, 2002). EGU plus a buffer zone representing the Atlantic inflow/outflow. Air-sea physical parameterizations account for the heat budget at the air-sea interface using sea surface temperature prognosed by the model. Six-hourly atmospheric fields from European Centre for Medium-Range Weather Forecast -ECMWF are used. The biharmonic horizontal eddy viscosity is 0.5×10 18 cm 4 s −1 and the biharmonic horizontal eddy diffusivity is 5 OSD 3, 2006 EGU Biogeochemical equations, without any recourse to boundary conditions inside the integration domain, are solved in insulating, conservative way through second-order finite difference method on the numerical B-grid, at the time step of 900 s (15 min). When instabilitites occur in the biogeochemistry, biological sources and sinks are set to zero and the calculations proceed after borrowing.
doi:10.5194/osd-3-503-2006 fatcat:l57yogr7zrgvjfipheabwuwun4