Results 61-80 of 136.
Bookmark and Share    
See detailEstimation of tidal boundary conditions and surface winds by assimilation of high-frequency radar surface currents in the German Bight
Barth, Alexander ULiege; Alvera Azcarate, Aïda ULiege; Gurgel, Klaus-Werner et al

Conference (2010)

Numerical ocean models are affected by errors of various origins: errors in the initial conditions, boundary conditions and atmospheric forcings, uncertainties in the turbulence parametrization and ... [more ▼]

Numerical ocean models are affected by errors of various origins: errors in the initial conditions, boundary conditions and atmospheric forcings, uncertainties in the turbulence parametrization and discretization errors. In data assimilation, observations are used to reduce the uncertainty in the model solution. Ensemble-based assimilation schemes are often implemented such that the expected error of the model solution is minimized. It is shown that the observations can also be used to obtain improved estimates of the, in general, poorly known boundary conditions and atmospheric forcings. An ensemble smoother scheme is presented to assimilate high-frequency (HF) radar surface currents to improve tidal boundary conditions and wind forcings of a circulation model of the German Bight. To create an ensemble of dynamically realistic tidal boundary conditions, a cost function is formulated which is directly related to the probability of each perturbation. This cost function ensures that the 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 using all observations within the assimilation period using the covariances of the ensemble simulation. The approach acts like a smoother scheme since past and future observations are taken into account. The final analysis is obtained by rerunning the model using the optimal perturbation of the boundary conditions. The analyzed model solution satisfies thus the model equations exactly and does not suffer from spurious adjustments often observed with sequential assimilation schemes. Model results are also compared to independent tide gauge data. The assimilation also reduces the model error compared to those sea level observations. The same scheme is also used to correct surface winds. Surface winds are crucial for accurately modeling the marine circulation in coastal waters. The method is validated directly by comparing the analyzed wind speed to in situ measurements and indirectly by assessing the impact of the corrected winds on sea surface temperature (SST) relative to satellite SST. [less ▲]

Detailed reference viewed: 25 (5 ULiège)
See detailEnsemble smoother for optimizing tidal boundary conditions and wind forcing by assimilation of High-Frequency Radar surface currents measurements of the German Bight
Barth, Alexander ULiege; Alvera Azcarate, Aïda ULiege; Staneva, Joanna et al

Conference (2010)

An ensemble smoother scheme is presented to assimilate HF radar surface currents to improve tidal boundary conditions and wind forcings of a circulation model of the German Bight. To create an ensemble of ... [more ▼]

An ensemble smoother scheme is presented to assimilate HF radar surface currents to improve tidal boundary conditions and wind forcings of a circulation model of the German Bight. To create an ensemble of dynamically realistic tidal boundary conditions, a cost function is formulated which is directly related to the probability of each perturbation. This cost function ensures that the 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 using all observations within the assimilation period using the covariances of the ensemble simulation. The approach acts like a smoother scheme since all observations are taken into account. Since the scheme aims to derive the optimal perturbation, it might be called Ensemble Perturbation Smoother. The final analysis is obtained by rerunning the model using the optimal perturbation to the boundary conditions. The analyzed model solution satisfies thus the model equations exactly and does not suffer from spurious adjustments often observed with sequential assimilation schemes. Model results are also compared to independent tide gage data. The assimilation did also reduce the model error compared to those sea level observations. The same scheme has also been used to correct surface winds. Surface winds are crucial for accurately modeling the marine circulation in coastal waters. The method is validated directly by comparing the analyzed wind speed to in situ measurements and indirectly by assessing the impact of the corrected winds on sea surface temperature (SST) relative to satellite SST. [less ▲]

Detailed reference viewed: 27 (1 ULiège)
See detailAssimilation of high-frequency radar surface currents measurements to optimize tidal boundary conditions and wind forcing
Barth, Alexander ULiege; Alvera Azcarate, Aïda ULiege; Gurgel, Klaus-Werner et al

Conference (2010)

An ensemble smoother scheme is presented to assimilate high-frequency (HF) radar surface currents to improve tidal boundary conditions and wind forcings of a circulation model of the German Bight. To ... [more ▼]

An ensemble smoother scheme is presented to assimilate high-frequency (HF) radar surface currents to improve tidal boundary conditions and wind forcings of a circulation model of the German Bight. To create an ensemble of dynamically realistic tidal boundary conditions, a cost function is formulated which is directly related to the probability of each perturbation. This cost function ensures that the 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 using all observations within the assimilation period using the covariances of the ensemble simulation. The approach acts like a smoother scheme since past and future observations are taken into account. The final analysis is obtained by rerunning the model using the optimal perturbation of the boundary conditions. The analyzed model solution satisfies thus the model equations exactly and does not suffer from spurious adjustments often observed with sequential assimilation schemes. Model results are also compared to independent tide gage data. The assimilation also reduces the model error compared to those sea level observations. The same scheme is also used to correct surface winds. Surface winds are crucial for accurately modeling the marine circulation in coastal waters. The method is validated directly by comparing the analyzed wind speed to in situ measurements and indirectly by assessing the impact of the corrected winds on sea surface temperature (SST) relative to satellite SST. [less ▲]

Detailed reference viewed: 19 (2 ULiège)
See detailEnsemble-based assimilation of high-frequency radar surface currents in regional ocean models
Barth, Alexander ULiege; Alvera Azcarate, Aïda ULiege; Beckers, Jean-Marie ULiege et al

Conference (2010)

The results of coastal ocean models depend critically on the accuracy of boundary and initial conditions and atmospheric forcing. The precision of coastal ocean models is limited among others by ... [more ▼]

The results of coastal ocean models depend critically on the accuracy of boundary and initial conditions and atmospheric forcing. The precision of coastal ocean models is limited among others by uncertainty in those forcing fields. Since high-frequency (HF) radar installations provide measurements over a relatively large area, the assimilation of these data has a high potential to reduce the errors in ocean models and to provide a dynamically consistent estimation of the ocean circulation. The assimilation of HF radar data is not without its own challenges: the spatial variation of the surface currents uncertainty, the high temporal resolution of HF radar data, the simultaneous presence of a wide range of processes with distinct spatial and temporal scales (tides and other surface gravity waves, mesoscale and wind-driven circulation), and the generally strong sensitivity of regional models to errors in the boundary conditions and atmospheric forcings. These processess are important aspects to consider in the application of data assimilation methods to HF radar measurements. The results of two data assimilation experiments on the West Florida Shelf (WFS) and the German Bight are presented. HF radar currents are assimilated in a nested West Florida Shelf based on an ensemble of model realizations with different wind forcings. The model is sequentially updated and a filter is implemented to reduce spurious surface-gravity waves. Results of the WFS model assimilating surface currents show an improvement of the model currents not only at the surface but also at depth compared to independent ADCP observations. This West Florida Shelf assimilation experiment does not include tides. Tides are not generated within the domain, but are rather propagated inside the domain through the boundary conditions. The potential of using HF radar data to reduce errors in tidal boundary conditions is shown in a model setup of the German Bight. For improving the modeled tidal variability it is not sufficient to update the model state without updating the boundary conditions. An ensemble smoother to improve the tidal boundary values is presented and validated with independent HF radar measurements and tide-gage data. The ensemble-scheme is also applied to improve the wind forcing by assimilation of surface currents. The improvement of the analyzed wind forcing is assessed by using in-situ wind measurements. [less ▲]

Detailed reference viewed: 35 (2 ULiège)
Full Text
See detailDIVA: new features
Beckers, Jean-Marie ULiege; Alvera Azcarate, Aïda ULiege; Barth, Alexander ULiege et al

Scientific conference (2009, October 23)

Detailed reference viewed: 21 (3 ULiège)
See detailHF radar observations of surface currents in the German Bight: descriptive analysis, model-data comparison and non-sequential ensemble data assimilation
Port, A.; Staneva, J.; Schulz-Stellenfleth, J. et al

Conference (2009, September)

Detailed reference viewed: 44 (3 ULiège)
See detailEnsemble smoother for optimizing tidal boundary conditions and bottom roughness by assimilation of High-Frequency Radar surface currents
Barth, Alexander ULiege; Alvera Azcarate, Aïda ULiege; Staneva, J. et al

Conference (2009, September)

High-Frequency (HF) radars measure the ocean currents at various spatial and temporal scales. These include tidal currents, wind-driven circulation, density-driven circulation and Stokes drift. Sequential ... [more ▼]

High-Frequency (HF) radars measure the ocean 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 perturbation. This cost function ensures that the 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-dimension General Estuarine Ocean Model (GETM). Optimized boundary values are obtained using all observations within the assimilation period using the covariances of the ensemble simulation. [less ▲]

Detailed reference viewed: 54 (10 ULiège)
Full Text
See detailWeekly satellite sea surface temperature around Corsica, a DINEOF analysis of AVHRR data (1998), foreseeing comparison with interpolated and modelled fields.
Sirjacobs, Damien ULiege; Lenartz, Fabian ULiege; Troupin, Charles ULiege et al

Poster (2009, January)

Providing wide coverage and high spatio-temporal resolution, SST satellite archives are valuable sources of information for sound understanding of the ocean dynamics, including validation of ... [more ▼]

Providing wide coverage and high spatio-temporal resolution, SST satellite archives are valuable sources of information for sound understanding of the ocean dynamics, including validation of hydrodynamical modelling studies. Yet original SST fields have also many gaps (clouds, retrieval problems), but they are known to exhibit strong spatial and temporal correlations for regions of similar dynamics. This is exploited by the parameter free statistical technique DINEOF (Data Interpolation with Empirical Orthogonal Functions) [Alvera-Azcárate et al. (2005) Ocean Modell.; Beckers et al. (2006) Ocean Sciences] to produce full weekly analysis of the variability of the sea surface temperature (SST) around Corsica and in the Ligurian Sea at weekly temporal resolution during the year 1998. A detection of outliers implemented in DINEOF analysis is tested for pointing out unusual or invalid SST data. This study is realised foreseeing a comparison of DINEOF weekly averaged reconstructed fields with those obtained by interpolating methods on the same dataset (Data Interpolating Variationnal Analysis and Optimal Interpolation schemes), and with outputs of an implementation of the GHER 3D model in this area. [less ▲]

Detailed reference viewed: 93 (24 ULiège)
See detailSurface circulation of the Caribbean Sea and Gulf of Mexico using 13 years of satellite altimetry data
Alvera Azcarate, Aïda ULiege; Barth, Alexander ULiege; Weisberg, R. H.

Conference (2009)

The surface circulation of the Caribbean Sea and Gulf of Mexico is studied using thirteen years of satellite altimetry data. In the Caribbean Sea, processes at several temporal scales are analyzed: the ... [more ▼]

The surface circulation of the Caribbean Sea and Gulf of Mexico is studied using thirteen years of satellite altimetry data. In the Caribbean Sea, processes at several temporal scales are analyzed: the Caribbean eddies and meanders characteristics, the annual cycle and its variability through time, and the interannual variability, with a cycle of about 4 years affecting the SSH slope across the current and hence the ntensity of the Caribbean Current. Our analyses suggest that this cycle is related to changes in the wind intensity, wind curl and El Niño Southern Oscillation. In the Gulf of Mexico, the variability of the Loop Current is studied. We analyze the timing of anticyclonic eddy detachment from the Loop Current, the relation between the size of these eddies and the period between detachments, and the intrusion of the Loop Current into the Gulf of Mexico. A series of extreme Loop Current intrusions into the Gulf of Mexico, when the current is observed as far as 92°W, are described. The frequency of such events appears to have increased in recent years, with only one event happening from 1992 to 2002 (in 1993) versus three from 2002 to 2006. [less ▲]

Detailed reference viewed: 57 (1 ULiège)
See detailReconstruction of missing data in satellite data sets using DINEOF with constraints to reduce spurious high-frequency variations in the temporal EOFs
Alvera Azcarate, Aïda ULiege; Barth, Alexander ULiege; Sirjacobs, Damien ULiege et al

Conference (2009)

DINEOF (Data Interpolating Empirical Orthogonal Functions) is a method to reconstruct missing data in geophysical data sets, such as gaps originated by the presence of clouds in infrared satellite sensors ... [more ▼]

DINEOF (Data Interpolating Empirical Orthogonal Functions) is a method to reconstruct missing data in geophysical data sets, such as gaps originated by the presence of clouds in infrared satellite sensors. Based on Empirical Orthogonal Functions (EOFs), DINEOF uses an iterative procedure to calculate the missing values. DINEOF has been compared to Optimal Interpolation, showing that more accurate results are achieved, with up to 30 times less computational time (tests made with sea surface temperature of the Adriatic Sea, and validated with in situ data). Another advantage of this technology is that there is no need for a priori knowledge of the reconstructed data set statistics (such as covariance or correlation length). The technique can be applied to a broad range of data (physical, biological, chemical), and to a variety of platforms (satellite data, in situ data...). Given the nature of the EOFs, it is not necessary that data sets are regularly distributed in time. Irregularly distributed data sets, however, may lead to discontinuities in the temporal EOFs calculated from them, and these discontinuities can affect in turn the quality of the DINEOF reconstruction. In satellite data, some images can present a large amount of cloud cover, and only a few pixels with valid data. EOF projection to such images can also lead to discontinuities in the temporal modes, as there might be an over-fitting to the scarce information present in those images. After briefly describe DINEOF and its applications, we present a study aiming to reduce these discontinuities by including a time constraint to the covariance matrix used in the EOF decomposition. The approach is tested with sea surface temperature data of the Black Sea, and the results are compared to independent data. [less ▲]

Detailed reference viewed: 57 (9 ULiège)
See detailModelling error of a hydrodynamic model of the Mediterranean Sea
Vandenbulcke, Luc ULiege; Rixen, M.; Beckers, Jean-Marie ULiege et al

Conference (2009)

Detailed reference viewed: 20 (2 ULiège)
See detailApplication of a 3-D Super Ensemble to ocean forecast
Lenartz, Fabian ULiege; Barth, Alexander ULiege; Beckers, Jean-Marie ULiege et al

Conference (2009)

Super Ensemble (SE) techniques have recently allowed improving the forecast of various important oceanographic parameters, such as the significant wave height, the speed of sound or the surface drift, by ... [more ▼]

Super Ensemble (SE) techniques have recently allowed improving the forecast of various important oceanographic parameters, such as the significant wave height, the speed of sound or the surface drift, by correcting the prediction at a single or multiple locations, where data were available during the whole training period. However, nowadays common observation systems, such as satellite imagery or drifters, do not always provide information at the exact same locations, hence it is necessary to generalize the approach in order to take benefit of every image or track available. In this study, we try and apply a SE, fed with remote sensing and gliders data, to 3-D hydrodynamic models. The basic idea on which rely the SE methods is that a certain combination of several model runs and possibly data could yield better results than just one single model, even if it has a higher temporal or spatial resolution. As the most efficient techniques are the ones using observations, they rapidly developed and increased in complexity by copying what had been done in the data assimilation community; getting from the simple ensemble mean of the model outputs to their linear combination based on a particle filter. In our present study, we have decided to use the Kalman filter (KF) as it alleviates the need of an a priori determination of the training period length, and does not require the run of a very large ensemble of members. In addition, we apply it in a 3-D framework in order to take benefit of the spatial information contained by each source of measurements. For example, satellite images of sea surface temperature (SST) are very useful to correct the value of this parameter, but depending on the structure of the water column, it can also give a precious guess of how warm or cold is the ocean at 20 m deep. In our experiment the domain of interest is the Ligurian Sea during the last week of September, when part of the set-up for the CalVal08 campaign (SiC Charles Trees) had already taken place. The data assimilated during the training of the filter are SST images from AVHRR, as well as temperature and salinity profiles from two Rutgers University gliders. The models used for the study are three nested models of NCOM, run without data assimilation. The two considered variables are the temperature and the salinity. As our method is designed to work in a multivariate way, salinity forecast can possibly be improved by observing temperature profiles. Statistics are computed for both the training and the testing periods with an independent set of data. In four test cases, we review the impact of both the nature of the assimilated data, and the formulation of the model covariance matrix. At the end, we show that, on the basis of previous model outputs from which we’ve drawn an estimate of the model covariance, RMS error of the forecast in the whole 3-D domain can be reduced by 30%, thanks to the only assimilation of satellite SST images. [less ▲]

Detailed reference viewed: 57 (5 ULiège)
See detailA web interface for griding arbitrarily distributed in situ data based on Data-Interpolating Variational Analysis (Diva)
Barth, Alexander ULiege; Alvera Azcarate, Aïda ULiege; Troupin, Charles ULiege et al

Conference (2009)

Spatial interpolation of observations on a regular grid is a common task in many ceanographic disciplines (and geosciences in general). It is often used to create climatological maps for physical ... [more ▼]

Spatial interpolation of observations on a regular grid is a common task in many ceanographic disciplines (and geosciences in general). It is often used to create climatological maps for physical, biological or chemical parameters representing e.g. monthly or seasonally averaged fields. Since instantaneous observations can not be directly related to a field representing an average, simple spatial interpolation of observations is in general not acceptable. Diva (Data-Interpolating Variational Analysis) is an analysis tool which takes the error in the observations and the typical spatial scale of the underlying field into account. Barriers due to the coastline and the topography in general and also currents estimates (if available) are used to propagate the information of a given observation spatially. Diva is a command-line driven application written in Fortran and Shell Scripts. The observations and parameters are specified by the user using text files. The analyzed field and the expected error variance are returned as NetCDF files. This form of interaction with Diva is very similar to other high-performance codes and is a familiar approach for ocean modelers. However it represents a steep learning curve for oceanographers from other disciplines not familiar with command-line applications and programming. To make Diva easier to use, a web interface has been developed (http://gher-diva.phys.ulg.ac.be). Installation and compilation of Diva is therefore not required. The user can directly upload his/her data in ASCII format and enter several parameters for the analysis. The analyzed field, location of the observations, and the error mask are presented as different layers using the Web Map Service protocol. They are visualized in the browser using the Javascript library OpenLayers allowing the user to interact with layers (for example zooming and panning). Finally, the results can be downloaded as a NetCDF file, Matlab file (also readable in Octave, an open source program similar to Matlab) and Keyhole Markup Language (KML) file for visualisation in applications such as Google Earth. [less ▲]

Detailed reference viewed: 72 (7 ULiège)
See detailDynamically constrained ensemble perturbations. Application to tides on the West Florida Shelf
Barth, Alexander ULiege; Alvera Azcarate, Aïda ULiege; Beckers, Jean-Marie ULiege et al

Conference (2009)

A method is presented to create an ensemble of perturbations that satisfies linear dynamical constraints. A cost function is formulated defining the probability of each perturbation. It is shown that the ... [more ▼]

A method is presented to create an ensemble of perturbations that satisfies linear dynamical constraints. A cost function is formulated defining the probability of each perturbation. It is shown that the perturbations created with this approach take the land-sea mask into account in a similar way as variational analysis techniques. The impact of the land-sea mask is illustrated with an idealized configuration of a barrier island. Perturbations with a spatially variable correlation length can be also created by this approach. The method is applied to a realistic configuration of the West Florida Shelf to create perturbations of the M2 tidal parameters for elevation and depth-averaged currents. The perturbations are weakly constrained to satisfy the linear shallow-water equations. Despite that the constraint is derived from an idealized assumption, it is shown that this approach is applicable to a non-linear and baroclinic model. The amplitude of spurious transient motions created by constrained perturbations of initial and boundary conditions is significantly lower compared to perturbing the variables independently or to using only the momentum equation to compute the velocity perturbations from the elevation [less ▲]

Detailed reference viewed: 27 (9 ULiège)
See detailA West Florida Shelf ROMS model nested in HYCOM: validation, application and data assimilation
Barth, Alexander ULiege; Alvera Azcarate, Aïda ULiege; Zheng, L. et al

Conference (2008, November)

A high-resolution regional model for the West Florida Shelf based on ROMS is nested into the GODAE Atlantic and Global HYCOM models. The benefit of using HYCOM boundary conditions instead of climatology ... [more ▼]

A high-resolution regional model for the West Florida Shelf based on ROMS is nested into the GODAE Atlantic and Global HYCOM models. The benefit of using HYCOM boundary conditions instead of climatology is established by comparing the nested model to altimetry, in situ temperature time series, andADCP and high-frequency (HF) radar currents. Jointly with the Fish and Wildlife Research Institute, the model is applied to assist in harmful algae bloom forecasts. The model currents are used to estimate the trajectories of water parcels where red tide organisms have been detected. This system has been able to simulate the onset of the 2005 red tide event. HF Radar Currents are assimilated to improve the shelf circulation. An ensemble simulation of the WFS ROMS model is carried out under different wind forcings in order to estimate the error covariance of the model state vector and the covariance between ocean currents and winds. Radial currents measured by HF Radar antennas near St. Petersburg and Venice, FL, are assimilated using this ensemble-based error covariance. Different assimilation techniques using a time-average ensemble, a filter to reduce surface gravity waves and an extended state vector including wind stress were tested. Results of WFS ROMS model assimilating surface currents show an improvement of the model currents not only at the surface but also at depth. This model is run daily, to provide forecasts for the next 3.5 days. The model predictions and validations are available at http://ocgweb.marine.usf.edu under "WFS Model Forecasts". [less ▲]

Detailed reference viewed: 134 (1 ULiège)