PLUSNet Coupled Ocean and Acoustic Modeling Data Assimilation

PLUSNet Coupled Ocean and Acoustic Modeling, Data Assimilation, Predictions and Adaptive Sampling Recommendations Pierre Lermusiaux, Patrick Haley, Wayne Leslie, Oleg Logutov Ding Wang, Henrik Schmidt and Donald Eickstedt HARVARD – MIT Modeling, Predictions and Adaptive Sampling Team MB 06 – 07 Planning PLUSNet Meeting APL, UW, Oct 25, 2006 1

Harvard PLUSNet Goals and Objectives 2. 5 Environmentally Adaptive Sensing and Network Control PLUSNet Goal: To provide (sub)-mesoscale environmental fields/picture to MIT and PLUSNet, using new multi-scale environmental data-driven forecasting systems and new HU-MIT physical-acoustical adaptive sampling schemes Specific objectives are to: (i) Research and develop a new nested sub-mesoscale (non)-hydrostatic ocean modeling system within coarser regional domains for improved acoustic predictions (ii) Investigate and carry out physical-acoustical-seabed estimation and data assimilation (iii) Evaluate oceanic sub-mesocale parameterizations and study selected sub-mesoscale/mesoscale interactions and their acoustical impacts (iv) Collaborate with other efforts sponsored by ONR and NRL (v) Lead the environmental PLUSNet scientific research, coordinating the HU and Scripps contributions (e. g. internal tide conversions, sub-mesoscale eddy mixing and atmospheric forcing) 2

PLUSNet Major MB 06 Accomplishments 1. Daily ocean environment data assimilation and prediction. 2. Daily acoustic-environment prediction and TL prediction. 3. Daily optimization and recommendations for Adaptive Rapid Environmental Assessment (AREA). 4. Recommendations for capturing fronts. 5. Ocean data management and checks for PLUSNet assets, relay to MBARI COOP. 6. Web-based distribution of all our outputs and results. 3

http: //ocean. deas. harvard. edu/PLUSNet 4

http: //ocean. deas. harvard. edu/PLUSNet 5

Example of Daily Summary Description Emails PLUSNet • Environmental and acoustic nowcasts and predictions up to 00 Z Aug 26 are available from: http: //oceans. deas. harvard. edu/PLUSNet • http: //oceans. deas. harvard. edu/PLUSNet/Aug 25/plusnet_product. html • Product files are in: http: //oceans. deas. harvard. edu/PLUSNet/Products • Upwelling event is ending and our forecasts indicate that relaxation conditions (warmer surface temperatures and no wind-driven currents) will start tomorrow. See for example tomorrow's afternoon surface sound speed and current vectors: http: //oceans. deas. harvard. edu/PLUSNet/Aug 25/PN/Forecast/ssp 000_aug 27_fcst_00 Z_PN. gif • Based on our studies of the past week of transmission loss estimates, relaxation conditions correspond to -more- mean loss at 100 Hz and 400 Hz than upwelling conditions. The average difference at 100 -400 Hz and at 10 km range is about 3 to 8 db more mean loss. It would be very interesting to find out if the measured mean acoustic performance (TL)at these frequencies over 7 to 15 km ranges is lower tomorrow afternoon than it has been in the past days of upwelling conditions. • Vertically averaged current velocities in the PLUSNet ops area are expected to be to the northwest, following the coastline. For example, for tomorrow morning: http: //oceans. deas. harvard. edu/PLUSNet/Aug 25/PN/Forecast/vave_PN_aug 26_fcst_12 Z. gif • Surface velocities are expected to be patchy. For drifting behaviors, we recommend to combine it with a vertical yoyo so as to be drifting in vertically averaged currents, which are forecast to the northwest. • Transmission loss predictions available from: http: //oceans. deas. harvard. edu/PLUSNet/Prolog/Aug 25/index. html 6

Major Accomplishment I (cont. ) PLUSNet I: Daily Ocean Environment DA and Predictions: Methodology Harvard Ocean Prediction System Error Subspace Statistical Estimation Free-surface PE, Generalized biological Uncertainty forecasts, Ensemble-based, Multivariate DA, Adaptive sampling, models, Coupled to acoustic models, Towards multi-model estimates XML schemes to check configuration 7

Major Accomplishment I (cont. ) PLUSNet Free-Surface Ocean Model (HOPS) • Tidal and atmospheric forcing • Twice-daily data assimilation Nested Ocean Modeling with Grid-computing in Two Domains • Monterey Bay/San Francisco Domain: 1. 5 km resolution • PLUSNet - Ano Nuevo Domain: 0. 5 km resolution 8

Major Accomplishment I (cont. ) PLUSNet I: Daily Ocean Environment DA and Predictions “Undersea Weather” 6 4 OPAREA Bathymetry and Analysis Radials Surface SSP and Currents Measured and Estimated Sound Speed Radial 4 30 m Depth SSP and Currents Measured and Estimated Averaged Currents 0 -200 m 100 m Depth SSP and Currents See: http: //ocean. deas. harvard. edu/PLUSNet for twice daily plots and hourly data files (including tides) 9

APL- HU Collaboration: Use and Evaluation of PLUSNet Surface Drift Predictions by the Seaglider Team Bob Miyamoto, Bruce Howe and APL collaborators utilized our HOPS ocean model predictions of surface currents to plan their drifting missions APL uses the HOPS currents to compute a drift prediction. The result (green curve) is compared to the observed glider drift while on the surface. Case 1: better agreement than dreamed! Note sharp east-west gradient captured by the nested ocean models including tides Case 2: Not good. Atmos. and ocean models seem in phase error in time (see |u| change) and space (see weaker currents to the east)10

Major Accomplishment II PLUSNet II: Daily Coupled Acoustic-Environment and TL Predictions: MB 06 is the first time this has been done Seabed Modeling (thanks to M. Porter et al) 6 4 OPAREA Bathymetry and Analysis Radials Sound Speed Radial 6 Forecast Surface Currents (used by PLUSNet assets for planning drifting missions) TL Prediction (issued on Aug 22 12 Z) along Radial 6, for a source at 5 m and a receiver at 75 m See: http: //ocean. deas. harvard. edu/PLUSNet for twice daily TL plots and hourly data files

Major Accomplishment II (cont. ) PLUSNet II: Daily Acoustic-Environment and TL Predictions (cont. ) Relaxation: B 6 D 3 8_25_06 Upwelling: B 6 D 1 8_22_06 SVP(m/s) TL filed (d. B) SVP(m/s) f=100 Hz, sz=40 m TL filed (d. B) 12

Major Accomplishment II (cont. ) PLUSNet II: Daily Acoustic-Environment and TL Predictions (cont. ) Relaxation: B 6 D 3 8_25_06 Upwelling: B 6 D 1 8_22_06 Relaxation conditions (warmer surface temperatures and no winddriven currents) correspond to more mean loss at 100 Hz and 400 Hz than upwelling conditions: 3 to 8 db more mean loss over 7 to 15 km. 13

Major Accomplishment III PLUSNet III: Daily Optimization and Recommendations for Adaptive Sampling and AREA The True Ocean MIT SVP Forecast From HOPS/ESSE SVP Nowcast Data Assimilation Smaller Objective: Find the optimal path so as to minimize Acoustic Modeling Acoustic Prediction Uncertainty 14

Major Accomplishment III (cont. ) PLUSNet III: Daily Optimization and Recommendations for Adaptive Sampling and AREA Suboptimal predetermined path Suboptimal yoyo control Suboptimal predetermined path (2 -way) Suboptimal on-board adaptive path 15

Major Accomplishment III (cont. ) PLUSNet III: Daily Optimization and Recommendations for Adaptive Sampling and AREA Posteriori TL std map after the suboptimal predetermined path A priori TL std map Smoothed TL std (d. B) 0 3 (m) 50 2. 5 100 2 150 1. 5 200 1 0. 5 250 2 4 6 (km) 8 10 12 0 Above TL uncertainty estimates (std) do not account for most internal wave effects. They account for impacts of mesoscale ocean uncertainties and tidal effects. 16

Major Accomplishment III (cont. ) PLUSNet III: Daily Optimization and Recommendations for Adaptive Sampling and AREA Example of email sent daily: 8/22/06 Starting point: lat=36. 9414, long=-122. 2232 1. Bearing=180 degree clockwise from the north, waypoints are: r: 0 1 2 3 4 5 6 7 8 9 10 km z: 0 50 75. 0000 100. 0000 50. 0000 100. 0000 m 2. Bearing=180 degree clockwise from the north, max range=10 km, optimal yoyo control parameters are(10, 0. 5): XXX. Initialize(5, 100, 0, 10, 0. 5)//upper bound: 5 m; lower bound: 100 m; initial depth: 0 m; points: 10; threshold: 0. 5. 3. Bearing=180 degree clockwise from the north ADP<double> DW_ADP; DW_ADP. Initialize( ); 17

Major Accomplishment III (cont. ) PLUSNet III: Daily Optimization and Recommendations for Adaptive Sampling and AREA 18

PLUSNet Major Accomplishment IV IV: Recommendations for capturing fronts 1. Make in-situ measurements crossing the fronts. 2. Make a horizontal yoyo control focusing on the fronts. 8_24_06 Morning 8_24_06 Afternoon 3. Tracking temperature gradients. w 8_24_06 Morning w 8_24_06 Afternoon 19

Major Accomplishment IV (cont. ) PLUSNet IV: Recommendations for capturing fronts: issued twice Example of one of the email sent: 8/23/06 In the morning: 1. run predetermined path: Waypoints are: (lat=36. 916774452289 long=-122. 190015784337), (36. 959751554578 -122. 262152715663), (36. 966729966627, -122. 255703125301), (36. 923748958313 , -122. 183563374699), depth=0 m 2. or run this: Starting point: lat= 36. 916774452289 , long =-122. 190015784337, depth=0 m. Ding. Wang_2 D_Gradient<double> DW_Gradient; DW_Gradient. Initialize(270. , -1); //the initial AUV direction is 270 degree clockwise from the north; -1: follow the opposite gradient direction let AUV hit the northwest side of PLUSNet box, and then directly come back. In the afternoon: 3. Waypoints are: (36. 927933661928 -122. 179691928916), (36. 970917013855, -122. 251833371084), (36. 977895425904, -122. 245383780723), (36. 934908167952 -122. 173239519277), depth=0 m. 4. Starting point: lat= 36. 927933661928 , long = -122. 179691928916, depth=0 m. Ding. Wang_2 D_Gradient<double> DW_Gradient; DW_Gradient. Initialize(306. 663991. , -1); //the initial AUV direction is 306 degree clockwise from the north; -1: follow the opposite gradient direction let AUV hit the northwest side of PLUSNet box, and then directly come back. 20

Major Accomplishment IV (cont. ) PLUSNet IV: Recommendations for capturing fronts: Realized at sea once 21

HU MB 06 Deliverables and Needs PLUSNet Deliverables • Substantial measured T, S data and data-driven ocean-acoustic model estimates, which we will continue to improve and distribute • Model estimates also available in relocated area (Monterey Bay) Needs • To write collaborative papers, we need quantitative evidence (plots, data) to support successful or failed usage of our estimates, e. g: • Marc Stewart’s email based on our SE currents predictions in S corner of domain : “Pierre, You'll be interested to know that our Restech reported substantial SE surface current in the S corner of the PLUSNet box today as she attempted to recover a Bluefin in an Avon. I personally doubted it, but now I know better! Thanks, Marc” • For FY 07, PLUSNet should create a real-time coop web-page where such feedback and other data can be posted (e. g. Warren Fox’s suggestion: Data transfer interface / client-server architecture) 22

Research-Work: FY 07 Focus in Blue 1. Multi-Scale and Non-Hydrostatic Nested Ocean Modeling • Research and develop relocatable submesoscale nested modeling capability: • Higher-resolution hydrostatic model (Mini. HOPS) • HOPS coupled with non-hydrostatic models (2 D to 3 D: Lamb, Smolarkiewicz or MIT-GCM) • Compare parameterizations of submesoscales and boundary layers, and evaluate with HOPS and ROMS (run at HU, continue to collaborate with Scripps) PLUSNet 2. Coupled Physical-Acoustical Data Assimilation in real-time • Integrate and optimize physicalacoustical DA software with Mini-HOPS and AREA • Initiate coupled physical-acousticalseabed estimation and DA • Couple mini-HOPS/ESSE with selected sonar performance prediction (End-2 End System) 3. Acoustical-Physical Nonlinear Adaptive Sampling with ESSE and AREA • Implement and progressively demonstrate in FY 05 -06 -07 automated adaptive environmental sampling, integrating mini-HOPS and ESSE with AREA/NAFCON • Continue ocean-TL predictions (done for first time during MB 06) • Compute TL on denser spatial grids and for varied (higher) frequencies (with APL) • Provide field and uncertainty estimates to NAFCON and coordinate with Penn State for real-time display (towards end-to-end TDAs) 23

Papers with PLUSNet support PLUSNet Published/Submitted Lermusiaux P. F. J. , 2006. Uncertainty Estimation and Prediction for the Interdisciplinary Ocean. Special issue of the Journal of Computational Physics on ``Uncertainty Quantification''. James Glimm and George Karniadakis, Editors. 176 -199. Lermusiaux P. F. J, 2006. Adaptive Sampling, Adaptive Data Assimilation and Adaptive Modeling. Special issue of the J. of Inverse Problems on ``Mathematical Issues and Challenges in Data Assimilation for Geophysical Systems: Interdisciplinary Perspectives''. Christopher K. R. T. Jones and Kayo Ide, Eds. Under review. Wang D. , P. F. J. Lermusiaux, P. J. Haley, W. G. Leslie and H. Schmidt, 2006. Adaptive Acoustical. Environmental Assessment for the Focused Acoustic Field-05 At-sea Exercise, Proceedings of IEEE/MTS Oceans'06 Conference, Boston, MA, September 18 -21, 2006, 6 pp. In Press. Yilmaz N. K. , C. Evangelinos, N. M. Patrikalakis, P. F. J. Lermusiaux, P. J. Haley, W. G. Leslie, A. R. Robinson, D. Wang and H. Schmidt, 2006. Path Planning Methods for Adaptive Sampling of Environmental and Acoustical Ocean Fields, Proceedings of IEEE/MTS Oceans'06 Conference, Boston, MA, September 18 -21, 2006, 6 pp In Press. In Preparation (working titles) Wang, Schmidt, Lermusiaux et al. AREA for PLUSNet: MB 06 and FAF 05. To be prepared for IEEE. Lermusiaux, Wang, Schmidt, Haley, Leslie, Fox? , Wilson? , et al. Coupled ocean-acoustic Predictions during MB 06 for PLUSNet. Short paper to be prepared to JASA Shipley, Stewart, Schmidt et al. PLUSNet-MB 06 team short paper? (to do asap? ) Seaglider team and HU team. Planning of drifting missions using seaglider and HOPS? Curcio et al. Ocean data from autonomous kayaks for physical oceanography and ocean modeling, JRL 24