Enhancing temporal correlations in EOF expansions for the reconstruction of missing data using DINEOF

DINEOF (Data Interpolating Empirical Orthogonal Functions) is an EOF-based technique for the reconstruction of missing data in geophysical fields, such as those produced by clouds in sea surface temperature satellite images. A technique to reduce spurious time variability in DINEOF reconstructions is presented. The reconstruction of these images within a long time series using DINEOF can lead to large discontinuities in the reconstruction. Filtering the temporal covariance matrix allows to reduce this spurious variability and therefore more realistic reconstructions are obtained. The approach is tested in a three years sea surface temperature data set over the Black Sea. The effect of the filter in the temporal EOFs is presented, as well as some examples of the improvement achieved with the filtering in the SST reconstruction, both compared to the DINEOF approach without filtering

Ensemble perturbation smoother for optimizing tidal boundary conditions by assimilation of High-Frequency radar surface currents - application to the German Bight

High-Frequency (HF) radars measure the ocean surface currents at various spatial and temporal scales. These include tidal currents, wind-driven circulation, density-driven circulation and Stokes drift. Sequential assimilation methods updating the model state have been proven successful to correct the density-driven currents by assimilation of observations such as sea surface height, sea surface temperature and in-situ profiles. However, the situation is different for tides in coastal models since these are not generated within the domain, but are rather propagated inside the domain through the boundary conditions. For improving the modeled tidal variability it is therefore not sufficient to update the model state via data assimilation without updating the boundary conditions. The optimization of boundary conditions to match observations inside the domain is traditionally achieved through variational assimilation methods. In this work we present an ensemble smoother to improve the tidal boundary values so that the model represents more closely the observed currents. To create an ensemble of dynamically realistic boundary conditions, a cost function is formulated which is directly related to the probability of each boundary condition perturbation. This cost function ensures that the boundary condition perturbations are spatially smooth and that the structure of the perturbations satisfies approximately the harmonic linearized shallow water equations. Based on those perturbations an ensemble simulation is carried out using the full three-dimensional General Estuarine Ocean Model (GETM). Optimized boundary values are obtained by assimilating all observations using the co-variances of the ensemble simulation

Study of the combined effects of data assimilation and grid nesting in ocean models - application to the Gulf of Lions

Modern operational ocean forecasting systems routinely use data assimilation techniques in order to takeobservations into account in the hydrodynamic model. Moreover, as end users require higher and higher resolution predictions, especially in coastal zones, it is now common to run nested models, where the coastal model gets its opensea boundary conditions from a low-resolution global model. This configuration is used in the “Mediterranean Forecasting System: Towards environmental predictions” (MFSTEP) project. A global model covering the whole Mediterranean Sea is run weekly, performing 1 week of hindcast and a 10-day forecast. Regional models, using different codes and covering different areas, then use this forecast to implement boundary conditions. Local models in turn use the regional model forecasts for their own boundary conditions. This nested system has proven to be a viable and efficient system to achieve high-resolution weekly forecasts. However, when observations are available in some coastal zone, it remains unclear whether it is better to assimilate them in the global or local model. We perform twin experiments and assimilate observations in the global or in the local model, or in both of them together. We show that, when interested in the local models forecast and provided the global model fields are approximately correct, the best results are obtained when assimilating observations in the local model

Half a century of satellite remote sensing of sea-surface temperature

Sea-surface temperature (SST) was one of the first ocean variables to be studied from earth observation satellites. Pioneering images from infrared scanning radiometers revealed the complexity of the surface temperature fields, but these were derived from radiance measurements at orbital heights and included the effects of the intervening atmosphere. Corrections for the effects of the atmosphere to make quantitative estimates of the SST became possible when radiometers with multiple infrared channels were deployed in 1979. At the same time, imaging microwave radiometers with SST capabilities were also flown. Since then, SST has been derived from infrared and microwave radiometers on polar orbiting satellites and from infrared radiometers on geostationary spacecraft. As the performances of satellite radiometers and SST retrieval algorithms improved, accurate, global, high resolution, frequently sampled SST fields became fundamental to many research and operational activities. Here we provide an overview of the physics of the derivation of SST and the history of the development of satellite instruments over half a century. As demonstrated accuracies increased, they stimulated scientific research into the oceans, the coupled ocean-atmosphere system and the climate. We provide brief overviews of the development of some applications, including the feasibility of generating Climate Data Records. We summarize the important role of the Group for High Resolution SST (GHRSST) in providing a forum for scientists and operational practitioners to discuss problems and results, and to help coordinate activities world-wide, including alignment of data formatting and protocols and research. The challenges of burgeoning data volumes, data distribution and analysis have benefited from simultaneous progress in computing power, high capacity storage, and communications over the Internet, so we summarize the development and current capabilities of data archives. We conclude with an outlook of developments anticipated in the next decade or so

Amplified surface temperature response of cold, deep lakes to inter-annual air temperature variability

Summer lake surface water temperatures (LSWTs) have previously been shown to respond more rapidly to climatic warming compared to local summer surface air temperatures (SATs). In a global- scale analysis, we explore the factors underpinning the observation of an amplified response of summer LSWT to SAT variability using 20 years of satellite-derived temperatures from 144 lakes. We demonstrate that the degree of amplification in inter-annual summer LSWT is variable, and is greater for cold lakes (e.g. high latitude and high altitude), which are characterised by a short warming season, and deep lakes, that exhibit long correlation timescales of temperature anomalies due to increased thermal inertia. Such lakes are more likely to display responses in excess of local inter-annual summer SAT variability. Climatic modification of LSWT has numerous consequences for water quality and lake ecosystems, so quantifying this amplified response at a global scale is important