## 6 Conclusions  
Summarizing the above, we can formulate the following main conclusions:  
1. The algorithm MatNERApor and its Matlab implementation were introduced for numerical modeling of the nonlinear response of porous saturated soil deposits to P- and SH-waves vertical propagation. The presented algorithm is based on the NERA algorithm, but develops significantly, expanding its applicability to the cases of a porous saturated soil structure, propagation of P-waves instead of SH-waves and occurrence of soil deposits in the water area. In addition, many shortcomings of the software implementation of the original algorithm have been overcome. We added the ability to run the calculation for a set of input seismic signals, soil profiles, shear and constrained moduli reduction curves with subsequent averaging of the output response spectra. The source code can be further developed.
2. The package was tested and verified using the records at two sites of the Kik-net network: the TCGH15 site is characterized by sandy soils lying on the bedrock, and the KSRH10 site is characterized by clayey soils. The effect of the water column on the reduction of vertical motion spectra was demonstrated using the records of local earthquakes obtained by ocean bottom seismographs in the Laptev Sea in 2019-2020. The results of the calculations showed good agreement with the data obtained from real seismic records, which justifies the correctness of the theoretical basis of the presented algorithm and its software implementation.
3. The MatNERApor package has significant prospects for practical use. This is especially actual for sites located in the water area within the zones of continuous or sparse permafrost, as well as the massive release of bubble gas from bottom sediments, such as the Arctic shelf of Eastern Siberia. In this case, it is especially important to consider the complex structure of marine soils, i.e., their saturation with water, gas mixture and ice.  
The MatNERApor package and the manual are available in GitHub repository: [[https://github.com/ArtWings/MatNERApor.git](https://github.com/ArtWings/MatNERApor.git)]([https://github.com/ArtWings/MatNERApor.git](https://github.com/ArtWings/MatNERApor.git)) (accessed on 3 May 2022).  
Conceptualization, A.A.K.; methodology, A.A.K.; software, A.A.K., E.A.R., M.A.N.; validation, A.A.K.; formal analysis, A.A.K.,; investigation, A.A.K., E.A.R., M.A.N.; data curation, A.A.K., L.L., I.P.S.; writing--original draft preparation, A.A.K.; writing--review and editing, all authors; visualization, A.A.K.; supervision, A.A.K., L.L., I.P.S.; project administration, A.A.K., K.A.R., L.L., I.P.S.; funding acquisition, A.A.K., O.S., K.A.R., I.P.S. All authors have read and agreed to the published version of the manuscript.  
This work was financially supported by the Russian Science Foundation, project no. 20-77-00034 (developing of the site response algorithm and its approach); Russian Foundation for Basic Research, project no. 20-05-00533 (collecting and processing the ocean bottom seismograph records, calculating spectral ratios of earthquake signals); the Russian State Assignment, project no. FMWE-2021-0004 (a review of publications on modeling the earthquake site response); the Ministry of Science and Higher Education of the Russian Federation, project no. 1021052204209-8 and no. FMWE-2021-0011 (discussion of the results in the broad interdisciplinary context, text editing); the Russian President Grant for young scientists, project no. MK-45.2022.1.5 (study of earthquake site response in context of its relation to endogenous geological hazards); the Russian Science Foundation, project no. 21-77-30001 (developing of the models of porous soil structure, typical for marine soils in the area of active methane seepage from seafloor); the Ministry of Science and High Education of the Russian Federation, project no. 075-15-2020-928, HSE (study of the seismic load effect on sub-bottom permafrost condition).  
Not applicable.  
The authors declare no conflict of interest.