The frequency-dependent behavior in subsurface and through-the-wall media is analyzed in this paper as well as the formation of the Brillouin precursor waveforms inherently related to this feature. The emergence of these forerunners is presented as a plausible form to explain classical impairments observed in imaging technologies. The evolution of mono- and multicycle rectangular and first derivative Gaussian pulses through two dispersive media—concrete blocks and soil—is analyzed using a frequency-domain technique and the Debye dielectric model to characterize the media, at operating frequencies below 3 GHz. The frequency-domain approach facilitated to check the influence of some parameters considered critical for the precursor emergence—operating frequency, input pulse configuration, and internal structure of the underlying medium—results in a versatile tool suitable for any kind of modulated input signal propagated through any dispersive medium. The internal multireflection model has been considered as the most suitable model to describe the transmission process underlying both subsurface and through-wall imaging technologies. Two different moisture contents have been considered for concrete as a parameter to determine the performance of through-wall imaging radar from the precursor formation perspective. The results reveal that precursor is a phenomenon to take into account for application demanding larger signal-to-noise ratios.
1. Introduction
The frequency dependence of the dielectric properties in some materials is an important issue affecting the performance of radar imaging applications mainly in terms of both signal-to-noise ratio and pulse duration alterations. The frequency-dependent features can lead to the appearance of the waves, namely, Brillouin and Sommerfield precursors [1–4]. Despite precursors being a well-known phenomenon, we usually do not find the imaging and the expected precursor linkage, especially ultrawideband system configurations for dispersive medium.
The main distinguishing properties of the precursor waves have awakened interest for their probable use in different applications in a broad range of frequency bands, such as radio standoff inspection, biosensing, subsurface (underwater/underground) communications or through-wall-tissue imaging [5–9]. For instance, the nonexponential peak amplitude decay of these precursor waves would improve the link range which also means an optimal penetration depth in a given media; walls made of concrete, therefore, present an ideal opportunity to detect and utilize the precursors to obtain high signal-to-noise ratio (SNR) and/or greater penetration depths.
The precursor formation in the transmitted or propagating signal, however, can produce undesired features such as pulse broadening and downshifting of the carrier frequency of the modulated pulses [1, 3]. The time spreading undergone by the signal through the dispersive media turns frequency dispersion into an extremely important impairment in receiver systems based on correlation processing [3].
Nevertheless, the appearance of the precursors depends on many factors, including the frequency window, enough frequency bandwidth, and the input signal configuration. In any case, neglecting or eliminating this forerunner component removes an important portion of the received signal power [4], and this fact can lead to under optimal receiver hardware performance as well as inaccurate inversion methods in imaging applications. The pulse broadening or the algebraically amplitude decay introduced in the received signal can be better understood by an exhaustive analysis of the precursor fields.
For the cases of through-wall imaging (TWI) and underground imaging (UGI), the frequency-dependent attenuation and phase constants of the walls constructed of concrete structures [7–11] can obscure the interior due to reduced SNR.
We apply a frequency-domain method based on considering the frequency response H(z,f) of the medium under study. This response is derived from the theoretically or experimentally achieved propagation constant, and it depends on the material thickness or distance travelled within the medium z. The frequency response can be seen as a frequency filter containing information about the effects of attenuation and phase for each frequency component of the potential signal travelling through it. We can then analyze the evolution of any input signal x(z,t) through the dispersive medium for any propagation depth z [3–6, 12–16].
Nevertheless, the transmission process undergone by the signal travelling through the material can strongly condition the output. For a medium consisting of slabs of materials rather than half-space medium, usually assumed in theoretical simulations, the double interface or boundary conditions present a finite thickness medium, and the underlying transmission process can give rise to internal successive reflections arising out on the material edges with a more than minor influence on the propagated field.
In order to extract the material frequency response to an incoming EM signal for TWI and UGI, an internal multireflection model is assumed thereby considering the different spurious components due to the internal reflections occurring inside the material considered as a finite slab in the TWI case or as a finite half-space or single slab in the UGI case. For the latter subsurface case, it would be plausible to consider a stack of slabs to model the gradually changing dielectric properties of the soil medium which can be seen as a superposition of different layers, from a superficial cover of turf, following sandy soil and dry soil till deep stratums of rocky soil.
For the purposes of this paper, we considered two pulses: the first derivative of the Gaussian pulse and a rectangular pulse. Monocycle and multicycles settings are presented to better illustrate the precursor performance. The dispersive media consisted of soil and concrete structures characterized by an extended Debye model for frequencies below 3 GHz [1, 10].
In Section 2 we introduce the main elements required to perform the analysis of the propagation through a dispersive medium: the theoretical formulation and the dielectric model as well as the frequency transfer of the medium. In Section 3 we present and discuss the results achieved for the simulations done for the soil and concrete models for this study. Finally, the paper is concluded in Section 4.
2. Theoretical Analysis of Dispersive Propagation
We assumed a sine-modulated signal x(t) as in (1), with f0 being the carrier frequency, and s(t) represents the baseband information-bearing signal of duration Tb:
(1)x(z=0,t)=s(tTb)·sin(2πf0t).
The modulated pulse x(t), or equivalently x(z=0,t), passes through the medium having a frequency response H(z,f), acting as a frequency filter or channel response, as indicated in (2). The final result, y(z,t), transformed to the time domain by an inverse Fourier Transform, is the signal affected by the dispersion phenomenon after travelling through the dispersive medium at a distance z, as expressed in (3):
(2)Y(z,f)=X(f)·H(z,f),(3)y(z,t)=F-1{X(z=0,f)·H(z,f)}.
Only a linear polarization is considered for this study. As can be noticed, all parameters controlling the pulse shape and the sinusoidal carrier frequency can easily be changed, within the constraints given by the bandwidth to be covered.
In (2), we observe that each frequency component X(fk) is differently affected upon propagation through the dispersive material, modeling the behavior of the different spectral components of the input pulse x(t). The material frequency response H(z,f) depends on the transmission process occurring through the medium. Several electromagnetic models have been developed trying to provide a realistic formulation of the transmission process of an EM wave through the material. When the material thickness dm is large compared to the wavelength and/or the material losses are high, the transmission and reflection coefficients, respectively, Tg and Rg of a plane wave on a dielectric slab can be approximated by the Fresnel’s model [17]. For this case, only one ray contribution is ensured to be received. When material thickness is not large enough, the internal successive reflections become a dominant phenomenon; the internal multireflection model (IMR) [18] seems a more accurate approach of the underlying transmission physics.
The transmission coefficient Tg is a frequency-dependent function and can be understood as the material frequency response H(z,f) also representing the boundary conditions due to the medium interfaces (points A and B in Figure 1), as well as the transmission through the medium to reach a given propagation distance z. According to the geometry explained in Figure 1, H(z,f) can be defined as in (4), analogous to the transmission coefficient Tg given for the IMR model [18]:
(4)H(f,z)=EBEA=Tam·Tma·e-jγm(f)z1+Γam·Γma·e-j2γm(f)z
with Tam, Tma, Γam, and Γma denoting the Fresnel transmission and reflection coefficients for the air-material (am) and material-air (ma) interfaces, with Tij=2εj(f)/(εi(f)+εj(f)) and Tij=1+Γij, i,j=a,m. The value of the propagation constant γm(f) required for (4) is derived from the theoretical εr(f) as indicated in
(5)γm(f)=α(f)+jβ(f)=j2πfcεrm′(f)-jεrm′′(f),
where εrm′(f) and εrm′′(f) are the real and imaginary parts of the complex relative permittivity, c is the speed of light in free space, and f is the frequency in Hz.
IMR transmission model for concrete and soil slabs.
From (2) and (5) we can infer that the assumed dielectric model becomes vital to achieve a correct estimation of signal propagation through the dispersive medium and consequently for an accurate evaluation of the precursor formation, as mentioned in [3–5]. For instance, the relative dielectric permittivity model also gives us an idea of the most suitable frequency window to achieve the formation of the precursors.
We have considered Debye models to characterize the soil and concrete media considered in this study. The complex relative dielectric permittivity for soil εr,soil(f) is given as a Rocard-Powles-Debye model [9] as indicated in
(6)εr,soil(f)=∑j=i=13aj(1-i·ωτj)·(1-i·ωτfj)
with constants a1=13.63, a2=0.33, and a3=0.293; rotational relaxation times τ1=1.82·10-6s, τ2=1.97·10-8s, τ3=2.36·10-10s; frictional relaxation times τf1=8.78·10-11s, τf2=1.12·10-10s, and τf3=1.00·10-15s; and with ω the frequency in rad/s.
The Debye model for the concrete structures is given in [10] as by (7):
(7)εr(ω)=ε∞+Δε1+j·(ω/ωr)2-j·(ω/ωr)·Δε1+(ω/ωr)2-j·σdcωε0,
where Δε=εs-ε∞ is the difference between the real part of the complex relative permittivity at DC and high frequency ωr is the relaxation angular frequency, and σdc is the DC electrical conductivity. The parameter values are summarized in Table 1 according to [10] for two moisture contents, 0.2% (dry concrete) and 12% (wet concrete).
Parameters for concrete Debye model as per [10].
Moisture content
εs
ε∞
ωr
σdc
0.2%
4.814
4.507
1.22·109
6.06·10-4
12%
12.84
7.42
1.64·109
20.6·10-3
In Figure 2 we show the frequency response H(z,f) for the two Debye media, soil, and concrete structures, as a function of the material thickness, or distance travelled within the medium, z/zd, normalized to the reference propagation depth zd=α-1(f0) meters, with α(f), the attenuation constant of the dispersive media; zd=1.0184m (soil), zd=2.9178m (dry concrete), and zd=0.1512m (wet concrete). The influence of the moisture content is not a determining parameter of the trend shown by H(z,f), and it mostly affects the lower frequencies.
Frequency response assuming: (a) soil and (b) concrete.
For simulation studies in this paper, we choose as input signals two different pulses: a rectangular pulse and first derivative of Gaussian (FDG) pulse, this latter as an approach of impulse case. The rectangular pulse is a good choice to show the pros and cons of the precursor formation. Even when widely used for UGI, in [5] we determined that the Gaussian pulse is not the best option to produce the emergence of the Brillouin and/or Sommerfield precursor with large amplitude levels due to the smoothness of its noninfinite leading edges. Because of this reason, we preferred to consider the first derivative Gaussian (FDG) that works like a modulated Gaussian pulse with a bandwidth and central frequency that are related to the time width of the Gaussian pulse, and it is related to larger amplitude level precursors due to its spread-out frequency spectrum. In case of sine carrier modulation of the baseband pulse p(t), central frequency of the Gaussian pulse would equal f0 corresponding to the carrier component. Other pulses have been analyzed as optimal [19], and they could also be applied further improving the results.
In (8) we indicate the settings used for the sequence of delta pulses and FDG:
(8)p(t)=-2(t-a)b2·exp{-(t-ab)2}
with (a,b)=(0.05·Tb,0.0375·Tb). The expression for the rectangular pulse just involves one parameter, the pulse width, or period, Tb:p(t)=∏(t/Tb), 0≤t≤Tb.
The input pulse configuration influences the formation and evolution of the resulting precursors [3–5]. In order to show this fact, we have considered two different settings, mono- and multicycles, by choosing Tb=1/f0 and Tb=10/f0 for both rectangular and FDG pulses. In Figure 3 we show the input pulses used for this study, modulated at a carrier frequency of f0=1.5 GHz. The resulting energy contents for both types of pulses are not equal.
Baseband rectangular and FDG pulses in time domain.
3. Dynamical Evolution within the Dispersive Medium for Drude Model
The dynamical evolution of the total field ET(z,t) is theoretically predicted using (1) for the Debye model and considering IMR structures of soil and concrete, for both rectangular and FDG pulses, modulated using a sine carrier frequency f0=1.5 GHz. From ET(z,t) we can analyze the peak amplitude decay law of the propagated field which helps to determine the penetration depth that the transmitted waveform can achieve into the given medium. The electric-field intensity and the frequency spectra are also analyzed.
In the results plotted, the propagation distance is normalized with respect to the reference distance zd that is defined as the propagation distance z at which the amplitude decay exp{-α(f)·z} of the signal travelling through a dispersive medium equals exp{-α(f0)·zd}=e-1, equivalently 1 Np, for a given frequency value, in our case f=f0 corresponding to the carrier component, with α(f) the attenuation constant of the medium.
3.1. Dynamical Evolution
In Figures 4, 5, and 6 we show the spatial and time evolution for ten propagation distances zd at f0=1.5 GHz assuming the IMR model for soil, dry concrete, and wet concrete (Figure 6), respectively. In order to emphasize the weaker values for the longer propagation distances, each signal has been normalized by its peak value.
Normalized electric field evolution through different propagation distances within a dielectric slab of soil for a rectangular pulse with f=1.5 GHz and Tb=10/f0.
Normalized electric field evolution through different propagation distances within a dielectric slab of dry concrete (0.2% moisture content) for a FDG pulse with f=1.5 GHz and Tb=10/f0.
Normalized electric field evolution through different propagation distances within a dielectric slab of wet concrete (12% moisture content) for a FDG pulse with f=1.5 GHz and Tb=1/f0.
For the rectangular pulse case with Tb=10/f0, shown in Figure 4, the Brillouin precursors are observable as a superimposed wave in the leading and trailing edges of the output rectangular pulse. For the monocycle setting cases of Figures 5 and 6, even that not so much graphically observable, we detect the appearance of the Brillouin precursor at the trailing edge of the rectangular and FDG pulses, explained by the broad lobule-shaped component indicated in Figure 5. For larger propagation distances, the Brillouin precursor becomes the predominant forerunner, and the remaining carrier frequency component is hardly observable.
For the wet concrete, the Brillouin precursor component is less noticeable, and the carrier frequency is also less attenuated, perhaps due to the larger moisture contents.
In Figures 7(a) and 7(b) we show the intensity of the normalized electric field at any given location (z,t):I(z,t)=|E(z,t)|2, estimated for a multicycle rectangular pulse and a monocycle FDG pulse through soil. In these plots, the amplitude levels present a range between 0 and 1, the blue pixels represent locations where the intensity is at minimum, and the red pixels correspond to locations of maximum intensity [19]. In Figure 7(a), we observe that the carrier cycles extinguish after a few propagation depths and the Brillouin precursors corresponding to the leading and trailing edges of the rectangular pulse evolve following a linear track. For the monocycle case in Figure 7(b), we observe just one track due to the fact of that both precursors merge after a short propagation distance.
Spatial evolution of the normalized electric field intensity distribution through a dielectric slab of soil at f0=1.5 GHz, IMR model: (a) rectangular pulse with Tb=10/f0 and (b) FDG pulse with Tb=1/f0. Amplitude levels range in the interval [0,1].
3.2. Peak Amplitude Decay
The maximum value of the evolving field ET(z=zk,t) at kth depth z=zk is plotted in Figures 8, 9, and 10 for ten normalized depths Δz=z/zd, for rectangular and FDG pulses with Tb=1/f0 and Tb=10/f0. These maxima values include the carrier and precursor components given that they have not been segregated. Figure 8 illustrates the peak decay amplitude assuming the soil model. The dry concrete case is shown in Figure 9 and wet concrete in Figure 10.
Peak amplitude as a function of propagation distance through soil for rectangular and FDG pulses of duration Tb=10/f0 and Tb=1/f0 at f0=1.5 GHz.
Peak amplitude as a function of propagation distance through dry concrete (0.2% moisture) for rectangular and FDG pulses of duration Tb=10/f0 and Tb=1/f0 at f0=1.5 GHz.
Peak amplitude as a function of propagation distance through wet concrete (12% moisture) for rectangular and FDG pulses of duration Tb=10/f0 and Tb=1/f0 at f0=1.5 GHz.
We conclude that the total propagated field ET(t,z) follows an algebraically decaying peak amplitude instead of an exponentially decay as given by the attenuation corresponding to the carrier frequency, with similar performance of both settings and mono- and multicycle, for the rectangular pulse; however, the pulse configuration takes a large influence for the FDG pulse even for shorter propagation distances.
We notice that for wet concrete the peak amplitude decay is slightly less attenuated, likely due to a larger presence of water.
3.3. Frequency Spectrum
In Figures 11(a) and 11(b) we show the normalized spectra for rectangular and FDG pulses after travelling through a propagation distance of 4·zd within soil, dry, and wet concrete, respectively. We observe that the FDG pulse case undergoes a larger alteration of the spectrum shape than that of a rectangular pulse; the increasing level of the low frequencies—larger as closer to zero—is an effect of the dispersion phenomenon more remarkable for the rectangular pulse.
Frequency spectra for rectangular and FDG pulses after propagating through z=4·zd through soil and concrete at f0=1.5 GHz, with pulse duration (a) Tb=10/f0 and (b) Tb=1/f0.
The percentage of downshifting with respect to the original frequency carrier f0 was measured as Δf0=100·(f0′-f0)/f0 and shown in Figures 12(a), 12(b), and 12(c), only for multicycle configuration and soil medium. The rectangular cases offer similar results in terms of frequency carrier downshifting; however, it becomes a critical impairment for the FDG pulse.
Frequency carrier shifting estimation: (a) soil, (b) dry concrete (0.2%), and (c) wet concrete (12%).
4. Conclusions
Generally speaking, the results presented in this paper indicate an undeniable effect of the frequency dispersion in subsurface and through-wall imaging applications, particularly for the soil and concrete structures considered here. The propagation model has been assumed as frequency dispersive and described by the internal multireflection model, and the output waveforms indicate a more than plausible malfunction of the receiver systems which are not usually built to consider the Brillouin precursor formation.
The method herein described provides a versatile framework to theoretically analyze the phenomenon of the Brillouin precursor formation and evolution through a given medium. This tool would be also useful to observe the influence of different variables such as the type of the input signal, the carrier frequency, the pulse width, the dielectric model, the bandwidth limiting, the multipath, the filtering, and so on [3–5]. All of them condition the precursor formation.
We demonstrated the importance of the precursor fields and the further research needed in this direction. Outcomes indicate that the propagation phenomenon, derived from the analysis of Brillouin precursor waveform, cannot be neglected, especially due to the large amount of energy carried by these components. In order to take advantage of the potential benefits of the precursor fields further research must be done, both in the experimental and theoretical field.
Unfortunately, there is a gap in experimental corroboration regarding the topic of precursors and microwave frequency band. In [5] we presented experimental results achieved by measuring the precursors arising from an ad hoc prepared soil layer. Even when limited to one type of soil and one thickness, the outcomes in [5] can be considered in good agreement with the model herein proposed, concluding from both works that the SNR can be enhanced if the precursor formation is considered, thereby improving the imaging for through-wall and underground scenarios.
In [5] we did not contemplate the influence of the underlying transmission model that is introduced here as a key to further understand the precursor formation and the related results for radar imaging purposes. The combination of larger SNR, that is, larger amplitude echoes, and the transmission model could be useful in the inversion analysis through correctly interpreting the backscattered signals as a part of the internal multireflection process.
In a pulse-dependent analysis, the processing of the received signal is attached to a priori configured pulse. Then, if precursors are formed during the transmission, it could derive a malfunction of the receiver system. Differently from a pulse-dependent enunciation, a frequency domain solution provides many degrees of freedom to achieve optimal processing results derived from the possibility of selecting any input pulse, frequency bandwidth windows, number of frequency components, and so on, once the transfer function H(z,f) is obtained, well experimentally from in situ scenario or from a theoretical setting. Different pulses can be analyzed in order to improve the SNR values so that weak echoes backscattered by plausible targets can be differentiated from the noise level.
Inversely, from the approach here described, an optimal pulse can be selected [19], in terms of penetration distance, amplitude decay, and time width spreading, for a specific scenario/medium in order to implement a pulse-dependent hardware.
In this paper we have defined a plausible theoretical scenario in order to show the potential advantages from considering the precursor formation. If the theoretical study is used as a way to design an optimal pulse, the practical performance of such waveform can be affected by a deviation from the assumed dielectric model and scenario parameters, as described in [3]. However, any other approach or theory can be affected by the same circumstances.
The sensitivity in Brillouin precursor formation with respect to uncertainty in the model is not as much critical as the effect of uncertainty in the dielectric properties. The precursor formation is more related to the dielectric properties of the underlying medium rather than the geometrical features of the propagation scenario. As aforementioned, the dielectric medium joint to the frequency bandwidth allocation and input pulse affects the amplitude and width of the precursor field; however, the geometrical model of the medium influences the multipath behavior and this influence is observable in the expression for the medium frequency response H(f,z), and therefore it must be taken into consideration for inverse techniques in radar imaging. In this paper we have considered a scenario model different from the ideal half-space medium, and as previously indicated any other scenarios can be object of a similar analysis.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by Galician Regional Government under Project EM2012/138. The authors would also like to thank the support received from the New Mexico State University (NMSU), Las Cruces, NM, USA.
OughstunK. E.20092New York, NY, USASpringerOughstunK. E.Dynamical evolution of the Brillouin precursor in Rocard-Powles-Debye model dielectrics2005535158215902-s2.0-1914437144910.1109/TAP.2005.846452AlejosA. V.DawoodM.MohammedH. U.Analysis of brillouin precursor propagation through foliage for digital sequences of pulses20118159632-s2.0-7865096279510.1109/LGRS.2010.2051213AlejosA. V.DawoodM.Estimation of power extinction factor in presence of brillouin precursor formation through dispersive media20112544554652-s2.0-7975147021410.1163/156939311794500313MohammedH. U. R.DawoodM.AlejosA. V.Experimental detection and characterization of Brillouin precursor through loamy soil at microwave frequencies20125024364452-s2.0-8485632841710.1109/TGRS.2011.2160983MohammedH.DawoodM.Vázquez AlejosA.Experimental detection of Brillouin precursors through tap water at microwave frequencies2010422516451647HoncharenkoW.BertoniH. L.Transmission and reflection characteristics at concrete block walls in the UHF bands proposed for future PCS19944222322392-s2.0-002836847010.1109/8.277217ThajudeenC.HoorfarA.AhmadF.DogaruT.Measured complex permittivity of walls with different hydration levels and the effect on power estimation of twri target returns2011301771992-s2.0-79956363503LandronO.FeuersteinM.RappaportT.In situ microwave reflection coefficient measurements for the rough exterior wall surfacesProceedings of the 43rd IEEE Vehicular Technology ConferenceMay 199377802-s2.0-0027206531OgunsolaA.ReggianiU.SandroliniL.Modelling shielding properties of concreteProceedings of the 17th International Zurich Symposium on Electromagnetic CompatibilityMarch 200634372-s2.0-33751048117OgunsolaA.ReggianiU.SandroliniL.Shielding properties of conductive concrete against transient electromagnetic disturbancesProceedings of the IEEE International Conference on Microwaves, Communications, Antennas and Electronics Systems (COMCAS '09)November 2009Tel-Aviv, Israel2-s2.0-7794984947810.1109/COMCAS.2009.5385975Raymond LuebbersJ.Tom UnoKen KumagaiComments on ‘pulse propagation in a linear, causally dispersive medium’1993814631639YoungJ. L.NelsonR. O.A summary and systematic analysis of FDTD algorithms for linearly dispersive media200143161772-s2.0-003524570810.1109/74.920019SamnS. W.Modeling dispersive dielectric media in FDTD: a systematic approach20055310336733732-s2.0-2764458705710.1109/TAP.2005.856318AlbaneseR.PennJ.MedinaR.Short-rise-time microwave pulse propagation through dispersive biological media19896914411446JordanP. M.PuriA.Digital signal propagation in dispersive media1999853127312822-s2.0-0343455253BalanisC. A.1989John Willey and SonsDemarestK. R.19971st Prentice HallAlejosA. V.DawoodM.MohammedH. U. R.Empirical pseudo-optimal waveform design for dispersive propagation through loamy soil2012599539572-s2.0-8485873248310.1109/LGRS.2012.2186431