Loading [MathJax]/jax/output/SVG/jax.js
Communication

Experimental test of temperature and moisture controls on the rate of microbial decomposition of soil organic matter: preliminary results

  • Soil organic matter (SOM) is a major reservoir of carbon derived from the biosphere that is returned to the atmosphere largely via microbial decomposition. The potential for feedbacks between climate change and SOM decomposition makes a full understanding of the controls on SOM decomposition rates essential to modeling future climate changes. We measured soil CO2 flux in a laboratory setting using pots containing a uniform mix of soil in which we varied both temperature and moisture. Following initial desiccation, a strong CO2 pulse was measured within two hours of rewetting and a return to equilibrium conditions obtained within 168 hours, with the magnitude of the initial pulse varying by soil temperature and moisture addition. At equilibrium conditions, no correlation was found between CO2 flux and temperature across all moisture levels, although a weak positive correlation (r2=0.1 to 0.2) was seen at moderate to high moisture levels. A much stronger correlation (r2 > 0.4) was found between CO2 flux and soil moisture across the full range of temperatures and at both low and high temperatures. Thus, we conclude that when all other variables were constrained, soil moisture fluctuations appeared to have greater impact than temperature variations on the rate of microbial decomposition of SOM. These preliminary results suggest directions for future research examining the relationships between soil moisture, temperature and CO2 flux for soils in which clay mineral and/or SOM composition are varied.

    Citation: Corey A. Palmer, Katherine P. Markstein, Lawrence H. Tanner. Experimental test of temperature and moisture controls on the rate of microbial decomposition of soil organic matter: preliminary results[J]. AIMS Geosciences, 2019, 5(4): 886-898. doi: 10.3934/geosci.2019.4.886

    Related Papers:

    [1] Gayane Semerjyan, Inesa Semerjyan, Mikayel Ginovyan, Nikolay Avtandilyan . Characterization and antibacterial/cytotoxic activity of silver nanoparticles synthesized from Dicranum scoparium moss extracts growing in Armenia. AIMS Biophysics, 2025, 12(1): 29-42. doi: 10.3934/biophy.2025003
    [2] Shen Helvig, Intan D. M. Azmi, Seyed M. Moghimi, Anan Yaghmur . Recent Advances in Cryo-TEM Imaging of Soft Lipid Nanoparticles. AIMS Biophysics, 2015, 2(2): 116-130. doi: 10.3934/biophy.2015.2.116
    [3] Christophe A. Monnier, David C. Thévenaz, Sandor Balog, Gina L. Fiore, Dimitri Vanhecke, Barbara Rothen-Rutishauser, Alke Petri-Fink . A guide to investigating colloidal nanoparticles by cryogenic transmission electron microscopy: pitfalls and benefits. AIMS Biophysics, 2015, 2(3): 245-258. doi: 10.3934/biophy.2015.3.245
    [4] O.S. Sorzano Carlos, Vargas Javier, Otón Joaquín, Abrishami Vahid, M. de la Rosa-Trevín José, del Riego Sandra, Fernández-Alderete Alejandro, Martínez-Rey Carlos, Marabini Roberto, M. Carazo José . Fast and accurate conversion of atomic models into electron density maps. AIMS Biophysics, 2015, 2(1): 8-20. doi: 10.3934/biophy.2015.1.8
    [5] Nicholas Spellmon, Xiaonan Sun, Wen Xue, Joshua Holcomb, Srinivas Chakravarthy, Weifeng Shang, Brian Edwards, Nualpun Sirinupong, Chunying Li, Zhe Yang . New open conformation of SMYD3 implicates conformational selection and allostery. AIMS Biophysics, 2017, 4(1): 1-18. doi: 10.3934/biophy.2017.1.1
    [6] Suleyman Yilmaz . Using elastic scattering to determination of diseases via urine samples. AIMS Biophysics, 2021, 8(4): 307-317. doi: 10.3934/biophy.2021024
    [7] Michelle de Medeiros Aires, Janine Treter, Antônio Nunes Filho, Igor Oliveira Nascimento, Alexandre José Macedo, Clodomiro Alves Júnior . Minimizing Pseudomonas aeruginosa adhesion to titanium surfaces by a plasma nitriding process. AIMS Biophysics, 2017, 4(1): 19-32. doi: 10.3934/biophy.2017.1.19
    [8] Ateeq Al-Zahrani, Natasha Cant, Vassilis Kargas, Tracy Rimington, Luba Aleksandrov, John R. Riordan, Robert C. Ford . Structure of the cystic fibrosis transmembrane conductance regulator in the inward-facing conformation revealed by single particle electron microscopy. AIMS Biophysics, 2015, 2(2): 131-152. doi: 10.3934/biophy.2015.2.131
    [9] Adam Redzej, Gabriel Waksman, Elena V Orlova . Structural studies of T4S systems by electron microscopy. AIMS Biophysics, 2015, 2(2): 184-199. doi: 10.3934/biophy.2015.2.184
    [10] Vittoria Raimondi, Alessandro Grinzato . A basic introduction to single particles cryo-electron microscopy. AIMS Biophysics, 2022, 9(1): 5-20. doi: 10.3934/biophy.2022002
  • Soil organic matter (SOM) is a major reservoir of carbon derived from the biosphere that is returned to the atmosphere largely via microbial decomposition. The potential for feedbacks between climate change and SOM decomposition makes a full understanding of the controls on SOM decomposition rates essential to modeling future climate changes. We measured soil CO2 flux in a laboratory setting using pots containing a uniform mix of soil in which we varied both temperature and moisture. Following initial desiccation, a strong CO2 pulse was measured within two hours of rewetting and a return to equilibrium conditions obtained within 168 hours, with the magnitude of the initial pulse varying by soil temperature and moisture addition. At equilibrium conditions, no correlation was found between CO2 flux and temperature across all moisture levels, although a weak positive correlation (r2=0.1 to 0.2) was seen at moderate to high moisture levels. A much stronger correlation (r2 > 0.4) was found between CO2 flux and soil moisture across the full range of temperatures and at both low and high temperatures. Thus, we conclude that when all other variables were constrained, soil moisture fluctuations appeared to have greater impact than temperature variations on the rate of microbial decomposition of SOM. These preliminary results suggest directions for future research examining the relationships between soil moisture, temperature and CO2 flux for soils in which clay mineral and/or SOM composition are varied.


    It is universally acknowledged that multi-agent cooperative control has a lot of applications, such as unmanned air vehicles, traffic control, animal groups and automated highway systems [1,2,3,4]. The cooperative control consensus that has received a lot of attention can be roughly divided into leader-follower consensus and leaderless consensus [5,6]. Most of the research aims have been to design a consensus protocol to make agents exchange local information with their neighbors, so that a cluster of agents are capable of achieving the consistent state.

    As a paradigmatic instance, a digital microprocessor is installed in each agent of the system, which is responsible for gathering information from its neighboring agents and updating the controller accordingly. The majority of research studies utilize continuous measurement signals, yet continuous communication in a constrained energy exchange network is not feasible. In order to avoid continuous communication, some studies have introduced sampling data control [7,8,9], and each agent transmits the corresponding data at the sampling instants. Control using periodic sampling data often requires estimation of an optimal sampling period since improper selection can result in either excessively frequent or infrequent sampling, both of which can be detrimental to system performance. Stochastic sampling can better solve this problem. Stochastic sampling has been gaining increasing attention due to its flexibility in terms of dynamically switching sampling periods between different values. In [10], the authors investigated the leaderless multi-agent consensus problem under stochastic sampling, where the sampling period was chosen randomly at a given value.

    Compared to the conventional time-triggered mechanism, the utilization of an event-triggered mechanism presents a significant advantage in terms of improving the efficiency of communication resource allocation. Research on event-triggered mechanisms has been extensively conducted in various fields, including cyber-physical systems [11], cyber-control systems [12] and multi-agent systems [13,14,15,16]. To solve the high-order multi-agent consensus problem, in [16], Wu et al. presented by estimating the state of neighbor agents, a novel event-triggered protocol. For further research, many scholars put forward dynamic event-triggered protocol [17,18,19,20,21]. In [17], a dynamic event-triggered protocol was proposed for individual agents, which established a distributed adaptive consensus protocol. This protocol involves updating the coupling strength to achieve consensus among the agents. A dynamic event-triggered protocol was proposed in [20] for the investigation of consensus in multi-agent systems. This protocol involved internal dynamic variables, and the elimination of Zeno behavior played an important role. In [21], Du et al. studied the multi-agent problem of leader-follower consensus based on a dynamic event-triggered mechanism. However, since event-triggered control has its own limitation, there is a need for continuous event detection. To loosen that constraint, many researchers designed an event-triggered protocol based on sampling data [22,23,24]. The authors presented research on event-triggered consensus strategies for multi-agent systems in [25] and [26]. Specifically, Su et al. investigated sampled data-based leader-follower multi-agent systems with input delays, with a focus on making sampling-based mechanisms for event detection more realistic. Meanwhile, He et al. focused on mean-square leaderless consensus for networked non-linear multi-agent systems and presented an efficiently distributed event-triggered mechanism that reduces communication costs and controller updates for random sampling-based systems. In [27], Ruan et al. studied the consensus problem with bounded external disturbances under an event-triggered scheme based on two independent dynamic thresholds in the context of leader-follower dynamics. The authors presented a nonlinear dynamic event-triggered control strategy for achieving prescribed-time synchronization in networks of piecewise smooth systems in [28].

    The above studies are mostly the ones based on the fixed topology of multiple agents. A large number of switching topologies, some of which are studied by establishing Markov models, can be observed in our real world [29,30,31]. In [29], Hu et al. investigated the multi-agent consensus of Markov jump systems based on event-triggered strategies. In [30,31], the scholars studied consensus issues in multi-agent systems with a Markov network structure. Nevertheless, the application of time-varying topology based on Markov processes has certain limitations, primarily due to the exponential distribution of jump times in Markov chains. Some researchers have focused on semi-Markovian jump topologies [32,33,34,35,36]. This is because the sojourn time in the semi-Markovian exchange topology is a generic continuous random variable. Its probability distribution is general. The H consensus problem for multi-agent systems in a semi-Markov switching topology with incomplete known transmission rates was investigated in [35]. In a related research study by Xie et al. [36], the consensus problem of multi-agent systems was investigated under an attack scenario in which both the semi-Markov switching topology and network were susceptible to attacks, with the possibility of recovery. In [37], the authors presented a study on achieving cluster synchronization in finite/fixed time for semi-Markovian switching TS fuzzy complex dynamical networks with discontinuous dynamic nodes.

    Building upon the aforementioned literature review, this research article is focused on the development and analysis of an event-triggered mechanism for the stochastic sampling of leader-follower multi-agent systems. The key contributions of this study can be summarized as follows:

    (1). This paper proposes a novel event-triggered methodology by utilizing stochastic sampling, which is capable of significantly reducing the frequency of control updates and communication overhead among agents. Furthermore, the proposed mechanism ensures the avoidance of Zeno behavior.

    (2). The dynamic system switching investigated in this paper is modeled by using the semi-Markov jump process. A sufficient condition for mean-square consensus is derived.

    (3). The sampling period of stochastic sampling is randomly selected from a finite set. Stochastic sampling differs from the periodic sampling and stochastic sampling in the literature [26,29].

    The subsequent sections of this article are organized as follows. Section Ⅱ presents the problem formulation and introductory suggestions. Section Ⅲ reports the major findings. The numerical tests in Section Ⅳ validate the accuracy of the theoretical conclusions. Section Ⅴ concludes with final remarks.

    Notations: The n-dimensional identity matrix is denoted by In. A zero matrix of appropriate dimension is represented by O. A positive (negative) definite matrix A is denoted by A>0 (A<0). The element implied by the symmetry of a matrix is denoted by . The function C([a,b],Rn) maps the interval [a,b] to a continuous vector-valued function Rn. The Euclidean norm of a vector is represented by ||. A superscript T and the symbol indicate matrix transposition and the Kronecker product, respectively. Let IN={1,2,,N} denote a finite index set.

    We consider a directed graph, denoted by G={V,E,A}, where V represents a set of vertices {1,2,,N}, E denotes a set of directed edges and A is a weighted adjacency matrix of size N×N. The elements of A, denoted by aij, are positive if there exists a directed edge going from vertex i to vertex j, and zero otherwise. We can refer to the set of neighbors of vertex i as Ni, which is defined as Ni={jV:(j,i)E}. The degree matrix D is defined as a diagonal matrix of size N×N, where the diagonal entries di represent the weighted degree of vertex i. Specifically, di=jNiaij. The Laplacian matrix of G is defined as L=DA=(lij)N×N. The diagonal entries of L are given by lii=jNilij, and the off-diagonal entries are defined as follows: lij=aij if (i,j) is an edge in G; otherwise, lij=0.

    Consider a multi-agent system that has a leader and N followers. Label them respectively as 0 and 1,2,3,,N. The dynamic equations for each agent are illustrated below:

    {˙xi(t)=A(r(t))xi(t)+B(r(t))ui(t),iIN,˙x0(t)=A(r(t))x0(t). (2.1)

    The control input and state for the ith state are denoted by ui(t)Rn and xiRn, respectively. A state of x0(t)Rn indicates the leader agent. The constant matrices A(r(t)) and B(r(t)) have the appropriate dimensions. The semi-Markov chain {r(t),t0} is a right-continuous process defined on the complete probability space (Ω,F,P), where its values belong to the finite state space D=IM and its generators are denoted by λ=(λmn(v))M×M. The transition probabililty is

    Pr{r(t+v)=nr(t)=m}={λmn(v)v+o(v)mn,1+λmn(v)v+o(v)m=n. (2.2)

    v is the time interval, which stands for the amount of time that passes between two successive jumps. The transition rate from mode m at time t to mode n at time t+v is denoted by λmn(v), and it satisfies λmn(v)0. o(v) can be defined as limv0+o(v)v=0. λmm(v)=mnλmn(v).

    Remark 1: The dwell time v is the time elapsed from the last jump of the system, and it is distinct from the time t. When the system jumps, v resets to 0 and the transition probability λmn(v) only depends on v.

    A consensus protocol for event-triggered consensus that is grounded in the stochastic sampling of data is presented. Assuming that the sampling time is 0=t0<t1<t2<<ts<, the sampling period is h=ts+1ts, in which h is selected from a random finite set h1,h2,,hl. The probability is described as Pr{h=hs}=πs, sIl, πs[0,1], and ls=1πs=1. For the sake of generality, we can set 0=h0<h1<h2<<hs<<hl,l>1.

    Assuming that the ith agent has a K-time event-triggered time of tik, {tik}k=0 represents the event-triggered time sequence of the ith agent, tik{ts,sN}, ti0=0. tik+1represents the next event-triggering time of the ith agent, which is determined by the following formula

    tik+1=mints>tk{ts:(ei(ts))TΦ(ei(ts))>σi(zi(ts))TΦ(zi(ts))}. (2.3)

    The threshold parameter is represented in this case by σi>0. A positive definite event-triggered matrix is the intended matrix of matrix Φ. ei(ts)=xi(tik)xi(ts) and zi(ts)=Nj=1aij(xi(tik)xj(tjk))+bi(xi(tik)x0(ts)), where ts[tik,tik+1). tik is the latest transmitted sampled data of its neighbors before tik, that is tjk=max{tjktjktik}, k=0,1,2,.

    Remark 2: The stochastic sampling sequence is 0=t0<t1<t2<<ts<. The sampling period is h=ts+1ts, where h is selected from a random finite set {h1,h2,,hl}, where 0=h0<h1<h2<<hl. This stochastic sampling differs from the periodic sampling and stochastic sampling in the literature [26,29]. They represent the special form of stochastic sampling when h is constant and l equals 2.

    Remark 3: In accordance with the event-triggered condition (2.3), the ith agent broadcasts the most recent sampled data to its neighbors. The sampling sequence includes the event-triggered sequence because the sampling period h=ts+1ts is stochastic in the set of {h1,h2,,hl},h>0. Zeno behavior is precisely precluded.

    Remark 4: It should be noted that in an attempt to decrease unnecessary communication between agents, a stochastic sampling static event-triggered protocol is proposed.

    The following consensus protocol should be taken into consideration in light of the discussion above:

    ui(t)=K(r(t))[Nj=1aij(xi(tik)xj(tjk))+bi(xi(tik)x0(ts))]. (2.4)

    Remark 5: The consensus protocol relies upon the stochastic sampling of event-triggered conditions as well as upon a semi-Markov switching system, where the feedback gain K(r(t)) depends on r(t), which will be given by a theorem later. When the ith agent and its neighbor agents satisfy the trigger conditions, the controller will be updated. With a zero-order holder, ui(t) remains constant between two successive event instants.

    Define

    ei(ts)=xi(tik)xi(ts),ej(ts)=xj(tjk)xj(ts),δi(ts)=xi(ts)x0(ts). (2.5)

    Submitting (2.5) into (2.4), we can obtain

    ui(t)=K(r(t))[Nj=1lij(ej(ts)+δj(ts))+bi(ei(ts)+δi(ts))]. (2.6)

    Define τ(t)=tts, where τ(t) is a piecewise linear function with a slope of ˙τ(t)=1 for all t[ts,ts+1), except at time ts. The control protocol (2.6) may be expressed as

    ui(t)=K(r(t))[Nj=1lij(ej(tτ(t))+δj(tτ(t)))+bi(ei(tτ(t))+δi(tτ(t)))]. (2.7)

    A new stochastic variable is introduced as follows:

    βs(t)={1ts1τ(t)<ts,s=1,2,,l0 otherwise.  (2.8)

    In this way, we can obtain

    Pr{βs(t)=1}=Pr{ts1τ(t)<ts}=li=sπihihi1hi=βs. (2.9)

    The Bernoulli distribution is satisfied by βs(t), given that E{βs(t)}=βs and E{βs(t)β2s}=βs(1βs) respectively. We can obtain the calculation of the following formula from the above study:

    ˙δ(t)=(INA(r(t)))δ(t)ls=1βs(t)(HB(r(t))K(r(t)))e(tτs(t))ls=1βs(t)(HB(r(t))K(r(t)))δ(tτs(t)), (2.10)

    where δ=diag{δ1,δ2,,δN}, H=L+B1, B1=diag{b1,b2,,bN}.

    According to the initial conditions of Equation (2.10), let δ(t)=ϕ(t),hlt0, ϕ(t)=[ϕT1(t),ϕT2(t),,ϕTN(t)] and ϕi(t)C([hl,0],Rn).

    Definition 1 [34]. Under semi-Markov switching topologies, the leader-follower consensus of multi-agent system (2.1) with the consensus protocol is said to be achieved if limtExi(t)x0(t)=0, iIN holds for any initial distribution r0D and any initial condition ϕ(t),t[hl,0].

    Assumption 1. The directed spanning tree in the network graph G has the leader's root.

    Lemma 1 [34]. For symmetric matrices R>0 and X and any scalar μ, the following inequality holds:

    XR1Xμ2R2μX.

    Theorem 1. Under Assumption 1 and utilizing the protocol given by (2.7), the constants provided are 0=h0<h1<<hs<<hl and πi[0,1], σi>0, iIN. Consensus can be achieved in the mean-square sense for the multi-agent system (2.10) by employing the stochastic sampling event-triggered strategy given by (2.3), provided that constant metrics are present, namely P(m)>0,Qs>0,Rs>0,Ws>0,mDandsIl. In such a way, the inequality below holds:

    (ΞFT(m)Σ(m,v)Ψ0X2(m))<0, (3.1)

    where

    Ξ=1+2+3+4+5+6.

    1=FT(m)(INP(m))ε1+εT1(INP(m))F(m)+Mn=1λmn(v)(εT1(INP(n))ε1).

    2=β1εT1(INQs)ε1β1εTl+1(INQs)εi+1+ls=2βs(εTl+s(INQs)εl+sεTl+s+1(INQs)εl+s+1).

    3=ls=1βs1hshs1((εTl+sεTl+s+1)×(IN(Rs+Ws))(εl+sεl+s+1)).

    4=εT1(INP(m))ε1.

    5=ls=1βsεTs+1Φεs+1.

    6=(εs+1+ε2l+s+1)T(HTΛHΦ)(εs+1+ε2l+s+1).

    F(m)=(INA(m))ε1ls=1βs(HB(m)K(m))ε2l+s+1ls=1βs(HB(m)K(m))εs+1.

    Φ=diag{Φ1,Φ2,,ΦN}, Σ(m,v)=λ(m,v)X1(m).

    Ψ=(ls=1βs(hshs1)(IN(Rs+Ws)))1.

    λ(m,v)=(λm1(v),λm2(v),,λmm1(v),λmm+1(v),,λmM(v)).

    X1(m)=diag{INP(m),,INP(m)}M1.

    X2(m)=diag{INP(1),INP(2),,INP(m1),INP(m+1),,INP(M)}M1.

    A(m)=A(r(t)=m),P(m)=P(r(t)=m) and βs is defined the same way as in (2.10). Define εs as a block matrix consisting of 3l+1 block elements. The s-th block element is an Nn×Nn identity matrix, denoted by INn, while all other block elements are zero matrices. Therefore, εs can be expressed as εs=[0,0,,INn,0,0,,0]RNn×(3l+1)Nn for s=1,2,,3l+1.

    Proof of Theorem 1. Think about the Lyapunov-Krasovskii functional presented as follows:

    V(t,δt,˙δt,r(t))=3i=1Vi(t,δt,˙δt,r(t)),t[ts,ts+1], (3.2)

    where

    V1(t,δt,˙δt,r(t))=δT(t)(INP(r(t)))δ(t), (3.3)
    V2(t,δt,˙δt)=ls=1βsths1thsδT(μ)×(INQs)δ(μ)dμ, (3.4)
    V3(t,δt,˙δt)=ls=1βshs1hstt+v˙δT(μ)(IN(Rs+Ws))˙δ(μ)dμdv. (3.5)

    Consider the weak infinitesimal generator

    V(t,zt)=limΔ0+1Δ{E{V(t+Δ,δt+Δ,˙δt+Δ,r(t+Δ))δt,r(t)=m}V(t,δt,˙δt,r(t))}. (3.6)

    Introduce y(t)=(δT(t),δT(tτ1(t)),,δT(tτl(t)),δT(th1),, δT(thl),eT(tτ1(t)),,eT(tτl(t))),y(t)R(3l+1)Nn, A(m)=A(r(t)=m),P(m)=P(r(t)=m),mD.

    Thus, we obtain

    E[V1(t,δt,˙δt,r(t))]=E[limΔ0+1Δ{E{V1(t+Δ,δt+Δ,˙δt+Δ,r(t+Δ)δt,r(t)=m)}V1(t,δi,˙δi,r(t))}]=E[yT(t)(FT(m)(INP(m))ε1+ετ1(INP(m))F(m)+Mn=1λmn(v)ετ1(INP(n))ε1)y(t)], (3.7)

    where

    F(m)=(INA(m))ε1ls=1βs(HB(m)K(m))ε2l+s+1ls=1βs(HB(m)K(m))εs+1.
    E[V2(t,δi,˙δt)]=E{yT(t)(β1δT(th0)(INQ1)δ(th0)β1δT(th1)(INQ1)δ(th1)+ls=2βs(δT(ths1)(INQs)δ(ths1)δT(ths)(INQs)δ(ths)))y(t)}, (3.8)

    and

    E[V3(t,δt,˙δt)]=E[ls=1βs(hshs1)˙δT(t)(IN(Rs+Ws))˙δ(t)ls=1βsths1ths˙δT(v)(IN(Rs+Ws))˙δ(v)dv]. (3.9)

    According to Jensen's inequality, it can be obtained that

    ls=1βsths1ths˙δT(ν)(IN(Rs+Ws))˙δ(ν)dνls=1βs1hshs1ths1ths˙δT(ν)dν×(IN(Rs+Ws))ths1ths˙δ(ν)dν. (3.10)

    Submit (2.10) and (3.10) into (3.9), and we can obtain

    E[V3(t,δt,˙δt,r(t))]yT(t)[ls=1βs(hshs1)×FT(m)(IN(Rs+Ws))F(m)ls=1βs1hshs1(εTl+sεTl+s+1)×(IN(Rs+Ws))(εTl+sεTl+s+1)]y(t). (3.11)

    From (2.3), we can obtain

    eT(tτs)Φe(tτs)zT(tτs)(ΛΦ)z(tτs)=(δ(tτs)+e(tτs))T(HTΛHΦ)×(δ(tτs)+e(tτs)), (3.12)

    where e(tτs)=col{e1(tτs),e2(tτs),,eN(tτs)}, Λ=diag{σ1,σ2,,σN}. Additionally, z(tτs)=col{z1(tτs),z2(tτs),,zN(tτs)}.

    Thus, combine (3.8–3.10) with (3.11), and we can obtain

    E[V(t,δt,˙δt,r(t))]E[yT(t)ls=1βs(ΞETs+1ΦEs+1+(Es+1+E2l+s+1)T(HTΛHΦ)(Es+1+E2l+s+1)]y(t), (3.13)

    and

    E[V(t,δt,˙δt,r(t))]E[yT(t)ls=1βs(ΞETs+1ΦEs+1+(Es+1+E2l+s+1)T(HTΛHΦ)(Es+1+E2l+s+1)]y(t). (3.14)

    Applying the Schur complement and (3.1) leads to the conclusion that

    E[V(t,δt,˙δt,r(t))]<0, (3.15)

    where

    Ξ=1+2+3+4+5+6.

    Hence, the consensus (3.1) of the multi-agent system can be attained in a mean squared sense under the event-triggered method given by (2.3).

    Keeping Theorem 1's results in mind, this method provides an efficient approach to design consensus controller gains.

    Theorem 2. Under Assumption 1, and by utilizing the protocol given by (2.7), we have the following constants 0=h0<h1<<hs<<hl,πi[0,1],σi>0,iIN and μ>0. The multi-agent system consensus (2.10) can be achieved in the mean-square sense with the stochastic sampled even-triggered strategy (2.3) if there exist the matrices ˆP(m)>0,ˆQs>0,ˆRs>0,ˆWs>0 and ˆΦ>0 and the matrices ˆK(m),mD and sIl satisfy the following inequality:

    (ˆΞ(hmh(m1))ˆFT(m)ˆΣ(m,v)ˆΨOˆX2(m))<0, (3.16)

    where ˆΞ=ˆ1+ˆ2+ˆ3+ˆ4+ˆ5+ˆ6.

    ˆ1=ˆFT(m)ε1+εT1ˆF(m)+εT1λmm(v)(INˆP(m))ε1.

    ˆ2=β1εT1(INˆQ1)ε1β1εTl+1(INˆQ1)εl+1 +ls=2βs(εTl+s(INˆQs)εl+sεTl+s+1(INˆQs)εl+s+1).

    ˆ3=ls=1βs1hshs1((εTl+sεTl+s+1)(IN(ˆRs+ˆWs))(εl+sεl+s+1)).

    ˆ4=εT1(INˆP(m))ε1.

    ˆ5=ls=1βsεTs+1(ˆΦ)εs+1.

    ˆ6=(εs+1+ε2l+s+1)T(HTΛHˆΦ)(εs+1+ε2l+s+1).

    ˆF(m)=(INA(m)ˆP(m))ε1ls=1βs(HB(m)ˆK(m))ε2l+s+1ls=1βs(HB(m)ˆK(m))εs+1.

    ˆΦ=diag{ˆΦ1,ˆΦ2,,ˆΦN},ˆΦi=ˆP(m)ΦiˆP(m),ˆP(m)=P1(m), ˆΦs=ˆP(m)ΦsˆP(m).

    ˆΨ=ls=1(βs(hshs1))1(μ2(IN(Rs+Ws))2μ(INˆP(m))).

    ˆRs=ˆP(m)RsˆP(m),ˆWs=ˆP(m)WsˆP(m),ˆΣ(m,v)=λ(m,v)ˆX1(m).

    λ(m,v)=(λm1(v),λm2(v),,λmm1(v),λmm+1(v),,λmM(v),0,,0).

    ˆX1(m)=diag{INˆP(m),,INˆP(m),O,,O}.

    ˆX2(m)=diag{INˆP(1),INˆP(2),,INˆP(m1),INˆP(m+1),,INˆP(M),O,,O}.

    In addition, the feedback gain is supplied by K(m)=ˆK(m)ˆP1(m) and the event-triggered parameter matrix is given by Φ(m)=ˆP1(m)ˆΦ(m)ˆP1(m).

    Proof of Theorem 2. Here we present the definitions of matrix variables K(m)=ˆK(m)ˆP1(m),ˆP(m)=P1(m) and ˆΦ(m)=ˆP(m)Φ(m)ˆP(m). We pre- and post-multiply both sides of (3.3) by the matrix diag{INP1(m),INP1(m),INP1(m),InN}, and both sides of (3.2) by the matrix diag{INP1(m),INP1(m)}, respectively.

    Lemma 1 enables one to derive the subsequent inequality.

    (INˆP(m))(ls=1β(hshs1)(IN(Rs+Ws)))1×(INˆP(m))

    μ2(ls=1β(hshs1)(IN(Rs+Ws)))2μ(INˆP(m));

    we can get

    (ˆΞ(hmh(m1))ˆFT(m)ˆΣ(m,v)ˆΨOˆX2(m))<0, (3.17)

    where ˆΞ=ˆ1+ˆ2+ˆ3+ˆ4+ˆ5+ˆ6,

    ˆ1=ˆFT(m)ε1+εT1ˆF(m)+εT1λmm(v)(INˆP(m))ε1,

    ˆ2=β1εT1(INˆQ1)ε1β1εTl+1(INˆQ1)εl+1+ls=2βs(εTl+s(INˆQs)εl+s εTl+s+1(INˆQs)εl+s+1),

    ˆ3=ls=1βs1hshs1((εTl+sεTl+s+1)(IN(ˆRs+ˆWs))(εl+sεl+s+1)),

    ˆ4=εT1(INˆP(m))ε1,

    ˆ5=ls=1βsεTs+1(ˆΦ)εs+1,

    ˆ6=(εs+1+ε2l+s+1)T(HTΛHˆΦ)(εs+1+ε2l+s+1),

    ˆF(m)=(INA(m)ˆP(m))ε1ls=1βs(HB(m)ˆK(m))ε2l+s+1ls=1βs(HB(m)ˆK(m))εs+1,

    ˆΦ=diag{ˆΦ1,ˆΦ2,,ˆΦN}, ˆΦi=ˆP(m)ΦiˆP(m), ˆP(m)=P1(m), ˆΦs=ˆP(m)ΦsˆP(m),

    ˆΨ=ls=1(βs(hshs1))1(μ2(IN(Rs+Ws))2μ(INˆP(m))),

    ˆRs=ˆP(m)RsˆP(m), ˆWs=ˆP(m)WsˆP(m), ˆΣ(m,v)=λ(m,v)ˆX1(m),

    λ(m,v)=(λm1(v),λm2(v),,λmm1(v),λmm+1(v),,λmM(v),0,,0),

    ˆX1(m)=diag{INˆP(m),,INˆP(m),O,,O},

    ˆX2(m)=diag{INˆP(1),INˆP(2),,INˆP(m1),

    INˆP(m+1),,INˆP(M),O,,O}.

    The proof is therefore complete.

    In Theorems 1 and 2, we establish sufficient conditions for achieving consensus in event-triggered semi-Markov jump multi-agent systems through stochastic sampling. But the sufficient conditions do not satisfy linear matrix inequality (LMI) conditions, because λ(m,v) is time-varying. As a result, the problems cannot be directly solved by using the LMI toolbox in MATLAB. Nevertheless, we can establish lower and upper bounds for the transition rate and apply the theorem presented below to overcome this issue.

    Theorem 3. Under Assumption 1, and by utilizing the protocol given by (2.7), we have the following constants 0=h0<h1<<hs<<hl,πi[0,1],σi>0,iIN and μ>0. By utilizing the stochastic sampled event-triggered strategy (2.3) and assuming the existence of positive matrices ˆP(m),ˆQs,ˆRs,ˆWs,ˆΦ, consensus of the multi-agent system (2.10) can be achieved in a mean square sense. This is subject to the condition that the matrices ˆK(m),mD and sIl satisfy the following inequality:

    (ˆΞ_(hmh(m1))ˆFT(m)ˆΣ_(m)ˆΨOˆX2)<0, (3.18)
    (¯ˆΞ(hmh(m1))ˆFT(m)¯ˆΣ(m)ˆΨOˆX2)<0, (3.19)

    where ˆΞ_=ˆ_1+ˆ2+ˆ3+ˆ4+ˆ5+ˆ6,

    ˆ_1=ˆFT(m)ε1+εT1ˆF(m)+εT1λ_mm(INˆP(m))ε1,

    ¯ˆΞ=¯ˆ1+ˆ2+ˆ3+ˆ4+ˆ5+ˆ6,

    ¯ˆ1=ˆFT(m)ε1+εT1ˆF(m)+εT1¯λmm(INˆP(m))ε1,

    ˆΣ_(m)=λ_(m)ˆX1(m), ¯ˆΣ(m)=ˉλ(m)ˆX1(m),

    λ_(m)=(λ_m1,λ_m2,,λ_mm1,λ_mm+1,,λ_mM,0,,0),

    ¯λ(m)=(¯λm1,¯λm2,,¯λmm1,¯λmm+1,,¯λmM,0,,0).

    The definitions of Theorem 2 are applicable to the remaining terms in the inequalities. By using the same strategy for proof as Theorem 2 of [35], the theorem may be simply constructed. Therefore, it is omitted here.

    Remark 6. It is worth mentioning that Theorem 3's conclusion is relatively conservative. To decrease conservativeness, the sojourn-time division method is used by dividing the sojourn time υ by J and denoting the pth segment as ¯λnm,p and λ_nm,p to represent the upper and lower bounds on the transmission probability, respectively. The conclusions drawn are relatively lenient.

    Corollary 1. Under Assumption 1 and the protocol (2.7), where 0=h0<h1<<hs<<hl, πi[0,1], σi>0,iIN and μ>0, the multi-agent system consensus (2.10) can be achieved in a mean-square sense with the stochastic sampled event-triggered strategy (2.3). This can be achieved if there exists a positive matrix ˆP(m), and positive matrices ˆQs,ˆRs,ˆWs,ˆΦ, along with the matrices ˆK(m),mD,sIl, which satisfy the following LMI:

    (ˆΞ_(m,p)(hmh(m1))ˆFT(m,p)ˆΣ_(m,p)ˆΨ(m,p)OˆX2(m,p))<0,
    (¯ˆΞ(m,p)(hmh(m1))ˆFT(m,p)¯ˆΣ(m,p)ˆΨ(m,p)OˆX2(m,p))<0,

    in which ˆΞ_(m,p),¯ˆΞ(m,p),ˆΣ_(m,p),¯ˆΣ(m,p),ˆF(m,p),ˆP(m,p), ˆX2(m,p),ˆΨ(m,p),ˆΥ1(m,p) and ˆΥ2(m,p) are similarly defined as in Theorem 3, with the exception that (m) is substituted by (m,p). Furthermore, the feedback gain is given by K(m,p)=ˆK(m,p)ˆP1(m,p). Moreover, the expression for the feedback gain is defined as K(m,p)=ˆK(m,p)ˆP1(m,p). The parameter matrix for the event-triggered strategy is denoted as Φ(m,p)=ˆP1(m,p)ˆΦ(m,p)ˆP1(m,p), where ˆΦ(m,p) is an estimated parameter and m is an element in the set D while p is an element in the index set IJ.

    Within this part, we will provide a numerical illustration to showcase the efficacy of the suggested design methodology. For consideration of a {semi-Markov jump multi-agent system}, which contains a leader and five followers, we assume that the model is described in formula (2.1). The coefficient matrices of the system equation are Ar, Br, Cr, r=1,2,3. A(1)=(14181128), A(2)=(16161523), A(3)=(13181120), B(1)=(172), B(2)=(612), B(3)=(68). The topology of the network is shown in Figure 1. The corresponding Laplacian matrix L and the leader adjacency matrix B can be derived in the manner shown as follows:

    L=(1100000000101000011000101), B=diag(1,1,0,0,0).

    Figure 1.  The leader and five followers topology.

    Let the event-trigged parameters be σ1=6.353,σ2=7.163,σ3=6.093,σ4=7.533 and σ5=6.312. The stochastic sampling period h takes values from the set {h1,h2}={0.1s,0.2s} with probabilities of occurrence π1=Pr{h=h1}=0.2 and θ2=Pr{h=h2}=0.8. With these values, we can obtain that ρ1=0.6, ρ2=0.4. By utilizing MATLAB's LMI toolbox, we can verify the feasibility of solutions to LMIs} (3.18, 3.19) for μ=4 in Theorem 3. The event-triggered parameter metrics are derived as follows: Φ1=(3.11590.01890.01893.1294), Φ2=(3.11530.01950.01963.1279), and Φ3=(3.11490.01950.01953.1279). The consensus feedback matrices are as follows: K(1)=(0.00580.0038), K(2)=(0.00270.0011), K(3)=(0.00480.0023). The initial states of the leader and followers were selected as follows: x0(0)=(10), x1(0)=(3.5479.553), x2(0)=(6.1545.902), x3(0)=(3.5947.611), x4(0)=(9.3417.841), x5(0)=(2.30112.24). The tracking errors between the leader and the followers are shown in Figures 2 and 3.

    Figure 2.  The system state e1.
    Figure 3.  The system state e2.

    Figure 4 illustrates the point in time at which an event is triggered. It indicates that the triggering of the event occurs at a lower frequency than the sampling rate.

    Figure 4.  Event-triggered instants for each agent.
    Figure 5.  The semi-Markov switching system state trajectory under the event-triggered protocol.

    Furthermore, the stochastic sampling period h is shown in Figure 6.

    Figure 6.  Stochastic sampling period h.

    We can compare our simulation example with that in [34]. Both our article and [34] share the same state equations. However, [34] adopted a static event-triggered protocol, while we added a process of stochastic sampling to its event-triggered protocol. Compared with the figure in [34], the multi-agent system controlled by stochastic sampling event-triggered control has a faster convergence rate and smaller steady-state error. Thus, this example validates the validity of Theorem 3.

    The paper presented a study on the mean-square consensus of a semi-Markov jump multi-agent system based on event-triggered stochastic sampling. We have proposed a novel approach to improve the efficiency of multi-agent systems for consensus control.

    The results of the study showed that the proposed approach was effective in achieving mean-square consensus in multi-agent systems. The use of event-triggering via stochastic sampling reduced the communication frequency and improved the computational efficiency of the system. The semi-Markov jump model provided a more accurate representation of the state transitions in the system.

    However, there are some limitations to this study. The numerical examples presented in the paper were relatively small, and it is unclear how the proposed approach would scale to larger multi-agent systems. Additionally, the study assumed perfect knowledge of the system parameters, which may not be the case in real-world scenarios.

    Future research can further investigate the robustness of the proposed approach against uncertainties and disturbances in the system. The scalability of the approach can also be explored in more details, and the approach can be tested on more complex multi-agent systems.

    The study presented in this article focuses on the use of multi-agent systems for the leader-follower consensus control topic. We have proposed a novel event-triggered stochastic sampling approach and investigated the use of a semi-Markov switching system architecture. We have also developed appropriate measures for mean-square consensus in multi-agent systems.

    The results of the numerical example presented in this study demonstrate the accuracy of the theoretical computations. The proposed approach has the potential to improve the efficiency of multi-agent systems for consensus control in various applications.

    The authors declare that they have not used artificial intelligence tools in the creation of this article.

    This work was supported by the College Students' Innovation and Entrepreneurship Program of the Ministry of Education (202211306068), Excellent Scientific Research and Innovation Team of Anhui Colleges (2022AH010098), Innovation and Entrepreneurship Training Program for College Students in Anhui Province (S202211306114, S202211306134), Quality Engineering Project of Chizhou University (2022XXSKC09), Chizhou University Introducing Doctoral Research Startup Project (CZ2022YJRC08), and Key Research Project of Chizhou University (CZ2021ZR03, CZ2023ZRZ04).

    The authors declare that there is no conflict of interest.



    [1] Jobbágy EG, Jackson RB (2000) The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecol Appl 10: 423-436. doi: 10.1890/1051-0761(2000)010[0423:TVDOSO]2.0.CO;2
    [2] Rustad LE, Hungtington TG, Boone RD (2000) Controls on soil respiration: implications for climate change. Biogeochem 48: 1-6. doi: 10.1023/A:1006255431298
    [3] Balesdent J, Basile-Doelsch I, Chadoeuf J, et al. (2018) Atmosphere-soil carbon transfer as a function of soil depth. Nature 559: 599-602. doi: 10.1038/s41586-018-0328-3
    [4] Malik AA, Puissant J, Buckeridge KM, et al. (2018) Land use driven change in soil pH affects microbial carbon cycling processes. Nat Commun 9: 3591. doi: 10.1038/s41467-018-05980-1
    [5] Trumbore SE, Chadwick OA, Amundson R (1996) Rapid exchange between soil carbon and atmospheric carbon dioxide driven by temperature change. Science 272: 393-396. doi: 10.1126/science.272.5260.393
    [6] Cao MK, Woodward FI (1998) Dynamic responses of terrestrial ecosystem carbon cycling to global climate change. Nature 393: 249-252. doi: 10.1038/30460
    [7] Knorr W, Prentice IC, House JI, et al. (2005) Long-term sensitivity of soil carbon turnover to warming. Nature 433: 298-301. doi: 10.1038/nature03226
    [8] Davidson EA, Janssens IA (2006) Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 440: 165-173. doi: 10.1038/nature04514
    [9] Crowther TW, Todd-Brown KEO, Rowe CW, et al. (2016) Quantifying global soil carbon losses in response to warming. Nature 540: 104-108. doi: 10.1038/nature20150
    [10] Teramoto M, Liang N, Takagi M, et al. (2016) Sustained acceleration of soil carbon decomposition observed in a 6-year warming experiment in a warm-temperate forest in southern Japan. Sci Rep 6: 35563. doi: 10.1038/srep35563
    [11] Hicks Pries CE, Castanha C, Porras RC, et al. (2017) The whole-soil carbon flux in response to warming. Science 355: 1420-1423. doi: 10.1126/science.aal1319
    [12] Melillo JM, Frey SD, DeAngelis KM, et al. (2017) Long-term pattern and magnitude of soil carbon feedback to the climate system in a warming world. Science 358: 101-105. doi: 10.1126/science.aan2874
    [13] Nottingham AT, Whitaker J, Ostle NJ, et al. (2019) Microbial responses to warming enhance soil carbon loss following translocation across a tropical forest elevation gradient. Ecol Lett.
    [14] Conant RT, Ryan MJ, Ågren GI, et al. (2011) Temperature and soil organic matter decomposition rates-synthesis of current knowledge and a way forward. Glob Change Biol 17: 3392-3404. doi: 10.1111/j.1365-2486.2011.02496.x
    [15] Cox PM, Betts RA, Jones CD, et al. (2000) Acceleration of global warming due to carbon-cycle feedbacks in a coupled climate model. Nature 408: 184-187. doi: 10.1038/35041539
    [16] Hanson PJ, Edwards NT, Garten CT, et al. (2000) Separating root and soil microbial contributions to soil respiration: a review of methods and observations. Biogeochem 48: 115-146. doi: 10.1023/A:1006244819642
    [17] O'Connell CS, Ruan L, Silver WL (2018) Drought drives rapid shifts in tropical rainforest soil biogeochemistry and greenhouse gas emissions. Nat Comm 9: 1348. doi: 10.1038/s41467-018-03352-3
    [18] Subke JA, Bahn M (2010) On the 'temperature sensitivity' of soil respiration: can we use the immeasurable to predict the unknown? Soil Biol Biochem 42: 1653-1656. doi: 10.1016/j.soilbio.2010.05.026
    [19] Shi Z, Crowell S, Luo Y, et al. (2018) Model structures amplify uncertainty in predicted soil carbon responses to climate change. Nat Comm 9: 2171. doi: 10.1038/s41467-018-04526-9
    [20] Potts M (1994) Desiccation tolerance for prokaryotes. Microbiol Mol Biol Rev 58: 755-805.
    [21] Doetterl S, Stevens A, Six J, et al. (2015) Soil carbon storage controlled by interactions between geochemistry and climate. Nat Geosci 8: 780-783. doi: 10.1038/ngeo2516
    [22] Doetterl S, Behre AA, Arnold C, et al. (2018) Links among warming, carbon and microbial dynamics mediated by soil mineral weathering. Nat Geosci 11: 589-593. doi: 10.1038/s41561-018-0168-7
    [23] Kuzyakov Y (2010) Priming effects: interactions between living and dead organic matter. Soil Biol Biochem 42: 1363-1371. doi: 10.1016/j.soilbio.2010.04.003
    [24] Kuzyakov Y, Friedel JK, Stahr K (2000) Review of mechanisms and quantification of priming effects. Soil Biol Biochem 32: 1485-1498. doi: 10.1016/S0038-0717(00)00084-5
    [25] Liang J, Zhou Z, Huo C, et al. (2018) More replenishment than priming loss of soil organic carbon with additional carbon input. Nat Comm 9: 3175. doi: 10.1038/s41467-018-05667-7
    [26] Birch HF (1958) The effect of soil drying on humus decomposition and nitrogen. Plant Soil 10: 9-31. doi: 10.1007/BF01343734
    [27] Schimel J, Balser TC, Wallenstein M (2007) Microbial stress-response physiology and its implications for ecosystem function. Ecology 88: 1386-1394. doi: 10.1890/06-0219
    [28] Wang L, Manzoni S, Ravi S, et al. (2015) Dynamic interactions of ecohydrological and biogeochemical processes in water-limited systems. Ecosphere 6: 133. doi: 10.1890/ES15-00122.1
    [29] Waring BG, Powers JS (2016) Unraveling the mechanisms underlying pulse dynamics of soil respiration in tropical dry forests. Environ Res Lett 11: 105005. doi: 10.1088/1748-9326/11/10/105005
    [30] Placella SA, Brodie EL, Firestone MK (2012) Rainfall-induced carbon dioxide pulses result from sequential resuscitation of phylogenetically clustered microbial groups. Proc Nat Acad Sci 109: 10931-10936. doi: 10.1073/pnas.1204306109
    [31] Gabriel C, Kellman L (2014) Investigating the role of moisture as an environmental constraint in the decomposition of shallow and deep mineral soil organic matter of a temperate coniferous soil. Soil Biol Biochem 68: 373-384. doi: 10.1016/j.soilbio.2013.10.009
    [32] Lu H, Liu S, Wang H, et al. (2017) Experimental throughfall reduction barely affects soil carbon dynamics in a warm-temperate oak forest, central China. Sci Rep 7: 15099. doi: 10.1038/s41598-017-15157-3
    [33] Bond-Lamberty B, Bailey VL, Chen M, et al. (2018) Globally rising soil heterotrophic respiration over recent decades. Nature 560: 80-83. doi: 10.1038/s41586-018-0358-x
    [34] Moyano FE, Vasilyeva NA, Bouckaert L, et al. (2012) The moisture response of soil heterotrophic respiration: interaction with soil properties. Biogeosciences 9: 1173-1182. doi: 10.5194/bg-9-1173-2012
    [35] Wang Q, He T, Liu J (2016) Litter input decreased the response of soil organic matter decomposition to warming in two subtropical forest soils. Sci Rep 6: 33814. doi: 10.1038/srep33814
    [36] Jung M, Reichstein M, Schwalm CR, et al. (2017) Compensatory water effects link yearly global land CO2 sink changes to temperature. Nature 541: 516-520. doi: 10.1038/nature20780
    [37] Craine JM, Fierer N, McLauchlan KK (2010) Widespread coupling between the rate and temperature sensitivity of organic matter decay. Nat Geosci 3: 854-857. doi: 10.1038/ngeo1009
    [38] Schmidt MWI, Torn MS, Abiven S, et al. (2011) Persistence of soil organic matter as an ecosystem property. Nature 478: 49-56. doi: 10.1038/nature10386
    [39] Barré P, Fernandez-Ugalde O, Virto I, et al. (2014) Impact of phyllosilicate mineralogy on organic carbon stabilization in soils: incomplete knowledge and exciting prospects. Geoderma 235-236: 382-395. doi: 10.1016/j.geoderma.2014.07.029
  • This article has been cited by:

    1. Jiangtao Dai, Ge Guo, A leader-following consensus of multi-agent systems with actuator saturation and semi-Markov switching topologies, 2024, 21, 1551-0018, 4908, 10.3934/mbe.2024217
  • Reader Comments
  • © 2019 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)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4202) PDF downloads(433) Cited by(1)

Figures and Tables

Figures(3)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog