
There is great interest to understand and reduce the massive carbon footprint of the concrete industry. Recent descriptions of microbes incidentally living inside concrete materials ("concrete endoliths") raised questions about how much carbon is either stored in or released from concrete by these microbes. We generated preliminary global estimates of how much organic carbon is stored within the living biomass of concrete endoliths (biomass-carbon) and much CO2 is released from respiring concrete endoliths. Between 2020–2022, we collected widely varying samples of Portland cement-based concrete from Lubbock, Texas. After quantifying endolith DNA from 25 concrete samples and estimating the current global mass of concrete, we calculated that the global concrete endolith biomass-carbon as low as 5191.9 metric tons (suggesting that endoliths are a negligible part of concrete's carbon profile) or as high as 1141542.3 tons (suggesting that concrete endoliths are a pool of carbon that could equal or offset some smaller sources of concrete-related carbon emissions). Additionally, we incubated concrete samples in air-tight microcosms and measured changes in the CO2 concentrations within those microcosms. Two out of the ten analyzed samples emitted small amounts of CO2 due to the endoliths. Thus, "concrete respiration" is possible, at least from concrete materials with abundant endolithic microbes. However, the remaining samples showed no reliable respiration signals, indicating that concrete structures often do not harbor enough metabolically active endoliths to cause CO2 emissions. These results are preliminary but show that endoliths may alter the carbon dynamics of solid concrete and, thus, the carbon footprint of the concrete industry.
Citation: Jordan Brown, Corona Chen, Deborah Carr. Endolithic microbes may alter the carbon profile of concrete[J]. AIMS Environmental Science, 2024, 11(2): 221-247. doi: 10.3934/environsci.2024011
[1] | Tahir Afzal, Muhammad Asim Afridi, Muhammad Naveed Jan . Integrating LSTM with Fama-French six factor model for predicting portfolio returns: Evidence from Shenzhen stock market China. Data Science in Finance and Economics, 2025, 5(2): 177-204. doi: 10.3934/DSFE.2025009 |
[2] | Xiaozheng Lin, Meiqing Wang, Choi-Hong Lai . A modification term for Black-Scholes model based on discrepancy calibrated with real market data. Data Science in Finance and Economics, 2021, 1(4): 313-326. doi: 10.3934/DSFE.2021017 |
[3] | Nitesha Dwarika . Asset pricing models in South Africa: A comparative of regression analysis and the Bayesian approach. Data Science in Finance and Economics, 2023, 3(1): 55-75. doi: 10.3934/DSFE.2023004 |
[4] | Yi Chen, Zhehao Huang . Measuring the effects of investor attention on China's stock returns. Data Science in Finance and Economics, 2021, 1(4): 327-344. doi: 10.3934/DSFE.2021018 |
[5] | Moses Khumalo, Hopolang Mashele, Modisane Seitshiro . Quantification of the stock market value at risk by using FIAPARCH, HYGARCH and FIGARCH models. Data Science in Finance and Economics, 2023, 3(4): 380-400. doi: 10.3934/DSFE.2023022 |
[6] | Zimei Huang, Zhenghui Li . What reflects investor sentiment? Empirical evidence from China. Data Science in Finance and Economics, 2021, 1(3): 235-252. doi: 10.3934/DSFE.2021013 |
[7] | Tsumbedzo Mashamba, Modisane Seitshiro, Isaac Takaidza . A comprehensive high pure momentum equity timing framework using the Kalman filter and ARIMA forecasting. Data Science in Finance and Economics, 2024, 4(4): 548-569. doi: 10.3934/DSFE.2024023 |
[8] | Kuo-Shing Chen . Interlinkages between Bitcoin, green financial assets, oil, and emerging stock markets. Data Science in Finance and Economics, 2024, 4(1): 160-187. doi: 10.3934/DSFE.2024006 |
[9] | Qian Shen, Yifan Zhang, Jiale Xiao, Xuhua Dong, Zifei Lin . Research of daily stock closing price prediction for new energy companies in China. Data Science in Finance and Economics, 2023, 3(1): 14-29. doi: 10.3934/DSFE.2023002 |
[10] | Katleho Makatjane, Ntebogang Moroke . Examining stylized facts and trends of FTSE/JSE TOP40: a parametric and Non-Parametric approach. Data Science in Finance and Economics, 2022, 2(3): 294-320. doi: 10.3934/DSFE.2022015 |
There is great interest to understand and reduce the massive carbon footprint of the concrete industry. Recent descriptions of microbes incidentally living inside concrete materials ("concrete endoliths") raised questions about how much carbon is either stored in or released from concrete by these microbes. We generated preliminary global estimates of how much organic carbon is stored within the living biomass of concrete endoliths (biomass-carbon) and much CO2 is released from respiring concrete endoliths. Between 2020–2022, we collected widely varying samples of Portland cement-based concrete from Lubbock, Texas. After quantifying endolith DNA from 25 concrete samples and estimating the current global mass of concrete, we calculated that the global concrete endolith biomass-carbon as low as 5191.9 metric tons (suggesting that endoliths are a negligible part of concrete's carbon profile) or as high as 1141542.3 tons (suggesting that concrete endoliths are a pool of carbon that could equal or offset some smaller sources of concrete-related carbon emissions). Additionally, we incubated concrete samples in air-tight microcosms and measured changes in the CO2 concentrations within those microcosms. Two out of the ten analyzed samples emitted small amounts of CO2 due to the endoliths. Thus, "concrete respiration" is possible, at least from concrete materials with abundant endolithic microbes. However, the remaining samples showed no reliable respiration signals, indicating that concrete structures often do not harbor enough metabolically active endoliths to cause CO2 emissions. These results are preliminary but show that endoliths may alter the carbon dynamics of solid concrete and, thus, the carbon footprint of the concrete industry.
This research focuses its interest on the mathematical modeling of the population dynamics of semelparous biological species. A species is called semelparous when it has a single reproductive episode before dying. Semelparity (sometimes called big-bang reproduction) occurs in very diverse biological species, see [1], including amphibians (e.g., Hyla frogs), arachnids (e.g., Pardosa licosidae spider, australian redback spider, desert spider, or black widow spider), fish (e.g., Pacific salmon, or sockeye salmon), insects (e.g., some butterflies, cicadas, or mayflies), mammals (e.g., some didelphids or dasyurid marsupials), mollusks (some squids or octopuses), reptiles (e.g., Labord chameleon, or some lizards), etc.
Methodologies for modeling dynamic biological systems, such as those based on population viability analysis (see [2,3]) or compartmental modeling (see [4,5]) usually require information about environmental variables, mortality rates, growth rates, etc. In practice, such information is difficult to obtain. In this work, we shall consider the methodology based on branching processes. These stochastic processes are appropriate mathematical models to describe the evolution of dynamical systems whose components, after a certain life period, reproduce and die in such a way that the transition from one to another state of the system is made according to a certain probability distribution. For theoretical concepts and applications about such types of processes, we refer the reader to some classical monographs[6,7,8], where several applications to cell kinetics, cell biology, chemotherapy, gene amplification, human evolution, and molecular biology are presented. See also the contributions, based on branching processes, by [9,10], in nuclear physics and complex contagion adoption dynamics, respectively.
In fact, branching processes are routinely used to describe the population dynamics of biological species with both asexual and sexual reproduction. We are especially interested in the mathematical modeling of the demographic dynamics of biological species with sexual reproduction. To this end, a fairly rich literature has emerged about discrete-time two-sex branching processes, see the surveys in [11,12] and the discussions therein. Most of these branching processes assume that all of the progenitor couples have a similar reproductive behavior, see [13,14,15]. It is also frequently assumed that mating and reproduction depend on the current number of progenitor couples existing in the population, see [16,17,18]. However, it is known that, due to various environmental factors, e.g., weather conditions, food supply, fertility parameters, or predators, in many biological species mating and reproduction occur in a non-predictable environment influenced by the current number of females and males in the population. For stochastic modeling about the demographic evolution of such species, two-sex branching processes had not been sufficiently developed. With such motivation, in [19], a new class of two-sex branching processes, which takes into account the possibility of various mating and reproduction strategies, both depending on the number of females and males in the population, was introduced. This class of processes is appropriate for the description of the demographic dynamics of semelparous species, which are characterized by having diverse behaviors in the mating and reproduction phases. In [20,21], some results about such a class of two-sex processes were established. The main purpose of this work is to continue this research line by providing new methodological (probabilistic and statistical) contributions.
The paper is organized as follows. In Section 2, the probability model is mathematically described and interpreted. In Section 3, some probabilistic results are provided. The probability distribution associated with the number of generations elapsed before the possible extinction of the population (time to the extinction) is determined. The class of processes under study is then used to mathematically model phenomena concerning to populate or repopulate habitats with endangered semelparous species. In Section 4, statistical results are derived. By considering maximum likelihood and Bayesian estimation methodologies, approximations for the main reproductive parameters involved in the probability model are proposed. To this purpose, information about the reproduction of the couples is incorporated. As illustration, a simulated example contextualized with a species of chameleons (Labord's chameleon) is presented. Concluding remarks and some questions for research are included in Section 5.
Let us consider the discrete-time two-sex branching process introduced in [19], denoted by {Xn}∞n=0, Xn=(Fn,Mn) representing the number of female and male individuals at generation n. The underlying probability model is described as follows, with N and N+ denoting the non-negative and the positive integers, respectively:
1) The mating phase is represented by a sequence of nm≥1 two arguments integer-valued functions {Ll}l∈Nm, Nm:={1,…,nm}. Each Ll is assumed to be non-decreasing and such that Ll(F,0)=Ll(0,M)=0, F,M∈N. At generation n, according to the lth mating strategy (function), Ll(Fn,Mn) couples female-male are formed.
2) The reproduction phase is modeled by a sequence of nr≥1 offspring probability distributions {Ph}h∈Nr, Nr:={1,…,nr}, Ph:={phk,s}(k,s)∈Sh, with Sh⊆N2 and phk,s being the probability for a given couple to produce exactly k females and s males, when Ph is the reproductive strategy.
3) In each generation, the mating and reproduction strategies are determined through suitable functions φm and φr, both defined on N2, taking values on Nm and Nr, respectively.
Initially, X0=x0∈N2+, and the number of couples originated Lφm(x0)(x0)>0. Given that Xn=x∈N2, it is then derived that, in the nth generation, Lφm(x) and Pφr(x) are the corresponding mating and reproductive strategies, respectively. Hence, at generation n+1,
Xn+1:=Lφm(x)(x)∑i=1(Fhn,i,Mhn,i), h=φr(x), n∈N, | (2.1) |
with Fhn,i and Mhn,i denoting, respectively, the number of female and male individuals originated by the ith couple at generation n, providing that the reproductive strategy h has been considered. For each h∈Nr, the random vectors (Fhn,i,Mhn,i), i=1,…,Lφm(x)(x), are assumed to be independent and identically distributed (i.i.d.) with offspring probability distribution Ph, i.e.,
P(Fhn,1=k,Mhn,1=s)=phk,s, (k,s)∈Sh. |
Remark 2.1 Functions Ll, φm, and φr should be flexible enough in order to fit the main features of the semelparous species we pretend to describe. Usually, such functions will depend of certain biological/ecological parameters of interest in the demographic dynamics of the species.
Remark 2.2 It is verified that {Xn}∞n=0 is a homogeneous Markov chain. In fact, given x0,…,xn,xn+1∈N2, taking into account that for each h∈Nr, independent of n, the random vectors (Fhn,i,Mhn,i), i=1,…,Lφm(x)(x), are i.i.d.,
P(Xn+1=xn+1∣X0=x0,…,Xn=xn)=P(Xn+1=xn+1∣Xn=xn)=P(Lφm(xn)(xn)∑i=1(Fφr(xn)n,i,Mφr(xn)n,i)=xn+1). |
Note that if, for some n≥1, Xn=(0,0), then Xn+j=(0,0), j≥1, thus (0,0) is an absorbing state. Consequently, in such a case, the population will become extinct. Let us denote by:
q(x0):=P(limn→∞Xn=0∣X0=x0), x0∈N2+, |
the extinction probability, associated with {Xn}∞n=0, when initially X0=x0. Assuming that, for each h∈Nr:
max{P(Fh1,1=0),P(Mh1,1=0)}>0, | (2.2) |
it has been proved in [20] the extinction-explosion property:
q(x0)+P(limn→∞g(Xn)=∞∣X0=x0)=1, | (2.3) |
where for F,M∈N,
g(F,M):=αF+βM, α≥0, β≥0, α+β>0. |
Also, by considering the rate:
mg(x):=g(x)−1E[g(Xn+1)∣Xn=x]=Ll(x)g(μh)g(x)−1, |
sufficient conditions for the extinction/survival of a biological population, described through the model (2.1), have been established in [20].
In particular, the following cases have a special biological significance: g(F,M)=F, g(F,M)=M, or g(F;M)=F+M, namely, the number of females, males, or total individuals in the population, respectively.
Remark 2.3 For each h∈Nr, let us denote by μh:=(μh1, μh2) and Σh=(σhij)i,j=1,2 the mean vector and the covariance matrix of (Fh1,1,Mh1,1), respectively, i.e.,
μhi:=∑(k1,k2)∈Shkiphk1,k2,i=1,2, | (2.4) |
σhij:=∑(k1,k2)∈Sh(ki−μhi)(kj−μhj)phk1,k2,i,j=1,2. | (2.5) |
Given x∈N2,
E[Xn+1∣Xn=x]=E[Ll(x)∑i=1(Fhn,i,Mhn,i)]=Ll(x)∑i=1E[(Fhn,i,Mhn,i)]=Ll(x)μh, |
and
Var[Xn+1∣Xn=x]=Var[Ll(x)∑i=1(Fhn,i,Mhn,i)]=Ll(x)∑i=1Var[(Fhn,i,Mhn,i)]=Ll(x)Σh, |
where l=φm(x) and h=φr(x).
The research about this class of two-sex branching processes will now continue investigating new probabilistic and statistical questions of biological/ecological interest.
From a probabilistic point of view, attention will be focused on two issues of special interest. First, assuming the possible extinction of the population in the habitat, we will look at the determination of the probability distribution associated with the number of generations elapsed before the possible extinction occurs. Second, we will investigate the application of the class of two-sex branching processes under study to the phenomenon of populating or repopulating habitats with endangered semelparous species. We will assume condition (2.2) and q(x0)∈(0,1). Note that, taking into account (2.3), the cases q(x0)=0 or q(x0)=1 mean the explosion or the extinction of the population, respectively.
Let us denote by Zn:=Lφm(Xn)(Xn) the number of couples formed at generation n. Clearly, if for some n≥1, Zn=0, then the population will be extinct. Given that X0=x0=(F0,M0)∈N2+, let us introduce the random variable:
T(x0):=sup{n≥0: Zn>0}, |
representing the number of generations elapsed before the extinction of the population occurs, when there were F0 females and M0 males in the population. Assuming that T(x0)<∞, the next result provides the probability distribution associated with T(x0), and also its main moments (mean and variance). In the case that the population goes extinct, such probability distribution shows how quickly the extinction will probably occur.
Let un(s):=E[sZn], 0≤s≤1, be the probability generating function of Zn. By simplicity, it will be denoted by u∗n(0):=un+1(0)−un(0), n∈N.
Theorem 3.1
(a) P(T(x0)=n∣T(x0)<∞)=q(x0)−1u∗n(0), n∈N.
(b) E[T(x0)∣T(x0)<∞]=∞∑n=1(1−q(x0)−1un(0)).
(c) Var[T(x0)∣T(x0)<∞]=q(x0)−1[∞∑n=1n2u∗n(0)−q(x0)−1(∞∑n=1(q(x0)−un(0))2].
Proof. First, note that P(T(x0)<∞)=q(x0).
(a) Using that Z0:=Lφm(x0)(x0)>0, it is derived that u0(0)=P(Z0=0)=0. Therefore,
P(T(x0)=0∣T(x0)<∞)=P(T(x0)<∞)−1P(T(x0)=0)=q(x0)−1P(Z1=0)=q(x0)−1u1(0)=q(x0)−1(u1(0)−u0(0))=q(x0)−1u∗0(0). |
Now, for n∈N+,
P(T(x0)=n∣T(x0)<∞)=q(x0)−1P(T(x0)=n)=q(x0)−1(P(T(x0)≤n)−P(T(x0)≤n−1))=q(x0)−1(P(Zn+1=0)−P(Zn=0))=q(x0)−1(un+1(0)−un(0))=q(x0)−1u∗n(0). |
(b) Using that T(x0) is a non-negative random variable, it is derived, see [22, pg. 84], that
E[T(x0)∣T(x0)<∞]=∞∑n=0(1−P(T(x0)≤n∣T(x0)<∞)). |
Hence, using that P(T(x0)≤n)=P(Zn+1=0)=un+1(0),
E[T(x0)∣T(x0)<∞]=∞∑n=0(1−P(T(x0)≤n∣T(x0)<∞))=∞∑n=0(1−q(x0)−1P(T(x0)≤n))=∞∑n=0(1−q(x0)−1un+1(0))=∞∑n=1(1−q(x0)−1un(0)). |
(c) The result is derived, from Theorem 3.1(b), using that:
Var[T(x0)∣T(x0)<∞]=E[T(x0)2∣T(x0)<∞]−E[T(x0)∣T(x0)<∞]2, |
and taking into account that:
E[T(x0)2∣T(x0)<∞]=∞∑n=1n2P(T(x0)=n∣T(x0)<∞)=q(x0)−1∞∑n=1n2(P(T(x0)≤n)−P(T(x0)≤n−1))=q(x0)−1∞∑n=1n2(P(Zn+1=0)−P(Zn=0))=q(x0)−1∞∑n=1n2(un+1(0)−un(0))=q(x0)−1∞∑n=1n2u∗n(0). | (3.1) |
Let us consider a certain habitat in which a semelparous species has become extinct, or it is in serious danger of extinction. The purpose is to populate or repopulate the habitat with such species. To this end, of ecological significance, several attempts to repopulate the species in the habitat will probably be necessary. Next, the class of two-sex branching processes {Xn}∞n=0, defined in (2.1), will be used as the mathematical model.
In fact, let {X(j)n}∞n=0, j∈N+, be independent processes where, for each j∈N+, {X(j)n}∞n=0 is a two-sex branching process, like the one defined in (2.1), which describes the population dynamics concerning the jth attempt of repopulating. Thus, vector X(j)n=(F(j)n,M(j)n) represents the number of female and male individuals in the habitat at generation n, in the jth attempt of repopulating. All processes {X(j)n}∞n=0, j∈N+, have the same mating and reproductive strategies, i.e., they have the same sequences {Ll}l∈Nm and {Ph}h∈Nr. Assume condition (2.2) holds and q(x0)∈(0,1).
Initially (attempt j=1), it is assumed that X(1)0=x0=(F0,M0)∈N2+, i.e., F0 females and M0 males of the species under consideration are introduced in the habitat. If, after a certain number of generations, the population becomes extinct, then it will be necessary to restart the repopulating, by introducing again (attempt j=2) F0 females and M0 males in the habitat, i.e., X(2)0=x0, and so on. This iterative procedure continues until a sufficient number of females and males are achieved in the habitat, so that the risk of extinction disappears (for simplicity, this fact will be referred to as implementation of the species in the habitat). It is deduced that
P(limn→∞X(j)n=0∣X(j)0=x0)=q(x0), j∈N+. |
Let T(j)(x0):=sup{n≥0: Z(j)n>0}, Z(j)n:=Lφm(Xn)(j))(Xn(j)), j∈N+, namely, the time to the extinction at the jth attempt of repopulation. It is derived that {T(j)(x0)}∞j=1 is a sequence of i.i.d. random variables with the probability distribution given in Theorem 3.1(a).
Let us denote by N(x0) the number of attempts until the implementation of the species occurs. It is then verified that N(x0) is distributed according to a geometric law with parameter 1−q(x0). Hence,
P(N(x0)=n)=q(x0)n(1−q(x0)), n∈N, | (3.2) |
and
E[N(x0)]=q(x0)(1−q(x0))−1 and Var[N(x0)]=q(x0)(1−q(x0))−2. | (3.3) |
Finally, let us consider the variable:
T∗(x0):=N(x0)∑j=1T(j)(x0), | (3.4) |
representing the total number of generations elapsed until the implementation of the species in the habitat occurs. The next result establishes the probability distribution of T∗(x0) and its main moments.
Theorem 3.2
(a) P(T∗(x0)=n)=q∗(x0)−1(δn,0+∞∑j=1q(x0)jPn(x0)(∗j)), n∈N.
(b) E[T∗(x0)]=q∗(x0)∞∑n=1(q(x0)−un(0)).
(c) Var[T∗(x0)]=q∗(x0)[(1+q∗(x0))∞∑n=1n2u∗n(0)−q(x0)−1(∞∑n=1(q(x0)−un(0)))2],
where
q∗(x0):=(1−q(x0))−1, δn,0:=1 if n=0 or 0 if n≠0, |
Pn(x0):=P(T(1)(x0)=n∣T(1)(x0)<∞), Pn(x0)(∗j):=∑i1+…+ij=nPi1(x0)...Pij(x0). |
Proof.
(a) Taking into account (3.2) and (3.4),
P(T∗(x0)=n)=∞∑j=0P(T∗(x0)=n∣N(x0)=j)P(N(x0)=j)=q∗(x0)−1(δn,0+∞∑j=1q(x0)jP(j∑l=1T(l)(x0)=n))=q∗(x0)−1(δn,0+∞∑j=1q(x0)jPn(x0)(∗j) ), n∈N. |
(b) From (3.3), (3.4), and Theorem 3.1(b),
E[T∗(x0)]=E[N(x0)]E[T(1)(x0)∣T(1)(x0)<∞]=q∗(x0)∞∑n=1(q(x0)−un(0)). |
(c) Using (3.3), (3.4), and Theorem 3.1(b), (c), we can obtain that:
Var[T∗(x0)]=E[N(x0)]Var[T(1)(x0)∣T(1)(x0)<∞]+Var[N(x0)]E[(T(1)(x0))2∣T(1)(x0)<∞]. |
Remark 3.1 By simplicity in the underlying probabilistic development, it has been assumed that in the proposed iterative procedure, for each j∈N+, X(j)0=x0=(F0,M0), that is, each iteration begins with F0 females and M0 males in the habitat. In practice, this assumption is not a problem. It implies that the extinction probability q(x0) is the same for all iterations. Therefore, it is deduced that N(x0) is distributed according to a geometric law with parameter 1−q(x0).
With the aim to check some possible changes in the demographic dynamics of the semelparous species, it is important to determine accurate approximations for the main statistical parameters (offspring means, variances, and covariances) specified in (2.4) and (2.5), respectively. The application of estimation methodologies based on population viability analysis requires having information on various variables related to the biological species under consideration. In practice, the information about such variables is difficult to obtain. In this scenario, estimation based in Bayesian methodology provides a reasonable solution. In fact, by considering a parametric statistical context about the reproductive strategies Ph, h∈Nr, estimates for such statistical parameters have been determined in [19]. More recently, by using approximate Bayesian computation techniques, estimates have also been proposed in [21]. The application of such techniques requires a large number of simulations from the mathematical model. Consequently, it involves a significant computational effort. In this section, we shall consider the more general nonparametric framework about the reproductive strategies Ph, h∈Nr, that is, no functional form is assumed for the probabilities phk,s, (k,s)∈Sh. By using maximum likelihood and Bayesian estimation methodologies, we will determine approximations for such probabilities and also for μhi, σhij, i,j=1,2.
Let {Xn}∞n=0 be the two-sex branching process described in the previous section. Assume, for a given n, the following observed information over time in a set of generations, denoted by Gn, up to the nth generation is reached:
Sn:={X0=(F0,M0), Zi,(k,s), (k,s)∈Sφr(Xi), i∈Gn}, | (4.1) |
where
Zi,(k,s):=Lφm(Xi)(Xi)∑j=1I{(Fφr(Xi)i,j,Mφr(Xi)i,j)=(k,s)}, | (4.2) |
with IA denoting the indicator function of A. The variable Zi,(k,s) represents the total number of couples at generation i, which produce exactly k females and s males, (k,s)∈⋃h∈NrSh.
For h∈Nr, let Ghn:={i∈Gn: φr(Xi)=h}, i.e., the set of observed generations of Gn where Ph has been the reproductive strategy. Clearly,
⋃h∈NrGhn=Gn,Ghn∩Gh′n=∅, h,h′∈Nr, h≠h′. |
In what follows, it will be assumed that Ghn≠∅, h∈Nr. From (4.2), let us introduce the variables:
Vhn,(k,s):=∑i∈GhnZi,(k,s) , h∈Nr. |
Note that Vhn,(k,s) represents the total number of couples in the observed generations that have originated exactly k females and s males, with Ph being the underlying reproductive strategy.
According to this estimation methodology, taking into account the sample information given in (4.1) and (4.2), we have to determine, for h∈Nr, the values of phk,s,(k,s)∈Sh that maximize the corresponding likelihood function.
Theorem 4.1 The maximum likelihood estimates for phk,s, μhi, and σhij, i,j=1,2, h∈Nr, are given by:
(a) ^phk,s=(Vhn)−1Vhn,(k,s), (k,s)∈Sh.
(b) ^μhi=(Vhn)−1∑(k1,k2)∈ShkiVhn,(k1,k2).
(c) ^σhij=(Vhn(1+Vhn))−1[Vhn∑(k1,k2)∈ShkikjVhn,(k1,k2)−∑(k1,k2),(l1,l2)∈ShkiljVhn,(k1,k2)Vhn,(l1,l2)],
where Vhn:=∑(k,s)∈ShVhn,(k,s) is assumed to be positive.
Proof. Taking into account (4.1), the corresponding likelihood function is given by:
L(P1,…,Pnr∣Sn)=nr∏h=1∏(k,s)∈Sh(phk,s)Vhn,(k,s). | (4.3) |
From (4.3), the log-likelihood function is then deduced:
ℓ(P1,…,Pnr∣Sn):=log(L(P1,…,Pnr∣Sn))=nr∑h=1∑(k,s)∈ShVhn,(k,s)log(phk,s). | (4.4) |
In order to determine the maximum likelihood estimates for phk,s, (k,s)∈Sh, it is necessary to obtain the values that maximize the likelihood function (4.3), or equivalently, the log-likelihood function (4.4), subject to the constraints:
∑(k,s)∈Shphk,s=1,phk,s≥0,h=1,…,nr. |
Consequently, using the Lagrange multipliers technique, it is required to determine the values of phk,s that maximize the function:
Ψ(P1,…,Pnr):=ℓ(P1,…,Pnr∣Sn)+λh(1−∑(k,s)∈Shphk,s), |
where λh, h=1,…,nr denotes the corresponding Lagrange multipliers. Now, it is obtained:
∂∂phk,s[ℓ(P1,…,Pnr∣Sn)+λh(1−∑(k,s)∈Shphk,s)]=(phk,s)−1Vhn,(k,s)−λh=0. | (4.5) |
From (4.5), it is deduced that λh=Vhn. Hence, the solutions of the equations given in (4.5) are the expressions given in Theorem 4.1(a). It can be checked that such solutions maximize function (4.4). Intuitively, notice that ^phk,s is the proportion of progenitor couples that generate exactly k females and s males when Ph is the underlying reproductive strategy.
Estimates given in Theorem 4.1(b), (c), are then determined using (2.4), (2.5), and Theorem 4.1(a).
First, in order to apply the Bayesian estimation methodology, it is necessary to choose a suitable class of prior densities π(P1,…,Pnr) for (P1,…,Pnr). From the mathematical expression given in (4.3), it is derived that an appropriate class of prior densities for (P1,…,Pnr) is the product of Dirichlet densities:
π(P1,…,Pnr)=nr∏h=1Dh∏(k,s)∈Sh(phk,s)τhk,s−1, | (4.6) |
where for each h∈Nr,
τh=(τhk,s, (k,s)∈Sh) is a vector of positive constants,
Dh=∏(k,s)∈ShΓ(τh∗)(Γ(τhk,s))−1, τh∗=∑(k,s)∈Shτhk,s, Γ(u):=∫∞0e−xxu−1dx, u>0.
Theorem 4.2 By considering the squared error loss function, the Bayesian estimates for phk,l, μhi, and σhij, i,j=1,2, h∈Nr, are given by:
(a) ~phk,s=(Wh∗)−1Whk,s, (k,s)∈Sh.
(b) ~μhi=(Wh∗)−1∑(k1,k2)∈ShkiWhk1,k2.
(c) ~σhij=(Wh∗(1+Wh∗))−1(Wh∗∑(k1,k2)∈ShkikjWhk1,k2−∑(k1,k2),(l1,l2)∈ShkiljWhk1,k2Whl1,l2),
where
Whk,s:=τhk,s+Vhn,(k,s),Wh∗:=∑(k,s)∈ShWhk,s. |
Proof. Taking into account (4.6), the posterior density of (P1,…,Pnr), incorporating the information provided by (4.1), is given by the product of Dirichlet densities:
π(P1,…,Pnr∣Sn)=nr∏h=1D∗h∏(k,s)∈Sh(phk,s)Whk,s−1, | (4.7) |
where for each h∈Nr,
Wh=(Whk,s, (k,s)∈Sh),D∗h=∏(k,s)∈ShΓ(Wh∗)(Γ(Whk,s))−1. |
Given (k,l)∈Sh, from (4.7) we deduce, as marginal posterior density for phk,s, the Beta density:
π(phk,s∣Sn)=Γ(Wh∗)Γ(Whk,s)Γ(Wh∗−Whk,s)(phk,s)Whk,s−1(1−phk,s)Wh∗−Whk,s−1, phk,s∈(0,1). | (4.8) |
Finally, from (4.8), using the squared error loss function, it is obtained as Bayesian estimate for phk,s:
~phk,s:=∫10phk,lπ(phk,s∣Sn)dphk,s=(Wh∗)−1Whk,s. | (4.9) |
From (2.4), (2.5), and (4.9), the expressions given in Theorem 4.2(b), (c), are deduced as Bayesian estimates for μhi and σhij, i,j=1,2, h∈Nr, respectively.
Remark 4.1 From (4.8), we can determine the highest posterior density (HPD) credibility sets for phk,s, namely, sets of the form:
J(Q):={phk,s: π(phk,s∣Sn)≥Q}, Q>0, |
where given a certain credibility coefficient 1−α, the constant Q is calculated taking into account that:
∫J(Q)π(phk,s∣Sn)dphk,s=1−α. |
To determine HPD credibility sets for μhi and for σhij, i,j=1,2, h∈Nr, it is necessary to approximate the corresponding posterior densities for such parameters. An appropriate procedure is based on the simulation, from π(phk,s∣Sn), of a sufficiently large number of values of phk,s. Then, from (2.4) and (2.5), the corresponding values for μhi and for σhij, i,j=1,2, h∈Nr, are determined. From these values, applying a suitable Gaussian kernel method, see e.g., [23], approximated posterior densities for such parameters can be obtained. Using such approximated densities, we can then determine the corresponding empirical HPD credibility sets.
Remark 4.2 In order to simulate data from model (2.1), to calculate the proposed estimates for the parameters, and to determine the corresponding HPD credibility sets, we have implemented some specific programs using the statistical software R, see [24].
Labord's chameleon (Furcifer labordi) is a native reptile of southwestern Madagascar, where it usually lives in dry deciduous forests. It is considered the shortest-lived tetrapod animal. It spends most of its life in the developing embryo phase (8 to 9 months) and, after that, it experiences a very rapid growth. It has a short lifespan (4 to 5 months), reaching its sexual maturity at an early age (2 months). There are no rigorous studies in the specialized scientific literature about their social organization or on their mating and reproduction strategies. The few studies that have been carried out are based on monitoring experiments through radio telemetry. From the information recorded, it has been detected that females exhibit high habitat fidelity, moving small cumulative and linear distances with low dispersion rates. Males move greater distances, in a less predictable manner, with higher dispersal rates than females. This species of chameleon constitutes a particular example of semelparous life (progenitors die shortly after reproducing). Their mating and reproductive strategies, highly conditioned by the number of females and males in the habitat, must adapt to these temporal limitations existing in intense competition and fighting between males. It has been suggested that they possess a wide range of different mating systems, generally polygamous matings. Males can mate with more than one female and females can mate with different males during the same ovarian cycle (the dynamics of male color change could also affect the choice of partner). The female lays a clutch of eggs and the progenitor male and female die. Some studies reported that females can lay between 6 and 8 eggs. Due to various random factors, mainly predators and environmental factors, a high percentage of eggs will not hatch. See [25,26,27,28] for more information about this species of chameleon.
Unfortunately, there is no real data available on the demographic dynamics of this reptile species. Next, taking into account the special characteristics of this species of chameleon, a simulated example is presented where mating and reproduction strategies close to reality are assumed. Let us consider a population of chamaleon of Labord with the following population dynamics:
1) Females and males form couples according to the nm=2 mating strategies:
L1(F,M)=min{F,M}, L2(F,M)=Fmin{1,M}, F,M∈N. |
According to L1, the females and males practice fidelity, and they are allowed to have at most one mate. According to L2, in each generation, a dominant male mates with each female. The other males do not participate in the mating process.
2) Since the number of eggs laid by a female is between 6 and 8, and it is known that a high percentage of eggs do not hatch, it will be assumed in our simulation the nr=2 reproductive strategies Ph={phk,s}(k,s)∈Sh, h=1,2, where
S1=S2={(0,0),(3,3),(4,3),(3,4),(4,4),(5,3),(3,5)}, |
p10,0=0.71, p13,3=0.08, p14,3=0.06, p13,4=0.05, p14,4=0.03, p15,3=0.03, p13,5=0.04, |
p20,0=0.71, p23,3=0.08, p24,3=0.05, p23,4=0.06, p24,4=0.03, p25,3=0.04, p23,5=0.03. |
From P1 and P2,
μ11=1.02, μ12=1.03, σ111=2.679, σ122=2.749, σ112=2.519, |
μ21=1.03, μ22=1.02, σ211=2.749, σ222=2.679, σ212=2.519. |
P1 slightly favors the birth of males, with a ratio of males/females of the means equal to 1.01. This ratio has a value of 0.99 for P2, which slightly favors the birth of females. From such reproductive strategies, it is deduced that:
max{P(Fh0,1=0),P(Mh0,1=0)}=0.71, h=1,2. |
Hence, condition (2.2) holds.
3) In each generation, it is assumed that the functions φm and φr are given by:
φm(F,M):=1⋅I{M(F+M)−1>K}+2⋅I{M(F+M)−1≤K}, φr(F,M):=1⋅I{F≤M}+2⋅I{F>M}, F,M∈N. |
K<1 is a suitable threshold for the proportion of males in the population.
Let {Xn}∞n=0 be the two-sex branching process described in (2.1) with the mating and reproduction strategies mentioned previously. By way of illustration, taking X0=(10,10) and K=0.8, we have developed a simulation for the first n=100 generations of such a chameleon Labord population, i.e., in this case, G100={1,2,…,100}, see Figure 1.
According to the obtained simulation, it is derived that P2 has been the reproductive strategy in 19 generations and P1 has been the reproductive strategy in 81 generations. In fact,
G2100={2,3,5,6,8,9,10,12,17,18,20,24,26,33,34,36,52,54,56},G1100=G100−G2100. |
Taking into account the special reproductive characteristics of this biological species, with a high probability of unsuccessful hatching, we have considered the prior density given in (4.7), where
τ1=τ2=(29, 3.5, 3.5, 3.5, 3.5, 3.5, 3.5). |
Using the information provided in the simulation, from Theorems 4.1 and 4.2, we have determined the corresponding maximum likelihood and Bayesian estimates for phk,l, μhi, and σhij, i,j,h=1,2, see Tables 1 and 2.
p10,0 | p13,3 | p14,3 | p13,4 | p14,4 | p15,3 | p13,5 | |
True values | 0.71 | 0.08 | 0.06 | 0.05 | 0.03 | 0.03 | 0.04 |
Maximum likelihood estimates | 0.712 | 0.079 | 0.061 | 0.049 | 0.028 | 0.029 | 0.041 |
Bayesian estimates | 0.712 | 0.079 | 0.061 | 0.049 | 0.028 | 0.029 | 0.041 |
μ11 | μ12 | σ111 | σ122 | σ112 | |||
True values | 1.02 | 1.03 | 2.679 | 2.749 | 2.519 | ||
Maximum likelihood estimates | 1.012 | 1.023 | 2.664 | 2.738 | 2.505 | ||
Bayesian estimates | 1.013 | 1.024 | 2.666 | 2.740 | 2.507 |
p20,0 | p23,3 | p24,3 | p23,4 | p24,4 | p25,3 | p23,5 | |
True values | 0.71 | 0.08 | 0.05 | 0.06 | 0.03 | 0.04 | 0.03 |
Maximum likelihood estimates | 0.706 | 0.072 | 0.051 | 0.068 | 0.028 | 0.049 | 0.026 |
Bayesian estimates | 0.700 | 0.072 | 0.052 | 0.068 | 0.030 | 0.049 | 0.028 |
μ21 | μ22 | σ211 | σ222 | σ212 | |||
True values | 1.03 | 1.02 | 2.749 | 2.679 | 2.519 | ||
Maximum likelihood estimates | 1.038 | 1.030 | 2.857 | 2.670 | 2.556 | ||
Bayesian estimates | 1.081 | 1.054 | 2.898 | 2.722 | 2.595 |
maxh=1,2{max(k,l)∈Sh|phk,l−^phk,l|}=0.009,maxh=1,2{maxi=1,2|μhi−^μhi|}=0.009,maxh=1,2{maxi,j=1,2|σhij−^σhij|}=0.108, |
maxh=1,2{max(k,l)∈Sh|phk,l−~phk,l|}=0.008,maxh=1,2{maxi=1,2|μhi−~μhi|}=0.051,maxh=1,2{maxi,j=1,2|σhij−~σhij|}=0.149, |
which shows good accuracy of the estimators.
As an illustration, Figures 2 and 3 show the evolution of the estimates for ph0,0 and ph3,5, in the successive generations belonging to Gh100, h=1,2, respectively. Similarly, Figures 4 and 5 show the evolution of the estimates for μ1i and μ2i, i=1,2, in the successive generations belonging to G1100 and G2100, respectively. The 95% HPD credibility sets are also included in all figures (the true values for the parameters are represented by horizontal lines). For a better visualization of the graphs, the generations belonging to G1100 and G2100 have been renumbered, from 1 to 81 and from 1 to 19, respectively.
In this research, we have continued the working line started in previous papers about mathematical modeling concerning the demographic dynamics of semelparous biological species through two-sex branching processes with multiple mating and reproduction strategies. It has been assumed that, in each generation, both phases (mating and reproduction) are influenced by the current number of females and males existing in the population (habitat). Several probabilistic and statistical contributions have been established. In fact, the probability distribution associated to the variable number of generations elapsed before the possible extinction of the population has been derived. Also, the aforementioned two-sex branching processes have been used to mathematically describe the phenomenon of populating or repopulating habitats with a semelparous species that has become extinct, or is in danger of extinction. By considering the most general non-parametric statistical setting on the reproductive strategies associated with the biological species, various inferential questions about the main parameters governing the reproduction phase have been determined. To this purpose, information concerning the reproduction of the couples has been included in the observed data sample. Using this information, estimates for the reproductive parameters have been proposed. To this end, maximum likelihood and Bayesian estimates, and the corresponding 95% HPD credibility sets, have been determined. To simulate data from the mathematical model and also to calculate the estimates and the HPD credibility sets, the necessary computing programs have been developed. By way of illustration, the proposed estimation methodologies have been applied through a simulated example with Labord's chamaleon populations. The simulation performed has showed the accuracy of the proposed estimates.
Some possible direcions for research are, for example, to extend this class of two-sex branching processes, including in the probability model the immigration of females, males, or couples, from external populations; considering mating and/or reproduction in random environments; or assuming multi-type populations. It is also necessary to explore other possible mathematical methodologies that allow modeling the evolution of semelparous biological species with sexual reproduction. In this regard, the survey in [29] provides interesting information to be considered.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors wish to express their gratitude to the referees. Their technical questions, suggestions, and corrections have improved the presentation and understanding of this paper. This research has been supported by the Ministerio de Ciencia e Innovación of Spain (grant PID2019-108211GB-I00/AEI/10.13039/501100011033).
The authors declare there is no conflict of interest.
[1] |
Head MJ, Steffen W, Fagerlind D, et al. (2021) The Great Acceleration is real and provides a quantitative basis for the proposed Anthropocene Series/Epoch. Episodes 45: 359–376. https://doi.org/10.18814/epiiugs/2021/021031 doi: 10.18814/epiiugs/2021/021031
![]() |
[2] |
Phillips S (2014) Review of Concrete and culture: a material history. J Soc Archit Hist 73: 417–419. https://doi.org/10.1525/jsah.2014.73.3.417 doi: 10.1525/jsah.2014.73.3.417
![]() |
[3] |
Elhacham E, Ben-Uri L, Grozovski J, et al. (2020) Global human-made mass exceeds all living biomass. Nature 588: 442–444. https://doi.org/10.1038/s41586-020-3010-5 doi: 10.1038/s41586-020-3010-5
![]() |
[4] |
Monteiro PJM, Miller SA, Horvath A (2017) Towards sustainable concrete. Nat Mater 16: 698–699. https://doi.org/10.1038/nmat4930 doi: 10.1038/nmat4930
![]() |
[5] | Meyer C (2005) Concrete as a green building material. Proceedings of ConMat '05 10. |
[6] |
Van Damme H (2018) Concrete material science: Past, present, and future innovations. Cem Concr Res 112: 5–24. https://doi.org/10.1016/j.cemconres.2018.05.002 doi: 10.1016/j.cemconres.2018.05.002
![]() |
[7] |
Miatto A, Schandl H, Fishman T, et al. (2016) Global patterns and trends for non-metallic minerals used for construction. J Ind Ecol 21: 924–937. https://doi.org/10.1111/jiec.12471 doi: 10.1111/jiec.12471
![]() |
[8] |
Liew KM, Sojobi AO, Zhang LW (2017) Green concrete: Prospects and challenges. Constr Build Mater 156: 1063–1095. https://doi.org/10.1016/j.conbuildmat.2017.09.008 doi: 10.1016/j.conbuildmat.2017.09.008
![]() |
[9] | Pilon BS, Tyner JS, Yoder DC, et al. (2019) The effect of pervious concrete on water quality parameters: A case study. MDPI Water 11. http://dx.doi.org/10.3390/w11020263 |
[10] | Halauniova A (2022) Good and bad concrete. Fugitive modern and the aesthetics of renovation in Poland. City https://doi.org/10.1080/13604813.2021.2019490 |
[11] | The Cement Sustainability Initiative (2016) Cement Industry Energy and CO2 Performance. World Business Council for Sustainable Development. Available from: https://www.wbcsd.org/Sector-Projects/Cement-Sustainability-Initiative/Resources/Cement-Industry-Energy-and-CO2-Performance |
[12] | Internation Energy Agency (2018) Technology roadmap: low-carbon transition in the cement industry. Available from: https://www.iea.org/reports/technology-roadmap-low-carbon-transition-in-the-cement-industry |
[13] |
Barcelo L, Kline J, Walenta G, et al. (2014) Cement and carbon emissions. Mater Struct 47: 1055–1065. https://doi.org/10.1617/s11527-013-0114-5 doi: 10.1617/s11527-013-0114-5
![]() |
[14] | Hammond G, Jones C (2011) Embodied Carbon. The inventory of carbon and energy (ICE). A Building Services Research Information Association guide. |
[15] |
Naqi A, Jang JG (2019) Recent progress in green cement technology utilizing low-carbon emission fuels and raw materials: a review. MDPI Sustainability 11: 537. https://doi.org/10.3390/su11020537 doi: 10.3390/su11020537
![]() |
[16] | Marinkovic SB (2013) Life cycle assessment of (LCA) aspects of concrete, In: Eco-efficient concrete. Pacheco-Torgal F, Jalali S, Labrincha J, et al. Eds., United Kingdom: Woodhead Publishing, 45–80. https://doi.org/10.1533/9780857098993.1.45 |
[17] | Xi F, Davis SJ, Ciais P, et al. (2016) Substaintial global carbon uptake by cement carbonation. Nat. Geosci. 9: 880–887. https://doi.org/10.1038/NGEO2840 |
[18] |
Possan E, Thomaz WA, Aleandri GA, et al. (2017) CO2 uptake potential due to concrete carbonation: A case study. Case Stud Constr Mater 6: 147–161. http://dx.doi.org/10.1016/j.cscm.2017.01.007 doi: 10.1016/j.cscm.2017.01.007
![]() |
[19] | Gajda J (2001) Absorption of atmospheric carbon dioxide by portland cement concrete. Portland Cement Association |
[20] |
Pade C, Guimaraes M (2007) The CO2 uptake of concrete in a 100 year perspective. Cem Concr Res 37: 1348–1356. https://doi.org/10.1016/j.cemconres.2007.06.009 doi: 10.1016/j.cemconres.2007.06.009
![]() |
[21] |
Maresca JA, Moser P, Schumacher T (2017) Analysis of bacterial communities in and on concrete. Mater Struct 50: 25. https://doi.org/10.1617/s11527-016-0929-y doi: 10.1617/s11527-016-0929-y
![]() |
[22] |
Golubic S, Friedmann E, Schneider J (1981) The lithobiontic ecological niche, with special reference to microorganisms. J Sediment Res 51: 475–478. https://doi.org/10.1306/212F7CB6-2B24-11D7-8648000102C1865D doi: 10.1306/212F7CB6-2B24-11D7-8648000102C1865D
![]() |
[23] |
Colwell FS, D'Hondt S (2013) Nature and extent of the deep biosphere. Rev Mineral Geochem 75: 547–574. https://doi.org/10.2138/rmg.2013.75.17 doi: 10.2138/rmg.2013.75.17
![]() |
[24] |
Santelli CM, Banergee N, Bach W, et al. (2010) Tapping the subsurface ocean crust biosphere: Low biomass and drilling-related contamination calls for improved quality controls. Geomicrobiol J 27: 158–169. https://doi.org/10.1080/01490450903456780 doi: 10.1080/01490450903456780
![]() |
[25] |
Kallmeyer J, Pockalny R, Adhikari RR, et al. (2012) Global distribution of microbial abundance and biomass in subseafloor sediment. PNAS 109: 16213–16216. https://doi.org/10.1073/pnas.1203849109 doi: 10.1073/pnas.1203849109
![]() |
[26] |
Fordyce AJ, Ainsworth TD, Leggat W (2021) Light capture, skeletal morphology, and the biomass of corals' boring endoliths. mSphere 6: e00060–21. https://doi.org/10.1128/mSphere.00060-21 doi: 10.1128/mSphere.00060-21
![]() |
[27] | Matthes-Sears U, Gerrath JA, Larson DW (1997) Abundance, biomass, and productivity of endolithic and epilithic lower plants on the temperatre-zone cliffs of the Niagara Escarpment, Canada. Int J Plant Sci 158: 451–460. |
[28] | Sun HJ, Friedmann EI (1999) Growth on geological time scales in the Antarctic cryptoendolithic microbial community. Geomicrobiol J 16: 193–202. |
[29] |
Tyler NA, Ziolkowski LA (2021) Endolithic microbial carbon cycling in East Antarctica. Astrobiology 21: 165–176. https://doi.org/10.1089/ast.2019.2109 doi: 10.1089/ast.2019.2109
![]() |
[30] |
Fang J, Zhang L (2011) Exploring the deep biosphere. Sci China Earth Sci 54: 157–165. https://doi.org/10.1007/s11430-010-4148-z doi: 10.1007/s11430-010-4148-z
![]() |
[31] |
Brown J, Chen C, Fernandez M, et al. (2023) Urban endoliths: incidental microbial communities occurring inside concrete. AIMS Microbiol 9: 277–312. https://doi.org/10.3934/microbiol.2023016 doi: 10.3934/microbiol.2023016
![]() |
[32] |
Kiledal EA, Keffer JL, Maresca JA (2021) Bacterial communities in concrete reflect its composite nature and change with weathering. mSystems 6: e01153–20. https://doi.org/10.1128/mSystems.01153-20 doi: 10.1128/mSystems.01153-20
![]() |
[33] |
Coombes MA, Naylor LA, Thompson RC, et al. (2011) Colonization and weathering of engineering materials by marine microorganisms: an SEM study. Earth Surf Proc Land 36: 582–593. https://doi.org/10.1002/esp.2076 doi: 10.1002/esp.2076
![]() |
[34] | Kiledal EA, Shaw M, Polson SW, et al. (2023) Metagenomic analyis of concrete bridge reveals a microbial communitiy dominated by halophilic Bacteria and Archaea. Microbiol Spectr 11. https://doi.org/10.1128/spectrum.05112-22 |
[35] |
Kaur P, Singh V, Arora A (2022) Microbial concrete—a sustainable solution for concrete construction. Appl Biochem Biotechnol 194: 1401–1416. https://doi.org/10.1007/s12010-021-03604-x doi: 10.1007/s12010-021-03604-x
![]() |
[36] |
Al-Mansour A, Chow CL, Feo L, et al. (2019) Green concrete: By-products utilization and advanced approaches. MDPI Sustainability 11: 5145. https://doi.org/10.3390/su11195145 doi: 10.3390/su11195145
![]() |
[37] | United States Climate Data (2022) Lubbock, Texas climate data. Available from: https://www.usclimatedata.com/climate/lubbock/texas/united-states/ustx2745 |
[38] | Kiledal EA, Maresca JA (2021) DNA extraction from concrete V.2. Available from: https://doi.org/10.17504/protocols.io.b2i5qcg6 |
[39] |
Semenov M, Blagodatskaya E, Stepanov A, et al. (2018) DNA-based determination of soil microbial biomass in alkaline and carbonaceous soils of semi-arid climate. J Arid Environ 150: 54–61. https://doi.org/10.1016/j.jaridenv.2017.11.013 doi: 10.1016/j.jaridenv.2017.11.013
![]() |
[40] |
Anderson TH, Martens R (2013) DNA determinations during growth of soil microbial biomasses. Soil Biol Biochem 57: 487–495. http://dx.doi.org/10.1016/j.soilbio.2012.09.031 doi: 10.1016/j.soilbio.2012.09.031
![]() |
[41] |
Fornasier F, Ascher J, Ceccherini MT, et al. (2014) A simplified rapid, low-cost and versatile DNA-based assessment of soil microbial biomass. Ecol Indic 45: 75–82. http://dx.doi.org/10.1016/j.ecolind.2014.03.028 doi: 10.1016/j.ecolind.2014.03.028
![]() |
[42] | Gong H, Du Q, Xie S, et al. (2021). Soil microbial DNA concentration is a powerful indicator for estimating soil microbial biomass C and N across arid and semi-arid regions in northern China. Appl Soil Ecol 103869. https://doi.org/10.1016/j.apsoil.2020.103869 |
[43] | United States Geological Survey (2023) Cement Production Worldwide from 1995 to 2022 (in Billion Metric Tons). Available from: https://www.statista.com/statistics/1087115/global-cement-production-volume. |
[44] |
Waters CN, Zalasiewicz J (2018) Concrete: the most abundant novel rock type of the Anthropocene. Encyclopedia Anthrop 1: 75–85. https://doi.org/10.1016/B978-0-12-809665-9.09775-5 doi: 10.1016/B978-0-12-809665-9.09775-5
![]() |
[45] |
Moriarty R, O'Brien TD (2013) Distribution of mesozooplankton biomass in the global ocean. Earth Syst Sci Data 5: 45–55. https://doi.org/10.5194/essd-5-45-2013 doi: 10.5194/essd-5-45-2013
![]() |
[46] |
Cho SU, Chae CU (2016) A study on life cycle CO2 emissions of low-carbon building in South Korea. MDPI Sustainability 8: 579. https://doi.org/10.3390/su8060579 doi: 10.3390/su8060579
![]() |
[47] |
Llatas C, Quinones R, Bizcocho N (2022) Environmental impact assessment of construction waste recycling versus disposal scenarios using an LCA-BIM tool during the design stage. MDPI Recycling 7: 82. https://doi.org/10.3390/recycling7060082 doi: 10.3390/recycling7060082
![]() |
[48] |
Wierzchos J, de los Ríos A, Ascaso C (2012) Microorganisms in desert rocks: the edge of life on Earth. Int Microbiol 15: 171–181. https://doi.org/10.2436/20.1501.01.170 doi: 10.2436/20.1501.01.170
![]() |
[49] |
Bone JR, Stafford R, Hall AE, et al. (2022) The intrinsic primary bioreceptivity of concrete in the coastal environment – A review. Dev Built Environ 10: 100078. https://doi.org/10.1016/j.dibe.2022.100078 doi: 10.1016/j.dibe.2022.100078
![]() |
[50] |
Marshall D, Meuller C, Clifford B (2020) Computational arrangement of demolition debris. Detritus 11: 3–18. https://doi.org/10.31025/2611-4135/2020.13967 doi: 10.31025/2611-4135/2020.13967
![]() |
[51] |
Convey P, Key RS, Key RJD (2010) The establishment of a new ecological guild of pollinating insects of sub-Antarctic South Georgia. Antarct Sci 22: 508–512. https://doi.org/10.1017/S095410201000057X doi: 10.1017/S095410201000057X
![]() |
[52] |
Lebouvier M, Lambret P, Garnier A, et al. (2020) Spotlight on the invasion of a carabid beetle on an oceanic island over a 105-year period. Sci Rep 10: 17103. https://doi.org/10.1038/s41598-020-72754-5 doi: 10.1038/s41598-020-72754-5
![]() |
[53] |
Gjørv OE (2011) Durability of Concrete Structures. Arab J Sci Eng 36: 151–172. https://doi.org/10.1007/s13369-010-0033-5 doi: 10.1007/s13369-010-0033-5
![]() |
[54] |
Giannantonio DJ, Kurth JC, Kurtis KE, et al. (2009) Molecular characterizations of microbial communities fouling painted and unpainted concrete structures. Int Biodeterior Biodegradation 63: 30–40. https://doi.org/10.1016/j.ibiod.2008.06.004 doi: 10.1016/j.ibiod.2008.06.004
![]() |
[55] |
Hayek M, Slagues M, Souch J, et al. (2023) How to improve the bioreceptivity of concrete infrastructure used in marine ecosystems? Literature review for mechanisms, key factors, and colonization effects. J Coast Res 39: 553–568. https://doi.org/10.2112/JCOASTRES-D-21-00158.1 doi: 10.2112/JCOASTRES-D-21-00158.1
![]() |
[56] |
Wiscovitch-Russo R, Singh H, Oldfield LM, et al. (2022) An optimized approach for processing of frozen lung and lavage samples for microbiome studies. PLoS One 17: e0265891. https://doi.org/10.1371/journal.pone.0265891 doi: 10.1371/journal.pone.0265891
![]() |
[57] |
Joergensen RG, Emmerling C (2006) Methods for evaluating human impact on soil microorganisms based on their activity, biomass, and diversity in agricultural soils. J Plant Nutr Soil Sci 169: 295–309. https://doi.org/10.1002/jpln.200521941 doi: 10.1002/jpln.200521941
![]() |
[58] |
Rotbart N, Borisover M, Bukhanovsky N, et al. (2017) Examination of residual chloroform interference in the measurement of microbial biomass C by fumigation-extraction. Soil Biol Biochem 111: 60–65. http://dx.doi.org/10.1016/j.soilbio.2017.03.018 doi: 10.1016/j.soilbio.2017.03.018
![]() |
[59] | Zelles L, Palojarvi A, Kandeler E, et al. (1997) Changes in soil microbial properties and phospholipid fatty acid fractions after chloroform fumigation. Soil Biol Biochem 29: 1325–1336. |
[60] | Ross DJ (1987) Soil microbial biomass estimated by the fumigation-incubation procedure: Seasonal fluctuations and influence of soil moisture content. Soil Biol Biochem 19: 397–404. |
[61] | Krausmann F, Wiedenhofer D, Lauk C, et al. (2017) Global socioeconomic material stocks rise 23-fold over the 20th century and require half of annual resource use. PNAS 114: 1880–1885. www.pnas.org/cgi/doi/10.1073/pnas.1613773114 |
[62] |
Muller DB, Liu G, Løvik AN, et al. (2013) Carbon emissions of infrastructure development. Environ Sci Technol 47: 11739–11746. http://dx.doi.org/10.1021/es402618m doi: 10.1021/es402618m
![]() |
[63] |
Krausmann F, Lauk C, Haas W, et al. (2018) From resource extraction to outflows of wastes and emissions: The socioeconomic metabolism of the global economy, 1900–2015. Glob Environ Change 52: 131–140. https://doi.org/10.1016/j.gloenvcha.2018.07.003 doi: 10.1016/j.gloenvcha.2018.07.003
![]() |
[64] |
Nedeljkovic M, Visser J, Savija B, et al. (2021) Use of fine recycled concrete aggregates in concrete: A critical review. J Build Eng 38: 102196. https://doi.org/10.1016/j.jobe.2021.102196 doi: 10.1016/j.jobe.2021.102196
![]() |
[65] |
McNeil K, Kang THK (2013) Recycled concrete aggregates: A review. Int J Concr Struct 7: 61–69. https://doi.org/10.1007/s40069-013-0032-5 doi: 10.1007/s40069-013-0032-5
![]() |
[66] |
Hashimoto S, Tanikawa H, Moriguchi Y (2007) Where will large amounts of materials accumulated within the economy go? – A material flow analysis of construction minerals for Japan. Waste Manage 27: 1725–1738. https://doi.org/10.1016/j.wasman.2006.10.009 doi: 10.1016/j.wasman.2006.10.009
![]() |
[67] | Carlson CA, Del Giorgio PA, Herndl GJ (2007) Microbes and the dissipation of energy and respiration: From cells to ecosystems. Oceanography 20: 89–100. |
[68] |
Konopka MC, Strovas TJ, Ojala DS, et al. (2011) Respiration response imaging for real-time detection of microbial function at the single-cell level. Appl Environ Microbiol 77: 67–72. https://doi.org/10.1128/AEM.01166-10 doi: 10.1128/AEM.01166-10
![]() |
[69] |
Ziolkowski LA, Mykytczuk NCS, Omelon CR, et al. (2013) Arctic gypsum endoliths: a biogeochemical characterization of a viable and active microbial community. Biogeosciences 10: 7661–7675. https://doi.org/10.5194/bg-10-7661-2013 doi: 10.5194/bg-10-7661-2013
![]() |
[70] | Davila AF, Hawes I, Araya JG, et al. (2015) In situ metabolism in halite endolithic microbial communities of the hyperarid Atacama Desert. Front Microbiol 6 https://doi.org/10.3389/fmicb.2015.01035 |
[71] | Mergelov N, Mueller CW, Prater I, et al. (2018) Alteration of rocks by endolithic organisms is one of the pathways for the beginning of soils on Earth. Sci Rep 8. https://doi.org/10.1038/s41598-018-21682-6 |
[72] | Sajjad W, Ilahi N, Kang S, et al. (2022) Endolithic microbes of rocks, their community, function and survival strategies. Int Biodeterior. Biodegradation 169. https://doi.org/10.1016/j.ibiod.2022.105387 |
[73] | Rochette P, Hutchinson GL (2005) Measurement of soil respiration in situ: Chamber techniques. In: Hatfield JL, Baker JM. Authors, Micrometerology in Agricultural Systems, Agron. Monogr. 47., Madison, WI: ASA, CSSA, and SSSA, 247–286. https://doi.org/10.2134/agronmonogr47.c12 |
[74] |
Midwood AJ, Millard P (2011) Challenges in measuring the δ13C of the soil surface CO2 efflux. Rapid Commun Mass Spectrom 25: 232–242. https://doi.org/10.1002/rcm.4857 doi: 10.1002/rcm.4857
![]() |
[75] | MacFayden A (1973) Inhibitory effects of carbon dioxide on microbial activity in soil. Pedobiologia 13: 140–149. |
[76] |
Shestak CJ, Busse MD (2005) Compaction alters physical but not biological indices of soil health. Soil Sci Soc Am J 69: 236–246. https://doi.org/10.2136/sssaj2005.0236 doi: 10.2136/sssaj2005.0236
![]() |
[77] |
Frouz J, Luděk B (2018) Flow of CO2 from soil may not correspond with CO2 concentration in soil. Sci Rep 8: 10099. https://doi.org/10.1038/s41598-018-28225-z doi: 10.1038/s41598-018-28225-z
![]() |
[78] | Neville AM (2011) Properties of concrete, 5th edition. Essex, United Kingdom: Pearson Education Limited. |
[79] |
Dauti D, Tengattini A, Dal Pont S, et al. (2018) Analysis of moisture migration in concrete at high temperature through in-situ neutron tomography. Cem Concr Res 111: 41–55. https://doi.org/10.1016/j.cemconres.2018.06.010 doi: 10.1016/j.cemconres.2018.06.010
![]() |
[80] |
Goordial J, Davila A, Lacelle D, et al. (2016) Nearing the cold-arid limits of microbial life in permafrost of an upper dry valley, Antarctica. ISME J 10: 1613–1624 https://doi.org/10.1038/ismej.2015.239 doi: 10.1038/ismej.2015.239
![]() |
[81] |
Levin GV, Straat PA (2016) The case for extant life on Mars and its possible detection by the Viking Labeled Release Experiment. Astrobiology 16: 798–810. https://doi.org/10.1089/ast.2015.1464 doi: 10.1089/ast.2015.1464
![]() |
[82] | Xie J, Li Y, Zhai C, et al. (2009) CO2 absorption by alkaline soils and its implication to the global carbon cycle. Environ Geol 56: 953–961. |
[83] | Tang B, Fan M, Yang Z, et al. (2023) A comparison study of aggregate carbonation and concrete carbonation for the enhancement of recycled aggregate pervious concrete. Constr Build Mater 371. https://doi.org/10.1016/j.conbuildmat.2023.130797 |
[84] |
Šantrůčková H, Bird MI, Elhottová D, et al. (2005) Heterotrophic fixation of CO2 in soil. Microb Ecol 49: 218–225. https://doi.org/10.1007/s00248-004-0164-x doi: 10.1007/s00248-004-0164-x
![]() |
[85] |
Kuzyakov Y (2006) Sources of CO2 efflux from soil and review of partitioning methods. Soil Biol Biochem 38: 425–448. https://doi.org/10.1016/j.soilbio.2005.08.020 doi: 10.1016/j.soilbio.2005.08.020
![]() |
[86] | Mehta PK (2010) Sustainable cements and concrete for the climate change era – A review. Proceedings of the Second International Conference on Sustainable Construction Materials and Technologies, Ancona, Italy. |
[87] |
Shi J, Xu Y (2006) Estimation and forecasting of concrete debris amount in China. Resour Conserv Recycl 49: 147–158. https://doi.org/10.1016/j.resconrec.2006.03.011 doi: 10.1016/j.resconrec.2006.03.011
![]() |
[88] | Bahmani H, Mostafaei H, Ghiassi B, et al. (2023) A comparitive study of calcium hydroxide, calcium oxide, calcined dolomite, and metasilicate as activators for slag-based HPC. Structures 58. https://doi.org/10.1016/j.istruc.2023.105653 |
[89] | Mostafaei H, Bahmani H, Mostofinejad D, et al. (2023) A novel development of HPC without cement: Mechanical properties and sustainability evalution. J Build Eng 76. https://doi.org/10.1016/j.jobe.2023.107262 |
[90] | Rowell A, Ghebrab T, Jeter R (2023) Bacterial treatment of recycled concrete aggregate. MDPI Recycling 8. https://doi.org/10.3390/recycling8050068 |
1. | Fangzhou Huang, Jiao Song, Nick J. Taylor, The impact of business conditions and commodity market on US stock returns: An asset pricing modelling experiment, 2022, 6, 2573-0134, 433, 10.3934/QFE.2022019 |
p10,0 | p13,3 | p14,3 | p13,4 | p14,4 | p15,3 | p13,5 | |
True values | 0.71 | 0.08 | 0.06 | 0.05 | 0.03 | 0.03 | 0.04 |
Maximum likelihood estimates | 0.712 | 0.079 | 0.061 | 0.049 | 0.028 | 0.029 | 0.041 |
Bayesian estimates | 0.712 | 0.079 | 0.061 | 0.049 | 0.028 | 0.029 | 0.041 |
μ11 | μ12 | σ111 | σ122 | σ112 | |||
True values | 1.02 | 1.03 | 2.679 | 2.749 | 2.519 | ||
Maximum likelihood estimates | 1.012 | 1.023 | 2.664 | 2.738 | 2.505 | ||
Bayesian estimates | 1.013 | 1.024 | 2.666 | 2.740 | 2.507 |
p20,0 | p23,3 | p24,3 | p23,4 | p24,4 | p25,3 | p23,5 | |
True values | 0.71 | 0.08 | 0.05 | 0.06 | 0.03 | 0.04 | 0.03 |
Maximum likelihood estimates | 0.706 | 0.072 | 0.051 | 0.068 | 0.028 | 0.049 | 0.026 |
Bayesian estimates | 0.700 | 0.072 | 0.052 | 0.068 | 0.030 | 0.049 | 0.028 |
μ21 | μ22 | σ211 | σ222 | σ212 | |||
True values | 1.03 | 1.02 | 2.749 | 2.679 | 2.519 | ||
Maximum likelihood estimates | 1.038 | 1.030 | 2.857 | 2.670 | 2.556 | ||
Bayesian estimates | 1.081 | 1.054 | 2.898 | 2.722 | 2.595 |
p10,0 | p13,3 | p14,3 | p13,4 | p14,4 | p15,3 | p13,5 | |
True values | 0.71 | 0.08 | 0.06 | 0.05 | 0.03 | 0.03 | 0.04 |
Maximum likelihood estimates | 0.712 | 0.079 | 0.061 | 0.049 | 0.028 | 0.029 | 0.041 |
Bayesian estimates | 0.712 | 0.079 | 0.061 | 0.049 | 0.028 | 0.029 | 0.041 |
μ11 | μ12 | σ111 | σ122 | σ112 | |||
True values | 1.02 | 1.03 | 2.679 | 2.749 | 2.519 | ||
Maximum likelihood estimates | 1.012 | 1.023 | 2.664 | 2.738 | 2.505 | ||
Bayesian estimates | 1.013 | 1.024 | 2.666 | 2.740 | 2.507 |
p20,0 | p23,3 | p24,3 | p23,4 | p24,4 | p25,3 | p23,5 | |
True values | 0.71 | 0.08 | 0.05 | 0.06 | 0.03 | 0.04 | 0.03 |
Maximum likelihood estimates | 0.706 | 0.072 | 0.051 | 0.068 | 0.028 | 0.049 | 0.026 |
Bayesian estimates | 0.700 | 0.072 | 0.052 | 0.068 | 0.030 | 0.049 | 0.028 |
μ21 | μ22 | σ211 | σ222 | σ212 | |||
True values | 1.03 | 1.02 | 2.749 | 2.679 | 2.519 | ||
Maximum likelihood estimates | 1.038 | 1.030 | 2.857 | 2.670 | 2.556 | ||
Bayesian estimates | 1.081 | 1.054 | 2.898 | 2.722 | 2.595 |