The goal of SSRS is to predict movements of soaring raptors (such as Golden Eagles) for a given wind conditions with the aim of determining potential negative conflict between soaring raptors and wind turbines. SSRS uses a stochastic agent-based model for predicting raptor movements through an orographic updraft field estimated using the spatially varying wind conditions and ground features (altitude, slope, aspect). SSRS can be applied to any rectangular region within the US without the need for any eagle-centric or atmosphere-related data collection efforts, using only the publicly available data sources. SSRS implements and extends the capability of the fluid-flow model from 'Brandes, D., & Ombalski, D. (2004). Modelling raptor migration pathways using a fluid-flow analogy. The Journal of Raptor Research, 38(3), 195-207.' More details can be found in the publication:
Rimple Sandhu, Charles Tripp, Eliot Quon, Regis Thedin, Michael Lawson, David Brandes, Christopher J. Farmer, Tricia A. Miller, Caroline Draxl, Paula Doubrawa, Lindy Williams, Adam E. Duerr, Melissa A. Braham, Todd Katzner,
Stochastic agent-based model for predicting turbine-scale raptor movements during updraft-subsidized directional flights <https://authors.elsevier.com/a/1eWBY15DJ-5ecl>_,
Ecological Modelling, Volume 466, 109876, 2022.
SSRS uses the following publically available data sources:
3D Elevation Program (3DEP) <https://www.usgs.gov/core-science-systems/ngp/3dep>_ dataset for terrain altitude, slope and aspect at a spatial resolution of 10 meters within the US.United States Wind Turbine Database (USWTDB) <https://eerscmap.usgs.gov/uswtdb/>_ for up-to-date turbine locations within the US.Wind ToolKit (WTK) <https://www.nrel.gov/grid/wind-toolkit.html>_ dataset for atmospheric conditions such as wind speed and direction at 1-hour temporal resolution and 2 km spatial resolution within the US.SSRS operates under three modes:
Clone the GitHub repository on local machine,
cd into the SSRS directory and using
conda virtual environments from
Anaconda <https://docs.anaconda.com/anaconda/install/index.html>_, do:
.. code-block:: bash
conda env create -f environment.yml
conda activate ssrs_env
pip install .For editable install (for development purpose), do
.. code-block:: bash
conda env create -f environment.yml
conda activate ssrs_env
pip install -e .For running conda environment ssrs_env in Jupyter Notebook (Comes with Anaconda),
.. code-block:: bash
conda install ipykernel
ipython kernel install --user --name=ssrs_envWithout Anaconda (requires python>=3.8 and pip>21.3):
.. code-block:: bash
pip install git+https://github.com/NREL/SSRS.git#egg=ssrsFor SSRS to access NREL's WTK dataset in the snapshot mode, need to get an API key from https://developer.nrel.gov/signup/ and save it in .hscfg file located in the working directory. Both examples/ and notebooks/ directories contain a sample .hscfg_need_api_key file. Make sure to rename this file as .hscfg after inserting your API key in this file.
The Jupyter notebooks in notebooks/ and python scripts in examples/ show the usage of SSRS for a given region. For instance, ssrs simulation in uniform mode for a 60 km by 50 km region in Wyoming and simulating 1000 eagles travelling north can be implemented using the following code:
.. _notebook: notebooks/sample_ssrs_uniform.ipynb
.. code-block:: python
config_wy_uniform = Config(
run_name='run_wy',
southwest_lonlat=(-106.21, 42.78),
region_width_km=(60., 50.),
resolution=100.,
sim_mode='uniform',
uniform_winddirn=270.,
uniform_windspeed=10.,
track_direction=0.,
track_count = 1000,
track_start_region=(5, 55, 1, 2)
)
sim = Simulator(config_wy_uniform)
sim.simulate_tracks()
sim.plot_terrain_elevation(show=True)
sim.plot_updrafts(show=True)
sim.plot_simulated_tracks(show=True)
sim.plot_presence_map(show=True)This will produce the following figures:
Ground elevation and turbine locations:
.. image:: docs/figs/elevation.png :width: 400 px :align: left :alt: Ground elevation and turbine locations
Orographic updrafts:
.. image:: docs/figs/s10d270_orograph.png :width: 400 px :align: left :alt: Orographic updrafts
1000 simulated tracks travelling towards north:
.. image:: docs/figs/s10d270_0_tracks.png :width: 400 px :align: right :alt:
Relative eagle presence density
.. image:: docs/figs/s10d270_0_presence.png :width: 400 px :align: right :alt: Relative eagle presence density
SSRS settings can be changed through a set of parameters defined using ssrs.Config attribute. The default setting can be viewed through following code:
.. code-block:: python
from ssrs import Config
print(Config())Here is a description of the parameters available to the users to vary:
.. code-block:: python
run_name: str = 'default' # name of this run, determines directory names
out_dir: str = os.path.join(os.path.abspath(os.path.curdir), 'output')
max_cores: int = 8 # maximum number of cores to use
sim_seed: int = -1 # random number seed
sim_mode: str = 'uniform' # snapshot, seasonal, uniform
print_verbose: bool = False # if want to print verboseParameters for setting up the region:
.. code-block:: python
southwest_lonlat: Tuple[float, float] = (-106.21, 42.78)
projected_crs: str = 'ESRI:102008' # ESRI, EPSG, PROJ4 or WKT string
region_width_km: Tuple[float, float] = (60., 50.)
resolution: int = 100. # desired terrain resolution (meters)Parameters for setting up the uniform mode:
.. code-block:: python
uniform_winddirn: float = 270. # northerly = 0., easterly = 90, westerly=270
uniform_windspeed: float = 10. # uniform wind speed in m/sParameters for setting up the snapshot mode:
.. code-block:: python
snapshot_datetime: Tuple[int, int, int, int] = (2010, 6, 17, 13)Parameters for setting up the seasonal mode:
.. code-block:: python
seasonal_start: Tuple[int, int] = (3, 20) # start of season (month, day)
seasonal_end: Tuple[int, int] = (5, 15) # end of season (month, day)
seasonal_timeofday: str = 'daytime' # morning, afternoon, evening, daytime
seasonal_count: int = 8 # number of seasonal updraft computationsParameters for importing data from NREL's WTK dataset:
.. code-block:: python
wtk_source: str = 'AWS' # 'EAGLE', 'AWS', 'EAGLE_LED'
wtk_orographic_height: int = 100 # WTK wind conditions at this height
wtk_thermal_height: int = 100 # WTK pressure, temperature, at this height
wtk_interp_type: str = 'linear' # 'nearest' 'linear' 'cubic'Parameters for simulating tracks:
.. code-block:: python
track_direction: str = 0. # movement direction measured clockwise from north
track_count: str = 1000 # number of simulated eagle tracks
track_start_region: Tuple[float, float, float, float] = (5, 55, 1, 2) # xmin, xmax, ymin, ymax [km]
track_start_type: str = 'random' # uniform, random
track_stochastic_nu: float = 1. # scaling of move probs, 0 = random walk
track_dirn_restrict: int = 3 # restrict within 45 deg of this previous movesParameters for plotting:
.. code-block:: python
fig_height: float = 6. # height of the figure window
fig_dpi: int = 200 # increase this to get finer plots
turbine_minimum_hubheight: float = 50. # for plotting turbine locations
turbine_mrkr_styles = ('1k', '2k', '3k', '4k',
'+k', 'xk', '*k', '.k', 'ok')
turbine_mrkr_size: float = 3. # marker size for plotting turbinesSandhu, Rimple, Tripp, Charles, Quon, Eliot, Thedin, Regis, Williams, Lindy, Doubrawa, Paula, Draxl, Caroline, and Lawson, Mike. SSRS (Stochastic Soaring Raptor Simulator). Computer Software. https://github.com/NREL/SSRS. USDOE Office of Energy Efficiency and Renewable Energy (EERE), Renewable Power Office. Wind Energy Technologies Office. 18 Oct. 2021. Web. doi:10.11578/dc.20210903.2.