Research article

Estimation and inversion across the spectrum of carbon cycle modeling

  • Received: 04 December 2017 Accepted: 02 May 2018 Published: 08 May 2018
  • Understanding of the carbon cycle is particularly important because of the role of carbon dioxide as a greenhouse gas. Carbon cycle models play an essential role in the interpretation of observational data. The analysis of the carbon cycle involves statistical estimation in various contexts. These include various types of model calibration, including the estimation of feedbacks. A range of inverse calculations are involved in estimating the spatial and/or temporal dependence of carbon dioxide sources and sinks, given observations of concentrations. The uncertainties in these estimates propagate into uncertainties in projections of future carbon cycle behavior. These disparate analyses are discussed in terms of a modeling spectrum that runs from empirical statistical models through to reductionist mechanistic models. The use of the modeling spectrum allows a comparison of di erent modeling approaches. Comparing di erent levels of modeling can provide a basis for assessing the extent to which estimation is being applied consistently.

    Citation: Ian Enting. Estimation and inversion across the spectrum of carbon cycle modeling[J]. AIMS Geosciences, 2018, 4(2): 126-143. doi: 10.3934/geosci.2018.2.126

    Related Papers:

    [1] Chunqiang Cui, Jin Chen, Shixun Lin . Metric and strong metric dimension in TI-power graphs of finite groups. AIMS Mathematics, 2025, 10(1): 705-720. doi: 10.3934/math.2025032
    [2] Saba Al-Kaseasbeh, Ahmad Erfanian . A bipartite graph associated to elements and cosets of subgroups of a finite group. AIMS Mathematics, 2021, 6(10): 10395-10404. doi: 10.3934/math.2021603
    [3] Adeel Farooq, Musawwar Hussain, Muhammad Yousaf, Ahmad N. Al-Kenani . A new algorithm to compute fuzzy subgroups of a finite group. AIMS Mathematics, 2023, 8(9): 20802-20814. doi: 10.3934/math.20231060
    [4] Madeleine Al Tahan, Sarka Hoskova-Mayerova, B. Davvaz, A. Sonea . On subpolygroup commutativity degree of finite polygroups. AIMS Mathematics, 2023, 8(10): 23786-23799. doi: 10.3934/math.20231211
    [5] Zehra Velioǧlu, Mukaddes Balçik . On the Tame automorphisms of differential polynomial algebras. AIMS Mathematics, 2020, 5(4): 3547-3555. doi: 10.3934/math.2020230
    [6] Aneeza Imtiaz, Hanan Alolaiyan, Umer Shuaib, Abdul Razaq, Jia-Bao Liu . Applications of conjunctive complex fuzzy subgroups to Sylow theory. AIMS Mathematics, 2024, 9(1): 38-54. doi: 10.3934/math.2024003
    [7] Aman Ullah, Muhammad Ibrahim, Tareq Saeed . Fuzzy cosets in AG-groups. AIMS Mathematics, 2022, 7(3): 3321-3344. doi: 10.3934/math.2022185
    [8] Huaguo Shi, Zhangjia Han, Pengfei Guo . A complete classification of weakly Dedekind groups. AIMS Mathematics, 2024, 9(4): 7955-7972. doi: 10.3934/math.2024387
    [9] Shitian Liu . Finite groups for which all proper subgroups have consecutive character degrees. AIMS Mathematics, 2023, 8(3): 5745-5762. doi: 10.3934/math.2023289
    [10] Fawaz Aseeri, Julian Kaspczyk . The conjugacy diameters of non-abelian finite p-groups with cyclic maximal subgroups. AIMS Mathematics, 2024, 9(5): 10734-10755. doi: 10.3934/math.2024524
  • Understanding of the carbon cycle is particularly important because of the role of carbon dioxide as a greenhouse gas. Carbon cycle models play an essential role in the interpretation of observational data. The analysis of the carbon cycle involves statistical estimation in various contexts. These include various types of model calibration, including the estimation of feedbacks. A range of inverse calculations are involved in estimating the spatial and/or temporal dependence of carbon dioxide sources and sinks, given observations of concentrations. The uncertainties in these estimates propagate into uncertainties in projections of future carbon cycle behavior. These disparate analyses are discussed in terms of a modeling spectrum that runs from empirical statistical models through to reductionist mechanistic models. The use of the modeling spectrum allows a comparison of di erent modeling approaches. Comparing di erent levels of modeling can provide a basis for assessing the extent to which estimation is being applied consistently.


    Associating graphs to algebraic structures and studying their algebraic properties employing the methods in graph theory has been an interesting topic for mathematicians in the last decades and continuously arousing people's wide attention. For instance, Cayley graphs of groups [14], power graphs of groups ([2,15,24]), zero-divisor graphs of semigroups [17], zero-divisor graphs of commutative rings [5], total graphs of rings [1], unit graphs of rings [6], zero-divisor graphs of modules [9]. Bosák initiated the study of the intersection graphs of semigroups in [10]. Later, Csákány and Pollák defined the intersection graph of subgroups of a finite group and studied the graphic properties in [16]. In the recent years, several interesting properties of the intersection graphs of subgroups of groups have been obtained in the literature. In 1975, Zelinka investigated the intersection graphs of subgroups of finite abelian groups [33]. Akbari, Heydari and Maghasedi obtained some properties of the intersection graphs of subgroups of solvable groups in [7]. Shen classified finite groups with disconnected intersection graphs of subgroups in [29]. Ma gave an upper bound for the diameter of intersection graphs of finite simple groups in [23]. Rajkumar and Devi defined the intersection graph of cyclic subgroups of groups and obtained the clique number of intersection graph of subgroups of some non-abelian groups [26,27]. Recall that, for a finite group G, the intersection graph of subgroups of G, denoted by Γ(G), is a graph with all non-trivial subgroups of G as its vertices and two distinct vertices H and K are adjacent in Γ(G) if and only if HK1.

    The thickness of a graph, as a topological invariant of a graph, was first defined by Tutte in [30]. The problem is the same as determining whether K9 is biplanar or not, that is, a union of two planar graphs and it was showed in [11]. Tutte generalized the problem by defining the concept of the thickness of a graph. It is an important research object in topological graph theory, and it also has important applications to VLSI design and network design (see [4,28]). Outerthickness seems to be studied first in Geller's unpublished manuscript, where it was shown that θ0(K7) is 3 by similar exhaustive search as in the case of the thickness of K9. It is known that determining the thickness of a graph is NP-hard [22]. Until now, there are only a few results on the thickness of graphs. The attention has been focused on finding upper bounds of thickness and outerthickness. In 1964, Beineke, Frank and Moon obtained the thickness of the complete bipartite graph [13]. Beineke and Harary obtained the thickness of the complete graph ([3,12]). Chartrand, Geller and Hedetniemi conjectured that planar graph has an edge partition into two outerplanar graphs in [21]. Guy and Nowwkowski obtained the outerthickness of the complete graph and complete bipartite graph in [19,20]. In [18], Ding et al. proved that every planar graph has an edge partition into two outerplanar graphs.

    For the planarity of intersection graph of subgroups of a finite group, Ahmadi and Taeri, Kayacan and Yaraneri independently classified all finite groups with planar intersection graphs in [8,25], respectively. By their main theorem, it is easy to obtain the planarity of intersection graph of subgroups of a finite abelian group.

    Theorem 1.1. [8,25] Let G be a finite abelian group. Then Γ(G) is planar if and only if G is one of the following groups: Zpα(α=2,3,4,5),Zpαq(α=1,2),Zpqr,ZpZp,Z2Z2p, where p,q,r are distinct primes.

    Motivated by the above theorem, the aim of this paper is to classify finite abelian groups whose thickness and outerthickness of intersection graphs of subgroups are 1 and 2, respectively. We also investigate the thickness and outerthickness of intersection graphs of subgroups of some finite non-abelian groups in the last section.

    Let us recall some basic definitions and notations in graph theory. Let Γ be a graph with vertex set V(Γ) and edge set E(Γ), and all graphs in this paper are simple graphs, that are, undirect, no loops and no multiple edges. We use Kn to denote the complete graph on n vertices and Km,n to denote the complete bipartite graph with one partition consisting of m vertices and another partition consisting of n vertices. The degree of a vertex v in Γ, denoted by degΓ(v), is the number of vertices incident to v. If degΓ(v)=k, then the vertex is denoted by vk. A subgraph Γ of Γ is provided V(Γ)V(Γ) and E(Γ)E(Γ). A clique of a graph Γ is a complete subgraph of Γ. A maximal clique of a graph is a clique such that adding any vertex into the clique cannot be a clique. The clique number of a graph Γ is the number of vertices in a maximal clique in Γ with maximum size and it is denoted by ω(Γ). A subgraph H of Γ is a spanning subgraph of Γ if V(H)=V(Γ). The union of graphs Γ1 and Γ2 is the graph Γ1Γ2 with vertex set V(Γ1)V(Γ2) and edge set E(Γ1)E(Γ2). Let Γ1+Γ2 be a graph obtained from the union Γ1Γ2 by adding edges between each vertex in Γ1 and each vertex in Γ2. If Γ2 has exactly one vertex v, we shall use Γ1+v instead of Γ1+Γ2. A graph is planar if it can be drawn in a plane such that no two edges intersect expect at their end vertices. An outerplanar graph is a planar graph that can be embedded in the plane without crossings in such a way that all the vertices lie in the boundary of the unbounded face of the embedding. Recall that the thickness of a graph Γ, denoted by θ(Γ), is the minimum number of planar subgraphs whose union is Γ. Similarly, the outerthickness θ0(Γ) is obtained where the planar subgraphs are replaced by outerplanar subgraphs in the previous definition. Obviously, the thickness of a planar graph is 1 and the outerthickness of an outerplanar graph is 1.

    The following lemmas are obvious and are used frequently later.

    Lemma 2.1. [32,Lemma 2.1] If Γ1 is a subgraph of Γ2, then θ(Γ1)θ(Γ2) and θ0(Γ1)θ0(Γ2).

    Lemma 2.2. [12,19,3] (1) θ(Kn)={n+76n9,103n=9,10.

    (2) θ0(Kn)={n+14n73n=7.

    where x denotes the least integer greater than or equal to x, x denotes the largest integer greater than or equal to x

    By Lemma 2.2, one can see that θ(Kn)3 for n9, θ0(Kn)3 for n7.

    From Euler's formula, a planar graph Γ has at most 3|V|6 edges and an outerpalnar graph Γ has at most 2|V|3 edges, furthermore, one can get the lower bound for the thickness and outerthickness of a graph as follows.

    Lemma 2.3. [13] If Γ=(V,E) is a graph with |V|=n and |E|=m, then θ(Γ)m3n6 and θ0(Γ)m2n3.

    Lemma 2.4. Let Γ be a graph. If θ(Γ)=k, then θ(Γ+vk)=k. In particular, if Γ is a complete graph, then θ(Γ+vt)=k(1t2k).

    Proof. If Γ is a planar graph, then Γ+v1 is clearly a planar graph. If Γ is not a planar graph, then we may assume that {Γ1,Γ2,,Γk} is a planar decomposition of Γ, where each Γi is a spanning subgraph of Γ, then Γi+v1 are planar graphs. Hence {Γ1+v1,Γ2+v1,,Γk+v1} is a decomposition of Γ+vk. Thus, θ(Γ+vk)k. By Lemma 2.1, θ(Γ+vk)θ(Γ)=k and thus θ(Γ+vk)=k. In particular, if Γ is a complete graph, then θ(Γ+vt)=k(1t2k).

    Lemma 2.5. [31] A graph is outerplanar if and only if it has no subgraph homeomorphic to K4 or K2,3.

    The aim of this section is to classify finite abelian groups whose intersection graphs have thickness and outerthickness 1 and 2. We first consider the case of finite cyclic groups.

    Theorem 3.1. Let G=Zn be a cyclic group of order n and Γ(G) is not an empty graph. Then the following statements are hold.

    (1) θ(Γ(G))=1 if and only if G is one following groups:

    Zpα(α=2,3,4,5),Zp1pα2(α=1,2),Zp1p2p3,

    where p1, p2, p3 and p are distinct primes.

    (2) θ(Γ(G))=2 if and only if G is one following groups:

    Zpα(α=6,7,8,9),Zp1pα2(α=3,4),Zp21pα2(α=2,3),Zp1p2p23,Zp1p2p3p4,

    where p1, p2, p3, p4 and p are distinct primes.

    Proof. Suppose that n=pα11pα22pαss with α1α2αs. Then we have V(Γ(Zn))=si=1(αi+1)2. We complete our proofs by discussing the number of prime divisors of n.

    Case 1. s5. It is clear that

    p1,p2,p3,p4,p1p2,p1p3,p1p4,p1p2p3,p1p2p3p4

    are nine non-trivial subgroups of G. Note that any pair of these nine subgroups has non-trivial intersection. So they form a complete graph K9, which is a subgraph of Γ(G). Thus, ω(Γ(G))9. By Lemma 2.1, θ(Γ(G))θ(K9)=3.

    Case 2. s=4. If α42, then G has at least nine non-trivial subgroups, namely,

    p1,p2,p3,p4,p1p2,p1p3,p1p4,p1p2p3,p1p2p3p4.

    Note that any pair of these nine subgroups has non-trivial intersection. So they form a complete graph K9, which is a subgraph of Γ(G). Thus, θ(Γ(G))θ(K9)=3 by Lemma 2.1. If α4=1, then GZp1p2p3p4. So, G has totally 14 non-trivial subgroups, namely,

    p1,p2,p3,p4,p1p2,p1p3,p1p4,p2p3,p2p4,
    p3p4,p1p2p3,p1p2p4,p1p3p4,p2p3p4.

    Note that p1,p2,p3, p1p2, p1p3, p2p3, p1p2p3 form K7 as a subgraph of Γ(G). So ω(Γ(Zp1p2p3p4))7 and V(Γ(Zp1p2p3p4))=14. By Lemma 2.1 and Lemma 2.2, 2=θ(K7)θ(Γ(Zp1p2p3p4))θ(K14)=3. A planar decomposition of Γ(Zp1p2p3p4) shown in Figure 1 deduces that θ(Γ(Zp1p2p3p4))=2.

    Figure 1.  A planar decomposition of Γ(Zp1p2p3p4).

    Case 3. s=3. If α33, then

    p1,p2,p3,p23,p1p2,p1p3,p2p3,p1p2p3,p1p23

    are nine non-trivial subgroups of G. Any pair of these nine subgroups has non-trivial intersection. So they form K9, which is a subgraph of Γ(G). By Lemma 2.1, θ(Γ(G))θ(K9)=3. If α2=α3=2, then

    p1,p2,p3,p23,p1p2,p1p3,p2p3,p1p2p3,p1p23

    are nine non-trivial subgroups of G. Any pair of these nine subgroups has non-trivial intersection. They also form K9 as a subgraph of Γ(G). So θ(Γ(G))θ(K9)=3 by Lemma 2.1.

    If α1=α2=1,α3=2, then G has totally ten non-trivial subgroups, namely,

    p1,p2,p3,p23,p1p2,p1p3,p1p23,p2p3,p2p23,p1p2p3.

    Note that p1, p2, p3, p1p2, p1p3, p2p3, p1p2p3 form K7 as a subgraph of Γ(G). Then ω(Γ(G))7 and V(Γ(G))=10. By Lemma 2.1 and Lemma 2.2, 2=θ(K7)θ(Γ(G))θ(K10)=3. We can give a planar decomposition of Γ(Zp1p2p23) as shown in Figure 2. Thus, θ(Γ(G)=2.

    Figure 2.  A planar decomposition of Γ(Zp1p2p23).

    If α1=α2=α3=1, by Theorem 1.1, Γ(Zp1p2p3) is a planar graph. So θ(Γ(Zp1p2p3))=1.

    Case 4. s=2. In this case, it easy to see that Γ(G)Kα1α21+(Kα1Kα2). Then ω(Γ(Zpα11pα22))=α1α21+α2. If α1α21+α29, then θ(Γ(G))θ(K9)=3. So α1α21+α28 and implies that α24.

    If α2=4, then α1=1. So, Γ(G)K3+(K4K1). By Lemma 2.1 and Lemma 2.2, 2=θ(K7)θ(Γ(Zp1p42))θ(K8)=2. This deduces that θ(Γ(Zp1p42))=2.

    If α2=3, then α12. If α1=2, then Γ(G)K5+(K3K2), By Lemma 2.1 and Lemma 2.2, 2=θ(K8)θ(Γ(Zp21p32))θ(K10)=3. A planar decomposition of Γ(Zp21p32) as shown in Figure 3 deduces that θ(Γ(Zp21p32))=2. If α1=1, then Γ(G)K2+(K3K1). By Lemma 2.1 and Lemma 2.2, 2=θ(K5)θ(Γ(Zp1p32))θ(K6)=2. So, θ(Γ(Zp1p32))=2.

    Figure 3.  A planar decomposition of Γ(Zp21p32).

    If α2=α1=2, then Γ(G)K3+(K2K2). By Lemma 2.1 and Lemma 2.2, 2=θ(K5)θ(Γ(Zp21p22))θ(K7)=2. Thus, θ(Γ(Zp21p22))=2.

    If α2=2 and α1=1 or α1=α2=1, then by Theorem 1.1, Γ(Zp1p2) and Γ(Zp1p22) are planar graphs. So, θ(Γ(Zp1p2))=θ(Γ(Zp1p22))=1.

    Case 5. s=1. In this case, Γ(Zpα)Kα1. By Lemma 2.2, it deduces that if α=2,3,4,5, then θ(Γ(Zpα))=1; and if α=6,7,8,9, then θ(Γ(Zpα))=2.

    This completes the proof.

    For the outerthickness of Γ(G) of a cyclic group G, we have following theorem.

    Theorem 3.2. Let G=Zn be a cyclic group of order n and Γ(G) is not an empty graph. Then the following statements are hold.

    (1) θ0(Γ(G))=1 if and only if G is one following groups: Zpα(α=2,3,4), Zp1pα2(α=1,2), Zp1p2p3, where p1, p2, p3 and p are primes.

    (2) θ0(Γ(G))=2 if and only if G is one following groups: Zpα(α=5,6,7), Zp21p22, Zp1p32, where p1, p2 and p are primes.

    Proof Suppose that n=pα11pα22pαss with α1α2αs. According to the proof of Theorem 3.1, if s4, then θ0(Γ(G))3. We proceed with three cases.

    Case 1. s=3. Again according to the proof of Theorem 3.1, if θ0(Γ(G))2, then GZp1p2p3. In this case, it is clear that Γ(G) contains no subgraph homeomorphic to K4 or K2,3. By Lemma 2.5, Γ(Zp1p2p3) is an outerplanar graph. So, θ0(Γ(Zp1p2p3))=1.

    Case 2. s=2. In this case, Γ(G)Kα1α21+(Kα1Kα2). Then ω(Γ(G))=α1α21+α2. If α1α21+α27, then θ0(Γ(G))3. So α1α21+α26. This deduces that α23.

    If α2=3, then α1=1. So Γ(G)K2+(K3K1). By Lemma 2.1 and Lemma 2.2, 2=θ0(K5)θ0(Γ(Zp1p32))θ0(K6)=2. Hence θ0(Γ(Zp1p32))=2.

    If α2=2,α1=2, then Γ(G)K3+(K2K2). By Lemma 2.1 and Lemma 2.2, 2=θ0(K5)θ0(Γ(Zp21p22))θ0(K7)=3. An outerplanar decomposition of Γ(Zp21p22) as shown in Figure 4 deduces that θ0(Γ(Zp21p22))=2.

    Figure 4.  An outerplanar decomposition of Γ(Zp21p22).

    If α2=2,α1=1, then Γ(G)K1+(K2K1). It is clear that Γ(G) contains no subgraph homeomorphic to K4 or K2,3. By Lemma 2.5, Γ(Zp1p22) is an outerplanar graph. Then, θ0(Γ(Zp1p22))=1.

    If α1=α2=1, then Γ(G)K2. So, Γ(Zp1p2) is an outerplanar graph and thus θ0(Γ(Zp1p2))=1.

    Case 3. s=1. In this case, Γ(Zpα)Kα1. By Lemma 2.2, it deduces that if α=2,3,4, then θ0(Γ(Zpα))=1; and if α=5,6,7, then θ0(Γ(Zpα))=2. The proof is complete.

    Now, we consider the case of finite non-cyclic groups. Let G be a finite Abelian group but not a cyclic group. Then by fundamental theorem of finite abelian group, we have GG1G2Gs, where Gi(i=1,2,,s) is a cyclic group of prime power order.

    Lemma 3.3. Let GG1G2Gs, where s2 and Gi(i=1,2,,s) is a cyclic group of prime power order. If s4, then θ(Γ(G))3.

    Proof. Set

    H1=(1,0,0,0,,0),H2=(1,0,0,0,,0),(0,1,0,0,,0),
    H3=(1,0,0,0,,0),(0,0,1,0,,0),H4=(1,0,0,0,,0),(0,0,0,1,,0),
    H5=(1,0,0,0,,0),(0,1,1,0,,0),H6=(1,0,0,0,,0),(0,1,0,1,,0),
    H7=(1,0,0,0,,0),(0,0,1,1,,0),H8=(1,0,0,0,,0),(0,1,1,1,,0),
    H9=(1,0,0,0,,0),(0,1,0,0,,0),(0,0,1,0,,0),(0,1,1,0,,0).

    Then they are nine non-trivial subgroups of G. Note that any pair of these nine subgroups has non-trivial intersection. So they form K9, which is a subgraph of Γ(G). Thus, θ(Γ(G))θ(K9)=3.

    Theorem 3.4. Let G be a finite non-cyclic abelian group. Then the following statements are hold.

    (1) θ(Γ(G))=1 if and only if G is one following groups: Z2Z2p, ZpZp, where p is a prime.

    (2) θ(Γ(G))=2 if and only if G is one following groups: Z3Z3q, Z5Z5q, Z4Z4, Z2Z8, Z2Z2Z2, where q is a prime.

    Proof. Let the number of distinct prime factors of order of |G| be s. By Lemma 3.3, we know that if θ(Γ(G))2, then s3. By assumption that G is not cyclic, we have three cases for our discussion.

    Case 1. GZpαZpβZpγ, p is a prime.

    Case 2. GZpαZpβZqγ, p,q are distinct primes.

    Case 3. GZpαZpβ, p is a prime.

    For the Case 1, if there is at least one in {α,β,γ} greater than 1, then we may assume γ2. Set

    A1=(0,0,1),(0,1,0),A2=(0,0,1),(1,0,0),A3=(1,1,0),(1,1,1),
    A4=(0,1,0),(1,0,0),A5=(0,1,0),(1,0,1),A6=(1,0,0),(0,1,1),
    A7=(1,1,0),(0,1,1),A8=(0,0,p),(0,1,0),A9=(0,0,p),(1,0,0).

    It is easy to see that they are nine non-trivial subgroups of G. Note that any pair of these nine subgroups has non-trivial intersection. So they form K9, which is a subgraph of Γ(G). Thus, θ(Γ(G))θ(K9)=3.

    Now, we consider the case of α=β=γ=1. If p3, we let

    B1=(0,0,1),(0,1,0),B2=(0,0,1),(1,0,0),B3=(0,0,1),(1,1,0),
    B4=(0,0,1),(1,1,0),B5=(0,1,0),(1,0,0),B6=(0,1,0),(1,0,1),
    B7=(0,1,0),(1,0,1),B8=(0,1,1),(1,0,0),B9=(0,1,1),(1,0,1).

    They also form K9 and θ(Γ(G))θ(K9)=3.

    If p=2, then GZ2Z2Z2. G has totally 14 non-trivial subgroups, namely,

    C1=(0,0,1),(0,1,0),C2=(0,0,1),(1,0,0),C3=(1,1,0),(1,1,1),
    C4=(0,1,0),(1,0,0),C5=(0,1,0),(1,0,1),C6=(1,0,0),(0,1,1),
    C7=(1,1,0),(0,1,1),C8=(0,0,1),C9=(1,0,0),C10=(1,1,0),
    C11=(0,1,0),C12=(1,0,1),C13=(0,1,1),C14=(1,1,1).

    It is easy to see that Γ(G)=K7+7v3. By Lemma 2.4, we get θ(Γ(G))=θ(K7)=2.

    For the Case 2, we may assume αβ. If β2, we let

    D1=(0,p,0),D2=(0,1,0),(0,0,1),D3=(0,1,0),(1,0,0),
    D4=(0,1,0),D5=(0,p,0),(0,0,1),D6=(0,p,0),(1,0,0),
    D7=(1,1,0),D8=(0,p,0),(0,0,1),(1,0,0),D9=(0,p,0),(0,0,1),(1,1,0).

    Note that DiDj=D1 for all ij. So, they form a K9 and θ(Γ(G))θ(K9)=3.

    If γ2, then we let

    E1=(0,0,q),E2=(0,0,1),(1,0,0),E3=(0,0,1),(1,1,0),
    E4=(0,0,q),(0,1,0),E5=(0,0,q),(1,0,0),E6=(0,0,q),(1,1,0),
    E7=(0,0,1),E8=(0,0,1),(0,1,0),E9=(0,0,q),(0,1,0),(1,0,0).

    Note that EiEj=E1 for all ij. So, they form a K9 and θ(Γ(G))θ(K9)=3.

    Now, let GZpZpZq. Then all non-trivial subgroups of G are:

    I=(0,0,1),(0,1,0),
    Mi=(0,0,1),(1,i,0),where0ip1,
    T=(0,1,0),(1,0,0),K=(0,1,0),
    Nj=(1,j,0),where0jp1.

    It is not difficult to see that Γ(G)=(Kp+3e)+(p+1)v2. If p7, then θ(Γ(G))θ(K9)=3. So, p5. By Lemma 2.2, if p=2, then θ(Γ(G))=1; and if p=3,5, then θ(Γ(G))=2.

    For the Case 3, GZpαZpβ. We may assume that βα. If α2 and β3, then we set

    F1=(0,pα),F2=(0,p),F3=(0,1),F4=(0,1),(1,0),F5=(0,p),(1,0),
    F6=(0,1),(p,0),F7=(p,1),F8=(0,p),(p,0),F9=(p,p).

    They are non-trivial subgroups of G. Also, F1=FiFj for all ij. It follows that they form K9 as a subgraph of Γ(G). Thus, θ(Γ(G))θ(K9)=3.

    If α=β=2, then GZp2Zp2. Now, G has totally p+1 subgroups of order p, i.e., (0,p), (p,0), (p,p), (p2p,p); G has totally p2+p+1 subgroups of order p2, i.e., Q=(0,p),(p,0), (1,0), (0,1), (1,1), (p21,1); (1,p), (1,2p), (1,3p), (1,p2p); G has totally p+1 subgroups of order p3, i.e., (0,1),(p,0), (0,p),(1,0), (0,p),(p,0),(1,1), (0,p),(p,0),(1,2), (0,p),(p,0),(1,p1). It is not difficult to see that Γ(Zp2Zp2)Kp+2+(p+1)Kp+1. So, ω(Γ(Zp2Zp2))=2p+3,V(Γ(Zp2Zp2))=p2+3p+3. If p3, then ω(Γ(G))9 and θ(Γ(G))3. So p=2. By Lemma 2.1 and Lemma 2.2, 2=θ(K7)θ(Γ(Z4Z4)θ(K13)=3. On the other hand, a planar decomposition of Γ(Z4Z4) shown in Figure 5 gives θ(Γ(Z4Z4))=2.

    Figure 5.  A planar decomposition of Γ(Z4Z4).

    If α=1, then GZpZpβ. All subgroups of G can be listed below.

    (0,1),(1,1),(p1,1),

    (0,p),(1,p),(p1,p),

    (0,p2),(1,p2),(p1,p2),

    (0,pβ1),(1,pβ1),(p1,pβ1).

    It is clear that Γ(G)Kβ1+(Kpβp+1pK1). So Kβp+βp is a subgraph of Γ(G). If β4, then ω(Γ(G))9 and θ(Γ(G))3. So, β3.

    If β=3, then p=2 and GZ2Z8. So Γ(G)K2+(K52K1) and θ(Γ(Z2Z8)θ(K7)=2. By Lemma 2.4, θ(Γ(Z2Z8))=θ(K7)=2.

    If β=2, then p5. If p=2, then GZ2Z4 and Γ(G)K1+(K32K1). So, θ(Γ(Z2Z4))θ(K4)=1. By Lemma 2.4, θ(Γ(Z2Z4))=θ(K4)=1. If p=3, then GZ3Z9 and Γ(G)K1+(K43K1). Thus, θ(Γ(Z3Z9))θ(K5)=2. By Lemma 2.4, θ(Γ(Z3Z9))=θ(K5)=2. If p=5, then GZ5Z25 and Γ(G)K1+(K65K1). Thus, θ(Γ(Z5Z25))θ(K7)=2. By Lemma 2.4, then θ(Γ(Z5Z25))=θ(K7)=2.

    If β=1, then GZpZp. So, Γ(G)(p+1)K1. By Theorem 1.1, θ(Γ(ZpZp))=1.

    This completes our proof.

    Theorem 3.5. Let G be a finite non-cyclic abelian group. Then the following statements are hold.

    (1) θ0(Γ(G))=1 if and only if GZpZp, where p is a prime.

    (2) θ0(Γ(G))=2 if and only if GZ2Z2p or Z3Z3p, where p is a prime.

    Proof. Suppose GG1G2G3Gs, where s2 and Gi(i=1,2,3,,s) are cyclic groups. According to the proof of Theorem 3.4, if θ0(Γ(G))2, then GZpZp or GZpZpZq or GZpZp2.

    Case 1. GZpZp. In this case, Γ(G)(p+1)K1. It is clear that Γ(G) contains no subgraph homeomorphic to K4 or K2,3. By Lemma 2.5, Γ(ZpZp) is an outerplanar graph, then θ0(Γ(ZpZp))=1.

    Case 2. GZpZpZq. Again by the proof of Theorem 3.4, we know that ω(Γ(ZpZpZq))=p+2. If p5, then θ0(Γ(G))θ0(K7)=3. So p=2, 3. If p=2, then GZ2Z2q. It is easy to see that Γ(Z2Z2q)=(K5e)+3v2. By Lemma 2.1 and Lemma 2.2, 2=θ0(K4)θ0(Γ(Z2Z2q))θ0(K5)=2. Thus, θ0(Γ(Z2Z2q))=2. If p=3, then GZ3Z3q. It is easy to see that Γ(Z3Z3q)=(K6e)+4v2. By Lemma 2.1 and Lemma2.2, 2=θ0(K5)θ0(Γ(Z3Z3q))θ0(K6)=2. So θ0(Γ(Z3Z3q))=2.

    Case 3. GZpZp2. In this case, Γ(ZpZp2)K1+(Kp+1+pK1). Thus, ω(Γ(G))=p+2. If ω(Γ(G))7, then θ0(G)3. Thus, p=2,3. If p=2, then GZ2Z4 and Γ(G)K1+(K32K1). By Lemma 2.4, θ0(Γ(Z2Z4))=θ0(K4)=2. If p=3, then GZ3Z9 and Γ(G)K1+(K43K1). By Lemma 2.4, θ0(Γ(Z3Z9))=θ0(K5)=2.

    The proof is complete.

    Combine Theorems 3.1–3.5, we get the following theorem.

    Theorem 3.6. Let G be a finite abelian group and let p1,p2,p3,p4 and p be primes. Then the following statements are hold.

    (1) θ(Γ(G))=1 if and only if G is one following groups:

    Zpα(α=2,3,4,5),Zp1pα2(α=1,2),Zp1p2p3,Z2Z2p,ZpZp.

    (2) θ(Γ(G))=2 if and only if G is one following groups:

    Zpα(α=6,7,8,9),Zp1pα2(α=3,4),Zp21pα2(α=2,3),Zp1p2p23,Zp1p2p3p4,
    Z3Z3p,Z5Z5p,Z4Z4,Z2Z8,Z2Z2Z2.

    (3) θ0(Γ(G))=1 if and only if G is one following groups:

    Zpα(α=2,3,4),Zp1pα2(α=1,2),Zp1p2p3,ZpZp.

    (4) θ0(Γ(G))=2 if and only if G is one following groups:

    Zpα(α=5,6,7),Zp21p22,Zp1p32,Z2Z2p,Z3Z3p.

    In this section, we investigate the thickness and outerthickness of intersection graphs of subgroups of some finite non-abelian groups. For a positive integer n, τ(n) denotes the number of positive divisor of n; σ(n) denotes the sum of all the positive divisors of n. Sn and An are denoted the symmetric group and alternating group of degree n acting on {1,2,3,,n}, respectively.

    For any integer n3, the dihedral group of order 2n is defined by

    Dn={a,b|an=b2=1,ab=ba1}.

    For any integer n2, the generalized quaternion group of order 4n is defined by

    Qn={a,b|a2n=b4=1,an=b2,b1ab=a1}.

    For any prime p and integer α3, the modular group of order pα is defined by

    Mpα={a,b|apα1=bp=1,bab1=apα2+1}.

    For any integer n4, the quasi-dihedral group of order 2n is defined by

    QD2n={a,b|a2α1=b2=1,bab1=a2α21}.

    Lemma 4.1. [27] The following statements are hold.

    (1) For n3, n=pα11pα22pαkk, V(Γ(Dn))=τ(n)+σ(n)2 and ω(Γ(Dn))=σ(n)n1+ki=1αi.

    (2) Let n2 be an integer and 2n=2α1pα22pαkk, where p2,,pk are distinct odd primes and αi1. Then V(Γ(Qn))=τ(2n)+σ(n)2 and ω(Γ(Qn))=σ(n)1+α1ki=2(αi+1).

    (3) For α4, V(Γ(QD2α))=α+3(2α11) and ω(Γ(QD2α))=α+2α13.

    Theorem 4.2. Let G=Sn. Then the following statements are hold.

    (1) If n=3, then θ(Γ(G))=1.

    (2) If n4, then θ(Γ(G))3.

    Proof. If G=S3, then it is easy to see that Γ(S3)4K1, which is a planar graph. So, θ(Γ(G))=1. If n4, then we set

    S1={(1),(12),(34),(12)(34)},

    S2={(1),(12),(13),(23),(123),(132)},

    S3={(1),(12),(14),(24),(124),(142)},

    S4={(1),(13),(14),(34),(134),(143)},

    S5={(1),(23),(24),(34),(234),(243)},

    S6={(1),(13),(24),(13)(24),(14)(23),(12)(34),(1234),(1432)},

    S7={(1),(14),(23),(13)(24),(14)(23),(12)(34),(1243),(1342)},

    S8={(1),(13),(24),(13)(24),(14)(23),(12)(34),(1324),(1423)},

    S9={(1),(123),(132),(124),(142),(134),(143),(234),(243),(13)(24),(14)(23),(12)(34)}.

    Then they are nine non-trivial subgroups of Sn and form a complete graph K9 as a subgraph of Γ(Sn). Thus, θ(Γ(G))θ(K9)3.

    Theorem 4.3. Let G=An. Then the following statements are hold.

    (1) θ(Γ(G))=1 if and only if n=4.

    (2) If n5, then θ(Γ(G))3.

    Proof. If GA4, then Γ(A4)4K1K1,3, which is a planar graph. So, θ(Γ(G))=1. If n5, we let

    L1={(1),(12345),(13524),(14253),(15432),(12)(35),(13)(45),(14)(23),(15)(24),(25)(34)},

    L2={(1),(12453),(13542),(14325),(15234),(12)(34),(13)(25),(14)(35),(15)(24),(23)(45)},

    L3={(1),(12354),(13425),(14532),(15243),(12)(34),(13)(45),(14)(25),(15)(23),(24)(35)},

    L4={(1),(12435),(13254),(14523),(15342),(12)(45),(13)(24),(14)(35),(15)(23),(25)(34)},

    L5={(1),(12534),(13245),(14352),(15423),(12)(45),(13)(25),(14)(23),(15)(34),(24)(35)},

    L6={(1),(12543),(13452),(14235),(15324),(12)(35),(13)(24),(14)(25),(15)(34),(23)(45)},

    L7={(1),(123),(124),(134),(234),(132),(142),(143),(243),(13)(24),(14)(23),(12)(34)},

    L8={(1),(123),(125),(135),(235),(132),(152),(153),(253),(12)(35),(13)(25),(15)(23)},

    L9={(1),(124),(125),(145),(245),(142),(152),(154),(254),(12)(45),(14)(25),(15)(24)}.

    Then they are nine non-trivial subgroups of An and form a complete graph K9 as a subgraph of Γ(An). Thus, θ(Γ(G))θ(K9)3.

    Theorem 4.4. Let G=Dn(n3) be the dihedral group of order 2n. Then the following statements are hold.

    (1) θ(Γ(G))=1 if and only if GD4 or Dp, where p is an odd prime.

    (2) θ(Γ(G))=2 if and only if GD6,D9, D25 or D10.

    Proof. For n3, let n=pα11pα22pαkk with α1α2αk. By Lemma 4.1, we know that ω(Γ(Dn))=σ(n)n1+ki=1αi. If k2 and αk2, then ω(Γ(Dn))p1+pk+p1pk11. Thus, if θ(Γ(G))2, then n=p1p2 or n=pα.

    If n=p1p2, then by Lemma 4.1, ω(Γ(Dn))=σ(n)n1+ki=1αi=p1+p2+1. If p1+p2+19, then θ(Γ(Dn))3. Thus, if θ(Γ(Dn))2, then p1+p2+18. Thus, p1=2, p2=3 or p1=2, p2=5.

    If p1=2, p2=3, then ω(Γ(D6))=6. So, θ(Γ(D6))θ(K6)=2. Note that degΓ(D6)(aib)=τ(6)2=2 and V(aib)=6. By Lemma 2.4, θ(Γ(D6))=θ(Γ(D6)6v2). On the other hand, V(Γ(D6)6v2)=τ(6)1+σ(6)61=8, then 2=θ(K6)θ(Γ(D6)6v2)θ(K8)=2. Hence, θ(Γ(D6))=2.

    If p1=2 and p2=5, then ω(Γ(D10))=8. So, θ(Γ(D10))θ(K8)=2. Note that degΓ(D10)(aib)=τ(10)2=2 and V(aib)=n=10. By Lemma 2.4, θ(Γ(D10))=θ(Γ(D10)10v2). As V(Γ(D10)10v2)=τ(10)1+σ(10)101=10, we have 2=θ(K8)θ(Γ(D10)10v2)θ(K10)=3. On the other hand, a planar decomposition of Γ(D10)10v2 as shown in Figure 6 deduces that θ(Γ(D10))=2.

    Figure 6.  A planar decomposition of Γ(D10)10v2.

    If n=pα, then by Lemma 4.1, ω(Γ(Dn))=pαpp1+α. If ω(Γ(Dn))9, then θ(Γ(Dn))3. So, α2.

    If α=1, then Γ(Dp)(p+1)K1 and thus θ(Γ(Dp))=1. If α=2, then ω(Γ(Dp2))=p+2. Note that degΓ(Dp2)(aib)=τ(p2)2=1, V(aib)=n=p2, and V(Γ(Dp2))=τ(n)+σ(n)2=1+p+p2+32=p2+p+2. So, Γ(Dp2)=Kp+2+p2v1. By Lemma 2.4, θ(Γ(Dp2))=θ(Kp+2). Since p+29 deduces θ(Γ(Dp2))3, we have p+28, that is, p=2,3,5. By Lemma 2.2, we get if p=2, then θ(Γ(D4))=1; and if p=3,5, then θ(Γ(D9))=θ(Γ(D25))=2.

    Theorem 4.5. Let G=Qn. Then the following statements are hold.

    (1) θ(Γ(G))=1 if and only if G=Q2.

    (2) θ(Γ(G))=2 if and only if GQ3 or Q5.

    Proof. Let n2 be an integer and 2n=2α1pα22pαkk, where pi are distinct primes and αi1. By Lemma 4.1, ω(Γ(Qn))=σ(n)1+α1ki=2(αi+1). If k2, then ω(Γ(Qn))1+p1+p2+p1p21+1×2×2=4+p1+p2+p1p2>4+2+3+2×3>9 and θ(Γ(Qn))>3. So, if θ(Γ(Qn))2, then we have n=pα.

    If n=2α, ω(Γ(Q2α))=σ(n)1+α1ki=2(αi+1)=2α+11+α,V(Γ(Q2α))=τ(2n)+σ(n)2=2α+11+α, then Γ(Q2α) is a complete graph. If 2α+11+α9, then θ(Γ(Qn))3. So, if θ(Γ(Qn))2, then 2α+11+α8 and thus α=1. In this case, θ(Γ(Q2))=1.

    If n=pα and p2, then ω(Γ(Qn))=σ(n)1+α1ki=2(αi+1)=1+p+p2++pα1+pα+1×(α+1)1=pα+11p1+α9, then θ(Γ(Qn))3. So, if θ(Γ(Qn))2, then pα+11p1+α8, then α=1.

    For α=1, ω(Γ(Qn))=σ(n)1+α1ki=2(αi+1)=p+2 and V(Γ(Qn))=τ(2n)+σ(n)2=p+3. So, θ(Kp+2)θ(Γ(Qn))θ(Kp+3). If ω(Γ(Qn))=p+29, then θ(Γ(Qn))3. So, if θ(Γ(Qn))2, then p+28, then p=3,5.

    By Lemma 2.1 and Lemma 2.2, 2=θ(K5)θ(Γ(Q3))θ(K6)=2 and 2=θ(K7)θ(Γ(Q5))θ(K8)=2, then θ(Γ(Q3))=θ(Γ(Q5))=2.

    Theorem 4.6. Let G=QD2α. Then θ(Γ(G))3.

    Proof. For α4, by Lemma 4.1, we have ω(Γ(QD2α))=α+2α139. Thus, θ(Γ(QD2α))θ(K9)=3.

    Theorem 4.7. Let GMpα. Then the following statements are hold.

    (1) θ(Γ(G))=1 if and only if G=M8.

    (2) θ(Γ(G))=2 if and only if GM27,M125 or M16.

    Proof. If pα=8, then Γ(M8)Γ(D4). So θ(Γ(M8))=1. If pα8, then Γ(Mpα)Γ(ZpZpα1). If θ(Γ(G))2, then by the proof of Theorem 3.5, we know α=3,4. If α=3, then Γ(Mp3)K1+(Kp+1pK1). By θ(Γ(G))2, we get p=3,5. By Lemma 2.4, θ(Γ(M27))=θ(K5)=2 and θ(Γ(M125))=θ(K7)=2. If α=4, then p=2. By Theorem 3.5, θ(Γ(M16))=θ(Γ(Z2Z8))=2.

    Theorem 4.8. Let G=Sn. Then the following statements are hold.

    (1) If n=3, then θ0(Γ(G))=1.

    (2) If n4, then θ0(Γ(G))3.

    Proof. If G=S3, then Γ(S3)4K1, which is an outerplanar graph. So θ0(Γ(S3))=1. By the proof of Theorem 4.2, if n4, then θ0(Γ(G))θ0(Γ(K9))=3.

    Theorem 4.9. Let GAn. Then the following statements are hold.

    (1) If n=4, then θ0(Γ(G))=1.

    (2) If n5, then θ0(Γ(G))3.

    Proof. Note that Γ(A4)4K1K1,3. It is clear that Γ(A4) contains no subgraph homeomorphic to K4 or K2,3. By Lemma 2.5, Γ(A4) is an outerplanar graph. So θ0(Γ(A4))=1. By the proof of Theorem 4.3, if n5, then θ0(Γ(G))θ0(K9)=3.

    Theorem 4.10. Let GDn. Then the following statements are hold.

    (1) θ0(Γ(G))=1 if and only if GDp, where p is a prime.

    (2) θ0(Γ(G))=2 if and only if G is one following groups: D4, D6, D9.

    Proof. For n3, let n=pα11pα22pαkk with α1α2αk. According to the proof of Theorem 4.4, if θ0(Γ(Dn))2, then n=p1p2 or pα.

    If n=p1p2, then ω(Γ(Dn))=σ(n)n1+ki=1αi=p1+p2+1. If p1+p2+17, then θ0(Γ(Dn))3. So, if θ0(Γ(Dn))2, then p1+p2+16, that is, p1=2, p2=3. As ω(Γ(D6))=2+3+1=6, we know θ0(Γ(D6))θ0(K6)=2. An outerplanar decomposition of Γ(D6)6v2 as shown in Figure 7 gives θ0(Γ(D6))=2.

    Figure 7.  An outerplanar decomposition of Γ(D66v2).

    If n=pα, then ω(Γ(Dn))=σ(n)n1+ki=1αi=p+p2++pα1+α=pαpp1+α7, then θ0(Γ(Dn))3. So, if θ0(Γ(Dn))2, then pαpp1+α6, so α=1,2.

    If α=1, then Γ(Dp)(p+1)K1. So, θ0(Γ(Dp))=1.

    If α=2, Γ(Dp2)=Kp+2+p2v1. By Lemma 2.4, θ0(Γ(Dp2))=θ0(Kp+2). As p+27, then θ0(Γ(Dp2))3. So, if θ0(Γ(Dp2))2, then p+26, then p=2,3. By Lemma 2.2, we get if p=2, then θ0(Γ(D4))=θ0(K4)=2; and if p=3, then θ0(Γ(D9))=θ0(K5)=2.

    Theorem 4.11. Let GQn. Then the following statements are hold.

    (1) θ0(Γ(G))2.

    (2) θ0(Γ(G))=2 if and only if GQ2 or Q3.

    Proof. According the proof of Theorem 4.5, if θ0(Γ(Qn))2, then n=pα.

    If n=2α, again by the proof of Theorem 4.5, we have Γ(Q2α)=K2α+11+α. If 2α+11+α7, then θ0(Γ(Qn))3. So θ0(Γ(Qn))2, then 2α+11+α6, so α=1. In this case, θ0(Γ(Q2))=2.

    If n=pα and p2, then ω(Γ(Qn))=pα+11p1+α7 deduces θ0(Γ(Qn))3. So, if θ0(Γ(Qn))2, then pα+11p1+α6 and α=1.

    In this case, ω(Γ(Qn))=σ(n)1+α1ki=2(αi+1)=p+2 and V(Γ(Qn))=τ(2n)+σ(n)2=p+3. So, θ0(Kp+2)θ0(Γ(Qn))θ0(Kp+3). If ω(Γ(Qn))=p+27, then θ0(Γ(Qn))3. So, if θ0(Γ(Qn))2, then p+26, so p=3. By Lemma 2.1 and Lemma 2.2, 2=θ0(K5)θ0(Γ(Q3))θ0(K6)=2, then θ0(Γ(Q3))=2.

    Theorem 4.12. Let GQD2α. Then θ0(Γ(G)3.

    Proof. For α4, by Lemma 4.1, we have ω(Γ(QD2α))=α+2α139. Thus, θ0(Γ(QD2α))θ0(K9)=3.

    Theorem 4.13. Let G=Mpα. Then the following statements are hold.

    (1) θ0(Γ(G))2.

    (2) θ0(Γ(G))=2 if and only if GM8 or M27.

    Proof. If pα=8, then Γ(M8)Γ(D4). By Theorem 4.4, θ0(Γ(M8))=2. If pα8, then Γ(Mpα)Γ(ZpZpα1). By the proof of Theorem 4.7, we get α=3,p=3. In this case, we have θ0(Γ(M27))=2.

    In this paper, we classify all finite abelian groups whose thickness and outerthickness of subgroup intersection graphs are 1 and 2, respectively. We also investigate the thickness and outerthickness of subgroup intersection graphs for some finite non-abelian groups. The method is finding a forbidden subgraph of Γ(G), like K9 and then determining the remaining cases. The key point is how to separate Γ(G) into two planar graphs. Next, one can consider to determine all finite abelian groups G whose Γ(G) has thickness or outerthickness 3. For a general finite group G, how to determine the thickness and outerthickness of Γ(G) is a challenge. Of course, investigating other properties of Γ(G) is a valuable exploration.

    We would like to thank the referees for a careful reading of the paper and for their valuable comments. This work was supported by the National Natural Science Foundation of China (Grant No. 11661013, 11961050), the Guangxi Natural Science Foundation (Grant No. 2020GXNSFAA159103, 2020GXNSFAA159053) and High-level talents for scientific research of Beibu Gulf University (2020KYQD07). The second author is supported by the Science and Technology Research Project of Jiangxi Education Department.

    The authors declared that they have no conflicts of interest to this work.

    [1] Field CB, Raupach MR (2004) The Global Carbon Cycle: Integrating Humans, Climate and the Natural World. Island Press, Washington DC.
    [2] Falkowski P, Scholes RJ, Boyle E, et al. (2000) The global carbon cycle: A test of our knowledge of the earth as a system. Sci 290: 291–296. doi: 10.1126/science.290.5490.291
    [3] Canadell JG, Mooney HA, Baldocchi DD, et al. (2000) Carbon metabolism of the terrestrial biosphere: a multitechnique approach for improved understanding. Ecosyst 3: 115–130. doi: 10.1007/s100210000014
    [4] Raupach MR, Rayner PJ, Barrett DJ et al. (2005) Model-data synthesis in terrestrial carbon observation: methods, data requirements and data uncertainty specifications. GCB 11: 378–397.
    [5] Enting IG, Rayner PJ, Ciais P (2012) Carbon cycle uncertainty in REgional Carbon Cycle Assessment and Processes (RECCAP). Biogeosci 9: 2889–2904. doi: 10.5194/bg-9-2889-2012
    [6] Karplus W (1977) The spectrum of mathematical modeling and systems simulation. Math Comp Simul 19: 3–10 doi: 10.1016/0378-4754(77)90034-9
    [7] Enting IG (1987) A modeling spectrum for carbon cycle studies. Math Comp Simul 29: 75–85. doi: 10.1016/0378-4754(87)90099-1
    [8] Enting IG (2010) Inverse problems and complexity in earth system science. In R. L. Dewar and F. Detering, editors, Complex Physical, Biophysical and Econophysical Systems. World Scientific, Singapore.
    [9] Enting IG (2002) Inverse Problems in Atmospheric Constituent Transport. CUP, Cambridge, UK.
    [10] Boschetti F, Grigg N, Enting IG (2010) Modelling = conditional prediction. Ecological Complexity.
    [11] Enting IG (2008) Assessing the information content in environmental modelling: A carbon cycle perspective. Entropy 10: 556–575. doi: 10.3390/e10040556
    [12] Rayner PJ, Brien DO (2001) The utility of remotely sensed CO2 concentration data in surface source inversions. Geophys Res Lett 28: 175–178. doi: 10.1029/2000GL011912
    [13] Evans SN, Stark PB (2002) Inverse problems as statistics. Inverse Probl 18: R55–R97. doi: 10.1088/0266-5611/18/4/201
    [14] Trudinger CM, Raupach MR, Rayner PJ, et al. (2008) OptIC project: an intercomparison of optimization techniques for parameter estimation in terrestrial biogeochemical models. J Geophys Res 112: G02027 .
    [15] Schneider SH (1983) The problem of pre-industrial CO2 concentration: An editorial. Clim Change 5: 311–313. doi: 10.1007/BF02423527
    [16] Laurmann JA, Spreiter JR (1983) The e ects of carbon cycle model error in calculating future atmospheric carbon dioxide levels. Clim Change 5: 145–181. doi: 10.1007/BF02423488
    [17] Gloor M, Sarmiento JL, Gruber N (2010) What can be learned about carbon cycle climate feedbacks from the CO2 airborne fraction? Atmos Chem Phys 10: 7739–7751. doi: 10.5194/acp-10-7739-2010
    [18] Oeschger H, Heimann M (1983) Uncertainties of predictions of future atmospheric CO2 concentrations. J Geophys Res 88C: 1258–1262.
    [19] Joos F, Roth R, Fuglestvedt JS, et al. (2013) Carbon dioxide and climate impulse response functions for the computation of greenhouse gas metrics: a multi-model analysis. Atmos Chem Phys 13: 2793–2825. doi: 10.5194/acp-13-2793-2013
    [20] Wigley TML (1991) A simple inverse carbon cycle model. Global Biogeochem Cycles 5: 373–382. doi: 10.1029/91GB02279
    [21] Joos F, Bruno M (1998) Long-term variability of the terrestrial and oceanic carbon sinks and the budgets of the carbon isotopes 13C and 14C. Global Biogeochem Cycles 12: 277–295. doi: 10.1029/98GB00746
    [22] Joos F, Prentice IC, Sitch S, et al. (2001) Global warming feedbacks on terrestrial carbon uptake under the Intergovernmental Panel on Climate Change (IPCC) scenarios. Global Biogeochem Cycles 15: 891–907. doi: 10.1029/2000GB001375
    [23] Meinshausen M, Raper SCB, Wigley TML (2011) Emulating coupled atmosphere-ocean and carbon cycle models with a simpler model: - part 1 model description and calibration MAGICC6. Atmos Chem Phys 11: 1417–1456. doi: 10.5194/acp-11-1417-2011
    [24] Enting IG (2007) Laplace transform analysis of the carbon cycle. Environ Modelling Software 22: 1488–1497. doi: 10.1016/j.envsoft.2006.06.018
    [25] Trudinger CM, Enting IG, Rayner PJ, et al. (2002) Kalman filter analysis of ice core data. 2 Double deconvolution of CO2 and 13C measurements. J Geophys Res 107: 4423.
    [26] Enting IG, Pearman GI (1987) Description of a one-dimensional carbon cycle model calibrated using techniques of constrained inversion. Tellus 39B: 459–476. doi: 10.1111/j.1600-0889.1987.tb00206.x
    [27] Young P, Parkinson S, Lees M (1996) Simplicity out of complexity in environmental modelling: Occam's razor revisited. J Appl Statist 23: 165–210. doi: 10.1080/02664769624206
    [28] Bodman RW (2011) Estimating Uncertainties in Future Global Warming using a Simple Climate Model. PhD thesis, University of Melbourne.
    [29] Bodman RW, Rayner PJ, Karoly DJ (2013) Uncertainty in temperature projections reduced using carbon cycle and climate observations. Nature Clim Change pages 725–729.
    [30] Sundquist ET (1985) Geological perspectives on carbon dioxide and the carbon cycle. In E. T. Sundquist and W. S. Broecker, editors, The Carbon Cycle and Atmospheric CO2: Natural Variations Archean to Present, Geophysical Monograph 32, pages 5–59. AGU, Washington.
    [31] Wigley TML, Raper SCB (1992) Implications for climate and sea level rise of the revised IPCC emission scenarios. Nature 357: 293–300. doi: 10.1038/357293a0
    [32] Willeit M, Ganopolski A, Dalmonech D, et al. (2014) Time-scale and state dependence of the carbon-cycle feedback to climate. Clim Dynamics 42: 1699–1713. doi: 10.1007/s00382-014-2102-z
    [33] Friedlingstein P, Dufresne JL, Cox PM, et al. (2003) How positive is the feedback between climate change and the carbon cycle? Tellus 55B: 692–700.
    [34] Rubino M, Etheridge DM, Trudinger CM, et al. (2016) Low atmospheric CO2 levels during the Little Ice Age due to cooling-induced terrestrial uptake. Nature Geosci 9: 691–694. doi: 10.1038/ngeo2769
    [35] Enting IG, Mansbridge JV (1987) The incompatibility of ice-core CO2 data with reconstructions of biotic CO2 sources. Tellus B 39B: 318–325. doi: 10.1111/j.1600-0889.1987.tb00102.x
    [36] Enting IG (1992) The incompatibility of ice-core CO2 data with reconstructions of biotic CO2sources (II). Tellus B 44: 23–32. doi: 10.3402/tellusb.v44i1.15418
    [37] Broecker WS, Peng TH, Engh R (1980) Modeling the carbon system. Radiocarbon 22: 565–598. doi: 10.1017/S0033822200009966
    [38] Enting IG, Mansbridge JV (1987) Inversion relations for the deconvolution of CO2 data from ice cores. Inverse Problems 3: L63–69. doi: 10.1088/0266-5611/3/4/001
    [39] Siegenthaler U, Oeschger H (1987) Biospheric CO2 emissions during the past 200 years reconstructed by deconvolution of ice core data. Tellus 39B: 140–154. doi: 10.1111/j.1600-0889.1987.tb00278.x
    [40] Enting IG (2010) Inversion of atmospheric CO2 concentrations. In D. Moreira and M. Vilhena, editors, Air Pollution and Turbulence: Modeling and Applications, chapter 11, pages 287–316. CRC Press (Taylor and Francis), Boca Raton, Florida.
    [41] Bolin B, Keeling CD (1963) Large-scale atmospheric mixing as deduced from the seasonal and meridional variations of carbon dioxide. J Geophys Res 68: 3899–3920. doi: 10.1029/JZ068i013p03899
    [42] Enting IG, Mansbridge JV (1989) Seasonal sources and sinks of atmospheric CO2: Direct inversion of filtered data. Tellus 41B: 111–126. doi: 10.1111/j.1600-0889.1989.tb00129.x
    [43] Tans PP, Conway TJ, Nakazawa T (1989) Latitudinal distribution of the sources and sinks of atmospheric carbon dioxide derived from surface observations and an atmospheric transport model. J Geophys Res 94D: 5151–5172.
    [44] Enting IG, Mansbridge JV (1991) Latitudinal distribution of sources and sinks of CO2: Results of an inversion study. Tellus 43B: 156–170.
    [45] Law RM (1999) CO2 sources from a mass-balance inversion: sensitivity to the surface constraint. Tellus 51B: 254–265.
    [46] Dargaville RJ, Simmonds I (2000). Calculating CO2 fluxes by data assimilation coupled to a three dimensional massbalance inversion. In P. Kasibhatla et al., editors, Inverse Methods on Global Biogeochemical Cycles, pages 255–264. AGU, Washington, DC.
    [47] Keeling CD, Piper SC, Heimann M (1989) A three-dimensional model of atmospheric CO2 transport based on observed winds. 4: Mean annual gradients and interannual variations. In D. H. Peterson, editor, Aspects of Climate Variability of the Pacific and Western Americas . Geophysical Monograph 55. AGU, Washington.
    [48] Tans PP, Fung IY, Takahashi T. Observational constraints on the global atmospheric CO2 budget. Sci 247: 1431–1438.
    [49] Fung IY, John J, Lerner J, et al. (1991) Three-dimensional model synthesis of the global methane cycle. J Geophys Res 96D: 13033–13065.
    [50] Enting IG, Trudinger CM, Francey RJ (1995). A synthesis inversion of the concentration and 13C of atmospheric CO2. Tellus 47B: 35–52.
    [51] Trudinger CM, Raupach MR, Rayner PJ, et al. (2008) Using the Kalman filter for parameter estimation in biogeochemical models. Environ.
    [52] Newsam GN, Enting IG (1988) Inverse problems in atmospheric constituent studies: I. Determination of surface sources under a di usive transport approximation. Inverse Problems 4: 1037–1054.
    [53] Griewank A (2000). Evaluating Derivatives: Principles and Techniques of Algorithmic Di erentiation. SIAM, Philadelphia.
    [54] Enting IG (2011) Tangents, adjoints and computational complexity in terrestrial carbon modelling. In W. McLean and A. J. Roberts, editors, Proceedings of the 15th Biennial Computational Techniques and Applications Conference, CTAC-2010, volume 52 of ANZIAM J., pages C806–C822, October.
    [55] R¨odenbeck C, Houweling S, Gloor M, et al. (2003) CO2 flux history 1982–2001 inferred from atmospheric data using a global inversion of atmospheric transport. Atmos Chem Phys 3: 1919– 1964. doi: 10.5194/acp-3-1919-2003
    [56] Kaminski T, Knorr W, Rayner PJ, et al. (2002) Assimilating atmospheric data into a terrestrial biosphere model: A case study of the seasonal cycle. Global Biogeochem Cycles 16: 1066.
    [57] Rayner PJ, Scholze M, Knorr W, et al. (2005) Two decades of terrestrial carbon fluxes from a carbon cycle data assimilation system (CCDAS). Global Biogeochemical Cycles, 19:GB2026.
    [58] Enting IG (2011) Seeking carbon-consistency in the climate-science-to-policy interface. Biogeochem 104: 59–67. doi: 10.1007/s10533-009-9351-7
    [59] Tans PP, Berry JA, Keeling RF (1993) Oceanic 13C/12C observations: A new window on ocean CO2 uptake. Global Biogeochem Cycles 7: 353–368. doi: 10.1029/93GB00053
    [60] Quay PD, Tilbrook B, Wong CS (1992) Oceanic uptake of fossil fuel CO2: Carbon-13 evidence. Sci 256: 74–79. doi: 10.1126/science.256.5053.74
    [61] Sarmiento JL, Sundquist ET (1992) Revised budget for the oceanic uptake of anthropogenic carbon dioxide. Nature 356: 589–593. doi: 10.1038/356589a0
    [62] Heimann M, Maier-Reimer E (1996) On the relations between the oceanic uptake of CO2 and its isotopes. Global Biogeochem Cycles 10: 89–110. doi: 10.1029/95GB03191
    [63] Francey RJ, Tans PP, Allison CE, et al. (1995) Changes in oceanic and terrestrial carbon uptake since 1982. Nature 373: 326–330. doi: 10.1038/373326a0
    [64] Keeling CD, Whorf TP, Wahlen M, et al. (1995) Interannual extremes in the rate of rise of atmospheric carbon dioxide since 1980. Nature 375: 666–670. doi: 10.1038/375666a0
    [65] Keeling RF, Najjar RP, Bender ML, et al. (1993) What atmospheric oxygen measurements can tell us about the global carbon cycle. Global Biogeochem Cycles 7:37–67. doi: 10.1029/92GB02733
    [66] Bender M, Ellis T, Tans P, et al. (1996) Variability in the O2/N2 ratio of southern hemisphere air, 1991–1994: Implications for the carbon cycle. Global Biogeochem Cycles 10: 9–21. doi: 10.1029/95GB03295
    [67] Manning AC, Keeling RF, Severinghaus JP (1999) Precise atmospheric oxygen measurements with a paramagnetic oxygen analyzer. Global Biogeochem Cycles 13: 1107–1115. doi: 10.1029/1999GB900054
    [68] Denning AS, Fung IY, Randall D (1995) Latitudinal gradient of atmospheric CO2 due to seasonal exchange with the land biota. Nature 376: 240–243. doi: 10.1038/376240a0
    [69] Friedlingstein P, Cox P, Betts R, et al. (2006) Climate-carbon cycle feedback analysis: Results from the C4MIP model intercomparison. J Clim 19: 3337–3353. doi: 10.1175/JCLI3800.1
    [70] Canadell JG, Ciais P, Gurney K, et al. (2011) An international e ort to quantify regional carbon fluxes. EOS Trans AGU 92: 81–82.
    [71] Le Qu´er´e C, Raupach MR, Canadell JG, et al. (2009) Trends in the sources and sinks of carbon dioxide. Nature Geosci 2: 831–836. doi: 10.1038/ngeo689
    [72] Le Qu´er´e C, Andrew RM, Canadell JG, et al. (2016) Global carbon budget 2016. Earth Syst Sci Data 8: 605–649. doi: 10.5194/essd-8-605-2016
  • This article has been cited by:

    1. Arezoo Beheshtipour, Seyyed Majid Jafarian Amiri, The Clique Number of the Intersection Graph of a Finite Group, 2023, 49, 1017-060X, 10.1007/s41980-023-00804-5
    2. Ling Zhu, Huadong Su, Finite Abelian groups with positive genus subgroup intersection graphs, 2024, 23, 0219-4988, 10.1142/S0219498824502505
  • Reader Comments
  • © 2018 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(4434) PDF downloads(915) Cited by(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog