DISTRIBUTION STATEMENT A: Approved for public release; distribution is unlimited. Eddy Resolving Global Ocean Prediction including Tides Alan J. Wallcraft NRL Code 7323 Stennis Space Center, MS 39529-5004 phone: (228) 688-4813 fax: (228) 688-4759 email: [email protected] Award Number: N0001409WX20491 http://www.hycom.org http://www7320.nrlssc.navy.mil/GLBhycom1-12/skill.html LONG-TERM GOALS Use the HYbrid Coordinate Ocean Model (HYCOM) with tides, dynamic sea ice, and data assimilation in an eddy-resolving, fully global ocean prediction system with 1/25° horizontal resolution that will run in real time at the Naval Oceanographic Office (NAVOCEANO) starting in 2012. The model will include shallow water and provide boundary conditions to finer resolution coastal models that may use HYCOM or a different model. OBJECTIVES To develop, evaluate, and investigate the dynamics of 1/25° global HYCOM (HYbrid Coordinate Ocean Model) with tides coupled to CICE (the Los Alamos sea ice model) with atmospheric forcing only, with data assimilation via NCODA (NRL Coupled Ocean Data Assimilation), and in forecast mode. Also to incorporate advances in dynamics and physics from the science community into the HYCOM established and maintained within the Navy. APPROACH Traditional ocean models use a single coordinate type to represent the vertical, but no single approach is optimal for the global ocean. Isopycnal (density tracking) layers are best in the deep stratified ocean, pressure levels (nearly fixed depths) provide high vertical resolution in the mixed layer, and σ- levels (terrain-following) are often the best choice in coastal regions. The generalized vertical coordinate in HYCOM allows a combination of all three types (and others), and it dynamically chooses the optimal distribution at every time step via the layered continuity equation. HYCOM use a C-grid, has scalable, portable computer codes that run efficiently on available DoD High Performance Computing (HPC) platforms, and has a data assimilation capability. Global HYCOM with 1/12° horizontal resolution at the equator (~7 km at mid-latitudes) is the ocean model component of the eddy-resolving nowcast/forecast system currently running in real time in the operational queue on the Cray XT5 at the Naval Oceanographic Office (NAVOCEANO). It provides nowcasts and forecasts of the three dimensional global ocean environment. HYCOM is coupled to the Los Alamos CICE sea-ice model (Hunke and Lipscomb, 2004) via the Earth System Modeling 1 Report Documentation Page Form Approved OMB No. 0704-0188 Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number. 1. REPORT DATE 3. DATES COVERED 2009 2. REPORT TYPE 00-00-2009 to 00-00-2009 4. TITLE AND SUBTITLE 5a. CONTRACT NUMBER Eddy Resolving Global Ocean Prediction including Tides 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 6. AUTHOR(S) 5d. PROJECT NUMBER 5e. TASK NUMBER 5f. WORK UNIT NUMBER 7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION Naval Research Laboratory,Code 7323,Stennis Space REPORT NUMBER Center,MS,39529-5004 9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSOR/MONITOR’S ACRONYM(S) 11. SPONSOR/MONITOR’S REPORT NUMBER(S) 12. DISTRIBUTION/AVAILABILITY STATEMENT Approved for public release; distribution unlimited 13. SUPPLEMENTARY NOTES 14. ABSTRACT 15. SUBJECT TERMS 16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF 18. NUMBER 19a. NAME OF ABSTRACT OF PAGES RESPONSIBLE PERSON a. REPORT b. ABSTRACT c. THIS PAGE Same as 8 unclassified unclassified unclassified Report (SAR) Standard Form 298 (Rev. 8-98) Prescribed by ANSI Std Z39-18 Framework (ESMF) (Hill et al., 2004), although currently for Arctic-only configurations. Coupling between the ocean and ice models more properly accounts for the momentum, heat and salt fluxes at the ocean/ice interface. The final component of the nowcast/forecast system is the NRL Coupled Ocean Data Assimilation (NCODA) which is a multivariate optimal interpolation scheme that assimilates surface observations from satellites, including altimeter and Multi-Channel Sea Surface Temperature (MCSST) data, sea ice concentration and also profile data such as XBTs (expendable bathythermographs), CTDs (conductivity temperature depth) and ARGO floats (Cummings, 2005). By combining these observations via data assimilation and using the dynamical interpolation skill of the model, the three dimensional ocean state can be accurately nowcast and forecast. The principal goal of this project is to perform the necessary R&D to prepare to provide a next- generation ocean nowcast/forecast system with real time depiction of the three-dimensional global ocean state at fine resolution (1/25° on the equator, 3.5 km at mid-latitudes, and 2 km in the Arctic). A major sub-goal of this effort is to test new capabilities in the existing 1/12° global HYCOM nowcast/forecast system and to transition some of these capabilities to NAVOCEANO in the 1/12° system, and others in the 1/25° global system. The new capabilities support (1) increased nowcast and forecast skill, the latter out to 30 days in many deep water regions, including regions of high Navy interest such as the Western Pacific and the Arabian Sea/Gulf of Oman, (2) boundary conditions for coastal models in very shallow, and (3) external and internal tides, the latter will initially be tested at 1/12°, to minimize computational cost, but will transition to NAVOCEANO only in the 1/25° system because at this resolution it will replace regional models with tides (all these will greatly benefit from the increase to 1/25° resolution). In addition to the NRL core tasking covered here, this effort will collaborate with a core team of similar size at FSU COAPS, with other parties interested in HYCOM development, and ONR field programs to test and validate the model in different regions and different regimes. Demonstrated advancements in HYCOM numerics and physics from all sources will be incorporated through this project. WORK COMPLETED HYCOM version 2.2.18 was made available on the HYCOM web site. This includes numerous improvements over version 2.1. The most significant of these is an improved vertical remapper, hybgen, which was implemented in FY08 but first used extensively this year. The first .04º global HYCOM simulations (3.5 km resolution at mid-latitudes) were run a total of 8.7 model years in collaboration with other projects, starting on 12 January 2009. This is the highest resolution ever for a global ocean model with high vertical resolution (2x finer horizontal grid spacing). It was run at NAVOCEANO using an HPC Challenge grant of computer time from the DoD HPC Modernization Program. For comparison, near twin .08º global HYCOM simulations were also integrated, including the same climatological 6-hourly atmospheric forcing (derived from the European Centre for Medium-Range Weather Forecasts 40-year reanalysis, ERA-40, with wind speed scaled using a monthly climatology from the QuikScat scatterometer; Kara et al., 2009) and bottom topography derived from the same .04º data set. However, the .08º topography was subjected to numerous hand edits in past years, mainly in shallow water and particularly to straits and sill depths, while the .04º global topography has received relatively few edits, mainly in FY09 to the Philippine and Indonesian Archipelagos. 2 The coupled CICE+HYCOM (sea-ice and ocean) model has been upgraded to the latest version of each model component (CICE 4.0 and HYCOM 2.2.18). It has been tested globally at 0.72 º resolution, but the primary focus this year was a 3.5 km Arctic cap domain that is an exact copy of the 1/12 º degree global tri-pole grid north of 40N. It is one-way off-line nested inside global HYCOM at 40N. Both non-assimilative and assimilative (NCODA) cases have been run. The latter to support validation of the 1/12 degree global system's sea ice capability RESULTS Figures 1 illustrates the difference that model improvements have made to the simulation of the Equatorial Pacific. The largest factor in the improvement is the use of higher-order vertical remapping in HYCOM version 2.2 vs 1st order remapping in version 2.1. The strength of the equatorial undercurrent is now more realistic and the thermocline is better defined. However, there is still too little upwelling in the Eastern Equatorial Pacific. Figure 2. Along-equator annual-mean sub-surface velocity (top) and temperature (bottom) in the Pacific Ocean from (left to right): (a) a 1/12° global simulation using HYCOM 2.1, (b) observations at the TAO array (Johnson et al, 2007), and (c) a 1/25° global simulation using HYCOM 2.2. 3 Figure 2. Zooms of the Gulf Stream region, all using the same color scale in cm. (a) Quasi- contemporaneous along-track RMS SSH variability from satellite altimeter data in four orbits overlaid on topographic contours (depth in m). The tracks are overlaid in the following order from top to bottom: (1) Envisat, (2) GFO, (3) Jason-1, and (4) Topex interleaved. (b,c) RMS SSH variability over model year 3 from non-assimilative (b) 1/12° and (c) 1/25° global HYCOM simulations. 4 Figure 3. Sea ice contentration in September 14, 2007 from two 3.5 km Arctic cap HYCOM+CICE cases atmospherically forced by NOGAPS : (a) non-assimilative (b) with NCODA assimilation. Heavy black line is the National Ice Center’s daily ice edge location. This data is not assimilated into the ice model. 5 Figure 2 compares the sea surface height variability in the Gulf Stream region from identical twin simulations at 1/12° and 1/25° after just three model years with along-track variability from four satellite altimeters. The Gulf Stream greatly benefits from higher resolution and the Stream is already penetrating further with more realistic variability at 1/25°. Figure 3 compares twin simulations of the coupled HYCOM and CICE Arctic cap system. The heavy black line is the National Ice Center’s daily ice edge location for this day. This data is not assimilated into the ice model, and so provides an independent indication of how well the model is performing. The non-assimilative simulation is already relatively accurate (particularly since the NIC’s ice edge probably represents about 15% concentration), but its area of high concentration is too high. Adding SSMI ice concentration via NCODA significantly improves the comparison to the ice edge, although there are areas where ice is underestimated. IMPACT/APPLICATIONS The 1/25° (3.5 km mid-latitude) resolution, first used in some FY09 global HYCOM simulations, is the highest so far for a global ocean model with high vertical resolution. A global ocean prediction system, based on 1/25° global HYCOM with tides, is planned for real-time operation starting in 2012. At this resolution, a global ocean prediction system can directly provide boundary conditions to nested relocatable models with ~1 km resolution anywhere in the world, a goal for operational ocean prediction at NAVOCEANO. Internal tides and other internal waves can have a large impact on acoustic propagation and transmission loss (Chin-Bing et al., 2003, Warn-Varnas et al., 2003, 2007), which in turn significantly impacts Navy anti-submarine warfare and surveillance capabilities. At present, regional and coastal models often include tidal forcing but internal tides are not included in their open boundary conditions. By including tidal forcing and assimilation in a fully 3-D global ocean model we will provide an internal tide capability everywhere, and allow nested models to include internal tides at their open boundaries. TRANSITIONS The 1/12° global HYCOM system is running in real time, and was transitioned to NAVOCEANO at the end of FY08 with operational testing is planned in FY10. It is receiving 6.4 SPAWAR funding (see below) for evaluation/validation. RELATED PROJECTS Partnering projects at NRL include 6.1 Global Remote Littoral Forcing via Deep Water Pathways, the 6.1 PhilEx DRI project Flow through the straits of the Philippine Archipelago simulated by global HYCOM and EAS NCOM, 6.1 Dynamics of the Indonesian Throughflow (ITF) and its remote impact, 6.2 Full Column Mixing for Numerical Ocean Models, 6.3 Battlespace Environments Institute – ESMF for Atmospheric-Ice-Ocean Coupling and Component Interoperability, 6.4 Large Scale Ocean Modeling, 6.4 Ocean Data Assimilation, and 6.4 Ice Modeling Assimilation from NPOESS. The computational effort is strongly supported by DoD HPC Challenge and NRL non-challenge grants of computer time. In FY09 1/25° and 1/12° global HYCOM ran under a new FY09-11 DoD HPC Challenge grant. 6 REFERENCES Chin-Bing, S.A., A. Warn-Varnas, D.B. King, K.G. Lamb, M. Teixeira, and J.A. Hawkins, 2003: Analysis of coupled oceanographic and acoustic soliton simulations in the Yellow Sea: a search for soliton-induced resonances. Mathematics and Computers in Simulation, 62, 11-20. Cummings, J.A., 2005. Operational multivariate ocean data assimilation. Quart. J. Royal Met. Soc., 131 (613), 3583-3604. Hill C., C. DeLuca, V. Balaji, M. Suarez, A. da Silva, 2004. The Architecture of the Earth System Modeling Framework. Computing in Science and Engineering, 6, 18-28. Hunke, E.C. and W.H. Lipscomb, 2004. CICE: the Los Alamos sea ice model documentation and software user's manual. http://climate.lanl.gov/Models/CICE Kara, A.B., A.J. Wallcraft, P.J. Martin, and R.L. Pauley, 2009. Optimizing surface winds using QuikSCAT measurements in the Mediterranean Sea during 2000–2006, J. Mar. Syst., doi:10.1016/j.jmarsys.2009.01.020. Warn-Varnas, A.C., S.A. Chin-Bing, D.B. King, Z. Hallock, and J.A. Hawkins, 2003. Ocean-acoustic solitary wave studies and predictions. Surveys in Geophysics, 24, 39-79. Warn-Varnas, A.C., S.A. Chin-Bing, D.B. King, J.A. Hawkins, K.G. Lamb, and J.F. Lynch, 2007. Winter PRIMER ocean-acoustic solitary wave modeling studies. IEEE J. Oceanic Engineering, 32(2), 436-452. PUBLICATIONS Arbic, B.K., A.J. Wallcraft and E.J. Metzger, 2009: Concurrent simulation of the eddying general circulation and tides in a global ocean model, Ocean Modelling. [submitted, refereed] Chassignet, E.P., H.E. Hurlburt, E.J. Metzger, O.M. Smedstad, J.A. Cummings, G.R. Halliwell, R. Bleck, R. Baraille, A.J. Wallcraft and C. Lozano, 2009: US GODAE Global Ocean Prediction with the HYbrid Coordinate Ocean Model (HYCOM), Oceanography, 22, 64-75. [published, refereed] Hurlburt, H.E., G.B. Brassington, Y. Drillet, M. Kamachi, M. Benkiran, R. Bourdalle-Badie, E.P. Chassignet, G.A. Jacobs, O. Le Galloudec, J.M. Lellouche, E.J. Metzger, O.M. Smedstad and A.J. Wallcraft, et al., 2009: High-Resolution Global and Basin-Scale Ocean Analyses and Forecasts, Oceanography, 22, 110-127. [published, refereed] Kara, A.B., R.W. Helber, T. Boyer and H.E. Hurlburt, 2009: Where in the global ocean threshold and curvature defined mixed layer depth differ? J Geophys. Res. [submitted, refereed] Kara, A.B., R.W. Helber, T. Boyer, A.J. Wallcraft and H.E. Hurlburt, 2009: Global mixed layer depth climatologies based on climatology and profile derived density fields. J Geophys. Res. [submitted, refereed] 7 Kara, A.B., A.J. Wallcraft, H.E. Hurlburt and W-Y Loh, 2009: Quantifying SST errors from an OGCM in relation to atmospheric forcing variables, Ocean Modelling, 29, 43-57. [published, refereed] Kara, A.B., A.J. Wallcraft, H.E. Hurlburt and W.Y. Loh, 2009: Which surface atmospheric variable drives the seasonal cycle of sea surface temperature over the global ocean?, J. Geophys. Res., 114, D05101, doi:10.1029/2008JD010420. [published, refereed] Metzger, E.J., H.E. Hurlburt, A.J. Wallcraft, J.F. Shriver, L.F. Smedstad, O.M. Smedstad, P. Thoppil and D.S. Franklin, 2008: Validtion Test Report for the Global Ocean Prediction System v3.0 - 1/12° HYCOM/NCODA: Phase 1, NRL Tech Report NRL/MR/7320--08-9148. [published] Posey, P.G., E.J. Metzger, A. J. Wallcraft, O.M. Smedstad and M.W. Phelps, 2009: Real-time data assimilation of ice concentration into a coupled ice/ocean forecast system. MTS/IEEE Conference Proceedings, Biloxi, MS, 26-29 Oct 2009. [submitted] Wallcraft A.J., A.B. Kara, H.E. Hurlburt, E.P. Chassignet and G.H. Halliwell, 2008: Value of bulk heat flux parameterizations for ocean SST prediction, J. Marine Sys., 74, doi:10.1016/j.marsys.2008.01.009. [published, refereed] Wallcraft, A.J., A.B. Kara, C.N. Barron, E.J. Metzger, R.L. Pauley and M.A. Bourassa, 2009: Comparisons of monthly mean 10 m wind speeds from satellites and NWP products over the global ocean, J. Geophys. Res., 114, D16109, doi:10.1029/2008JD011696. [published, refereed] Wallcraft, A.J., E.J. Metzger and O.M. Smedstad, 2009: Eddy Resolving Global Ocean Prediction including Tides, Proceedings of the 2009 HPCMP User's Group Conference, San Diego, CA. [submitted] Wallcraft, A.J., E.J. Metzger and S.N. Carroll, 2009: Software Design Description for the HYbrid Coordinate Ocean Model (HYCOM) Version 2.2, NRL Report NRL/MR/7320---0-9166. [published] 8