# Introduction

I got interested in layered sedimentation from viewing a video and decided it would be interesting to code this into R. More on that in due course, but my first step was to code a syatem with one sediment “type”.

# Procedure

The following code drops sediment particles at x=1, and lets them roll downhill until they reach the bottom or a ledge. It draws numbers at the sedimented particles’ final positions. Since the numbers start at 1, the values are like inverse ages.

# Discussion and conclusions

Reading the numbers on the graph as inverse age, one can see an interesting age structure through the sediments.

1. Viewed along diagonals (i.e. along a given depth below the sediment surface), age increase by 1 time unit with every lateral step away from the source.

2. Viewed as a vertical slice (i.e. a core), the time step per unit depth depends on distance from sediment source.

3. Viewed along Z levels (i.e. what might be a surface trace if something sliced the sediment pile horizontally), the time step per unit X distance depends on Z.

I wonder if these patterns (or the code) are of interest to geologists?

# Resources

This website is written in Jekyll, and the source is available on GitHub.