2016 Narrabeen Storm EWS Performance

This repository investigates whether the storm impacts (i.e. Sallenger, 2000) of the June 2016 Narrabeen Storm could have been forecasted in advance.

Repository and analysis format

This repository follows the Cookiecutter Data Science structure where possible. The analysis is done in python (look at the /src/ folder) with some interactive,exploratory notebooks located at /notebooks.

Development is conducted using a gitflow approach. The master branch stores the officialrelease history and the develop branch serves as an integration branch for features. Other hotfix and feature branches should be created and merged as necessary.

How to start?

Getting software requirements

The following requirements are needed to run various bits:

• Anacond: Used for processing and analysing data. The Anaconda distribution is used for managing environments and is available for Windows, Mac and Linux. Jupyter notebooks are used for exploratory analyis and communication.
• QGIS: Used for looking at raw LIDAR pre/post storm surveys and extracting dune crests/toes
• rclone: Data is not tracked by this repository, but is backed up to a remote Chris Leaman working directory located on the WRL coastal drive. Rclone is used to sync local and remote copies. Ensure rclone.exe is located on your PATH environment.
• gnuMake: A list of commands for processing data is provided in the ./Makefile. Use gnuMake to launch these commands. Ensure make.exe is located on your PATH environment.
• git

Getting the repository

Clone the repository onto into your local environment:

git clone http://git.wrl.unsw.edu.au:3000/chrisl/nsw-2016-storm-impact.git
cd nsw-2016-storm-impact

Getting the python environment set up

Commands for setting up the python environment are provided in the Makefile. Simply run the following commands in the repo root directory:

make venv-init
make venv-activate
make venv-requirements-install

You can see what these commands are actually running by inspecting the Makefile.

Pull data

The actual raw, interim and processed data are not tracked by the repository as part of good git practices. A copy of the raw data is stored on the WRL Coastal J:\ drive and can be copied using the following command.

make pull-data

If you have updated the data and want to copy it back to the J:\ drive, use the following command. Note that it is probably not a good idea to modify data stored in ./data/raw.

make push-data

View notebooks

Jupyter notebooks have been set up to help explore the data. Once you have set up your environment and pulled the data, this is probably a good place to start as you. To run the notebook, use the following command and navigate to the ./notebooks folder once the jupyter interface opens in your web browser.

make notebook

Notebooks use nbstripout to remove output before committing.

Available data

Raw, interim and processed data used in this analysis is kept in the /data/ folder. Data is not tracked in the repository due to size constraints, but stored locally. A mirror is kept of the coastal folder J drive which you can use to push/pull to, using rclone. In order to get the data, run make pull-data.

List of data:

• /data/raw/processed_shorelines: This data was recieved from Tom Beuzen in October 2018. It consists of pre/poststorm profiles at every 100 m sections along beaches ranging from Dee Why to Nambucca . Profiles are based on raw aerial LIDAR and were processed by Mitch Harley. Tides and waves (10 m contour and reverse shoaled deepwater) for each individual 100 m section is also provided.
• /data/raw/raw_lidar: This is the raw pre/post storm aerial LIDAR which was taken for the June 2016 storm. .las files are the raw files which have been processed into .tiff files using PDAL. Note that these files have not been corrected for systematic errors, so actual elevations should be taken from the processed_shorelines folder. Obtained November 2018 from Mitch Harley from the black external HDD labeled "UNSW LIDAR".
• /data/raw/profile_features: Dune toe and crest locations based on prestorm LIDAR. Refer to /notebooks/qgis.qgz as this shows how they were manually extracted. Note that the shapefiles only show the location (lat/lon) of the dune crest and toe. For actual elevations, these locations need to related to the processed shorelines.
• /data/raw/profile_features_tom_beuzen: This mat file contains dune toes and crests that Tom Beuzen picked out for each profile. This is used as a basis for the toe/crest locations, but is overridden from data contained in /data/raw/profile_features_chris_leaman.
• /data/raw/profile_features_chris_leaman: An excel file containing manually selected dune toes, crests, berms and impacts by Chris Leaman. The values in this file should take preceedence over values picked by Tom Beuzen.

Notebooks

• /notebooks/01_exploration.ipynb: Shows how to import processed shorelines, waves and tides. An interactive widget plots the location and cross sections.
• /notebooks/qgis.qgz: A QGIS file which is used to explore the aerial LIDAR data in /data/raw/raw_lidar. Byexamining the pre-strom lidar, dune crest and dune toe lines are manually extracted. These are stored in the /data/profile_features/.

TODO

• Mitch updated the raw profiles.mat to include more information about the survey time. Our data scripts should be updated to parse this new information and include it in our dataframes.
• Setup precomit hook for automatic code formatting using black. Low priority as can run black using the command make format.
• Raw tide WL's are interpolated based on location from tide gauges. This probably isn't the most accurate method, but should have a small effect since surge elevation was low during this event. Need to assess the effect of this method.
• Estimate max TWL from elevation where pre storm and post storm profiles are the same. Need to think more about this as runup impacting dune toe will move the dune face back, incorrectly raising the observed twl. Perhaps this estimation of max TWL is only useful for the swash regime.
• Implement bayesian change detection algorithm to help detect dune crests and toes from profiles. Probably low priority at the moment since we are doing manual detection.
• Implement dune impact calculations as per Palmsten & Holman. Calculation should be done in a new dataframe.
• Implement data/interim/*.csv file checking using py.test. Check for correct columns, number of nans etc. Testing of code is probably a lower priority than just checking the interim data files at the moment. Some functions which should be tested are the slope functions in forecast_twl.py, as these can be tricky with different profiles.
• Investigate using modin to help speed up analysis.
• Need to think about how relative imports are handled, see here. Maybe the click CLI interface should be moved to the ./src/ folder and it can import all the other packages?
• Simplify runup_models in Stockdon06 - we should really only have one function for each runup model. Need to make it work with individual values or entire dataframe. Use np.maskedarray

Misc

If Qt5Agg error comes up when trying to plot in PyCharm: https://stackoverflow.com/a/42231526