coastDat-2 TRIM-NP-2d-PCA_Residual_Currents

Simulated 2D residual velocity fields in the inner German Bight were subjected to Principal Component Analysis (PCA). Residual currents were obtained from coastDat2 barotropic 2D simulations with the hydrodynamic model TRIM-NP V2.1.22 in barotropic 2D mode on a Cartesian grid (1.6km spatial resolution) stored on an hourly basis for the years 1948 - 2012 (doi:10.1594/WDCC/coastDat-2_TRIM-NP-2d) and later extended until August 2015. The present analysis refers to the period Jan 1958 - Aug 2015. The spatial domain considered is the region to the east of 6 degrees east and to the south of 55.6 degrees north. All grid nodes with a bathymetry of less than 10m were excluded. Residual velocities were calculated in two different ways: 1.) as 25h means, 2.) as monthly means. Both types of residual current data are available from * RESIDUAL_CURRENTS_195801_201508 The directory contains sub-directories for years and months. Daily residual currents for the 13th of September 1974, for instance, are stored in * RESIDUAL_CURRENTS_195801_201508/YEAR_1974/MONTH_09/TRIM2D_1974_09_13_means.nc while monthly mean residual currents for September 1974 are stored in: * RESIDUAL_CURRENTS_195801_201508/YEAR_1974/TRIM2D_1974_09_means.nc All current fields provided were interpolated from the original Cartesian model grid to a more convenient regular geographical grid (116x76 nodes). Mean residual currents are stored in: * mean_residual_currents.nc This data set contains residual velocities both on original Cartesian grid nodes and interpolated to the geographical grid. An example plot is provided: * mean_residual_currents.png For PCA, two residual velocity components from each of 12133 Cartesian grid nodes were combined into one data vector (length 2x12133), referring to 21061 daily or 692 monthly time levels. Results of two independent PCAs for either daily or monthly mean fields are stored in: * PCA_daily_residual_currents.nc * PCA_monthly_residual_currents.nc Files contain three leading Principal Components (PCs) and corresponding Emipirical Orthogonal Functions (EOFs). Again EOFs were also interpolated to a regular geographical grid. PC time series are also stored in plain ASCII format: * PCs_daily.txt * PCs_monthly.txt For monthly fields the number N of variables (N=2x12133) is much larger than the number T of time levels (T=692). Therefore, to reduce computational demands, the roles of time and space were formally interchanged. Having conducted the PCA the EOFs were then transformed back to the original spatial coordinates (cf. Section 12.2.6 in von Storch and Zwiers (1999), Statistical Analysis in Climate Research, Cambridge University Press). A much larger number of time levels made even this approach prohibitive for the full set of daily data. Therefore, PCAs were performed for six sub-periods (1958-1965, 1966-1975, 1976-1985, 1986-1995, 1996-2005, 2006-2015(Aug)) independently. EOFs obtained from these six sub-periods were then averaged to obtain EOFs representative for the whole period. Corresponding PCs were calculated by projecting daily fields onto these average EOFs. IMPORTANT: In contrast with PCA of monthly data, the PCA of daily data INVOLVES SOME APPROXIMATIONS! EOFs on the original nodes were normalized to have unit lengths. The following figures, * daily_EOF1.png * daily_EOF2.png * daily_EOF3.png show the first three EOFs obtained from daily data, assuming that corresponding PCs have the value of one standard deviation. The following two plots, * monthly_EOF1.png * monthly_EOF2.png show the leading EOFs for monthly mean data. EOF3 is omitted as it represents just a very small percentage of overall variance (1.7%).