revised tutorials

Author: baharmon

content/tutorials/earthworks/basics.ipynb

@@ -10,7 +10,7 @@
    "id": "8b05bf2d-4a12-459a-9b34-b1ade3769de3",
    "metadata": {},
    "source": [
-    "# Basic Earthworks\n",
+    "# Earthworks: Basics\n",
     "\n",
     "Brendan Harmon  \n",
     "2025-07-02\n",
@@ -172,7 +172,7 @@
     "[r.earthworks](https://grass.osgeo.org/grass-stable/manuals/addons/r.earthworks.html).\n",
     "These coordinates represents a local topographic maxima, i.e. a high\n",
     "point. The algorithm will fill locally up to that point, then the slopes\n",
-    "will falloff at a given rate of decay. Set `coordinates` to a pair of x-\n",
+    "will fall off at a given rate of decay. Set `coordinates` to a pair of x-\n",
     "and y-coordinates. Use the `z` parameter to set a z-coordinate for the\n",
     "top of the peak. Optionally use the `flat` parameter to create a plateau\n",
     "at the top of the peak. Using the default linear growth and decay\n",

content/tutorials/earthworks/basics.qmd

@@ -1,5 +1,5 @@
 ---
-title: "Basic earthworks"
+title: "Earthworks: Basics"
 author: "Brendan Harmon"
 date: 2025-07-02
 date-modified: today
@@ -38,7 +38,7 @@ jupyter: python3
 
 Learn the basics of terrain modeling with 
 {{< meta links.r_earthworks >}}.
-This tool uses cut and fill operations to tranform topography. 
+This tool uses cut and fill operations to transform topography. 
 Cut operations subtract from a topographic surface, 
 while fill operations add to the topography. 
 Transformations are based on proposed local topographic extrema,
@@ -81,6 +81,7 @@ with a Cartesian (XY) coordinate system.
 ## Command line
 
 ```{bash}
+grass --tmp-project XY
 ```
 
 ## Python
@@ -180,7 +181,7 @@ with {{< meta links.r_earthworks >}}.
 These coordinates represents a local topographic maxima, 
 i.e. a high point.
 The algorithm will fill locally up to that point,
-then the slopes will falloff at a given rate of decay.
+then the slopes will fall off at a given rate of decay.
 Set `coordinates` to a pair of x- and y-coordinates.
 Use the `z` parameter to set a z-coordinate for the top of the peak. 
 Optionally use the `flat` parameter to create a plateau at the top of the peak.

content/tutorials/earthworks/earthworks.qmd

@@ -1,6 +1,6 @@
 ---
-title: "Terrain modeling with r.earthworks in GRASS"
-author: "Brnedan Harmon"
+title: "Earthworks"
+author: "Brendan Harmon"
 date: 2025-07-07
 date-modified: today
 image: images/earthworks_01.webp
@@ -22,9 +22,7 @@ execute:
 
 ![Random fill operation with r.earthworks](images/earthworks_01.webp)
 
-<!-- Replace with new image -->
-
-Learn how model terrain with {{< meta links.r_earthworks >}}. With this tool you can generate new terrain or modify exisiting terrain using cut and fill operations inspired by earthworking. The first tutorial on the basics of terrain modeling introduces key concepts such as cut and fill operations. The subsequent tutorials introduce more advanced concepts such as relative operations and growth and decay functions. These tutorials demonstrate how to generate terrain, model natural landforms, and model anthropogenic landforms - i.e. earthworks - such as roads, dams, and levees. Read our [paper](https://raw.githubusercontent.com/baharmon/r.earthworks/main/paper/paper.pdf) to learn more. 
+Learn how model terrain with {{< meta links.r_earthworks >}}. With this tool you can generate new terrain or modify existing terrain using cut and fill operations inspired by earthworking. The first tutorial on the basics of terrain modeling introduces key concepts such as cut and fill operations. The subsequent tutorials introduce more advanced concepts such as relative operations and growth and decay functions. These tutorials demonstrate how to generate terrain, model natural landforms, and model anthropogenic landforms - i.e. earthworks - such as roads, dams, and levees. Read our [paper](https://raw.githubusercontent.com/baharmon/r.earthworks/main/paper/paper.pdf) to learn more. 
 
 <!-- Update link when paper published -->
 
@@ -55,7 +53,7 @@ Learn the basics of terrain modeling with r.earthworks.
 2. **[Gully modeling](./gullies.qmd)**  
 Model gullies using relative cut operations.
 
-3. **[Levee improvement & ridge restoration](./levees.qmd)**  
+3. **[Coastal infrastructure](./levees.qmd)**  
 Model coastal infrastructure projects including levees and restored alluvial ridges.
 
 4. **[Terrain synthesis](./synthesis.qmd)**  

content/tutorials/earthworks/gullies.ipynb

@@ -10,7 +10,7 @@
    "id": "5c102531-d33c-4a5a-94c3-94480ad449c5",
    "metadata": {},
    "source": [
-    "# Gully modeling\n",
+    "# Earthworks: Gully Modeling\n",
     "\n",
     "Brendan Harmon  \n",
     "2025-07-02\n",
@@ -28,13 +28,13 @@
     "and volumetric change calculations. While absolute operations calculate\n",
     "cut and fill based on a vertical datum of zero and can be used to model\n",
     "features at a given elevation, relative operations calculate cut and\n",
-    "fill relative to the exisiting topography and can be used to model\n",
-    "features that follow the exisiting terrain.\n",
+    "fill relative to the existing topography and can be used to model\n",
+    "features that follow the existing terrain.\n",
     "\n",
     "Use a relative earthworking operation with an exponential function to\n",
     "carve deep gullies into a landscape. In this tutorial all of the data\n",
     "will be created procedurally including the terrain and the streams.\n",
-    "Start by generating fractal terrain. Then model flow acculumation\n",
+    "Start by generating fractal terrain. Then model flow accumulation\n",
     "through the fractal terrain and extract its stream network. Finally use\n",
     "the vector map of streams as input for a relative cut operation to\n",
     "transform valleys into deep gullies. This tutorial covers:\n",

content/tutorials/earthworks/gullies.qmd

@@ -1,5 +1,5 @@
 ---
-title: "Gully modeling"
+title: "Earthworks: Gully Modeling"
 author: "Brendan Harmon"
 date: 2025-07-02
 date-modified: today
@@ -48,9 +48,9 @@ While absolute operations calculate cut and fill
 based on a vertical datum of zero
 and can be used to model features at a given elevation, 
 relative operations calculate cut and fill
-relative to the exisiting topography
+relative to the existing topography
 and can be used to model features 
-that follow the exisiting terrain.
+that follow the existing terrain.
 
 Use a relative earthworking operation 
 with an exponential function to 
@@ -59,7 +59,7 @@ In this tutorial
 all of the data will be created procedurally
 including the terrain and the streams.
 Start by generating fractal terrain. 
-Then model flow acculumation 
+Then model flow accumulation 
 through the fractal terrain
 and extract its stream network. 
 Finally use the vector map of streams 
@@ -95,6 +95,7 @@ with a Cartesian (XY) coordinate system.
 ## Command line
 
 ```{bash}
+grass --tmp-project XY
 ```
 
 ## Python

content/tutorials/earthworks/levees.ipynb

@@ -1,5 +1,5 @@
 ---
-title: "Modeling coastal infrastructure"
+title: "Earthworks: Modeling Coastal Infrastructure"
 author: "Brendan Harmon"
 date: 2025-07-06
 date-modified: today
@@ -70,13 +70,18 @@ Learn how to work with notebooks in the tutorial
 
 ## Project
 
-Start a GRASS session in the Bayou L'Ours project. 
+Download and unarchive the 
+[Bayou L'Ours Dataset](https://doi.org/10.5281/zenodo.15870441). 
+Then start a GRASS session in the Bayou L'Ours project. 
+Since the dataset is approximately 120MB, 
+it may take a couple minutes to download.
 
 ::: {.panel-tabset group="language"}
 
 ## Command line
 
 ```{bash}
+grass ~/grassdata/bayou_lours/PERMANENT
 ```
 
 ## Python
@@ -87,6 +92,8 @@ import os
 import sys
 import subprocess
 from pathlib import Path
+import urllib.request
+from zipfile import ZipFile
 
 # Find GRASS Python packages
 sys.path.append(
@@ -100,10 +107,26 @@ sys.path.append(
 import grass.script as gs
 import grass.jupyter as gj
 
-# Download dataset
+# Set GRASS database
+gisdbase = os.path.join(Path.home(), "grassdata")
 
-# Start GRASS in this project
-session = gj.init("~/grassdata/bayou_lours")
+# Download dataset
+url = "https://zenodo.org/records/15870442/files/bayou_lours.zip?download=1"
+filepath = os.path.join(gisdbase, "bayou_lours.zip")
+try:
+    urllib.request.urlretrieve(url, filepath)
+except Exception as e:
+    print(f"Error downloading file: {e}")
+
+# Unarchive dataset
+with ZipFile(filepath, 'r') as archive:
+    archive.extractall(gisdbase)
+
+# Delete archive
+os.remove(filepath)
+
+# Start GRASS in Bayou L'Ours project
+session = gj.init(Path(gisdbase, "bayou_lours"))
 ```
 
 :::
@@ -239,7 +262,7 @@ m.save("images/levees_02.webp")
 
 ## Flood Simulation
 
-Use {{< meta links.r_lake >}} to model innudation from storm surge. 
+Use {{< meta links.r_lake >}} to model inundation from storm surge. 
 Set the water level to one meter and the coordinates to `1109000, 102750`.
 The simulation will show that the degraded ridge
 would be overtopped by a meter of surge.
@@ -330,17 +353,22 @@ the borrow needed to build it.
 The restored ridge should have
 a crest elevation of 1.5 meters,
 a crest width of 15.2 meters,
-and 20% slopes with a ratio of 5H:1V.
+and 20% slopes with a ratio of 5H:1V,
+i.e. a horizontal run of 5 meters
+for a vertical rise of 1 meter.
 The soil to rebuild the ridge
-should be taken from a borrow area
+should be taken onsite 
+from an adjacent borrow area
 that would be offset by 7.6 meters 
 from the foot of the ridge. 
 
 ## Borrow Area
-
-Model the borrow area
-by offsetting the proposed ridge
-and then performing a cut operation.
+Start by modeling the borrow area - 
+an excavation for sourcing material for construction -
+needed to restore the ridge.
+To model the borrow area,
+offset the proposed ridge
+and then perform a cut operation.
 First use {{< meta links.v_transform >}}
 to shift the vector map of the proposed ridgeline 
 36 meters to the north.
@@ -405,13 +433,12 @@ Set `elevation` to the earthworks raster,
 `rate` to 0.2,
 and `flat` to 7.6. 
 
-
 ::: {.panel-tabset group="language"}
 
 ## Command line
 
 ```{bash}
-r.earthworks elevation=earthworks arthworks=earthworks operation=fill lines=ridges z=1.5 rate=0.2 flat=7.6 --overwrite
+r.earthworks elevation=earthworks earthworks=earthworks operation=fill lines=ridges z=1.5 rate=0.2 flat=7.6 --overwrite
 d.rast map=earthworks
 d.legend raster=earthworks color=white at=5,95,1,3
 ```
@@ -502,6 +529,9 @@ Use a fill operation with
 {{< meta links.r_earthworks >}}
 to model gap closures
 with sheet pile structures.
+These structures are walls 
+built from rows of interlocking vertical segments 
+of piles that are driven into the ground. 
 Set `elevation` to the earthworks raster,
 `lines` to the gaps vector,
 `z` to -0.6,
@@ -574,7 +604,7 @@ m.save("images/levees_08.webp")
 
 Test proposed coastal infrastructure with simulated storm surge.
 Use {{< meta links.r_lake >}} 
-to model innudation from storm surge over the earthworks raster.
+to model inundation from storm surge over the earthworks raster.
 Set the water level to one meter and the coordinates to `1109000, 102750`.
 The simulation will show that the restored ridge with closed gaps
 would be not overtopped by a meter of surge.

content/tutorials/earthworks/levees.qmd

@@ -1,5 +1,5 @@
 ---
-title: "Terrain synthesis"
+title: "Earthworks: Terrain Synthesis"
 author: "Brendan Harmon"
 date: 2025-06-30
 date-modified: today
@@ -87,6 +87,7 @@ with a Cartesian (XY) coordinate system.
 ## Command line
 
 ```{bash}
+grass --tmp-project XY
 ```
 
 ## Python