CPTU identification of regular, sensitive, and organic clays towards evaluating preconsolidation stress profiles

Shehab S Agaiby; Paul W Mayne; Shehab S Agaiby; Paul W Mayne

doi:10.3934/geosci.2021032

AIMS Geosciences

2021, Volume 7, Issue 4: 553-573. doi: 10.3934/geosci.2021032

Previous Article Next Article

Research article

CPTU identification of regular, sensitive, and organic clays towards evaluating preconsolidation stress profiles

Shehab S Agaiby ^{1
,
,},
Paul W Mayne ²

1.
Department of Public Works, Faculty of Engineering, Cairo University, Giza, Egypt
2.
Geosystems Engineering, Georgia Institute of Technology, Atlanta, GA 30332–0355 USA

Received: 05 August 2021 Accepted: 15 September 2021 Published: 13 October 2021

Soil classification by piezocone penetration tests (CPTU) is mainly accomplished using empirical soil behavior charts (SBT). While commonly-used SBT methods work well to separate fine-grained soils from granular coarse-grained soils, in many instances, the groupings often fail to properly identify different categories of clays, specifically: (a) "regular" clays that are inorganic and insensitive, (b) sensitive and quick clays; and (c) organic soils. Herein, a simple means of screening and sorting these three clay types is shown using three analytical CPTU expressions for evaluating the preconsolidation stress profile from net cone resistance, excess porewater pressure, and effective cone resistance. A number of case studies are utilized to convey the methodology.

Keywords:

Citation: Shehab S Agaiby, Paul W Mayne. CPTU identification of regular, sensitive, and organic clays towards evaluating preconsolidation stress profiles[J]. AIMS Geosciences, 2021, 7(4): 553-573. doi: 10.3934/geosci.2021032

Related Papers:

[1]	Qi Luo, Ryan Weightman, Sean T. McQuade, Mateo Díaz, Emmanuel Trélat, William Barbour, Dan Work, Samitha Samaranayake, Benedetto Piccoli . Optimization of vaccination for COVID-19 in the midst of a pandemic. Networks and Heterogeneous Media, 2022, 17(3): 443-466. doi: 10.3934/nhm.2022016
[2]	Prateek Kunwar, Oleksandr Markovichenko, Monique Chyba, Yuriy Mileyko, Alice Koniges, Thomas Lee . A study of computational and conceptual complexities of compartment and agent based models. Networks and Heterogeneous Media, 2022, 17(3): 359-384. doi: 10.3934/nhm.2022011
[3]	Xia Li, Chuntian Wang, Hao Li, Andrea L. Bertozzi . A martingale formulation for stochastic compartmental susceptible-infected-recovered (SIR) models to analyze finite size effects in COVID-19 case studies. Networks and Heterogeneous Media, 2022, 17(3): 311-331. doi: 10.3934/nhm.2022009
[4]	Eiman, Kamal Shah, Muhammad Sarwar, Thabet Abdeljawad . On rotavirus infectious disease model using piecewise modified $ABC$ fractional order derivative. Networks and Heterogeneous Media, 2024, 19(1): 214-234. doi: 10.3934/nhm.2024010
[5]	Ryan Weightman, Temitope Akinode, Benedetto Piccoli . Optimal control of pandemics via a sociodemographic model of non-pharmaceutical interventions. Networks and Heterogeneous Media, 2024, 19(2): 500-525. doi: 10.3934/nhm.2024022
[6]	Guillaume Cantin, Cristiana J. Silva, Arnaud Banos . Mathematical analysis of a hybrid model: Impacts of individual behaviors on the spreading of an epidemic. Networks and Heterogeneous Media, 2022, 17(3): 333-357. doi: 10.3934/nhm.2022010
[7]	Xavier Blanc, Claude Le Bris, Frédéric Legoll, Tony Lelièvre . Beyond multiscale and multiphysics: Multimaths for model coupling. Networks and Heterogeneous Media, 2010, 5(3): 423-460. doi: 10.3934/nhm.2010.5.423
[8]	Richard Carney, Monique Chyba, Taylor Klotz . Using hybrid automata to model mitigation of global disease spread via travel restriction. Networks and Heterogeneous Media, 2024, 19(1): 324-354. doi: 10.3934/nhm.2024015
[9]	Yimamu Maimaiti, Zunyou Lv, Ahmadjan Muhammadhaji, Wang Zhang . Analyzing vegetation pattern formation through a time-ordered fractional vegetation-sand model: A spatiotemporal dynamic approach. Networks and Heterogeneous Media, 2024, 19(3): 1286-1308. doi: 10.3934/nhm.2024055
[10]	Xia Li, Andrea L. Bertozzi, P. Jeffrey Brantingham, Yevgeniy Vorobeychik . Optimal policy for control of epidemics with constrained time intervals and region-based interactions. Networks and Heterogeneous Media, 2024, 19(2): 867-886. doi: 10.3934/nhm.2024039

Abstract

1. Introduction

In this paper we propose an epidemiological compartmental model where the efficacy time of vaccinations, i.e., the time a dosed individual needs to become immune, plays a key role. Furthermore, we also account for the concurrent effect of different vaccines, differing, for instance, in the time they need to provide immunization. The role of age can also be accounted for, letting vaccines' efficacy and immunization times depend on age, as well as on the choice of the vaccine type. These features are well known to be relevant in the present Covid-19 pandemic.

Compartmental models are a formidable tool in the description of a variety of real situations. The techniques above suggest a general framework able to introduce a specific dynamic evolution within compartments. In these models, each individual is considered to be of one compartment at any given time, its evolution consisting in passing from one compartment to another one, based on the structure and on the parameters of the model. In other words, the global dynamics consist in individuals entering the system (e.g., newborns), others leaving the system (e.g., casualties) and, during the evolution, passing from one compartment to another (e.g., getting ill, being vaccinated, recovering, $\ldots$ ). Intra-compartmental dynamics allows both to group the movements among specific compartments into a smooth evolution, and to detail the evolution within given compartments.

In the context of epidemiological models describing the spread of infectious diseases, the SIR model, named after its three compartments (Susceptible, Infected and Recovered), is the traditional prototype. This epidemic model dates back to 1927, see [16], less than a decade after the 1918 influenza pandemic. In what is probably its simplest form [22,§10.2], a prototype SIR model reads

$\begin{equation} \begin{array}{c@{\, }c@{\, }c@{\, }c@{\, }c} \boxed{ S } & \to & \boxed{ I } & \to & \boxed{ R } \\ & & \downarrow \end{array} \qquad\qquad \left\{ \begin{array}{rcl} \dot S & = & - \rho \, I \, S \\ \dot I & = & \rho \, I \, S - ( \vartheta + \mu) I \\ \dot R & = & \vartheta\, I \end{array} \right. \end{equation}$

(1)

where $t$ is time, $\rho$ and $\vartheta$ describe the transmission of the disease and the speed of recovery. Mortality of infected individuals is measured by $\mu$ .

The popularity of (1) is due partly to its simplicity, which allows to describe the disease behavior by estimating a small number of parameters, and partly to its being amenable to a variety of extensions. The dynamics of the SIR epidemic model are widely considered, here we recall for instance [1,4,21,25], or [14,Chapter 6], [22,Chapter 10], and [23,§1.5.1].

The literature on SIR-type models, typically containing additional compartments, is indeed huge. For instance, SEIR and SEIRS-type models contain also the compartment E (Exposed) where infected individuals spend a latent/incubation period prior to become infective, see [12,13,18,19] and the references therein. Vaccination campaigns and effects are described in SVIR-type models which consider individuals that get vaccinated, see [9,17,20] or also the different approach in [5].

Recall also, for instance, the SIHR model proposed in [6], where infected individuals are either Infective (I) or Hospitalized (H), the former ones spreading the disease, while the latter ones being isolated, typically hospitalized or in quarantine, thus taking into account lockdown effects.

Since 2020, with the Covid-19 pandemic, the development and use of these models boomed. The need to explicitly introduce the immunization time of a vaccine and the concurrent use of different vaccines lead us to introduce a dynamics within compartments. This machinery allows to fix a priori the time an individual spends in a compartment, as in the case of the immunization time of a vaccine, see §2, allowing also this time to be age dependent, as in §3. This general framework is sufficiently flexible to account for the simultaneous adoption of different vaccines, as in model (4).

In other instances, it might be appropriate to smoothen the change of status related to the change of compartment, an example being the SEIR model as modified in §4. Indeed, pass from (15) to (16) allows to account for a somewhat continuous evolution from exposed to infective. In all these examples, key statistics, such as the basic reproduction number, are naturally extended to these new frameworks keeping their original meaning.

Worth mentioning is the recent SIDARTHE model from [11], consisting of $8$ compartments. Indeed, while simplicity is a peculiarity of the SIR model, it is also a limit. The need to describe more complex dynamics often leads to more elaborate models. The framework we present below, at the expense of a few analytic technicalities, also allows to regroup similar compartments together obtaining a conceptual simplification, as in the extension to the SIDARTHE model proposed in §5.

In §6, the SIHR model is modified to comprehend spatial (geographic) movements. The introduction of the age structure, where advisable, is in general possible, as shown in §7. These two steps, namely the introduction of space and age structures, are indeed doable in all the examples discussed.

Remark that the present framework indeed provides an extension to the current compartmental modeling habits. In fact, suitable choices of various functions or parameters allow to trivialize intra–compartmental dynamics, thus recovering known models.

Pandemic dynamics can clearly be described through many other tools. For instance, the recent work [2] proposes a multiscale approach and discusses the current kinetic literature on the subject. Cellular automata are used for instance in [10] in a probabilistic setting, still with a compartmental structure, see [10,Figure 1 (b)]. The literature also offers essays correlating pandemics to other human activities, primarily to economy, see for instance [8].

Figure 1.

Solutions to (3)–(7)–(8) in the $4$ cases $T_* = 7, \, 21, \, 35, \, 49$

DownLoad: Full-Size Img PowerPoint

In the effective application of the models discussed below, a key issue is parameters' estimation. We defer for instance to [24] and to [27] for two entirely different approaches.

From the analytic point of view, a rigorous abstract approach to these classes of models is possible but at the cost of a quite intricate PDE based formalism, see a first attempt in [7]. Therefore, the sequel is devoted to show that intra-compartmental dynamics can both refine usual models, and also take into account new features, not captured by standard compartmental models.

2. Vaccination needs time $T_*$ to become effective

The model presented in [20], here slightly modified, amounts to this extension of (1):

$\begin{equation} \begin{array}{ccccc} & & \boxed{ V } \\ & \nearrow & \downarrow & \searrow \\ \boxed{ S } & \to & \boxed{ I } & \to & \boxed{ R } \\ & & \downarrow \end{array} \qquad\qquad \left\{ \begin{array}{rcl} \dot S & = & - \rho_S \, I \, S - p (t, S) \\ \dot V & = & p (t, S) - \rho_V \, I \, V - \vartheta_V \, V \\ \dot I & = & (\rho_S \, S + \rho_V \, V) I - ( \vartheta_I + \mu) I \\ \dot R & = & \vartheta_I\, I + \vartheta_V \, V \, . \end{array} \right. \end{equation}$

(2)

Here, $V$ is the number of vaccinated individuals, $p (t, S)$ measures the speed with which vaccinations of at most $S$ individuals take place at time $t$ , while $\mu$ is the mortality rate, $\rho_S, \rho_V$ measure infectivity, $\vartheta_I$ and $\vartheta_V$ are the recovery rates.

Before introducing the immunization time in (2), we remark that setting $p (t, S) = \alpha\, S$ and introducing the usual mortality terms for the $S$ , $V$ and $R$ populations, we recover precisely the case considered in [20].

The effect of doses, as is well known, is not immediate. On the contrary, the time $T_*$ that dosed individuals need to wait to get immunization plays a key role. Whenever a vaccine needs two (or more) doses, $T_*$ measures the time from the first dose to full immunization. Therefore, we propose the following model:

$\begin{equation} \begin{array}{@{}c@{\, }c@{\, }c@{\, }c@{\, }c@{}} & & \boxed{\begin{array}{*{20}{c}} V (0) \to V (T_*) \end{array}} \\ & \nearrow & \downarrow & \searrow \\ \boxed{ S } & & \to\hfill\boxed{ I }\hfill\to & & \boxed{ R } \\ & & \downarrow \end{array} \; \left\{ \begin{array}{r@{\;}c@{\;}l@{}} \dot S & = & - \rho_S \, I \, S - p (t, S) \\ \partial_t V + \partial_\tau V & = & - \rho_V \, I \, V \\ \dot I & = & (\rho_S \, S + \int_0^{T_*} \! \rho_V \, V) I - \vartheta \, I - \mu \, I \\ \dot R & = & \vartheta \, I + V (t, T_*) \\ V (t, 0) & = & p\left(t, S (t)\right) \, . \end{array} \right. \end{equation}$

(3)

Here, the number of individuals at time $t$ that were dosed at time $t-\tau$ is $V = V (t, \tau)$ , defined for $t \in {\mathbb{R}}_+$ and $\tau \in [0, T_*]$ . In other words, $\tau$ is the time since (the first) vaccination occurred. Note moreover that it is very reasonable to assume that $\rho_V$ is $\tau$ dependent, since $\rho_V(\tau)$ measures the ease with which vaccinated but not yet immunized individuals may get infected after time $\tau$ from the (first) dose. Note that also in the former model (2) the term $\rho_V \, I \, V$ is present and it implies that vaccinated individuals actually can get infected, but therein independently of the time from vaccination.

In connection with the present Covid-19 pandemic, we remark that model (3) can be easily extended to the case of different, say $m$ , vaccines having different efficacy times $T_*^1, \ldots, T_*^m$ . Indeed, denote by $V_\ell$ , for $\ell = 1, \ldots, m$ , the number of individuals dosed with the $\ell$ –th vaccine. We have

$\begin{equation} \begin{array}{@{}c@{\, }c@{\, }c@{\, }c@{\, }c@{}} & & \boxed{ \begin{array}{c} V_\ell (0) \to V_\ell (T_*^\ell) \\ \ell = 1, \ldots, m \end{array} } \\ & \nearrow & \downarrow & \searrow \\ \boxed{ S } & & \to\hfill\boxed{ I }\hfill\to & & \boxed{ R } \\ & & \downarrow \end{array} \left\{ \begin{array}{l} \dot S = - \rho_S \, I \, S - \sum\limits_{\ell = 1}^m p_\ell (t, S) \\ \partial_t V_\ell + \partial_\tau V_\ell = - \rho_{V_\ell} \, I \, V_\ell \\ \dot I = \left( \rho_S \, S + \sum\limits_{\ell = 1}^m \int_0^{T^\ell_*} \! \rho_{V_\ell} \, V_\ell \right) I - ( \vartheta + \mu) I \\ \dot R = \vartheta \, I + \sum\limits_{\ell = 1}^m V_\ell (t, T^\ell_*) \\ V_\ell (t, 0) = p_\ell \left(t, S (t)\right) \, . \end{array} \right. \end{equation}$

(4)

Note that, as is to be expected, according to (4) the total number of individuals is affected only by mortality, in the sense that

$\dfrac{\mathinner{\mathrm{d}{{\; }}}}{\mathinner{\mathrm{d}{{t}}}} \left( S (t) + \sum\limits_{\ell = 1}^m \int_0^{T_*^\ell} V_\ell (t, \tau) \, \mathinner{\mathrm{d}{\tau}} + I (t) + R (t) \right) = - \mu \, I (t) \, .$

A parameter often used to describe the tendency of the pandemic is the (time dependent) basic reproduction number $\mathcal{R}_o (t)$ . Within (4), it is explicitly computed:

$\begin{equation} \mathcal{R}_o (t) = \dfrac{ \rho_S \, S (t) + \sum\limits_{\ell = 1}^m \int_0^{T_*^\ell} \rho_{V_\ell} (\tau) \, V_\ell (t, \tau) \mathinner{\mathrm{d}{\tau}}} { \vartheta+\mu} \, , \end{equation}$

(5)

so that the increasing in the infected individuals at time $t$ is equivalent to $\mathcal{R}_o (t)>1$ :

$\begin{equation} \dfrac{\mathinner{\mathrm{d}{{\; }}}}{\mathinner{\mathrm{d}{{t}}}} I (t) \gtrless 0 \quad \iff \quad \mathcal{R}_o (t) \gtrless 1 \, . \end{equation}$

(6)

The above equivalence holds also in the case of the parameters $\rho_S$ , $\rho_{V_\ell}$ , $\vartheta$ and $\mu$ being time dependent.

At first sight, the numerator in the right hand side of (5) is misleading, for it apparently implies that higher values of $V_\ell$ lead to an increase in $\mathcal{R}_o (t)$ . On the contrary, increasing $V_\ell$ causes a decrease in $S$ and, even more important, individuals remain in the $V_\ell$ compartment only for a finite time $T_*^\ell$ , then passing to the $R$ compartment and thus reducing the value of $\mathcal{R}_o$ . Moreover, it is reasonable to assume that $\rho_S \gg \rho_{V_\ell}$ .

A direct comparison between (2) and (3), though suggestive, is inevitably highly arbitrary. Clearly, $V (t)$ in (2) has to be compared to $\int_0^{T_*} V (t, \tau) \mathinner{\mathrm{d}{\tau}}$ in (3). However the key mechanism with which $V$ individuals enter the $R$ population in (2) can hardly be related to that transforming $V$ into $R$ in (3). In the former case, at any time $t$ , the prescribed portion $\vartheta_V V (t)$ of the $V$ population enters $R$ . On the contrary, in (3), at time $t$ all $V (t, T_*)$ individuals after time $T_*$ from vaccination enter $R$ . Formally, this difference is evident noting that the parameter $\vartheta_V$ in (2) has no counterpart in (3) while, on the other hand, the time $T_*$ in (3) has no analogue in (2).

Thus, below, to investigate the role of $T_*$ , we compare the evolutions of instances of (3) differing only in the time necessary for vaccination to ensure immunization. Referring to (3), we choose the following parameters¹

¹ $\chi_{ {A}}$ is the characteristic function of the set $A$ : $\chi_{ {A}} (x) = 1 \iff x \in A$ and $\chi_{ {A}} (x) = 0 \iff x \not\in A$ .

$\begin{equation} \begin{array}{rcl} \rho_S & = & 5\times 10^{-4} \\ \vartheta & = & 1\times 10^{-3} \\ \mu & = & 1\times 10^{-4} \end{array} \qquad \begin{array}{rcl} \rho_V (\tau) & = & \rho_S \, \sqrt{1-\dfrac{\tau}{T_*}} \\ p (t, S) & = & 0.5 \; \chi_{ {[10, +\infty[}} (t) \;\; \chi_{ { {\mathbb{R}}_+}} (S) \end{array} \end{equation}$

(7)

and initial data

$\begin{equation} S_o = 99.9 \, , \quad V_o (\tau) \equiv 0 \, , \quad I_o = 0.1 \, , \quad R_o = 0 \, . \end{equation}$

(8)

We consider $8$ sample choices for the time required by vaccination to become effective and compute the corresponding number of deaths, the results are in Table 1. Clearly, a shorter value of $T_*$ is preferable. We provide in Figure 1 diagrams of the solutions corresponding to $T_* = 7, \, 21, \, 35, \, 49$ . Note first that the dynamics of the $S$ compartment is hardly affected by the changes in $T_*$ , the key differences being in the $I$ and $R$ populations.

Table 1.

Times necessary for the vaccination to provide immunity and corresponding casualties according to model (3)–(7)–(8). The initial total population is $100$

$T_*$ (days)	1	7	14	21	28	35	42	49
Deaths:	0.28	0.32	0.37	0.43	0.49	0.56	0.63	0.70

| Show Table

DownLoad: CSV

The relevant differences in the diagram of the map $t \to \int_0^{T_*} V (t, \tau)\mathinner{\mathrm{d}{\tau}}$ heavily depend on the different lengths of the time interval where $V$ is integrated. Moreover, higher values of $T_*$ results in more members of the $V$ population get infected, coherently with (3). This explains the rather horizontal plateau in the $\int V$ diagram for $T_* = 7$ , which is replaced by a clearly decreasing profile when $T_* = 35$ and $T_* = 49$ .

In Figure 1, the differences in the $I$ diagrams are evident and clear: higher values of $T_*$ enhance the spreading of the virus, causing a higher number of infected individuals and, hence, of casualties.

The evolutions in the $R$ populations are similar. Note however that they are slightly translated one with respect to the other, due to the delay with which dosed individuals are immunized and, hence, enter the $R$ population according to (3). It is mostly this delay that explains the differences in the casualties, as reported in Table 1.

The flexibility introduced by dynamics internal to the dosed ( $V$ ) population allows also to take into consideration variations in this evolution. Assume, for instance, that vaccines are unavailable between times, say, $30$ and $120$ , so that we replace $p$ in (7) with

$\begin{equation} p (t, S) = 0.5 \; \chi_{ {[10, 30] \cup [120, +\infty[}} (t) \;\; \chi_{ { {\mathbb{R}}_+}} (S) \, . \end{equation}$

(9)

Then, the number of casualties sharply grows with respect to the number of casualties in (3)–(7)–(8), as clearly shown when comparing Table 1 with Table 2.

Table 2.

Times necessary for the vaccination to provide immunity and corresponding casualties, according to model (3)–(7)–(8), in the case vaccinations are suspended as detailed in (9). The initial total population is $100$

$T_*$	1	7	14	21	28	35	42	49
Deaths:	1.11	1.18	1.25	1.32	1.38	1.43	1.48	1.53

| Show Table

DownLoad: CSV

The diagrams in Figure 2 also confirm the negative effect in a suspension of the vaccination campaign.

Figure 2.

Diagrams of the solutions to (3)–(7)–(8) with a suspension in the vaccination campaign as detailed in (9) in the $4$ cases $T_* = 7, \, 21, \, 35, \, 49$

DownLoad: Full-Size Img PowerPoint

As a further example, we consider the problem of choosing between two different vaccines. We integrate (4) with $m = 2$ and supplement the data and parameters in (7)–(8)

$\begin{equation} \begin{array}{@{}r@{\, }c@{\, }l@{\quad}r@{\, }c@{\, }l@{\quad}r@{\, }c@{\, }l@{}} T_*^1 & = & 7 & \rho_{V_1} (\tau) & = & \rho_S \sqrt{1-\dfrac{\tau}{T_*^1}} & p_1 (t, S) & = & \omega \, \chi_{ {\mathopen[20, +\infty\mathclose[}} (t) \; \chi_{ { {\mathbb{R}}_+}} (S) \\ T_*^2 & = & 35 & \rho_{V_2} (\tau) & = & 0.2 \; \rho_S \sqrt{1-\dfrac{\tau}{T_*^2}} & p_2 (t, S) & = & (0.5 - \omega) \, \chi_{ {\mathopen[20, +\infty\mathclose[}} (t) \; \chi_{ { {\mathbb{R}}_+}} (S) \, , \end{array} \end{equation}$

(10)

with the parameter $\omega \in [0, \, 0.5]$ determining the amount of vaccines of the two types dosed at each time. These choices correspond to the situation where two different vaccines are available. The first one is effective after only $T_*^1 = 7$ days, while the latter needs $T_*^2 = 35$ days to provide full immunization. On the other hand, the second vaccine provides a far better protection immediately after the first shot, as it follows from $\rho_{V_1} = 5 \, \rho_{V_2}$ in (10). The resulting integrations in Figure 3 show that the two rather different strategies resulting from $\omega = 0.1$ and $\omega = 0.4$ actually lead to similar results for what concerns the number of casualties. In the two cases, however, the global costs of the two campaigns might be significantly different, due to the different costs of the vaccines and of their distribution.

Figure 3.

Diagrams of the solutions to (4)–(7)–(8)–(10). On the left with $\omega = 0.1$ and, on the right, with $\omega = 0.4$

DownLoad: Full-Size Img PowerPoint

3. Intra–compartmental dynamics with age structure

The insertion of an intra–compartmental dynamics is not limited to ODE models. As it is well known, age differences often have a role in the spreading of diseases. Therefore, we introduce the following extension of (3):

$\begin{equation} \left\{ \begin{array}{@{\, }r@{\;}c@{\;}l@{}} \partial_t S + \partial_a S & = & - \int_{ {\mathbb{R}}_+} \rho (t, a, \alpha) \, I (t, \alpha) \mathinner{\mathrm{d}{\alpha}} \, S - p (t, a, S) \\ \partial_t V + \partial_a V + \partial_\tau V & = & - \rho_V \, I \, V \\ \partial_t I + \partial_a I & = & \int_{ {\mathbb{R}}_+} \rho (t, a, \alpha) \, I (t, \alpha) \mathinner{\mathrm{d}{\alpha}} S \\ & & + \int_{ {\mathbb{R}}_+} \int_0^{T_*} \rho_V (t, a, \tau, \alpha) \, V (t, \alpha, \tau)\mathinner{\mathrm{d}{\alpha}} \mathinner{\mathrm{d}{\tau}} I - \vartheta \, I - \mu \, I \\ \partial_t R + \partial_a R & = & \vartheta \, I + V (t, a, T_*) \, , \end{array} \right. \end{equation}$

(11)

Vaccination entering also as boundary datum in the evolution of the $V$ population:

$\begin{equation} V (t, a, 0) = p\left(t, a, S (t, a)\right) \, . \end{equation}$

(12)

In (11)–(12), we use the same symbols as in (3), allowing the various parameters to depend also on age. For instance $\rho = \rho (t, a, \alpha)$ describes the ease with which $I$ individuals of age $\alpha$ infect members of the $S$ population of age $a$ at time $t$ , and similarly for $\rho_V$ .

System (11) needs to be supplied with initial and boundary data, such as

$\begin{equation} \left\{ \begin{array}{rcl} S (0, a) & = & S_o (a) \\ V (0, a, \tau) & = & V_o (a, \tau) \\ I (0, a) & = & I_o (a) \\ R (0, a) & = & R_o (a) \end{array} \right. \qquad \left\{ \begin{array}{rcl} S (t, 0) & = & S_b (t) \\ V (t, 0, \tau) & = & 0 \\ I (t, 0) & = & 0 \\ R (t, 0) & = & 0 \end{array} \right. \end{equation}$

(13)

where, as it is realistic, we set to $0$ the birth rate of vaccinated, infected and recovered individuals.

Note that setting all data and parameters in (11)–(12)–(13) constant in $a$ , we essentially reobtain the previous model (3), essentially meaning that $S, V, I, R$ in (3) have to be replaced by their integrals over $a \in {\mathbb{R}}_+$ .

Extending the analogous property valid for system (3), note that the total number of individuals varies not only due to mortality but also due to the natality $S_b (t)$ , indeed

$\begin{eqnarray*} \dfrac{\mathinner{\mathrm{d}{{\; }}}}{\mathinner{\mathrm{d}{{t}}}} \! \int_{ {\mathbb{R}}_+} \!\! \left( S (t, a) + \int_0^{T_*} V (t, a, \tau) \mathinner{\mathrm{d}{\tau}} + I (t, a) + R (t, a) \right) \! \mathinner{\mathrm{d}{{a}}} = {-}\!\! \int_{ {\mathbb{R}}_+} \!\! \mu (t, a) \, I (t, a) \mathinner{\mathrm{d}{{a}}} + S_b (t) \, . \end{eqnarray*}$

The present model realistically allows infection to be propagated across different ages. Therefore, we can only have an index $\mathcal{R}_o (t)$ averaged over all ages. Indeed, defining

$\mathcal{R}_o (t) = \dfrac{ \int_{ {\mathbb{R}}_+} \int_{ {\mathbb{R}}_+} \! \left( \rho (t, \alpha, a) \, S (t, \alpha) + \int_0^{T_*} \rho_V (t, a, \tau, \alpha) \, V (t, \alpha, \tau)\mathinner{\mathrm{d}{\tau}} \right) \mathinner{\mathrm{d}{\alpha}} I (t, a) \mathinner{\mathrm{d}{{a}}} }{ \int_{ {\mathbb{R}}_+} \left( \vartheta (t, a) + \mu (t, a)\right)I (t, a) \, \mathinner{\mathrm{d}{{a}}}}$

we have the following extension of (6):

$\begin{equation} \dfrac{\mathinner{\mathrm{d}{{\; }}}}{\mathinner{\mathrm{d}{{t}}}} \int_{ {\mathbb{R}}_+} I (t, a) \mathinner{\mathrm{d}{{a}}} \gtrless 0 \iff \mathcal{R}_o (t) \gtrless 1 \, . \end{equation}$

(14)

Remark that also model (11) can easily be extended to the case of different vaccines having different full efficacy times through the introduction of different vaccinated populations $V_1, \ldots, V_M$ , now with $V_\ell = V_\ell (t, a, \tau)$ , for $\ell = 1, \ldots, m$ . The corresponding extensions of the above {formul\ae} for the variation of the total population and for the definition of $\mathcal{R}_o (t)$ are straightforward.

The use of model (11)–(12)–(13) allows to tackle the key issue of optimizing the use of the available vaccines on the basis of their different efficacy on different age classes.

4. Exposed are not immediately infective

The description provided by (1) is often too approximate. For instance, it might be necessary to distinguish between exposed ( $E$ ) and infective ( $I$ ) individuals, extending (1) to

$\begin{equation} \begin{array}{c@{\, }c@{\, }c@{\, }c@{\, }c@{\, }c@{\, }c@{\, }c} \boxed{ S } & \to & \boxed{ E } & \to & \boxed{ I } & \to & \boxed{ R } \\ & & \downarrow & & \downarrow \end{array} \qquad \qquad \left\{ \begin{array}{rcl} \dot S & = & - \rho \, I \, S \\ \dot E & = & \rho \, I \, S - ( \vartheta_E + \mu_E + \kappa) E \\ \dot I & = & \kappa\, E - ( \vartheta_I + \mu_I) I \\ \dot R & = & \vartheta_E \, E + \vartheta_I \, I \, , \end{array} \right. \end{equation}$

(15)

where we introduced the term $\kappa\, E$ representing the speed at which exposed individuals become infective. Refer for instance to [12,13,18] for various results about developments of SEIR models.

Assume now that individuals infected at time $t$ typically become infective at a later time. It is then natural, see also [15,16], to introduce the number $I (t, \tau)$ of individuals at time $t$ that were infected at time $t-\tau$ , so that (1) becomes

$\begin{equation} \begin{array}{c@{\, }c@{\, }c@{\, }c@{\, }c} \boxed{ S } & \to & \boxed{ I (0) \to I (\tau) \to } & \to & \boxed{ R } \\ & & \downarrow \end{array} \quad \left\{ \begin{array}{@{\, }rcl@{}} \dot S & = & - \int_{ {\mathbb{R}}_+} \rho (\tau) \, I (t, \tau) \, \mathinner{\mathrm{d}{\tau}} \, S \\ \partial_t I + \partial_\tau I & = & - \vartheta (\tau)\, I - \mu (\tau) \, I \\ \dot R & = & \int \vartheta (\tau) \, I (t, \tau) \, \mathinner{\mathrm{d}{\tau}} \\ I (t, 0) & = & \int_{ {\mathbb{R}}_+} \rho (\tau) \, I (t, \tau) \, \mathinner{\mathrm{d}{\tau}} \, S (t) . \end{array} \right. \end{equation}$

(16)

Observe that the last line above is in fact a boundary condition, prescribing how many individuals get infected at time $t$ .

Remark that the dependence of $\rho$ , $\vartheta$ and $\mu$ on $\tau$ allows to account for different infectivity, recovery and mortality rate at different stages of the infection. For instance, if an individual infected at time $t$ becomes infective at time $t+t_i$ , then we set $\rho (\tau) = 0$ for $\tau \in [0, t_i]$ .

Observe that, as is to be expected, the variation in the total number of individuals resulting from (16) is due exclusively to mortality. Indeed

$\begin{equation} \dfrac{\mathinner{\mathrm{d}{{\; }}}}{\mathinner{\mathrm{d}{{t}}}} \left( S (t) + \int_{ {\mathbb{R}}_+}I (t, \tau) \, \mathinner{\mathrm{d}{\tau}} + R (t) \right) = - \int_{ {\mathbb{R}}_+} \mu (\tau) \, I (t, \tau) \, \mathinner{\mathrm{d}{\tau}} \, . \end{equation}$

(17)

Within the framework of (16), an explicit expression for the basic reproduction number $\mathcal{R}_o (t)$ at time $t$ is available:

$\mathcal{R}_o (t) = \dfrac{\int_{ {\mathbb{R}}_+} \rho (\tau) \, I (t, \tau) \mathinner{\mathrm{d}{\tau}} \, S (t)} {\int_{ {\mathbb{R}}_+}\left( \vartheta (\tau) + \mu (\tau)\right) I (t, \tau) \mathinner{\mathrm{d}{\tau}}}$

and it is immediate to see that its exceeding unity is equivalent to the increase of the total number of infected individuals:

$\begin{equation} \dfrac{\mathinner{\mathrm{d}{{\; }}}}{\mathinner{\mathrm{d}{{t}}}} \left(\int_{ {\mathbb{R}}_+} I (t, \tau) \mathinner{\mathrm{d}{\tau}}\right) \gtrless 0 \quad \iff \quad \mathcal{R}_o (t) \gtrless 1 \, . \end{equation}$

(18)

It goes without saying that suitable choices of initial data and parameters in (16) allow to recover within (16) the solutions of (1). Here we only provide a quick example of a comparison among solutions to (1), (15) and (16), with data and parameters in (1) being an average of those in (15), while the choices in (16) are an interpolation of those in (15), as detailed in (19).

$\begin{equation} \begin{array}{r@{}c@{}l@{\qquad\qquad}r@{}c@{}l@{\qquad\qquad}r@{}c@{}l@{}} & (1) & & & (15) & & & (16) & \\ \rho & = & 0.1 & \rho & = & 0.1 & \rho (\tau) & = & 0.1 \, (1-e^{-0.2\, \tau}) \\ \vartheta & = & 3 & ( \vartheta_E, \vartheta_I) & = & (7, 1) & \vartheta (\tau) & = & \begin{cases} 7 & \tau < 2 \\ 1 & \tau \geq 2 \end{cases} \\ \mu & = & 0.5 &(\mu_E, \mu_I) & = & (0.1, 0.9) & \mu (\tau) & = & 0.9 - 0.8 \, e^{-2\tau} \\ & & & \kappa & = & 5 \\ S_o & = & 10 & S_o & = & 10 & S_o & = & 10 \\ I_o & = & 6. & (E_o, I_o) & = & (1, 5) & I_o (\tau) & = & \begin{cases} 0.5 & \tau \in [0, 2[ \\ 2.5 & \tau \in [2, 4] \end{cases} \\ R_o & = & 0 & R_o & = & 0 & R_o & = & 0 \, . \end{array} \end{equation}$

(19)

Figure 4 shows sample integrations of (1), (15) and (16). In the situation considered, according to mortality, the mixed ODE–PDE system (16) can be seen in some senses in the middle between the other two classical purely ODE based compartmental models (1) and (15). Intra-compartmental dynamics can thus recover the dynamics described through standard compartmental models.

Figure 4.

Above, the integrations of (1) and (15), below on the left that of (16) (19). The rightmost diagram on the second line displays the total number of living individuals in the three cases, showing that, with respect to mortality, the ODE–PDE model (16) can be seen in some senses in the middle between the ODE models (1) and (15)

DownLoad: Full-Size Img PowerPoint

On the other hand, other scenarios can hardly be recovered through the use of only standard compartments. Indeed, for instance, it can be reasonable to assume that infected individuals are most infective in given time intervals after infection where, in general, these intervals depend on the specific disease under consideration. Such a situation, within Model (16), is easily described and we consider, for instance, the following choices:

$\begin{equation} \rho (\tau) = \bar\rho \; t \; e^{- \alpha (t-\bar T)^2} \mbox{ with } \rho = 0.1 \mbox{ and } \begin{array}{@{\mbox{Case }}lr@{\, }c@{\, }l@{\qquad}r@{\, }c@{\, }l} (i) & \alpha & = & 1.0 \, , & \bar T & = & 10.0 \, , \\ (ii) & \alpha & = & 1.0 \, , & \bar T & = & 4.0 \, , \\ (iii) & \alpha & = & 0.25 \, , & \bar T & = & 4.0 \, . \end{array} \end{equation}$

(20)

The sample choices (20) differ in the time $\bar T$ at which $I$ individual are most infective and in the coefficient $\alpha$ ruling the width of the infective time interval.

The results are summarized in Figure 5. The differences in the evolutions prescribed by (16) corresponding to the choices (20) are very intuitive, Case $(ii)$ being the reference situation. In Case $(i)$ , infected individuals mostly either die or recover before being infective. In Case $(iii)$ , the wider distribution of $\rho$ helps the spreading of the infection.

Figure 5.

Above, from left to right, the integrations of Case $(i)$ , Case $(ii)$ and Case $(iii)$ in (20) with parameters and data as prescribed in (19). Below, the corresponding choices of the $\rho$ function as detailed in (20). The differences in the displayed evolutions are due to the intra–compartmental dynamics in the $I$ population

DownLoad: Full-Size Img PowerPoint

5. Intra-compartmental dynamics in the SIDARTHE model

We now consider the SIDARTHE model introduced in [11]. It consists of $8$ populations, which we briefly recall in Table 3, leaving the detailed explanation to [11].

Table 3.

Populations in model (21) from [11]

$S$	Susceptible	healthy	can be infected
$I$	Infected	asymptomatic	infective	undetected
$D$	Diagnosed	asymptomatic	infective	detected
$A$	Ailing	symptomatic	infective	undetected
$R$	Recognized	symptomatic	infective	detected
$T$	Threatened	acutely symptomatic	infected	detected
$H$	Healed	healthy	immune
$E$	Extinct

| Show Table

DownLoad: CSV

Its block diagram and the corresponding system of ordinary differential equations read:

$\begin{equation} \begin{array}{@{\!\!\!\!}cccccc@{\;}c@{}} \boxed{ S } & \to & \boxed{ I } & \to & \boxed{ D } \\ & \swarrow & \downarrow & & \downarrow & \searrow \\ \boxed{ H } & \leftarrow & \boxed{ A } & \to & \boxed{ R } & & \searrow\quad \\ \uparrow & & &\searrow & \downarrow & \searrow & \downarrow \\ \uparrow & & \boxed{ E } &\leftarrow & \boxed{ T } & \to & \downarrow \\ \uparrow & & & & & & \downarrow \\ \quad\leftarrow & \leftarrow & \leftarrow & \leftarrow & \leftarrow & \leftarrow & \;\leftarrow\quad \end{array} \left\{ \begin{array}{@{}r@{\;}c@{\;}l@{}} \dot S & = & - \alpha \, I \, S - \gamma \, A \, S - \beta \, D \, S - \delta \, R \, S \\ \dot I & = & \alpha \, I \, S + \gamma \, A \, S + \beta \, D \, S + \delta \, R \, S - ( {\varepsilon}+\zeta+\lambda)I \\ \dot D & = & {\varepsilon}\, I - (\eta + \rho)D \\ \dot A & = & \zeta\, I - ( \vartheta + \mu + \kappa)A \\ \dot R & = & \eta\, D + \vartheta\, A - (\nu + \xi) R \\ \dot T & = & \mu \, A + \nu \, R - (\sigma + \tau) T \\ \dot H & = & \lambda \, I + \rho \, D + \kappa \, A + \xi \, R + \sigma \, T \\ \dot E & = & \tau \, T \, . \end{array} \right. \end{equation}$

(21)

We refer to [11] for the specific meaning of each constant and for the motivations of this model.

Here, we note that both the evolutions $I\to A$ and $D \to R \to T$ are due to the development and worsening of the symptoms. It is therefore natural to single out these parts, namely

(22)

Introduce now the populations $\mathcal{I} = \mathcal{I} (t, s)$ and $\mathcal{R} = \mathcal{R} (t, s)$ . The former comprehends $I$ and $A$ , while the latter consists of $D$ , $R$ and $T$ , as it stems comparing the scheme (22) above with (23) below. The evolution of the symptoms is then described through the dependence on the $s$ variable, in completely independent ways in the two different compartments.

$\begin{equation} \begin{array}{ccccccc@{}} \boxed{ S } & \to & \boxed{ \mathcal{I} (0) \to \mathcal{I} (s) \to } & \to & \boxed{ \mathcal{R} (0) \to \mathcal{R} (s) \to } & \to & \boxed{ E } \\ & & \;\;\searrow & & \swarrow\;\; \\ & & \ \qquad\qquad \longrightarrow & \boxed{ H } & \longleftarrow \ \qquad\qquad \end{array} \end{equation}$

(23)

$\begin{equation} \left\{ \begin{array}{r@{\;}c@{\;}l@{}} \dot S & = & -\int \alpha(s) \, \mathcal{I} (t, s) \mathinner{\mathrm{d}{{s}}} \, S -\int \delta(s) \, \mathcal{R} (t, s) \mathinner{\mathrm{d}{{s}}} \, S \\ \partial_t \mathcal{I} \; + \;\partial_s \mathcal{I} & = & \int \alpha(s) \, \mathcal{I} (t, s) \mathinner{\mathrm{d}{{s}}} \, S + \int \delta(s) \, \mathcal{R} (t, s) \mathinner{\mathrm{d}{{s}}} \, S \\ & & \qquad - \int {\varepsilon} (s) \, \mathcal{I} (t, s) \mathinner{\mathrm{d}{{s}}} - \int \vartheta (s) \, \mathcal{I} (t, s) \mathinner{\mathrm{d}{{s}}} \\ \partial_t \mathcal{R} + \partial_s \mathcal{R} & = & \int \vartheta (s) \, \mathcal{I} (t, s) {s} - \int \eta (s)\, \mathcal{R} (t, s) \mathinner{\mathrm{d}{{s}}} \\ \dot H & = & \int {\varepsilon} (s) \, I \mathinner{\mathrm{d}{{s}}} \\ \dot E & = & \int \eta (s) \, \mathcal{R} (t, s) \mathinner{\mathrm{d}{{s}}} \, , \end{array} \right. \end{equation}$

(24)

All integrals being computed over ${\mathbb{R}}_+$ . Careful piecewise constant choices of the various functions $\alpha$ , $\delta$ , ${\varepsilon}$ , $\vartheta$ and $\eta$ allow to recover [11,Formulæ (1)–(8)]. However, system (24) is far more flexible, although it actually consists of only $3$ differential equations, since $H$ and $E$ are easily found once $S$ , $\mathcal{I}$ and $\mathcal{R}$ are available.

Note also that in model (24) the explicit expression of the basic reproduction number is immediately at hand. Indeed, defining

$\mathcal{R}_o (t) = \dfrac{\int \left(\alpha(s) \, \mathcal{I} (t, s) + \delta(s) \, \mathcal{R} (t, s) \right) \mathinner{\mathrm{d}{{s}}} \, S (t)} {\int \left( {\varepsilon} (s) + \vartheta (s)\right) \mathcal{I} (t, s) \mathinner{\mathrm{d}{{s}}}}$

the key property (18) still holds.

6. Movements as intra–compartmental dynamics

As a further example of intra-compartmental dynamics we consider movements. Indeed, as a prototype for the description of the role of quarantine in the spreading of a pandemic, consider the model [6,Formula (9)]

$\begin{equation} \begin{array}{c@{\, }c@{\, }c@{\, }c@{\, }c} \boxed{ S } & \to & \boxed{ I } & \searrow \\ & & \downarrow & & \boxed{ R } \\ & & \boxed{ H } & \nearrow \end{array} \qquad\quad \left\{ \begin{array}{r@{\, }c@{\, }l} \partial_t S + \mu_S \, S & = & - \rho \, I \, S \\ \partial_t I + \mu_I \, I & = & \rho \, I \, S - \kappa \, I - \vartheta \, I \\ \partial_t H + \mu_H \, H & = & \kappa\, I - \eta \, H \\ \partial_t R + \mu_R \, R & = & \vartheta \, I + \eta \, H \, , \end{array} \right. \end{equation}$

(25)

where we used the standard notation for the $S$ , $I$ and $R$ populations, while $H$ stands for the individuals that are infected but are hospitalized, or in quarantine. Here, $\rho$ describes, as usual, infectivity; $\kappa$ the hospitalization rate; $\eta$ , respectively $\vartheta$ , the recovery rate of hospitalized, respectively infected, individuals. The mortality terms $\mu_S\, S$ and $\mu_R\, R$ might as well be neglected on short time intervals.

In this connection, it can be relevant to take into consideration events that provoked relevant gatherings of crowds, a well known example being the Atalanta vs. Valencia football match played in Milan (Italy) on February 19th, 2020, see [26]. To this aim, following [6], we introduce a geographical movement as follows

$\begin{equation} \begin{array}{@{\!\!\!}c@{\, }c@{\, }c@{\, }c@{\, }c} \boxed{ S (x) \, \circlearrowleft } & \to & \boxed{ I (x) \, \circlearrowleft } & \searrow \\ & & \downarrow & & \boxed{ R (x) \, \circlearrowleft } \\ & & \boxed{ H } & \nearrow \end{array} \; \left\{ \begin{array}{@{}r@{\, }c@{\, }l@{}} \partial_t S + \mathinner{\mathop{{\rm{div}}}} (v_S \, S) +\mu_S \, S & = & - \rho \, I \, S \\ \partial_t I + \mathinner{\mathop{{\rm{div}}}} (v_I \, I) +\mu_I \, I & = & \rho \, I \, S - \kappa \, I - \vartheta \, I \\ \partial_t H + \mu_H \, H & = & \kappa\, I - \eta \, H \\ \partial_t R + \mathinner{\mathop{{\rm{div}}}} (v_R \, R) + \mu_R \, R & = & \vartheta \, I + \eta \, H \, . \end{array} \right. \end{equation}$

(26)

Above, all populations are also space dependent, so that, for instance, $S = S (t, x)$ , $x$ being the geographical coordinate in ${\mathbb{R}}^2$ . The various movements are described by the $2d$ vectors $v_S$ , $v_I$ and $v_R$ , while $H$ individuals are not assumed to be moving. We defer to [6,Section 6] for further details and for a sample integration of (26), while its well posedness is proven in [7].

7. An age dependent SIS model with immunization time

Introducing age structure in the SIS model [15], see also [3,Formula 10.1], we get

$\begin{equation} \boxed{ S } \; \longleftrightarrow \; \boxed{ I } \qquad\qquad\qquad \left\{ \begin{array}{r@{\, }c@{\, }l} \partial_t S + \partial_a S & = & -\rho \, I \, S + \gamma\, I \\ \partial_t I + \partial_a I & = & \rho \, I \, S - \gamma\, I \, . \end{array} \right. \end{equation}$

(27)

In this model, when infected individuals recover, they get back to being susceptible at a rate governed by the, possibly age dependent, parameter $\gamma$ . However, it might be reasonable to assume that, after infection, immunization lasts for a, possibly age dependent, time interval. Therefore, we are lead to introduce an $R$ population of individuals that recovered and are immune, but only for a predetermined age dependent time $T_* = T_* (a)$ :

$\begin{equation} \!\!\!\! \begin{array}{@{}rcl} \boxed{ S } \quad & \longrightarrow \qquad \longrightarrow & \quad \boxed{ I } \\ \nwarrow & &\swarrow \\ &\boxed{ R(0) \to R(T_*) } \end{array} \left\{ \begin{array}{@{\, }r@{\, }c@{\, }l@{}} \partial_t S + \partial_a S & = & -\rho \, I \, S + R \left(t, a, T_* (a)\right) \\ \partial_t I + \partial_a I & = & \rho \, I \, S - \gamma\, I \\ \partial_t R + \partial_a R + \partial_\tau R & = & 0 \\ R (t, a, 0) & = & \gamma \, I (t, a)\, . \end{array} \right. \end{equation}$

(28)

Note however that the very simple nature of the third equation allows to rewrite (28) as a system of $2$ PDEs with "delay" in both the $t$ and $a$ variables:

$\begin{equation} \left\{ \begin{array}{r@{\, }c@{\, }l} \partial_t S + \partial_a S & = & -\rho \, I \, S + \gamma \, I \left(t - T_* (a), a - T_* (a)\right) \\ \partial_t I + \partial_a I & = & \rho \, I \, S - \gamma\, I (t, a)\, . \end{array} \right. \end{equation}$

(29)

The form (28) is more prone to extensions than (29), since it allows to specify various evolutions of the $R$ population. Indeed, from the modeling point of view, systems (28) can easily be extended to comprise other factors, such as a $\tau$ –dependent rate of re-infection of the $R$ individuals, for instance.

From the analytic point of view, we remark that (28) provides a first minimal example of a "junction" in $2$ space dimensions. Indeed, the line $\tau = T_* (a)$ splits the plane of the "space" variables $(a, \tau)$ .

References

[1]	Agaiby S (2018) Advancements in the interpretation of seismic piezocone tests in clays and other geomaterials. School of Civil & Environmental Engineering, Georgia Institute of Technology, Atlanta, GA USA, 925.
[2]	Robertson PK (1990) Soil classification using the cone penetration test. Can Geotech J 27: 151-158. doi: 10.1139/t90-014
[3]	Lunne T, Robertson PK, Powell JJM (1997) Cone Penetration Testing in Geotechnical Practice, EF Spon/CRC Press, London, 352.
[4]	Eslami A, Fellenius BH (1997) Pile capacity by direct CPT and CPTU methods applied to 102 case histories. Can Geotech J 34: 886-904. doi: 10.1139/t97-056
[5]	Schneider JA, Hotstream JN, Mayne PW, et al. (2012) Comparing CPTU Q-F and Q-Δu₂/σ_vo' soil classification charts. Geotechnique Lett 2: 209-215. doi: 10.1680/geolett.12.00044
[6]	Jefferies M, Been K (2015) Soil Liquefaction: A Critical State Approach, Second Edition, Taylor & Francis Group, London, 712.
[7]	Robertson PK (2009) Interpretation of cone penetration tests—a unified approach. Can Geotech J 46: 1337-1345. doi: 10.1139/T09-065
[8]	Shahri AA, Malehmir A, Juhlin C (2015) Soil classification based on piezocone data. Eng Geo 189: 32-47. doi: 10.1016/j.enggeo.2015.01.022
[9]	Valsson SM (2016) Detecting quick clay with CPTU. 17th Nord Geotech Meet.
[10]	Sandven R, Gylland A, Montafia A, et al. (2016) In-situ detection of sensitive clays—Part I: selected test methods. 17th Nord Geotech Meet, Reykjavik.
[11]	Sandven R, Gylland A, Montafia A, et al. (2016) In-situ detection of sensitive clays—Part Ⅱ: Results. 17th Nord Geotech Meet. Reykjavik, Iceland: Icelandic Geotechnical Society.
[12]	Gylland AS, Sandven R, Montafia A, et al. (2017) CPTU classification diagrams for identification of sensitive clays. Landslides in Sensitive Clays, Springer Series on Advances in Natural & Technological Hazards Research, Cham, Switzerland, 57-66.
[13]	DeGroot DJ, Landon ME, Poirier SE (2019) Geology and engineering properties of sensitive Boston Blue Clay at Newbury, Massachusetts. AIMS Geosci 5: 412-447. doi: 10.3934/geosci.2019.3.412
[14]	Mayne PW, Benoît J (2020) Analytical CPTU Models Applied to Sensitive Clay at Dover, New Hampshire. J Geotech Geoenviron Eng 146: 04020130. doi: 10.1061/(ASCE)GT.1943-5606.0002378
[15]	Long M (2008) Design parameters from in-situ tests in soft ground. Geotechnical and Geophysical Site Characterization, Taylor & Francis, London, 90-116.
[16]	Coutinho RQ, Bello MI (2014) Geotechnical characterization of Suape soft clays, Brazil. Soils Rocks 37: 257-276.
[17]	Mlynarek Z, Wierzbicki J, Gogolik S, et al. (2014) Shear strength and deformation parameters of peat and gyttja from CPTu, SDMT, and VST tests, 5th Intl Workshop CPTu DMT Soft Clays Organic Soils, 193-209.
[18]	Nejaim PF, Jannuzzi GMF, Danziger FAB (2016) Soil behavior type of the Sarapuí Ⅱ test site. Geotechnical & Geophysical Site Characterization 5, Gold Coast, Australian Geomechanics Society, 1009-1014. Available from: www.usucger.org.
[19]	Zawrzykraj P, Rydelek P, Bąkowska A (2017) Geoengineering properties of Eemian peats from central Poland in the light of static cone penetration and dilatometer tests. Eng Geol 226: 290-300. doi: 10.1016/j.enggeo.2017.07.001
[20]	Mayne PW, Agaiby S (2019) Profiling yield tresses and identification of soft organic clays using piezocone tests, Proceedings XVI Pan American Conference on Soil Mechanics & Geotechnical Engineering, Paper 0149, Cancun, Mexican Society of Geotechnical Engineering (SMIG). Available from: www.issmge.org.
[21]	Mayne PW, Agaiby SS, Dasenbrock D (2020) Piezocone identification of organic clays and peats, GeoCongress 2020: Modeling, Geomaterials, and Site Characterization, (Minneapolis, GSP 317), ASCE, Reston, VA, 541-549.
[22]	Fellenius BH, Eslami A (2000) Soil profile interpreted from CPTU data. Geotech Eng Conf Year 2000 Geotech, 18.
[23]	Mayne PW (1991) Determination of OCR in clays by piezocone tests using cavity expansion and critical state concepts. Soils Found 31: 65-76. doi: 10.3208/sandf1972.31.2_65
[24]	Chen BY, Mayne PW (1994) Profiling the Overconsolidation Ratio of Clays by Piezocone Tests, Report No. GIT-CEE/GEO-94-1 submitted to National Science Foundation by Georgia Institute of Technology, Atlanta, 280.
[25]	Burns SE, Mayne PW (1998) Monotonic and dilatory porewater pressures during piezocone dissipation tests in clay. Can Geotech J 35: 1063-1073. doi: 10.1139/t98-062
[26]	Agaiby SS, Mayne PW (2018) Interpretation of piezocone penetration and dissipation tests in sensitive Leda Clay at Gloucester Test Site. Can Geotech J 55: 1781-1794. doi: 10.1139/cgj-2017-0388
[27]	Mayne PW, Greig J, Agaiby S (2018) Evaluating CPTu in sensitive Haney clay using a modified SCE-CSSM solution. 71st Can Geotech Conf GeoEdmonton, Paper ID No. 279, Canadian Geotechnical Society. Available from: www.cgs.ca.
[28]	Mayne PW, Paniagua P, L'heureux JS, et al. (2019) Analytical CPTu model for sensitive clay at Tiller-Flotten site, Norway, XVⅡ ECSMGE: Geotechnical Engineering Foundation of the Future, Paper 0153, Reykjavik, Icelandic Geotechnical Society. Available from: www.issmge.org.
[29]	Di Buò B, D'Ignazio M, Selãnpaã J, et al. (2019) Yield stress evaluation of Finnish clays based on analytical CPTU models. Can Geotech J 57: 1623-1638. doi: 10.1139/cgj-2019-0427
[30]	Agaiby SS, Mayne PW (2018) Evaluating undrained rigidity index of clays from piezocone data. Cone Penetration Testing (Delft), CRC Press/Balkema, 65-72.
[31]	Senneset K, Sandven R, Janbu N (1989) Evaluation of soil parameters from piezocone tests. Transp Res Rec, 24-37.
[32]	Mayne PW (2007) In-situ test calibrations for evaluating soil parameters, Characterization & Engineering Properties of Natural Soils, Taylor & Francis, London, 1602-1652.
[33]	Mayne PW (2007) NCHRP Synthesis 368: Cone Penetration Testing. Transportation Research Board, National Academies Press, Washington, DC, 118. Available from: www.trb.org.
[34]	Ouyang Z, Mayne PW (2018) Effective friction angle of clays and silts from piezocone. Can Geotech J 55: 1230-1247. doi: 10.1139/cgj-2017-0451
[35]	Ouyang Z, Mayne PW (2019) Modified NTH method for assessing effective friction angle of normally consolidated and overconsolidated clays from piezocone tests. ASCE J Geotech Geoenviron Eng 145.
[36]	Houlsby GT, The CI (1988) Analysis of the piezocone in clay. Penetration Testing 1988, Balkema, Rotterdam, 777-783.
[37]	Lunne T, Long M, Forsberg CF (2003) Characterization and engineering properties of Onsøy clay. Charact Eng Prop Nat Soils, 395-427.
[38]	Lunne T, Randolph M, Sjursen MA, et al. (2006) Shear strength parameters determined by in-situ tests for deep water soft soils. NGI-COFS Report 20041618-1. Joint Industry Project by the Norwegian Geotechnical Institute, Oslo and Centre for Offshore Foundation Systems, Perth: 558.
[39]	Gundersen A, Hansen R, Lunne T, et al. (2019) Characterization and engineering properties of the NGTS Onsøy soft clay site. AIMS Geosci 5: 665-703. doi: 10.3934/geosci.2019.3.665
[40]	Chung SG, Ryu CK, Min SC, et al. (2012) Geotechnical characterization of Busan clay. KSCE J Civ Eng 16: 341-350. doi: 10.1007/s12205-012-1433-8
[41]	Chung SG, Kweon HJ (2013) Oil-operated fixed-piston sampler and its applicability. J Geotech Geoenviron Eng 139: 134-142. doi: 10.1061/(ASCE)GT.1943-5606.0000730
[42]	Pineda JA, McConnell A, Kelly RB (2014) Performance of an innovative direct push piston sampler in soft clay. Proc 3rd Symp Cone Penetration Test, 279-288
[43]	Pineda JA, Suwal LP, Kelly RB, et al. (2016) Characterization of Ballina clay. Géotechnique 66: 556-577. doi: 10.1680/jgeot.15.P.181
[44]	Pineda JA, Kelly RB, Suwal L, et al. (2019) The Ballina soft soil field testing facility. AIMS Geosci 5: 509-534. doi: 10.3934/geosci.2019.3.509
[45]	Hight DW, Paul MA, Barras BF, et al. (2003) The characterization of the Bothkennar clay. Characterization and Engineering Properties of Natural Soils, Swets & Zeitlinger, Lisse, 543-597.
[46]	Mayne PW (2008) Piezocone profiling of clays for maritime site investigations. Geotechnics in Maritime Engineering, Polish Committee on Geotechnics, 333-350.
[47]	Getchell A, Santamaria A, Benoît J (2014) Geotechnical Test Embankment on soft marine clay in Newington—Dover, MS Thesis, Civil Engineering Dept, Univ of New Hampshire. Durham, NH: 103. Available from: https://scholars.unh.edu/thesis/825.
[48]	Locat A (2012) Rupture progressive et étalements dans les argiles sensible. PhD Dissertation, Université Laval, Quebec, 216.
[49]	Locat A, Locat P, Demers D, et al. (2017) The Saint-Jude landslide of 10 May 2010, Quebec, Canada: Investigation and characterization of the landslide and its failure mechanism. Can Geotech J 54: 1357-1374. doi: 10.1139/cgj-2017-0085
[50]	Locat A, Locat P, Michaud H, et al. (2019) Geotechnical characterization of the Saint-Jude clay, Quebec, Canada. AIMS Geosci 5: 273-302. doi: 10.3934/geosci.2019.2.273
[51]	Paniagua P, L'Heureux JS, Carroll R, et al. (2017) Evaluation of sample disturbance of three Norwegian clays. 19th ICSMGE Secr Seoul. Available from: www.issmge.org.
[52]	Lehtonen V (2015) Modelling undrained shear strength and pore pressure based on an effective stress soil model in Limit Equilibrium Method, Tampereen teknillinen yliopisto. Julkaisu-Tampere University of Technology. Publication, 213.
[53]	Di Buò B, D'Ignazio M, Selãnpaã J, et al. (2016) Preliminary results from a study aiming to improve ground investigation data. 17th Nord Geotech Meet 1: 25-28.
[54]	Wang B, Brooks GR, Hunter JAM (2015) Geotechnical data from a large landslide site at Quyon, Report 7904, Quebec Geological Survey of Canada, 54.
[55]	Wang B, Brooks GR, Hunter JAM (2015) Geotechnical investigations of a large landslide site at Quyon, Québec. 68th Can Geotech Conf.
[56]	Lafleur J, Silvestri V, Asselin R, et al. (1988) Behavior of a test excavation in soft Champlain Sea clay. Can Geotech J 25: 705-715. doi: 10.1139/t88-081
[57]	Chiasson P, Lafleur J, Soulié M, et al. (1995) Characterizing spatial variability of a clay by geostatistics. Can Geotech J 32: 1-10. doi: 10.1139/t95-001
[58]	Sandven R, Montafia A, Gylland A, et al. (2015) Detection of brittle materials. Summary report with recommendations. Final report. NIFS Report no. 126/2015,150.
[59]	Helle TE, Long M, Nordal S (2018) Interpreting improved geotechnical properties from RCPTUs in KCl-treated quick clays, Cone Penetration Testing 2018, CRC Press, Taylor & Francis Group, London, 339-345.
[60]	Edil TB, Wang X (2000) Shear strength and K₀ of peats and organic soils. Geotechnics of High Water Content Materials, American Society for Testing & Materials, West Conshohocken, PA, 209-225.
[61]	Den Haan EJ, Kruse GAM (2007) Characterization and engineering properties of Dutch peats, Characterization & Engineering Properties of Natural Soils, Taylor & Francis Group, London, 2101-2133.
[62]	Mesri G, Ajlouni M (2007) Engineering properties of fibrous peats. J Geotech Geoenviron Eng 133: 850-966. doi: 10.1061/(ASCE)1090-0241(2007)133:7(850)
[63]	Jannuzzi GMF, Danziger FAB, Martins ISM (2015) Geological-geotechnical characterization of Sarapuí Ⅱ clay. Eng Geol 190: 77-86. doi: 10.1016/j.enggeo.2015.03.001
[64]	Larsson R, Westerberg B, Albing D, et al. (2007) Sulfidjord: geoteknisk klassificering och odraneråd skjuvhållfastthet. SGI Report 69, Swedish Geotechnical Institute, Linköping, 138.
[65]	Lamb RA, Chow LC, Bentler JG (2018) US Highway 14 embankment over soft soils—success with ground improvement and modern instrumentation. 66th Annu Geotech Eng Conf, 117-126.
[66]	Chow LC, Bentler JG, Lamb RA (2019) Primary and post-surcharge secondary settlements of a highway embankment constructed over highly organic soils: a case history. Geo Congr 2019 Embankments Dams Slopes, 109-118.
[67]	McCabe BA (2002) Experimental investigations of driven pile group behaviour in Belfast soft clay, Doctoral dissertation, Trinity College Dublin, 415.
[68]	Lehane BM (2003) Vertically loaded shallow foundation on soft clayey silt. Proc Inst Civil Eng Geotech Eng 156: 17-26. doi: 10.1680/geng.2003.156.1.17
[69]	Westerberg B, Andersson M (2017) Sulfidjord—kompressionsegenskaper och sättningar. En studie av provbankar i Lampen och andra bankar. SGI Publikation 41, Swedish Geotechnical Institute, Linköping, 238.
[70]	Merani JM, Hunt CE, Donahue JL, et al. (2016) CPT interpretation in highly organic soils and soft clay soils. Geo Chicago 2016, 412-421.
[71]	Mayne PW (1987) Determining preconsolidation stress and penetration pore pressures from DMT contact pressures. Geotech Test J 10: 146-150. doi: 10.1520/GTJ10947J
[72]	Andersson M (2012) Kompressionsegenskaper hos sulfidjordar: En fält-och laboratoriestudie av provbankar, Luleå tekniska universitet, Sweden, 336.
[73]	Westerberg B, Andersson M, Winter MG, et al. (2015) Compression properties of an organic clay. XVI Eur Conf Soil Mech Geotech Eng, 3091–3096.
[74]	Baroni M, Almeida MSS (2017) Compressibility and stress history of very soft organic clays. Proce Inst Civ Eng 170: 148-160. doi: 10.1680/jgeen.16.00146
[75]	Mayne PW (2017) Stress history of soils from cone penetration tests. Soils Rocks 40: 203-218. doi: 10.28927/SR.403203
[76]	Agaiby SS, Mayne PW (2019) CPT evaluation of yield stress profiles in soils. J Geotech Geoenviron Eng 145: 04019104. doi: 10.1061/(ASCE)GT.1943-5606.0002164

This article has been cited by:

1.	Rinaldo M. Colombo, Mauro Garavello, Infectious Disease Spreading Fought by Multiple Vaccines Having a Prescribed Time Effect, 2023, 71, 0001-5342, 10.1007/s10441-022-09452-4
2.	Lihong Sun, Qiang He, Yueyang Teng, Qi Zhao, Xin Yan, Xingwei Wang, A complex network-based vaccination strategy for infectious diseases, 2023, 136, 15684946, 110081, 10.1016/j.asoc.2023.110081
3.	R.M. Colombo, M. Garavello, M. Tandy, On the coupling of well posed differential models, 2023, 232, 0362546X, 113290, 10.1016/j.na.2023.113290
4.	Rinaldo M. Colombo, Mauro Garavello, Francesca Marcellini, Elena Rossi, 2024, Chapter 28, 978-3-031-55263-2, 325, 10.1007/978-3-031-55264-9_28

Reader Comments

Your name:*

Email:*
© 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)