The R script IRFnnhs.R, which efficiently estimates impulse response functions for environmental systems that are nonlinear, nonstationary, or heterogeneous, based on their input and output time series. Scripts and results for a series of benchmark tests are also provided, to accompany Kirchner, J.W., Impulse response functions for heterogeneous, nonstationary, and nonlinear systems, estimated by deconvolution and demixing of noisy time series, _Sensors_, 22(9), 3291, https://doi.org/10.3390/s22093291, 2022.