apply_taphonomy function

model taphonomic effects