
Typical tree health assessment methods are destructive. Non-destructive Ground Penetrating Radar (GPR) technique can provide a diagnostic tool for assessing the health status of live tree trunks based on internal dielectric permittivity distribution. Typical GPR acquisition technique—common-offset—is not helpful in providing robust and high-resolution quantitative results. In the current work we evaluate the capabilities of GPR tomography on locating tree-decays in a number of different tree species, imaging the interval structure of a healthy tree and quantitative estimation of moisture content (MC) based on distribution of dielectric permittivity, directly related to MC. The measurements described in this work were made on the trunks of live trees of different species in different conditions: a "healthy" English oak (Quercus robur), a "dry" Siberian fir (Picea obovata), a Horse chestnut (Castanea dentata) and a European aspen (Populus tremula) with rot inside. The results of the suggested approach were confirmed by resistography. Different parts of the trunk (bark, core, sapwood), as well as healthy and affected areas differ in moisture content, so the method of GPR tomography allowed us to see both the structure of the trunk of a healthy tree, and the presence and dimensions of defects.
Citation: Maria Sudakova, Evgeniya Terentyeva, Alexey Kalashnikov. Assessment of health status of tree trunks using ground penetrating radar tomography[J]. AIMS Geosciences, 2021, 7(2): 162-179. doi: 10.3934/geosci.2021010
[1] | Yue Feng, Yujie Liu, Ruishu Wang, Shangyou Zhang . A conforming discontinuous Galerkin finite element method on rectangular partitions. Electronic Research Archive, 2021, 29(3): 2375-2389. doi: 10.3934/era.2020120 |
[2] | Chunmei Wang . Simplified weak Galerkin finite element methods for biharmonic equations on non-convex polytopal meshes. Electronic Research Archive, 2025, 33(3): 1523-1540. doi: 10.3934/era.2025072 |
[3] | Leilei Wei, Xiaojing Wei, Bo Tang . Numerical analysis of variable-order fractional KdV-Burgers-Kuramoto equation. Electronic Research Archive, 2022, 30(4): 1263-1281. doi: 10.3934/era.2022066 |
[4] | Guanrong Li, Yanping Chen, Yunqing Huang . A hybridized weak Galerkin finite element scheme for general second-order elliptic problems. Electronic Research Archive, 2020, 28(2): 821-836. doi: 10.3934/era.2020042 |
[5] | Victor Ginting . An adjoint-based a posteriori analysis of numerical approximation of Richards equation. Electronic Research Archive, 2021, 29(5): 3405-3427. doi: 10.3934/era.2021045 |
[6] | Jun Pan, Yuelong Tang . Two-grid H1-Galerkin mixed finite elements combined with L1 scheme for nonlinear time fractional parabolic equations. Electronic Research Archive, 2023, 31(12): 7207-7223. doi: 10.3934/era.2023365 |
[7] | Hongze Zhu, Chenguang Zhou, Nana Sun . A weak Galerkin method for nonlinear stochastic parabolic partial differential equations with additive noise. Electronic Research Archive, 2022, 30(6): 2321-2334. doi: 10.3934/era.2022118 |
[8] | Zexuan Liu, Zhiyuan Sun, Jerry Zhijian Yang . A numerical study of superconvergence of the discontinuous Galerkin method by patch reconstruction. Electronic Research Archive, 2020, 28(4): 1487-1501. doi: 10.3934/era.2020078 |
[9] | Suayip Toprakseven, Seza Dinibutun . A weak Galerkin finite element method for parabolic singularly perturbed convection-diffusion equations on layer-adapted meshes. Electronic Research Archive, 2024, 32(8): 5033-5066. doi: 10.3934/era.2024232 |
[10] | Bin Wang, Lin Mu . Viscosity robust weak Galerkin finite element methods for Stokes problems. Electronic Research Archive, 2021, 29(1): 1881-1895. doi: 10.3934/era.2020096 |
Typical tree health assessment methods are destructive. Non-destructive Ground Penetrating Radar (GPR) technique can provide a diagnostic tool for assessing the health status of live tree trunks based on internal dielectric permittivity distribution. Typical GPR acquisition technique—common-offset—is not helpful in providing robust and high-resolution quantitative results. In the current work we evaluate the capabilities of GPR tomography on locating tree-decays in a number of different tree species, imaging the interval structure of a healthy tree and quantitative estimation of moisture content (MC) based on distribution of dielectric permittivity, directly related to MC. The measurements described in this work were made on the trunks of live trees of different species in different conditions: a "healthy" English oak (Quercus robur), a "dry" Siberian fir (Picea obovata), a Horse chestnut (Castanea dentata) and a European aspen (Populus tremula) with rot inside. The results of the suggested approach were confirmed by resistography. Different parts of the trunk (bark, core, sapwood), as well as healthy and affected areas differ in moisture content, so the method of GPR tomography allowed us to see both the structure of the trunk of a healthy tree, and the presence and dimensions of defects.
For simplicity, we consider Poisson equation with a Dirichlet boundary condition as our model problem.
−Δu=f,inΩ, | (1) |
u=g,on∂Ω, | (2) |
where
Using integration by parts, we can get the variational form: find
(∇u,∇v)=(f,v),∀v∈H10(Ω). | (3) |
Various finite element methods have been introduced to solve the Poisson equations (1)-(2), such as the Galerkin finite element methods (FEMs)[2, 3], the mixed FEMs [15] and the finite volume methods (FVMs) [6], etc. The FVMs emphasis on the local conservation property and discretize equations by asking the solution satisfying the flux conservation on a dual mesh consisting of control volumes. The mixed FEMs is another category method that based on the variable
The classical conforming finite element method obtains numerical approximate results by constructing a finite-dimensional subspace of
(∇uh,∇vh)=(f,vh),∀vh∈V0h, | (4) |
where
One obvious disadvantage of discontinuous finite element methods is their rather complex formulations which are often necessary to ensure connections of discontinuous solutions across element boundaries. For example, the IPDG methods add parameter depending interior penalty terms. Besides additional programming complexity, one often has difficulties in finding optimal values for the penalty parameters and corresponding efficient solvers. Most recently, Zhang and Ye [21] developed a discontinuous finite element method that has an ultra simple weak formulation on triangular/tetrahedal meshes. The corresponding numerical scheme can be written as: find
(∇wuh,∇wvh)=(f,vh),∀vh∈V0h, | (5) |
where
Following the work in [21, 22], we propose a new conforming DG finite element method on rectangular partitions in this work. It can be obtained from the conforming formulation simply by replacing
In this paper, we keep the same finite element space as DG method, replace the boundary function with the average of the inner function, and use the weak gradient arising from local Raviart-Thomas (RT) elements [5] to approximate the classic gradient. Moreover, the derivation process in this paper is based on rectangular RT elements [16]. Error estimates of optimal order are established for the corresponding conforming DG approximation in both a discrete
The rest of this paper is organized as follows: In Section 2, we shall present the conforming DG finite element scheme for the Poisson equation on rectangular partitions. Section 3 is devoted to a discussion of the stability and solvability of the new method. In Section 4, we shall prepare ourselves for error estimates by deriving some identities. Error estimates of optimal order in
Throughout this paper, we adopt the standard definition of Sobolev space
H10(Ω)={v∈H1(Ω):v|∂Ω=0}, |
and the space
H(div,Ω)={q∈[L2(Ω)]d:∇⋅q∈L2(Ω)}. |
Assume that the domain
For any interior edge
{v}=12(v|∂T1+v|∂T2),[[v]]=v|∂T1n1+v|∂T2n2, | (6) |
where
{v}=v|eand[[v]]=v|en. | (7) |
We define a discontinuous finite element space
Vh={v∈L2(Ω):v|T∈Qk(T),∀T∈Th}, | (8) |
and its subspace
V0h={v∈Vh:v=0on∂Ω}, | (9) |
where
Definition 2.1. For a given
(∇dv,q)T:=−(v,∇⋅q)T+⟨{v},q⋅n⟩∂T,∀q∈RTk(T), | (10) |
where
The weak gradient operator
(∇dv)|T=∇d(v|T). |
We introduce the following bilinear form:
a(v,w)=(∇dv,∇dw), |
the conforming DG algorithm to solve the problems (1) - (2) is given by
Conforming DG algorithm 1. Find
a(uh,vh)=(f,vh),∀vh∈V0h, | (11) |
where
We will prove the existence and uniqueness of the solution of equation (11). Firstly, we present the following two useful inequalities to derive the forthcoming analysis.
Lemma 3.1 (trace inequality). Let
‖φ‖2e≤C(h−1T‖φ‖2T+hT‖∇φ‖2T), | (12) |
where
Lemma 3.2 (inverse inequality). Let
‖∇φ‖T≤C(n)h−1T‖φ‖T,∀T∈Th. | (13) |
Then, we define the following semi-norms in the discontinuous finite element space
|||v|||2=a(v,v)=∑T∈Th‖∇dv‖2T, | (14) |
‖v‖21,h=∑T∈Th‖∇v‖2T+∑e∈E0hh−1e‖[[v]]‖2e. | (15) |
We have the equivalence between the semi-norms
Lemma 3.3. For any
C1‖v‖1,h≤|||v|||≤C2‖v‖1,h, | (16) |
where
Proof. It follows from the definition of
‖∇dv‖2T1=(∇dv,∇dv)T1=−(v,∇⋅∇dv)T1+⟨{v}n,∇dv⟩∂T1=(∇v,∇dv)T1−⟨(v−{v})n,∇dv⟩∂T1≤‖∇v‖T1‖∇dv‖T1+‖(v−{v})n‖∂T1‖∇dv‖∂T1≤‖∇dv‖T1(‖∇v‖T1+h−12T1‖(v−{v})n‖∂T1). | (17) |
For any
(v−{v})|en1=v|∂T1n1−12(v|∂T1+v|∂T2)n1=12(v|∂T1n1+v|∂T2n2)=12[[v]]e. |
Then we can get
‖(v−{v})n‖2∂T1≤12∑e∈∂T1‖[[v]]‖2e. | (18) |
Substituting (18) into (17) gives
‖∇dv‖2T1≤C2‖∇dv‖T1(‖∇v‖T1+∑e∈∂T1h−12e‖[[v]]‖e), |
this completes the proof of the right-hand of (16).
To prove the left-hand of (16), we consider the subspace of
D(k,T):={q∈RTk(T):q⋅n=0on∂T}. |
Note that
‖∇v‖T=supq∈D(k,T)(∇v,q)T‖q‖T. | (19) |
Using the integration by parts, Cauchy-Schwarz inequality, the definition of
(∇v,q)T=−(v,∇⋅q)T+⟨v,q⋅n⟩∂T=(∇dv,q)T−⟨{v},q⋅n⟩∂T=(∇dv,q)T≤‖∇dv‖T⋅‖q‖T, |
where we have used the fact that
‖∇v‖T≤‖∇dv‖T. | (20) |
We define the space
‖[[v]]‖e=supq∈De(k,T)⟨[[v]],q⋅n⟩e‖q⋅n‖e. | (21) |
Following the integration by parts and the definition of
(∇dv,q)T=(∇v,q)T−⟨v,q⋅n⟩e+⟨{v},q⋅n⟩e. |
Together with (20), we obtain
|⟨[[v]],q⋅n⟩e|=2|(∇dv,q)T−(∇v,q)T|≤2|(∇dv,q)T|+2|(∇v,q)T|≤C(‖∇dv‖T‖q‖T+‖∇v‖T‖q‖T)≤C‖∇dv‖T‖q‖T. |
Substituting the above inequality into (21), by the scaling argument [13], for such
‖[[v]]‖e≤C‖∇dv‖T‖q‖T‖q⋅n‖e≤Ch12‖∇dv‖T. | (22) |
Combining (20) and (22) gives a proof of the left-hand of (16).
Lemma 3.4. The semi-norm
Proof. We shall only verify the positivity property for
The above two lemmas imply the well posedness of the scheme (11). We prove the existence and uniqueness of solution of the conforming DG method in Theorem 3.1.
Theorem 3.1. The conforming DG scheme (11) has and only has one solution.
Proof. To prove the scheme (11) is uniquely solvable, it suffices to verify that the homogeneous equation has zero as its unique solution. To this end, let
a(uh,uh)=0, |
which leads to
In this section, we will derive an error equation which will be used for the error estimates. For any
(∇⋅q,v)T=(∇⋅Πhq,v)T,∀v∈Qk(T). | (23) |
For any
‖Πh(∇w)−∇w‖≤Chk‖w‖1+k. | (24) |
Moreover, it is easy to verify the following property holds true.
Lemma 4.1. For any
∑T∈Th(−∇⋅q,v)T=∑T∈Th(Πhq,∇dv)T,∀v∈V0h. | (25) |
Proof.
∑T∈Th⟨{v},Πhq⋅n⟩∂T=0. | (26) |
By the definition of
∑T∈Th(−∇⋅q,v)T=∑T∈Th(−∇⋅Πhq,v)T=∑T∈Th(−∇⋅Πhq,v)T+∑T∈Th⟨{v},Πhq⋅n⟩∂T=∑T∈Th(Πhq,∇dv)T. |
This completes the proof of the lemma.
Before establishing the error equation, we define a continuous finite element subspace of
˜Vh={v∈H1(Ω):v|T∈Qk(T),∀T∈Th}. | (27) |
so as a subspace of
˜V0h:={v∈˜Vh:v|∂Ω=0}. | (28) |
Lemma 4.2. For any
∇dv=∇v. |
Proof. By the definition of
(∇dv,q)T=−(v,∇⋅q)T+⟨{v},q⋅n⟩∂T=−(v,∇⋅q)T+⟨v,q⋅n⟩∂T=(∇v,q)T, |
which gives
(∇dv−∇v,q)T=0,∀q∈RTk(T). |
Letting
Let
‖Ihu−u‖≤Chk+1‖u‖k+1, | (29) |
‖∇Ihu−∇u‖≤Chk‖u‖k+1. | (30) |
It is obvious that
Lemma 4.3. Denote
a(eh,vh)=lu(vh), | (31) |
where
lu(vh)=∑T∈Th(∇Ihu−Πh∇u,∇dvh). | (32) |
Proof. Since
∑T∈Th(∇dIhu,∇dvh)T=∑T∈Th(∇Ihu,∇dvh)T=∑T∈Th(∇Ihu−Πh∇u+Πh∇u,∇dvh)T=∑T∈Th(∇Ihu−Πh∇u,∇dvh)T+∑T∈Th(Πh∇u,∇dvh)T=lu(vh)−∑T∈Th(∇⋅∇u,vh)T=lu(vh)+(f,vh). |
By the definition of the scheme (11), we have
∑T∈Th(∇dIhu−∇duh,∇dvh)T=lu(vh). |
This completes the proof of the lemma.
The goal of this section is to derive the error estimates in
Theorem 5.1. Let
|||eh|||≤Chk|u|k+1. | (33) |
Proof. Letting
|||eh|||2=lu(eh). | (34) |
From the Cauchy-Schwarz inequality, the triangle inequality, the definition of
lu(vh)=∑T∈Th(∇Ihu−Πh(∇u),∇dvh)T≤∑T∈Th‖∇Ihu−Πh(∇u)‖T‖∇dvh‖T≤(∑T∈Th‖∇Ihu−Πh(∇u)‖2T)12(∑T∈Th‖∇dvh‖2T)12=(∑T∈Th‖∇Ihu−∇u+∇u−Πh(∇u)‖2T)12|||vh|||≤(∑T∈Th‖∇Ihu−∇u‖2T+‖∇u−Πh(∇u)‖2T)12|||vh|||≤Chk|u|k+1|||vh|||. |
Then, we have
lu(eh)≤Chk|u|k+1|||eh|||. | (35) |
Substituting (35) to (34), we obtain
|||eh|||2≤Chk|u|k+1|||eh|||, |
which completes the proof of the lemma.
It is obvious that
(∇˜uh,∇v)=(f,v),∀v∈˜V0h. | (36) |
For any
(∇duh−∇˜uh,∇v)=0,∀v∈˜V0h. | (37) |
In the rest of this section, we derive an optimal order error estimate for the conforming DG approximation (11) in
−∇⋅(∇Φ)=uh−˜uh,inΩ. | (38) |
Assume that the dual problem satisfies
‖Φ‖2≤C‖uh−˜uh‖. | (39) |
In the following of this paper, we note
Theorem 5.2. Assume
‖u−uh‖≤Chk+1|u|k+1. | (40) |
Proof. First, we shall derive the optimal order for
a(Φh,v)=(εh,v),∀v∈V0h. | (41) |
Since
(∇duh−∇˜uh,∇IhΦ)=0,∇dIhΦ=∇IhΦ, |
which gives
(∇duh−∇˜uh,∇dIhΦ)=0. | (42) |
Setting
‖εh‖2=a(Φh,εh)=∑T∈Th(∇dΦh,∇dεh)T=∑T∈Th(∇d(Φh−IhΦ),∇duh−∇˜uh)T≤|||Φh−IhΦ|||(|||uh−Ihu|||+‖∇(Ihu−˜uh)‖). |
Then, by the Cauchy-Schwarz inequality, (33) and (39), we obtain
‖εh‖2≤Ch|Φ|2hk|u|k+1≤Chk+1|u|k+1‖εh‖, |
which gives
‖εh‖≤Chk+1|u|k+1. | (43) |
Combining the error estimate of finite element solution, the triangle inequality and (43) yields (40), which completes the proof of the theorem.
In this section, we shall present some numerical results for the conforming discontinuous Galerkin method analyzed in the previous sections.
We solve the following Poisson equation on the unit square domain
−Δu=2π2sin(πx)sin(πy)in Ω | (44) |
u=0on ∂Ω. | (45) |
The exact solution of the above problem is
We first use the
level | rate | rate | |||
by |
|||||
6 | 0.1996E-02 | 1.97 | 0.8887E-02 | 1.98 | 1024 |
7 | 0.5013E-03 | 1.99 | 0.2228E-02 | 2.00 | 4096 |
8 | 0.1255E-03 | 2.00 | 0.5574E-03 | 2.00 | 16384 |
by |
|||||
6 | 0.2427E-02 | 1.97 | 0.1027E+00 | 1.02 | 3072 |
7 | 0.6100E-03 | 1.99 | 0.5105E-01 | 1.01 | 12288 |
8 | 0.1527E-03 | 2.00 | 0.2546E-01 | 1.00 | 49152 |
by |
|||||
5 | 0.1533E-03 | 3.00 | 0.2042E-01 | 2.03 | 1536 |
6 | 0.1915E-04 | 3.00 | 0.5061E-02 | 2.01 | 6144 |
7 | 0.2394E-05 | 3.00 | 0.1260E-02 | 2.01 | 24576 |
by |
|||||
5 | 0.7959E-05 | 4.00 | 0.1965E-02 | 3.00 | 2560 |
6 | 0.4971E-06 | 4.00 | 0.2451E-03 | 3.00 | 10240 |
7 | 0.3140E-07 | 3.98 | 0.3059E-04 | 3.00 | 40960 |
by |
|||||
4 | 0.1055E-04 | 4.97 | 0.1421E-02 | 4.05 | 960 |
5 | 0.3314E-06 | 4.99 | 0.8735E-04 | 4.02 | 3840 |
6 | 0.1057E-07 | 4.97 | 0.5417E-05 | 4.01 | 15360 |
by |
|||||
2 | 0.2835E-02 | 6.24 | 0.1450E+00 | 5.49 | 84 |
3 | 0.4532E-04 | 5.97 | 0.4718E-02 | 4.94 | 336 |
4 | 0.7115E-06 | 5.99 | 0.1478E-03 | 5.00 | 1344 |
The same test case is also computed using the
level | rate | rate | |||
by |
|||||
6 | 0.4006E-03 | 1.99 | 0.2389E-02 | 1.99 | 4096 |
7 | 0.1003E-03 | 2.00 | 0.5982E-03 | 2.00 | 16384 |
8 | 0.2510E-04 | 2.00 | 0.1496E-03 | 2.00 | 65536 |
by |
|||||
6 | 0.2360E-04 | 2.99 | 0.3186E-02 | 1.99 | 9216 |
7 | 0.2953E-05 | 3.00 | 0.7976E-03 | 2.00 | 36864 |
8 | 0.3692E-06 | 3.00 | 0.1995E-03 | 2.00 | 147456 |
by |
|||||
5 | 0.1413E-04 | 4.08 | 0.1650E-02 | 2.97 | 4096 |
6 | 0.8676E-06 | 4.03 | 0.2072E-03 | 2.99 | 16384 |
7 | 0.5398E-07 | 4.01 | 0.2593E-04 | 3.00 | 65536 |
by |
|||||
3 | 0.2226E-02 | 4.59 | 0.5414E-01 | 3.52 | 400 |
4 | 0.9610E-04 | 4.53 | 0.3723E-02 | 3.86 | 1600 |
5 | 0.3279E-05 | 4.87 | 0.2392E-03 | 3.96 | 6400 |
To test the superconvergence of
−Δu+u=fin Ωu=0on ∂Ω, |
where
u=(x−x2)(y−y3). | (46) |
Uniform square grids as shown in Figure 1 are used for numerical computation. The numerical results are listed in Table 3. Surprising, for this problem, the
level | rate | rate | |||
by |
|||||
3 | 0.8265E-02 | 1.06 | 0.4577E-01 | 1.14 | 16 |
4 | 0.2772E-02 | 1.58 | 0.1732E-01 | 1.40 | 64 |
5 | 0.7965E-03 | 1.80 | 0.6331E-02 | 1.45 | 256 |
6 | 0.2142E-03 | 1.90 | 0.2290E-02 | 1.47 | 1024 |
7 | 0.5564E-04 | 1.94 | 0.8213E-03 | 1.48 | 4096 |
8 | 0.1419E-04 | 1.97 | 0.2928E-03 | 1.49 | 16384 |
To test further the superconvergence of
−∇(a∇u)=fin Ωu=0on ∂Ω, |
where
u=(x−x3)(y2−y3). | (47) |
Uniform square grids as shown in Figure 1 are used for computation. The numerical results are listed in Table 4. Surprising, again, the
level | rate | rate | |||
by |
|||||
3 | 0.4929E-02 | 0.97 | 0.5371E-01 | 0.80 | 16 |
4 | 0.1917E-02 | 1.36 | 0.2401E-01 | 1.16 | 64 |
5 | 0.6004E-03 | 1.67 | 0.9407E-02 | 1.35 | 256 |
6 | 0.1682E-03 | 1.84 | 0.3507E-02 | 1.42 | 1024 |
7 | 0.4457E-04 | 1.92 | 0.1275E-02 | 1.46 | 4096 |
8 | 0.1148E-04 | 1.96 | 0.4576E-03 | 1.48 | 16384 |
In this paper, we establish a new numerical approximation scheme based on the rectangular partition to solve second order elliptic equation. We derived the numerical scheme and then proved the optimal order of convergence of the error estimates in
[1] | Johnstone D, Moore G, Tausz M, et al. (2010) The measurement of wood decay in landscape trees. Arboric Urban For 36: 121-127. |
[2] | Sambuelli L, Socco LV, Godio A, et al. (2003) Ultrasonic, electric and radar measurements for living trees assessment. B Geofis Teor Appl 44: 253-279. |
[3] |
Rabe C, Ferner D, Fink S, et al. (2004) Detection of decay in trees with stress waves and interpretation of acoustic tomograms. Arboricultural J 28: 3-19. doi: 10.1080/03071375.2004.9747399
![]() |
[4] | Göcke L, Rust S, Weihs U, et al. (2007) Combining sonic and electrical impedance tomography for the non-destructive testing of trees. Proceedings of 15th international symposium on nondestructive testing of wood. Forest products society. Madison, 31-42. |
[5] |
Leong EC, Burcham DC, Yok-King F (2012) A purposeful classification of tree decay detection tools. Arboricultural J Int J Urban For 34: 91-115. doi: 10.1080/03071375.2012.701430
![]() |
[6] |
Mattheck CG, Breloer H (1994) Field guide for visual tree assessment (VTA). Arboricultural J Int J Urban For 18: 1-23. doi: 10.1080/03071375.1994.9746995
![]() |
[7] | Hayes E (2001) Evaluating Tree Defects, Safe trees, 2 Eds., Rochester, 34. |
[8] | Luley CL (2005) Wood decay fungi common to living urban trees in the Northeast and Central United States, Naples: Urban Forestry, 61. |
[9] |
Rinn F, Schweingruber FH, Schiir E (1996) RESISTOGRAPH and X-ray density charts of wood comparative evaluation of drill resistance profiles and X-ray density charts of different wood species. Holzforschung 50: 303-311. doi: 10.1515/hfsg.1996.50.4.303
![]() |
[10] |
Martin T, Günther T (2013) Complex Resistivity Tomography (CRT) for fungus detection on standing oak trees. Eur J For Res 132 : 765-776. doi: 10.1007/s10342-013-0711-4
![]() |
[11] | Gilbert EA, Smiley ET (2004) Picus sonic tomography for the quantification of decay in white oak (Quercus alba) and hickory (Carya spp.). J Arboric 30: 277-280. |
[12] | Wessoly L (1995) Fracture diagnosis of trees. Part 3: Boring is no way for reliable fracture diagnosis. Stadt Und Grün 9: 635-640. |
[13] |
Hagrey SA (2007) Geophysical imaging of root-zone, trunk, and moisture heterogeneity. J Exp Bot 58: 839-854. doi: 10.1093/jxb/erl237
![]() |
[14] | Bassuk N, Grabosky J, Mucciardi A, et al. (2011) Ground-penetrating radar accurately locates tree roots in two soil media under pavement. Arboric Urban For 37: 160-166. |
[15] |
Zhu SP, Huang CL, Su Y, et al. (2014) 3D Ground penetrating radar to detect tree roots and estimate root biomass in the field. Remote Sens 6: 5754-5773. doi: 10.3390/rs6065754
![]() |
[16] |
Satriani A, Loperte A, Proto M, et al. (2010) Building damage caused by tree roots: laboratory experiments of GPR and ERT surveys. Adv Geosci 24: 133-137. doi: 10.5194/adgeo-24-133-2010
![]() |
[17] | Lorenzo H, Pérez-Gracia V, Novo A, et al. (2010) Forestry applications of ground-penetrating radar. For Syst 19: 5-17. |
[18] | Muller W (2003) Timber girder inspection using ground penetrating radar. OR Insight 45: 809-812. |
[19] | Pyakurel S (2009) 2D and 3D GPR imaging of wood and fiber reinforced polymer composites. West Virginia University, Morgantown |
[20] |
Halabe UB, Agrawal S, Gopalakrishnan B (2009) Non-destructive evaluation of wooden logs using ground penetrating radar. Nondestr Test Eval 24: 329-346. doi: 10.1080/10589750802474344
![]() |
[21] | Colla C (2010) GPR of a timber structural element. Proceedings of 13th international conference on ground penetrating radar, Lecce. |
[22] |
Rodriguez-Abad I, et al. (2010) Non-destructive methodologies for the evaluation of moisture content in sawn timber structures: ground-penetrating radar and ultrasound techniques. Near Surf Geophys 8: 475-482. doi: 10.3997/1873-0604.2010048
![]() |
[23] | Jezova J, Lambot S (2019) Influence of bark surface roughness on tree trunk radar inspection. Ground Penetrating Radar 2: 1-25. Available from: https://doi.org/10.26376/GPR2019001. |
[24] | Nicolotti G, Socco LV, Martinis R, et al. (2003) Application and comparison of three tomographic techniques for detection of decay in trees. J Arboric 29: 66-78. |
[25] | Sudakova MS, Kalashnikov AU, Terentieva EB (2017) GPR tomography as applied to delineation of voids. Construction Unique Build Struct 60: 57. |
[26] |
Sihvola A (2000) Mixing rules with complex dielectric coefficients. Subsurf Sens Technol Appl 1: 393-415. doi: 10.1023/A:1026511515005
![]() |
[27] |
Huisman JA, Hubbard SS, Redman JD, et al. (2003) Measuring soil water content with ground penetrating radar: a review. Vadose Zone J 2: 476-491. doi: 10.2136/vzj2003.4760
![]() |
[28] | Hartley I, Hamza MF (2016) Wood: Moisture Content, Hygroscopicity, and Sorption. In: Saleem Hashmi (Editor in chief), Reference Module in Materials Science and Materials Engineering, Oxford: Elsevier. |
[29] | Torgovnikov G (1993) Dielectric properties of wood and wood-based materials. Berlin, Heidelberg, New-York: Springer-Verlag, 196. |
[30] |
Sahin H, Ay N (2004) Dielectric properties of hardwood species at microwave frequencies. J Wood Sci 50: 375-380. doi: 10.1007/s10086-003-0575-1
![]() |
[31] | Sbartaï, ZM, Lataste JF (2008) Evaluation non-destructive de l'humidité du bois par mesures radar et résistivité électrique. Rapport interne, Université bordeaux 1, US2B, Ghymac. |
[32] | Jol HM (2009) Ground penetrating radar theory and applications, Oxford: Elsevier, 523. |
[33] | Kabir MF, Khalid KB, Daud WM, et al. (1997) Dielectric properties of rubber wood at microwave frequencies measured with an open-ended coaxial line. Wood Fiber Sci 29: 319-324. |
[34] | Afzal MT, Colpits B, Galik K (2003) Dielectric Properties of Softwood Species Measured with an Open-ended Coaxial Probe. In proceedings 8th International IUFRO Wood Drying Conference, Brasvo, Romania, August, 24-29, 2003; 110-115. |
[35] |
Koubaa A, Perré P, Hutcheon RM, et al. (2008) Complex Dielectric Properties of the Sapwood of Aspen, White Birch, Yellow Birch, and Sugar Maple. Drying Technol 26: 568-578. doi: 10.1080/07373930801944762
![]() |
[36] |
Topp GC, Davis JL, Annan AP (1980) Electromagnetic determination of soil water content: measurements in coaxial transmission lines. Water Resour Res 16: 574-582. doi: 10.1029/WR016i003p00574
![]() |
[37] | Lin G, Lv JX, Wen J (2014) Research on the relationship between moisture content and the dielectric constant of the tree trunk by the radar wave. Comput Modell New Technol 18: 1171-1175. |
[38] |
Razafindratsima S, Sbartaï ZM, Demontoux F (2017) Permittivity measurement of wood material over a wide range of moisture content. Wood Sci Technol 51: 1421-1431. doi: 10.1007/s00226-017-0935-4
![]() |
[39] |
Acuna L, Basterra LA, Casado M, et al. (2011) Application of resistograph to obtain the density and to differentiate wood species. Mater Construcc 61: 451-464. doi: 10.3989/mc.2010.57610
![]() |
1. | Xiu Ye, Shangyou Zhang, A weak divergence CDG method for the biharmonic equation on triangular and tetrahedral meshes, 2022, 178, 01689274, 155, 10.1016/j.apnum.2022.03.017 | |
2. | Jun Zhou, Da Xu, Wenlin Qiu, Leijie Qiao, An accurate, robust, and efficient weak Galerkin finite element scheme with graded meshes for the time-fractional quasi-linear diffusion equation, 2022, 124, 08981221, 188, 10.1016/j.camwa.2022.08.022 | |
3. | Xiu Ye, Shangyou Zhang, A conforming discontinuous Galerkin finite element method for the Stokes problem on polytopal meshes, 2021, 93, 0271-2091, 1913, 10.1002/fld.4959 | |
4. | Xiu Ye, Shangyou Zhang, Constructing a CDG Finite Element with Order Two Superconvergence on Rectangular Meshes, 2023, 2096-6385, 10.1007/s42967-023-00330-5 | |
5. | Yan Yang, Xiu Ye, Shangyou Zhang, A pressure-robust stabilizer-free WG finite element method for the Stokes equations on simplicial grids, 2024, 32, 2688-1594, 3413, 10.3934/era.2024158 | |
6. | Xiu Ye, Shangyou Zhang, A superconvergent CDG finite element for the Poisson equation on polytopal meshes, 2023, 0044-2267, 10.1002/zamm.202300521 | |
7. | Xiu Ye, Shangyou Zhang, Two-Order Superconvergent CDG Finite Element Method for the Heat Equation on Triangular and Tetrahedral Meshes, 2024, 2096-6385, 10.1007/s42967-024-00444-4 | |
8. | Xiu Ye, Shangyou Zhang, Order two superconvergence of the CDG finite elements for non-self adjoint and indefinite elliptic equations, 2024, 50, 1019-7168, 10.1007/s10444-023-10100-9 | |
9. | Fuchang Huo, Weilong Mo, Yulin Zhang, A locking-free conforming discontinuous Galerkin finite element method for linear elasticity problems, 2025, 465, 03770427, 116582, 10.1016/j.cam.2025.116582 |
level | rate | rate | |||
by |
|||||
6 | 0.1996E-02 | 1.97 | 0.8887E-02 | 1.98 | 1024 |
7 | 0.5013E-03 | 1.99 | 0.2228E-02 | 2.00 | 4096 |
8 | 0.1255E-03 | 2.00 | 0.5574E-03 | 2.00 | 16384 |
by |
|||||
6 | 0.2427E-02 | 1.97 | 0.1027E+00 | 1.02 | 3072 |
7 | 0.6100E-03 | 1.99 | 0.5105E-01 | 1.01 | 12288 |
8 | 0.1527E-03 | 2.00 | 0.2546E-01 | 1.00 | 49152 |
by |
|||||
5 | 0.1533E-03 | 3.00 | 0.2042E-01 | 2.03 | 1536 |
6 | 0.1915E-04 | 3.00 | 0.5061E-02 | 2.01 | 6144 |
7 | 0.2394E-05 | 3.00 | 0.1260E-02 | 2.01 | 24576 |
by |
|||||
5 | 0.7959E-05 | 4.00 | 0.1965E-02 | 3.00 | 2560 |
6 | 0.4971E-06 | 4.00 | 0.2451E-03 | 3.00 | 10240 |
7 | 0.3140E-07 | 3.98 | 0.3059E-04 | 3.00 | 40960 |
by |
|||||
4 | 0.1055E-04 | 4.97 | 0.1421E-02 | 4.05 | 960 |
5 | 0.3314E-06 | 4.99 | 0.8735E-04 | 4.02 | 3840 |
6 | 0.1057E-07 | 4.97 | 0.5417E-05 | 4.01 | 15360 |
by |
|||||
2 | 0.2835E-02 | 6.24 | 0.1450E+00 | 5.49 | 84 |
3 | 0.4532E-04 | 5.97 | 0.4718E-02 | 4.94 | 336 |
4 | 0.7115E-06 | 5.99 | 0.1478E-03 | 5.00 | 1344 |
level | rate | rate | |||
by |
|||||
6 | 0.4006E-03 | 1.99 | 0.2389E-02 | 1.99 | 4096 |
7 | 0.1003E-03 | 2.00 | 0.5982E-03 | 2.00 | 16384 |
8 | 0.2510E-04 | 2.00 | 0.1496E-03 | 2.00 | 65536 |
by |
|||||
6 | 0.2360E-04 | 2.99 | 0.3186E-02 | 1.99 | 9216 |
7 | 0.2953E-05 | 3.00 | 0.7976E-03 | 2.00 | 36864 |
8 | 0.3692E-06 | 3.00 | 0.1995E-03 | 2.00 | 147456 |
by |
|||||
5 | 0.1413E-04 | 4.08 | 0.1650E-02 | 2.97 | 4096 |
6 | 0.8676E-06 | 4.03 | 0.2072E-03 | 2.99 | 16384 |
7 | 0.5398E-07 | 4.01 | 0.2593E-04 | 3.00 | 65536 |
by |
|||||
3 | 0.2226E-02 | 4.59 | 0.5414E-01 | 3.52 | 400 |
4 | 0.9610E-04 | 4.53 | 0.3723E-02 | 3.86 | 1600 |
5 | 0.3279E-05 | 4.87 | 0.2392E-03 | 3.96 | 6400 |
level | rate | rate | |||
by |
|||||
3 | 0.8265E-02 | 1.06 | 0.4577E-01 | 1.14 | 16 |
4 | 0.2772E-02 | 1.58 | 0.1732E-01 | 1.40 | 64 |
5 | 0.7965E-03 | 1.80 | 0.6331E-02 | 1.45 | 256 |
6 | 0.2142E-03 | 1.90 | 0.2290E-02 | 1.47 | 1024 |
7 | 0.5564E-04 | 1.94 | 0.8213E-03 | 1.48 | 4096 |
8 | 0.1419E-04 | 1.97 | 0.2928E-03 | 1.49 | 16384 |
level | rate | rate | |||
by |
|||||
3 | 0.4929E-02 | 0.97 | 0.5371E-01 | 0.80 | 16 |
4 | 0.1917E-02 | 1.36 | 0.2401E-01 | 1.16 | 64 |
5 | 0.6004E-03 | 1.67 | 0.9407E-02 | 1.35 | 256 |
6 | 0.1682E-03 | 1.84 | 0.3507E-02 | 1.42 | 1024 |
7 | 0.4457E-04 | 1.92 | 0.1275E-02 | 1.46 | 4096 |
8 | 0.1148E-04 | 1.96 | 0.4576E-03 | 1.48 | 16384 |
level | rate | rate | |||
by |
|||||
6 | 0.1996E-02 | 1.97 | 0.8887E-02 | 1.98 | 1024 |
7 | 0.5013E-03 | 1.99 | 0.2228E-02 | 2.00 | 4096 |
8 | 0.1255E-03 | 2.00 | 0.5574E-03 | 2.00 | 16384 |
by |
|||||
6 | 0.2427E-02 | 1.97 | 0.1027E+00 | 1.02 | 3072 |
7 | 0.6100E-03 | 1.99 | 0.5105E-01 | 1.01 | 12288 |
8 | 0.1527E-03 | 2.00 | 0.2546E-01 | 1.00 | 49152 |
by |
|||||
5 | 0.1533E-03 | 3.00 | 0.2042E-01 | 2.03 | 1536 |
6 | 0.1915E-04 | 3.00 | 0.5061E-02 | 2.01 | 6144 |
7 | 0.2394E-05 | 3.00 | 0.1260E-02 | 2.01 | 24576 |
by |
|||||
5 | 0.7959E-05 | 4.00 | 0.1965E-02 | 3.00 | 2560 |
6 | 0.4971E-06 | 4.00 | 0.2451E-03 | 3.00 | 10240 |
7 | 0.3140E-07 | 3.98 | 0.3059E-04 | 3.00 | 40960 |
by |
|||||
4 | 0.1055E-04 | 4.97 | 0.1421E-02 | 4.05 | 960 |
5 | 0.3314E-06 | 4.99 | 0.8735E-04 | 4.02 | 3840 |
6 | 0.1057E-07 | 4.97 | 0.5417E-05 | 4.01 | 15360 |
by |
|||||
2 | 0.2835E-02 | 6.24 | 0.1450E+00 | 5.49 | 84 |
3 | 0.4532E-04 | 5.97 | 0.4718E-02 | 4.94 | 336 |
4 | 0.7115E-06 | 5.99 | 0.1478E-03 | 5.00 | 1344 |
level | rate | rate | |||
by |
|||||
6 | 0.4006E-03 | 1.99 | 0.2389E-02 | 1.99 | 4096 |
7 | 0.1003E-03 | 2.00 | 0.5982E-03 | 2.00 | 16384 |
8 | 0.2510E-04 | 2.00 | 0.1496E-03 | 2.00 | 65536 |
by |
|||||
6 | 0.2360E-04 | 2.99 | 0.3186E-02 | 1.99 | 9216 |
7 | 0.2953E-05 | 3.00 | 0.7976E-03 | 2.00 | 36864 |
8 | 0.3692E-06 | 3.00 | 0.1995E-03 | 2.00 | 147456 |
by |
|||||
5 | 0.1413E-04 | 4.08 | 0.1650E-02 | 2.97 | 4096 |
6 | 0.8676E-06 | 4.03 | 0.2072E-03 | 2.99 | 16384 |
7 | 0.5398E-07 | 4.01 | 0.2593E-04 | 3.00 | 65536 |
by |
|||||
3 | 0.2226E-02 | 4.59 | 0.5414E-01 | 3.52 | 400 |
4 | 0.9610E-04 | 4.53 | 0.3723E-02 | 3.86 | 1600 |
5 | 0.3279E-05 | 4.87 | 0.2392E-03 | 3.96 | 6400 |
level | rate | rate | |||
by |
|||||
3 | 0.8265E-02 | 1.06 | 0.4577E-01 | 1.14 | 16 |
4 | 0.2772E-02 | 1.58 | 0.1732E-01 | 1.40 | 64 |
5 | 0.7965E-03 | 1.80 | 0.6331E-02 | 1.45 | 256 |
6 | 0.2142E-03 | 1.90 | 0.2290E-02 | 1.47 | 1024 |
7 | 0.5564E-04 | 1.94 | 0.8213E-03 | 1.48 | 4096 |
8 | 0.1419E-04 | 1.97 | 0.2928E-03 | 1.49 | 16384 |
level | rate | rate | |||
by |
|||||
3 | 0.4929E-02 | 0.97 | 0.5371E-01 | 0.80 | 16 |
4 | 0.1917E-02 | 1.36 | 0.2401E-01 | 1.16 | 64 |
5 | 0.6004E-03 | 1.67 | 0.9407E-02 | 1.35 | 256 |
6 | 0.1682E-03 | 1.84 | 0.3507E-02 | 1.42 | 1024 |
7 | 0.4457E-04 | 1.92 | 0.1275E-02 | 1.46 | 4096 |
8 | 0.1148E-04 | 1.96 | 0.4576E-03 | 1.48 | 16384 |