Simulate Wind Speed Time Series

Source: R/simulate_wind_speed.R

Generates synthetic wind speed climate series for multiple stations using stochastic atmospheric dynamics.

simulate_wind_speed(
  stations,
  time_index,
  ar_coeff = 0.6,
  seasonal_strength = 1,
  noise_sd = 1.5,
  extreme_event_prob = 0.01,
  extreme_multiplier = 2,
  min_ws = 0,
  max_ws = 40,
  seed = NULL
)

Arguments

stations

data.frame from create_stations()

time_index

data.frame from generate_time_index()

ar_coeff

Numeric. AR(1) persistence coefficient. Default = 0.6

seasonal_strength

Numeric. Controls seasonal wind variability. Default = 1

noise_sd

Numeric. Standard deviation of stochastic turbulence. Default = 1.5

extreme_event_prob

Numeric. Probability of extreme wind events. Default = 0.01

extreme_multiplier

Numeric. Multiplier applied during extreme events. Default = 2

min_ws

Numeric. Minimum allowable wind speed. Default = 0

max_ws

Numeric. Maximum allowable wind speed. Default = 40

seed

Optional numeric seed.

Value

data.frame containing:

Station

Station name

LON

Longitude

LAT

Latitude

ELEV

Elevation

DATE

Simulation timestamp

Year

Calendar year

Month

Calendar month

Season

Climatological season

WindSpeed

Simulated wind speed (m/s)

Extreme_Wind

Extreme wind event indicator

Wind_Anomaly

Wind speed anomaly

Details

The simulation incorporates:

  • seasonal circulation variability,

  • coastal enhancement,

  • elevation acceleration,

  • temporal persistence,

  • stochastic turbulence,

  • extreme wind events,

  • climate-zone effects.

Wind speed is simulated in m/s.

Examples

stations <- create_stations(
  n = 5,
  seed = 123
)
#> Generating synthetic station network...
#> Generated 5 synthetic stations within bounding box.
#> Deriving climate-aware station attributes...

time_index <- generate_time_index(
  start_date = "2000-01-01",
  end_date = "2005-12-31",
  frequency = "monthly"
)
#> Generated 72 time steps at monthly resolution.

wind <- simulate_wind_speed(
  stations,
  time_index
)
#> Wind speed simulation complete for 5 stations.

head(wind)
#>      Station       LON      LAT  ELEV CLIMATE_ZONE       DATE Year Month
#> 1  Station_1 -2.062112 4.818895 765.5      Coastal 2000-01-16 2000     1
#> 6  Station_1 -2.062112 4.818895 765.5      Coastal 2000-02-15 2000     2
#> 11 Station_1 -2.062112 4.818895 765.5      Coastal 2000-03-16 2000     3
#> 16 Station_1 -2.062112 4.818895 765.5      Coastal 2000-04-15 2000     4
#> 21 Station_1 -2.062112 4.818895 765.5      Coastal 2000-05-16 2000     5
#> 26 Station_1 -2.062112 4.818895 765.5      Coastal 2000-06-15 2000     6
#>     Season WindSpeed Extreme_Wind Wind_Anomaly
#> 1      Dry     17.75            0         4.51
#> 6      Dry     12.00            0        -0.53
#> 11 Pre-Wet     13.92            0        -0.02
#> 16 Pre-Wet     15.09            0         1.93
#> 21 Pre-Wet     12.63            0         0.71
#> 26     Wet     12.87            0         1.83