Super Point Definition

In order to account for all the possible energy and to avoid the shadowing effect due to land presence, the first step consists on defining a super-point (Cagigal et al., 2021), which aggregated the energy from a number of points surrounding the study site.

This virtual super-point has been constructed using the directional spectra from the CAWCR (Centre for Australian Weather and Climate Research) hindcast (Durrant et al., 2014; Smith et al., 2020), which provides spectral data at 10° resolution for the global grid, and 0.5° for the Pacific Region, in which our study is located.

# common
import warnings
warnings.filterwarnings('ignore')
import os
import os.path as op
import sys

# pip
import xarray as xr
import numpy as np

# DEV: bluemath
sys.path.insert(0, op.join(op.abspath(''), '..', '..', '..', '..', '..'))

# bluemath modules
from bluemath.superpoint.spectra import stations_superposition, bulkparams_partitions
from bluemath.superpoint.plotting.stations import Plot_stations
from bluemath.superpoint.plotting.spectra import Plot_stations_spectra, Plot_spectrum, Plot_seasons_spectra
from bluemath.superpoint.plotting.bulk_parameters import Plot_bulk_parameters, Plot_partitions

from bluemath.plotting.nearshore import Plot_obstructions_CAWCR

Database and site parameters

# database
p_data = r'/media/administrador/HD2/SamoaTonga/data'
site = 'Samoa'
p_site = op.join(p_data, 'Samoa')

# spectra stations
p_stat = op.join(p_site, 'stations')

# CAWCR obstructions
p_resources = r'/media/administrador/HD2/SamoaTonga/resources'
p_obstr_cawcr = op.join(p_resources, 'ww3_grid_inputs', 'Obstructions_CAWCR.nc')

# deliverable folder
p_deliv = op.join(p_site, 'd04_seasonal_forecast_swells')

# spectra data
p_spec = op.join(p_deliv, 'spec')

# figures
p_figs = op.join(p_deliv, 'figs')


# output files
p_superpoint = op.join(p_spec, 'super_point_superposition.nc')
p_chunks = op.join(p_spec, 'chunks')
p_partitions = op.join(p_spec, 'partitions_stats.nc')
p_bulk = op.join(p_spec, 'bulk_params.nc')

# make output folders
for p in [p_deliv, p_spec, p_figs, p_chunks]:
    if not op.isdir(p):
        os.makedirs(p)
        

# SUPERPOINT parameters
deg_sup         = 15  # degrees of superposition
stations_id     = [2193, 2195, 2187, 2178, 2170, 2168, 2166, 2175, 2183, 2191]
stations_sector = [(337.5,22.5), (22.5,60), (60,90), (90,120), (120,157.5), 
                   (157.5,202.5), (202.5,240), (240,270), (270,300), (300,337.5)]
stations_pos    = [9, 0, 1, 8, 2, 7, 3, 6, 5, 1]  # positions of the stations for plotting
st_wind         = 2195
lat_p           = -13.76
lon_p           = -172.07

# WAVESPECTRA parameters
part0 = True           # False if we dont want a sea partition (that will change wind to 0)
wcut  = 0.00000000001  # wind cut
msw   = 8              # max number of swells
agef  = 1.7            # age Factor
model = 'csiro'

Motivation

CAWCR obstructions

Below, a plot of the obstructions included in the CAWCR global wave hindcast is included.

• Grey cells correspond to land
  

• Red and purple cells correspond to obstructions in the X and Y direction respectively.
  
  The misinterpretation of island settings puts into evidence the need of performing a detailed downscaling of the wave dynamics.
  

# load CAWCR obstructions
obstr_cawcr = xr.open_dataset(p_obstr_cawcr)

# plot CAWCR obstructions
Plot_obstructions_CAWCR(
    obstr_cawcr, 
    area = [lon_p -1 + 360, lon_p + 0.75 + 360, lat_p - 0.75, lat_p + 0.85], 
    figsize = [18, 9],
);

Constructing the “Super Point”

Location of the selected stations used for the Super Point construction. Each station will represent one directional sector of the Super-Point spectra

# plot Super Point stations
fig = Plot_stations(
    p_stat, stations_id, 
    lon_p = lon_p, lat_p = lat_p,
    extra_area = 2, 
    figsize = [12, 8],
)

# store figure
p_f = op.join(p_figs, 'stations_samoa.png')
fig.savefig(p_f, dpi = 600, transparent = True)

# Super Point stations superposition
super_point = stations_superposition(p_stat, stations_id, stations_sector, deg_sup, st_wind)
super_point.to_netcdf(p_superpoint)  

#super_point = xr.open_dataset(p_superpoint)

# remove wind so no sea partition is created
if not part0:
    super_point.Wspeed.values = np.zeros(np.shape(super_point.Wspeed.values))

print(super_point)

Station: 2193
Station: 2195
Station: 2187
Station: 2178
Station: 2170
Station: 2168
Station: 2166
Station: 2175
Station: 2183
Station: 2191
<xarray.Dataset>
Dimensions:  (dir: 24, freq: 29, time: 366477)
Coordinates:
  * time     (time) datetime64[ns] 1979-01-01 1979-01-01T01:00:00 ... 2020-10-01
  * dir      (dir) float32 262.5 247.5 232.5 217.5 ... 322.5 307.5 292.5 277.5
  * freq     (freq) float32 0.035 0.0385 0.04235 ... 0.4171 0.4589 0.5047
Data variables:
    efth     (time, freq, dir) float64 0.0 0.0 0.0 ... 2.09e-08 1.097e-09
    Wspeed   (time) float32 4.388 4.388 4.762 5.175 ... 6.471 6.407 6.339 6.268
    Wdir     (time) float32 330.2 330.2 321.8 319.2 ... 85.15 87.38 90.47 93.54
    Depth    (time) float32 4.946e+03 4.946e+03 ... 4.946e+03 4.946e+03

Here is an example for a specific time:

time_index = 5000

This is the spectr at each of the different stations at the same time:

The shadowing effect of the island is clear for some directions, which emphasized the importance of constructing the Super Point in order to capture all the energy arrivign towards our study site.

Plot_stations_spectra(
    p_stat, stations_id, stations_pos,
    time_index,
    gs_rows = 2, gs_cols = 5,
    figsize = [20, 14],
);

This is the Super-Point at that time:

Plot_spectrum(super_point, time_index);

Mean Super-Point spectrum for the period 1979-2020:

Plot_spectrum(super_point, average=True);

Mean spectrum by season:

Plot_seasons_spectra(super_point);

Obtention of bulk parameters and partitions

Once the Super-Point is defined, we obtain the representative partitions and bulk parameters by means of the library wavespectra (https://github.com/metocean/wavespectra)

Nevertheless, bulk parameters and/or partitions will not be used for constructing the seasonal forecast

# obtain superpoint partitons and bulk parameters
bulk_params, stats_part = bulkparams_partitions(
    p_chunks, super_point, chunks=3, 
    wcut=wcut, msw=msw, agef=agef,
)

# store partitions and bulk parameters
bulk_params.to_netcdf(p_bulk)
stats_part.to_netcdf(p_partitions)

#bulk_params = xr.open_dataset(p_bulk)
#stats_part = xr.open_dataset(p_partitions)

chunk 1/3 done.
chunk 2/3 done.
chunk 3/3 done.

Bulk parameters:

Plot_bulk_parameters(bulk_params);

Partitions:

Plot_partitions(stats_part, num_fig=1);