impulseResponse {astrochron} | R Documentation |

## Impulse response function calculation

### Description

Calculate the analytical response function from an impulse forcing using the 1-D advection diffusion model proposed in Schink and Guinasso (1975) to model the bioturbation impact on climate time series.

### Usage

```
impulseResponse(G, ML = NULL, v = NULL, nt = 500, genplot = FALSE,
verbose = FALSE)
```

### Arguments

`G` |
Bioturbation parameter. G = D/ML/v |

`ML` |
Mix layer depth (cm) |

`v` |
Sedimentation rate (cm/kyr) |

`nt` |
Number of steps after the signal is deposited. |

`genplot` |
Generate summary plots? (T or F) |

`verbose` |
Verbose output? (T or F) |

### Value

fc Impulse response function

### References

Guinasso, N.L. and Schinck, D.R., 1975, *Quantitative estimates of biological mixing rates in abyssal sediments*, J. Geophys. Res., 80, 3032-3043.

Goreau, T.J., 1977, *Quantitative effects of sediment mixing on stratigraphy and biogeochemistry: a signal theory approach*, Nature, 256, 730-732.

Goreau, T.J., 1980, *Frequency sensitivity of the deep-sea climatic record*, Nature, 287, 620-622.

Liu, H., Meyers, S.R., and Marcott, S.A., 2021, *Unmixing dee-sea paleoclimate records: A study on bioturbation effects through convolution and deconvolution*, Earth and Planetary Science Letters, 564, 116883.

### Examples

```
G <- 4
ML <- 10
v <- 1
# take a look at the IRF
impulseResponse(G=4, ML = 10, v = 1, genplot = TRUE)
```

*astrochron*version 1.3 Index]