The observation, computation and study of “Lagrangian Coherent Structures” (LCS) in turbulent geophysical flows have been active areas of research in fluid mechanics for the last 30 years. Growing evidence for the existence of LCSs in geophysical flows (e.g., eddies, oscillating jets, chaotic mixing) and other fluid flows (e.g., separation prole at the surface of an airfoil, entrainment and detrainment by a vortex) generates an increasing interest for the extraction and understanding of these structures as well as their properties. In parallel, realistic ocean modeling with dense data assimilation has developed in the past decades and is now able to provide accurate nowcasts and predictions of ocean flow fields to study coherent structures. Robust numerical methods and sufficiently fast hardware are now available to compute real-time forecasts of oceanographic states and render associated coherent structures. It is therefore natural to expect the direct predictions of LCSs based on these advanced models. The impact of uncertainties on the coherent structures is becoming an increasingly important question for practical applications. The transfer of these uncertainties from the ocean state to the LCSs is an unexplored but intriguing scientific problem. These two questions are the motivation and focus of this presentation. Using the classic formalism of continuous-discrete estimation [1], the spatially discretized dynamics of the ocean state vector x and observations are described by (1a) dx =M(x; t) + d yok (1b) = H(xk; tk) + k where M and H are the model and measurement model operator, respectively. The stochastic forcings d and k are Wiener/Brownian motion processes,   N(0;Q(t)), and white Gaussian sequences, k  N(0;Rk), respectively. In other words, Efd(t)d T (t)g := Q(t) dt. The initial conditions are also uncertain and x(t0) is random with a prior PDF, p(x(t0)), i.e. x(t0) = bx0 + n(0) with n(0) random. Of course, vectors and operators in Eqs. (1a-b) are multivariate which impacts the PDFs: e.g. their moments are also multivariate. The estimation problem at time t consists of combining all available information on x(t), the dynamics and data (Eqs. 1a-b), their prior distributions and the initial conditions p(x(t0)). Defining the set of all observations prior to time t by yt

Physical and biogeochemical ocean dynamics can be intermittent and highly variable, and involve interactions on multiple scales. In general, the oceanic fields, processes and interactions that matter thus vary in time and space. For efficient forecasting, the structures and parameters of models must evolve and respond dynamically to new data injected into the executing prediction system. The conceptual basis of this adaptive modeling and corresponding computational scheme is the subject of this presentation. Specifically, we discuss the process of adaptive modeling for coupled physical and biogeochemical ocean models. The adaptivity is introduced within an interdisciplinary prediction system. Model-data misfits and data assimilation schemes are used to provide feedback from measurements to applications and modify the runtime behavior of the prediction system. Illustrative examples in Massachusetts Bay and Monterey Bay are presented to highlight ongoing progress.

Ocean science and ocean acoustics today are engaged in coupled interdisciplinary research on both fundamental dynamics and applications. In this context interdisciplinary data assimilation, which melds observations and fundamental dynamical models for field and parameter estimation is emerging as a novel and powerful methodology, but computational demands present challenging constraints which need to be overcome. These ideas are developed within the concept of an interdisciplinary system for assessing sonar system performance. An end-to-end system, which couples meteorology-physical oceanography-geoacoustics-ocean acoustics-bottom-noise-target-sonar data and models, is used to estimate uncertainties and their transfers and feedbacks. The approach to interdisciplinary data assimilation for this system importantly involves a full, interdisciplinary state vector and error covariance matrix. An idealized end-to-end system example is presented based upon the Shelfbreak PRIMER experiment in the Middle Atlantic Bight. Uncertainties in the physics are transferred to the acoustics and to a passive sonar using fully coupled physical and acoustical data assimilation.

Oceanographic conditions in the Gulf of Procchio, along the northern Elba coast, are influenced by the circulation in the Corsica channel and the southeastern Ligurian Sea. In order to support ocean prediction by nested models, an initial 4-day CTD survey provided initial ocean conditions. The purposes of the forecasts were threefold: i) in support of AUV exercises; ii) as an experiment in the development of rapid environmental assessment (REA) methodology; and, iii) as a rigorous real time test of a distributed ocean ocean prediction system technology. The Harvard Ocean Prediction System (HOPS) was set up around Elba in a very high resolution domain (225 m horizontally) which was two-way nested in a high resolution domain (675 m) in the channel between Italy and Corsica. The HOPS channel domain was physically interfaced with a one-way nest to the CU-POM model run in a larger Ligurian Sea domain. Eleven nowcasts and 2-3 day forecasts were issued during the period 26 September to 10 October, 2000 for the channel domain and for a Procchio Bay operational sub-domain of the Elba domain.

After initialization with the NRV Alliance, CTD survey data adaptive sampling patterns for nightly excursions of the Alliance were designed on the basis of forecasts to obtain data for assimilation which would most efficiently maintain the structures and variability of the flow in future dynamical forecasts. Images of satellite sea surface temperature were regularly processed and used for track planning and also for model verification. Rapid environmental assessment (REA) techniques were used for data processing and transmission from ship to shore and vice versa for model results. ADCP data validated well the flow in the channel. Additionally and importantly, the direction and strength of the flow in Procchio Bay were correctly forecast by dynamics supported only by external observations. CU-POM model hydrographic and geostrophic flow data was assimilated successfully on boundary strips of the HOPS domain. Flow fields with/without CU-POM nesting were qualitatively similar and a quantitative analysis of differences is under study. A significant result was the demonstration of a powerful and efficient distributed ocean observing and prediction system with in situ data collected in the Ligurian Sea, satellite data collected at SACLANTCEN, forecast modeling at Harvard University and the University of Colorado, and adaptive sampling tracks designed at Harvard. The distributed system functioned smoothly and effectively and coped with the adverse six-hour time difference between Massachusetts and Italy.

Data-driven forecasts and simulations for Massachusetts Bay based on in situ observations collected during August – September 1998 and on coupled four-dimensional (4-D) physical and biogeochemical models are carried out, evaluated, and studied. The real-time forecasting and adaptive sampling took place from August 17 to October 5, 1998. Simultaneous synoptic physical and biogeochemical data sets were obtained over a range of scales. For the real-time forecasts, the physical model was initialized using hydrographic data from August 1998 and the new biogeochemical model using historical data. The models were forced with real-time meteorological fields and the physical data were assimilated. The resulting interdisciplinary forecasts were robust and the Bay-scale biogeochemical variability was qualitatively well represented. For the postcruise simulations, the August – September 1998 biogeochemical data are utilized. Extensive comparisons of the coupled model fields with data allowed significant improvements of the biogeochemical model. All physical and biogeochemical data are assimilated using an optimal interpolation scheme. Within this scheme, an approximate biogeochemical balance and dynamical adjustments are utilized to derive the non-observed ecosystem variables from the observed ones. Several processes occurring in the lower trophic levels of Massachusetts Bay during the summer – autumn period over different spatial and temporal scales are described. The coupled dynamics is found to be more vigorous and diverse than previously thought to be the case in this period. For the biogeochemical dynamics, multiscale patchiness occurs. The locations of the patches are mainly defined by physical processes, but their strengths are mainly controlled by biogeochemical processes. The fluxes of nutrients into the euphotic zone are episodic and induced in part by atmospheric forcing. The quasi-weekly passage of storms gradually deepened the mixed layer and often altered the Bay-scale circulation and induced internal submesoscale variability. The physical variability increased the transfer of biogeochemical materials between the surface and deeper layers and modulated the biological processes.