arxiv: v1 [cs.ni] 15 Feb 2016

Transcription

1 1 Energy-Effcent Tme Synchronzaton Based on Asynchronous Source Cloc Frequency Recovery and Reverse Two-Way Message Exchanges n Wreless Sensor Networs Kyeong Soo Km, Member, IEEE, Sanghyu Lee, and Eng Gee Lm, Member, IEEE arxv: v1 [cs.ni] 15 Feb 016 Abstract We consder energy-effcent tme synchronzaton n a wreless sensor networ where a head node (.e., a gateway between wred and wreless networs and a center of data fuson s equpped wth a powerful processor and suppled power from outlet, and sensor nodes (.e., nodes measurng data and connected only through wreless channels are lmted n processng and battery-powered. It s ths asymmetry that our study focuses on; unle most exstng schemes to save the power of all networ nodes, we concentrate on battery-powered sensor nodes n mnmzng energy consumpton for tme synchronzaton. We present a tme synchronzaton scheme based on asynchronous source cloc frequency recovery and reverse two-way message exchanges combned wth measuremenata report messages, where we mnmze the number of message transmssons from sensor nodes, and carry out the performance analyss of the estmaton of both measurement tme and cloc frequency wth lower bounds for the latter. Smulaton results verfy that the proposed scheme outperforms the schemes based on conventonal two-way message exchanges wth and wthout cloc frequency recovery n terms of the accuracy of measurement tme estmaton and the number of message transmssons and receptons at sensor nodes as an ndrect measure of energy effcency. Index Terms Tme synchronzaton, energy effcency, source cloc frequency recovery, two-way message exchanges, wreless sensor networs. I. INTRODUCTION A common tme frame among networ nodes n a wreless sensor networ (WSN s crtcal for ther carryng out mportant operatons le fusng data from dfferent sensor nodes, tme-based channel sharng and meda access control (MAC protocols, and coordnated sleep wae-up node schedulng mechansms [1]. In a typcal WSN, a head/master node (.e., a base staton that serves as a gateway between wred and wreless networs and a center for fuson of data from dstrbuted sensors s equpped wth a powerful processor and suppled power from outlet, whle sensor/slave nodes (.e., nodes measurng data wth sensors and connected to a WSN only through wreless channels are lmted n processng and battery-powered. It s ths asymmetry that our study focuses on; unle most exstng schemes tryng to save the power of all networ nodes (e.g., [] and [3], we concentrate on battery-powered sensor nodes, whch are Ths wor was supported n part by X an Jaotong-Lverpool Unversty Centre for Smart Grd and Informaton Convergence (CeSGIC and Research Development Fund (under Grant RDF An earler verson of ths paper was presented n part at ICISCA 015, Kuala Lumpur, Malaysa, June 015, and ISAE 015, Busan, Korea, October 015. K. S. Km, S. Lee, and E. G. Lm are wth the Department of Electrcal and Electronc Engneerng, X an Jaotong-Lverpool Unversty, Suzhou 1513, Jangsu Provnce, P. R. Chna (e-mal: {Kyeongsoo.Km, Sanghyu.Lee, Enggee.Lm}@xjtlu.edu.cn. many n number, n mnmzng energy consumpton for tme synchronzaton. In ths paper we mean by tme synchronzaton a process of establshng a common tme frame wth whch the nodes n a networ can operate one another, whether ther clocs are synchronzed or not. By cloc synchronzaton, on the other hand, we mean a process of synchronzng each node cloc to that of a common reference node, typcally through a logcal cloc that s a functon of a physcal cloc, whch s one way of achevng tme synchronzaton. Note that wthout synchronzng node clocs, we can stll provde a common tme frame for the operaton of networ nodes. For nstance, n mult-hop extenson of the reference broadcast synchronzaton (RBS algorthm [4], there are multple broadcasomans ndependently mantanng ther own clocs. In ths case, the common tme frame among multple broadcasomans s provded through tme converson at gateway nodes belongng to neghborng domans. In the lterature, tme synchronzaton s formulated as the problem of estmatng cloc parameters, often ncludng node dstances, n a parwse (dstrbuted or a networ-wde (global manner [5] [8]. Few wors, however, have focused on the mode of operaton (.e., tme vs. cloc synchronzaton dscussed above and the way of mplementaton related wth energy effcency. To acheve better performance, several jont estmaton algorthms for both cloc parameters and par-wse dstances have been proposed (e.g., [7] and [6]. Most of them are, however, based on centralzed and global offlne algorthms, and the estmaton of parameters and the way of mplementng cloc/tme synchronzaton together wth the delvery of the results of estmaton are not explctly dscussed. Whle jont estmaton algorthms usually assume that all sensor nodes are avalable n the begnnng, recursve and onlne operaton of synchronzaton schemes are crtcal for WSN applcatons because there are nodes who jon later or resume ther operatons n the mddle. In ths paper we present an energy-effcent tme synchronzaton scheme based on asynchronous source cloc frequency recovery (SCFR [9] and reverse two-way message exchanges combned wth measurng data report messages, where we mnmze the number of message transmssons and receptons at sensor nodes, especally the number of message transmssons notng that the energy for message transmsson s typcally hgher than that for message recepton [10]. In the proposed tme synchronzaton scheme, only the frequency of a sensor node cloc s synchronzed to that of the reference cloc at the head node, but not ts cloc offset; the proposed scheme s based on the dea of the separaton of the cloc

2 frequency 1 estmaton/compensaton at sensor nodes and the cloc offset and delay estmaton at the head node. For the cloc frequency recovery, each sensor node passvely lstens to any messages wth tmestamps ether broadcasted (e.g., beacons or uncasted (e.g., control messages to a specfc node from the head node and carres out asynchronous source cloc frequency recovery descrbed n [9], whch s bascally one-way cloc frequency estmaton. For the cloc offset and delay estmaton, a smple two-way message exchange procedure [11], [1] s used but n a reverse drecton where the head node ntates the procedure and eeps trac of the offsets between ts reference cloc and sensor node clocs; also, nstead of dedcated, perodcal synchronzaton message exchanges, we embed the synchronzaton Response messages of the two-way message exchange procedure n the measuremenata report messages from sensor nodes n order to mnmze the number of message transmssons. In ths way we can move most of tme synchronzaton operatons to the head node and reduce the complexty and thereby power consumpton of sensor nodes for tme synchronzaton. To carry out usual WSN operatons wth the proposed tme synchronzaton scheme, the head node translates tmestamp values based on cloc offset nformaton before transmttng control messages to sensor nodes. The major contrbutons of ths wor are three-fold: Frst, t provdes a new energy-effcent tme synchronzaton scheme for asymmetrc WSNs to mnmze the number of two-way message exchanges by cleverly combnng one-way cloc sew estmaton/compensaton and reverse two-way message exchanges. Note that the two-way message exchanges cannot be avoded n estmatng both cloc offset and cloc sew because t s mpossble to separate the effects of cloc offset and propagaton delay wth one-way message exchanges alone [1], [13]. The proposed scheme does not only mnmze the number of two-way message exchanges but also enable the use of smple, low-complexty one-way cloc sew estmaton/compensaton algorthms to avod the hgh complexty of advanced estmaton algorthms tang nto nusance parameters le mean and varance of the propagaton delay. Secondly, we formally descrbe the operatons of and analyze the tme synchronzaton performance of the proposed scheme: By separately modelng tmes of a hardware cloc and a logcal cloc at a networ node, we can capture the dynamc nature of the proposed tme synchronzaton scheme n the analyss, whch estmates cloc parameters of a counterpart and updates ts own logcal cloc recursvely. Thrdly, we also analyze the performance of measurement tme estmaton based on both conventonal and reverse two-way message exchanges. As for cloc sew estmaton/compensaton, we carry out a comparatve analyss of the performance of one-way and twoway maxmum lelhood ( and -le estmators and derve Cramér-Rao lower bounds and lower bounds for and -le estmators, respectvely. The rest of the paper s organzed as follows: In Secton II, we descrbe the proposed tme synchronzaton scheme wth 1 We use the terms cloc frequency and cloc sew nterchangeably n ths paper. (a (b Fg. 1. Comparson of two-way message exchanges: (a Conventonal twoway message exchanges as n the tme-sync protocol for sensor networs (TPSN [1]; (b reverse two-way message exchanges of the proposed scheme shown wth optonal bundlng of measurements n a Report/Response message. hardware and logcal cloc models; we also analyze the effect of cloc sew on measurement tme estmaton n both conventonal and reverse two-way message exchanges, present jont and -le one-way cloc sew estmators wth ther performance bounds, and dscuss an extenson to multhop tme synchronzaton through gateway nodes. Secton III presents the results of smulatons for a comparatve analyss of the performance of one-way and two-way cloc sew estmators for asynchronous SCFR and the nvestgaton of the tme synchronzaton performance and the energy effcency of the proposed scheme compared to the schemes based on conventonal two-way message exchanges wth and wthout SCFR. Secton IV revews the related wor n comparson to our wor before concludng ths paper n Secton V. II. ENERGY-EFFICIENT TIME SYNCHRONIZATION FOR ASYMMETRIC WSNS The major dea of the proposed scheme s to allow ndependent, unsynchronzed slave clocs at sensor nodes but runnng at the same frequency as the reference cloc at a head node through the asynchronous SCFR descrbed n [9], whch needs only the recepton of messages wth tmestamps at sensor nodes. The cloc offset, on the other hand, s estmated at the head node based on the reverse two-way message exchanges. Fg. 1 llustrates ths dea n comparson to conventonal twoway message exchanges. Frst, the proposed scheme shown n Fg. 1 (b does not have perodc, dedcated two-way message exchanges wth synchronzaton messages le Request and Response shown n Fg. 1 (a; nstead, the Request and Response messages

3 3 Fg.. Payload structure of a "Report/Response" message of the proposed scheme wth optonal bundlng of measuremenata and tmestamps. are embedded n the most recent tmestamped message ether broadcasted or uncasted to a specfc node from the head node and a measuremenata report message from a sensor node, respectvely. When there are no strct tmng requrements for the processng of measuremenata, the measuremenata and ther correspondng tmestamps can be optonally bundled together n a Report/Response message, whose payload structure s shown n Fg., n order to further reduce the number of message transmssons; n ths case the cloc offset estmated based on the tmestamp of the last measuremenata #N s collectvely appled to all bundled measuremenata n estmatng ther occurrences wth respect to the reference cloc at the head node. Secondly, the drecton of two-way message exchanges of the proposed scheme s reversed,.e., t s the master (.e., the head node that sends the Request messages, not the slave (.e., the sensor node, unle the conventonal twoway messages exchanges; as a result, the master nows the current status of the slave cloc, but the slave does not. So the nformaton of slave clocs (.e., cloc offsets wth respect to the reference cloc s centrally managed at the head node. For operatons le coordnated sleep wae-up node schedulng, before sendng a control message, the head node frst adjusts the tme for future operaton based on the cloc offset of the recpent sensor node. In ths way, even though sensor nodes n the networ have clocs wth dfferent cloc offsets, ther operatons can be coordnated based on the common reference cloc at the head node. A. Hardware and Logcal Cloc Models We consder an asymmetrc WSN where one head (master node and N sensor (slave nodes all equpped wth ndependent hardware clocs based on quartz crystal oscllators. For smplcty, we tae tme t of the head node cloc as a global reference and descrbe tmes of slave hardware clocs as functons of t. We use the frst-order affne cloc model [6] to descrbe tme T of the hardware cloc at the th sensor node as follows: For [0, 1,..., N 1], T (t = (1 + ɛ t + θ, (1 where (1+ɛ R + and θ R are cloc frequency rato and cloc offset, respectvely. Note that ɛ s called a cloc sew n the lterature, whch s defned as a normalzed cloc frequency dfference between hardware clocs, and ts typcal value for clocs based on quartz crystal oscllators s of the order of tens of ppm (.e., ɛ 1. WSN operatons at a sensor node s based on a logcal cloc, whose tme s agan a functon of the tme of ts physcal cloc and taes nto account the adjustments by an adopted tme synchronzaton scheme (e.g., offset adjustment by twoway message exchanges and frequency adjustment by SCFR. (a (b Fg. 3. Tmng dagrams for the analyss of the measurement tme estmaton error: (a Conventonal two-way message exchanges; (b reverse two-way message exchanges of the proposed scheme. Specfcally, tme T of the logcal cloc at the th sensor node can be modeled as a pecewse lnear functon as follows : For t <t t +1 (=0, 1,..., T (T (t = T T (t + T (t T (t 1 + ˆɛ, ˆθ,, ( where t s the reference tme when a th synchronzaton occurs, and ˆɛ, and ˆθ, are the estmated cloc sew and offset from the th synchronzaton. If the synchronzaton s only for frequency, we set ˆθ, to 0 n (; f the synchronzaton s only for offset, on the other hand, we set ˆɛ, to 0 n (. It should be noted that T n ( s the functon of T (t, but not t tself, because the reference tme t s not nown at a sensor node and the only avalable tme s from a local hardware cloc (.e., T (t. B. Effect of Cloc Sew on Measurement Tme Estmaton Here we compare the effect of cloc sew on the measurement tme estmaton n both conventonal and reverse twoway message exchange procedures by an approxmate analyss of the best-case performance under a determnstc delay (.e., no random component and no bundled measurements. Fg. 3 shows tmng dagrams for the analyss of the measurement tme estmaton error, where we assume the same amount of tme dfference T m between the measurement occurrence and the recepton of the last tme synchronzaton message.e., the Response message for the conventonal two-way message exchanges and a beacon ( Request message for the reverse two-way message exchanges. Because MAC-layer tmestampng of messages can remove all the sources of uncertantes except the propagaton delay [] Ths model fts for WSN operatons requrng tmes of dscrete events only. More complcated (often nonlnear models, however, s better suted for applcatons le playbac of multmeda streamng where an analog or dgtal phase-loced loop (PLL s used to generate a cloc sgnal [9].

4 4 and the varaton n propagaton delays s neglgble n snglehop transmssons, the cloc offset estmaton error ˆθ n the two-way message exchanges wthout cloc sew compensaton (.e., based on the hardware cloc model n (1 s gven by [1] ˆθ = (T 4 T 1ɛ = d ɛ, (3 where d denotes one-way propagaton delay. The measurement tme estmaton error for the conventonal two-way message exchanges shown n Fg. 3 (a, therefore, can be expressed as ˆt Conv. m = dɛ + T m ɛ = (d + T m ɛ. (4 If T m d, the measurement tme error can be approxmated as ˆt Conv. m T m ɛ. (5 When the cloc sew s compensated, (5 becomes ˆt Conv. m T m ˆɛ, (6 where ˆɛ s the cloc sew estmaton error. For the reverse two-way message exchanges shown n Fg. 3 (b, because the Request message of the reverse twoway message exchanges s embedded n the measuremenata Report message from a sensor node as shown n Fg., there s no tme dfference between the measurement tme and the end of two-way message exchange (.e., T 4. In ths case the two-way message exchange procedure s the only source of error,.e., ˆt Rev. m = (T 4 T 1 ˆɛ = (d + T m ˆɛ. (7 Agan, f T m d, the measurement tme estmaton error can be approxmated as ˆt Rev. m T m ˆɛ. (8 From (6 and (8, we can see that, when beng solated from the effect of random component of delay, the proposed scheme can reduce the effect of cloc sew on the measurement tme estmaton error by a factor of two for T m d. Because the tme dfference T m between the recepton of the last beacon (or any tmestamped message from the head node and the measurement occurrence at the sensor node can be a larger value n practce, however, t s stll mportant to compensate the cloc sew at sensor nodes n the proposed scheme. C. Asynchronous SCFR at Sensor Nodes: One-Way Cloc Sew Estmaton One of the essental components of the proposed tme synchronzaton scheme s the recovery of the reference cloc frequency at sensor nodes based on one-way message dssemnaton from the head node. Once we estmate the cloc sew ɛ n (1, the reference cloc frequency can be recovered through the cloc sew compensaton of the logcal cloc model n (. In case of two-way message exchanges, both jont estmaton of cloc offset and sew wth a nown fxed porton of delay and separate -le estmaton of cloc sew are studed n [14]. As for one-way message dssemnaton, here we derve both jont and separate -le estmators n a smlar manner but based on the problem formulaton n [9]. Let ( (=0, 1,... be the reference tme for the th message s departure from the head node; note that ( also denotes the value of the tmestamp carred by the th message. From the hardware cloc model n (1, we can obtan t a, (, the arrval tme of the th message wth respect to the th sensor node s hardware cloc, as follows: t a, ( = T ( ( + d( = (1 + ɛ ( + θ + d(, (9 where d( denotes a one-way paceelay from the head node to the th sensor node n terms of the th sensor node s hardware cloc. 3 For the observaton of tmestamps from (9, we can obtan the jont estmators (Es of cloc offset and sew as stated n Proposton 1. Proposton 1: For a whte Gaussan delay d( wth nown mean d and varance σ, the jont estmators (Es for cloc offset ˆθ ( and sew ˆR ( n (9 are gven by ˆθ ( = t d t a, t a, t d d, (10 ˆR ( = t a, t a, t d, (11 where the notatons x and xy denote the average values of x(j (.e., x(j/ and x(jy(j (.e., x(jy(j/, respectvely. Also ˆθ ( and ˆR ( are effcent estmators [16, p. 34] whch are unbased and attan the Cramér-Rao lower bounds (CRLBs gven by Var (ˆθ ( Var ˆR ( σ t d {t d }, (1 σ {t d }. (13 Proof: The proof of Proposton 1 s presented n Appendx A. Note that, even though the mean and the varance of whte Gaussan delay are assumed to be nown for the dervaton of ˆR the jont Es, the resultng cloc sew estmator ( the only one needed n the proposed scheme does noepend on t. Compared to the jont Es for two-way message exchanges derved n [14], (10 and (11 tae smpler expressons, but they are stll complcated and not sutable for recursve mplementaton. Because we are only nterested n the estmaton of the cloc sew equvalently the rato of cloc frequences at sensor nodes, we can formulate the estmaton problem where the cloc sew s the only parameter to estmate. In [9], the problem of asynchronous SCFR s formulated as a lnear regresson through the orgn (RTO 3 Ths defnton of the one-way delay (.e., n terms of the sensor nodes hardware cloc s dfferent from that n [14] and [15], where the delay s defned n terms of the head node cloc, and maes smpler the dervaton of the one-way estmators n ths paper.

5 5 model by subtractng both sdes of (9 wth ther ntal values: For = 1,,..., t a, ( = (1 + ɛ ( + d(, (14 where t a, ( t a, ( t a, (0, ( ( (0, and d( d( d(0. Note that d( now represents a nose process wth a zero mean for a statonary delay model. Note that t a, ( s n (14 are not ndependent one another due to d(. The dervaton of an E based on all observatons of t a, ( s, therefore, s not straghtforward. In [9], two practcal estmators are proposed n ths regard,.e., the recursve least squares (RLS and the cumulatve rato (CR, whch are not based on any assumpton on the delay dstrbuton. Of the two estmators, the CR estmator s best suted for battery-powered sensor nodes n the proposed tme synchronzaton scheme due to ts lower complexty. If we defne R as the rato of the th sensor node hardware cloc frequency to that of the reference cloc (.e., 1+ɛ, the CR estmator ˆR CR ( s gven by From (14, we can see that follows: ˆR CR ( = t a, ( (. (15 ˆRCR ( can be rewrtten as ˆR CR ( = R + d( t s (, (16 where ts nose component becomes zero as tme goes to nfnty rrespectve of ts statstcal characterstcs. Unle the RLS estmator, there are no desgn parameter values or ntal values to set for the CR estmator. Note that for a Gaussan delay model, the CR estmator becomes an unbased estmator achevng ts lower bound as stated n Proposton. Proposton : For a whte Gaussan delay d( wth zero mean and varance σ, the CR estmator n (16 becomes an unbased estmator whch attans the lower bound 4 gven by Var ˆRCR ( σ (. Proof: The proof of Proposton s presented n Appendx B. D. Extenson to Mult-Hop Tme Synchronzaton Fg. 4 shows how the proposed scheme can be extended to a herarchcal structure for networ-wde, mult-hop synchronzaton through smple pacet-relayng or more advanced tme-translatng gateway nodes, the latter of whch act as both head nodes (for the nodes n the lower herarchy and normal sensor nodes (for the gateway node n the hgher herarchy: For example, consder the message transmsson from the sensor node S to the head node through the two gateway nodes G1 and G as shown n Fg. 4. In case of pacet relay, G smply passes the receved tmestamped message from S to G1 and G1 to the head node wthout any change of the 4 Ths lower bound s n fact the CRLB for the case when the observatons are lmted to the tmestamps from the frst and last messages. Fg. 4. Extenson of the proposed tme synchronzaton scheme to a herarchcal structure for networ-wde, mult-hop synchronzaton through pacet-relayng or tme-translatng gateway nodes. value of tmestamp n the mddle. In case of tme translaton, because G acts as a head node for S, t frst translates the value of tmestamp based on the nformaton on the cloc offset of S. Then G relays the message from S to G1 wth translated tmestamp value (wth respect to ts own cloc. From G1 s pont of vew, G s just one of sensor nodes t manages. Agan, based on the nformaton on the cloc offset of G, G1 translates the value of cloc offset wth respect to ts own cloc and relays the message to the head node. Fnally, the head node receves the message from S, whch s just relayed by G1, and translates the tmestamp value based on the nformaton on the cloc offset of G1 t manages. In ths way the head node can obtan the measuremenata and ts occurrence n tme reported by S wth respect to ts own reference cloc. It s clear that, compared to tme-translatng gateway nodes, pacet-relayng gateway nodes could be much smpler because there s no addtonal functonalty (.e., estmaton and management of the offset nformaton and translaton of tmestamp values except pacet relayng. A downsde, however, s that they ntroduce larger paceelays resultng from queueng and MAC operatons, whch could deterorate the performance of tme synchronzaton. III. SIMULATION RESULTS We carry out a seres of smulaton experments to nvestgate the performance of the proposed tme synchronzaton scheme n terms of the accuracy of measurement tme estmaton and the number of message transmssons and receptons at a sensor node n comparson wth those based on the conventonal two-way message exchanges. We also analyze the performance of one-way cloc sew estmators ndependently of the proposed tme synchronzaton scheme to choose the best one for battery-powered sensor nodes. For the smulatons, we consder a smple WSN wth one head node and one sensor node that are deployed 100 m from each other because the tme synchronzaton of a sensor node n the proposed scheme can be carred out ndependently of that of other sensor nodes. As n [7], we set the cloc

6 B. Performance of Measurement Tme Estmaton and Energy Effcency For the analyss of the performance of the proposed tme synchronzaton scheme, we run smulatons for three dfferent values of synchronzaton nterval (SI5.e., 100 s, 1 s and 10 ms to nvestgate the effect of the tme dfference Tm between the last tme synchronzaton message from the head node and the measurement occurrence at the sensor node. The propagaton delay s modeled as an..d. Gaussan process wth a standard devaton of 1 ns. Durng the observaton nterval of 1 h, total 100 measurements are made where ther correspondng data arrvals are modeled as a Posson process. For SCFR, we use the CR estmator for the proposed scheme based on the observaton n Secton III-A and the GLE for the two-way scheme for comparson. The results of SCFR at the sensor node and measurement tme estmaton at the head node (.e., the estmaton of the 5 Ths s the nterval of beacons for the proposed scheme and the nterval of two-way message exchanges for others. Measurement Tme Error [s] 4E-11 Proposed (w/ CR Two-Way (w/ GLE E-11 -E-11-4E-11 Proposed (w/ CR Two-Way (w/ GLE 4E-03 E-03 Two-Way -E-03-4E-03-6E-03-8E-03-1E Tme [s] Measurement Tme Error [s] Frequency Dfference [ppm] (a 4E-11 Proposed (w/ CR Two-Way (w/ GLE E-11 -E-11-4E-11 Proposed (w/ CR Two-Way (w/ GLE 4E-05 E-05 Two-Way -E-05-4E-05-6E-05-8E-05-1E Tme [s] (b Frequency Dfference [ppm] A. Performance of One-Way Cloc Sew Estmaton Frst, we analyze the performance of one-way cloc sew estmators dscussed n Secton II-C, together wth the twoway -le estmator for Gaussan delay model (GLE proposed n [14] for comparson. Fg. 5 shows the mean square error (MSE of both one-way and two-way cloc sew estmators wth Gaussan delays wth standard devaton of 1 ns and 1 µs, whch are calculated over 10,000 smulaton runs. Fg. 6 also shows the results wth frst-order autoregressve (AR(1 random delays wth correlaton coeffcent (ρ of 0.6 and standard devaton of 1 µs and 1 ms; the results n ths case demonstrate the performance of the estmators under correlated delays, whch reflect the case of mult-hop extenson through pacet-relayng gateway nodes as dscussed n Secton II-D. The results for both Gaussan and AR(1 delays show that the jont E provdes the best performance as expected because t can use all observed tmestamp values. The performance gap between the jont E and other one-way estmators (.e., CR and RLS, however, becomes narrower when the delay does not follow the Gaussan dstrbuton and has correlaton. The estmaton performance of the two-way scheme (.e., GLE, on the other hand, s the worst gven the number of messages exchanged, whch ndcates that the two-way message exchange procedure s not energy-effcent n estmatng the cloc sew only. As dscussed n Secton II-C, due to ts lower complexty and robustness to delay characterstcs, we use the CR estmator at battery-powered sensor nodes n the nvestgaton of the performance of the proposed tme synchronzaton scheme n the followng sectons. Measurement Tme Error [s] frequency rato R and the cloc offset θ to ppm and 1 s, respectvely. We model propagaton delay, whch taes nto account tmestamp generaton and recepton nose as well, wth an ndependent and dentcally dstrbuted (..d. Gaussan process unless specfed otherwse. Frequency Dfference [ppm] 6 4E-11 Proposed (w/ CR Two-Way (w/ GLE E-11 -E-11-4E-11 Proposed (w/ CR Two-Way (w/ GLE 4E-07 E-07 Two-Way -E-07-4E-07-6E-07-8E-07-1E Tme [s] (c Fg. 7. Performance of SCFR at a sensor node and measurement tme estmaton at a head node for SI of (a 100 s, (b 1 s, and (c 10 ms. measurement tme at the sensor node wth respect to the reference cloc at the head node are shown n Fg. 7, and ther MSEs calculated over samples from one smulaton run wth the number of message transmssons (NTX and receptons (NRX at the sensor node are summarzed n Table I. As shown by the approxmate analyss n Secton II-B, the results show that, wthout SCFR (.e., Two-Way, the measurement tme estmaton errors hghly depend on the synchronzaton nterval. The use of SCFR (.e., Proposed (w/ CR and Two-Way (w/ GLE, on the other hand, greatly reduces ths dependency, and the resultng measurement tme estmaton errors (.e., the square root of MSE of measurement

7 RLS CR Jont E GLE (Two-Way LB for CR CRLB for Jont E LB for GLE RLS CR Jont E GLE (Two-Way LB for CR CRLB for Jont E LB for GLE MSE MSE Number of Messages (a Number of Messages Fg. 5. MSE of the estmated cloc sews from one-way and two-way cloc sew estmators for Gaussan random delays wth standard devaton (σ of (a 1 ns and (b 1 µs. (b RLS CR Jont E GLE (Two-Way RLS CR Jont E GLE (Two-Way MSE 10-1 MSE Number of Messages (a Number of Messages Fg. 6. MSE of the estmated cloc sews from one-way and two-way cloc sew estmators for AR(1 random delays wth correlaton coeffcent (ρ of 0.6 and standard devaton of (a 1 µs and (b 1 ms. (b TABLE I RESULTS OF TIME SYNCHRONIZATION SIMULATIONS FOR DIFFERENT SIS Synchronzaton Scheme Sew Estmaton MSE 1 Measurement Tme Estmaton MSE 1 N TX N RX SI = 100 s E E Proposed SI = 1 s E E SI = 10 ms E E Two-Way wth GLE Two-Way SI = 100 s 1.901E E SI = 1 s E E SI = 10 ms 9.099E E SI = 100 s E SI = 1 s N/A E SI = 10 ms E For the samples obtaned from 360 s (.e., one tenth of the observaton perod to avod the effect of a transent perod. tme estmaton are of the order of the nose standard devaton (.e., 1 ns for all values of the synchronzaton nterval, whch means that the synchronzaton nterval can be ncreased wthout much affectng the tme synchronzaton performance to save energy. Table I also shows that the number of message transmssons for the proposed scheme s the same as the number of measurements rrespectve of the synchronzaton nterval, whle the number of message receptons ncreases as the synchronzaton nterval decreases for all the schemes consdered; n case of the synchronzaton nterval of 1 s, the number of message transmssons s 100 for the proposed scheme and 3,700 (100 for the measuremenata and 3,600 for the synchronzaton messages for the schemes based on the conventonal two-way messages exchanges. C. Effect of Bundlng of Measurement Data We also nvestgate the effect of the bundlng of measurement data descrbed n Secton II wth fxed synchronzaton nterval of 1 s. Fg. 8 shows the measurement tme estmaton errors for dfferent numbers of bundled measurements (N BM n one Report/Response message, and agan ther MSEs calculated over samples from one smulaton run wth the

8 8 TABLE II RESULTS OF TIME SYNCHRONIZATION SIMULATIONS FOR DIFFERENT NUMBERS OF BUNDLED MEASUREMENTS (N BM WITH SI OF 1 S. Synchronzaton Scheme Measurement Tme Estmaton MSE 1 N TX N RX N BM = E Proposed N BM = E N BM = E N BM = E For the samples obtaned from 360 s (.e., one tenth of the observaton perod to avod the effect of a transent perod. Measurement Tme Error [s].0e E E E-09 N BM =1 N BM = N BM =5 N BM = Tme [s] Fg. 8. Measurement tme estmaton errors for dfferent numbers of bundled measurements (N BM wth SI of 1 s. number of message transmssons and receptons at the sensor node are summarzed n Table II. It s evdent from the fgure and the table that the measurement tme estmaton errors are not much affected by the number of bundled measurements after the ntal transent perod (.e., after around 500 s. In fact, the table shows that the MSE of measurement tme estmaton slghtly decreases as the number of bundled measurements ncreases; ths s because all the measuremenata reported n a Report/Response message share the same estmated value of the cloc offset as the last measurement n the bundle, whch n fact reduces the fluctuatons n the measurement tme estmaton errors as observed n Fg. 8. The bundlng of measuremenata further reduces the number of message transmssons: For nstance, f we bundle 10 measuremenata n a Report/Response message (.e., N BM =10, we can reduce the number of message transmssons from 100 to 10, whch s less than 0.3 percentage of the number of message transmssons for the schemes based on the conventonal two-way message exchanges (.e., 3,700. Note that we need to be careful n nterpretng the number of message transmssons as an ndrect measure of energy consumpton n ths case. Because the bundlng ncreases the length of message payload, the total number of bts transmtted for the case wth the bundlng s not the same as that for the case wthout the bundlng gven the number of messages. Stll, however, the measurement bundlng s an effectve opton to save the energy consumpton because we can save the energy for the transmsson of the overhead of a message ncludng a preamble, a frame header, and a frame chec sequence by reducng the number of message transmssons. Also note that the measurement bundlng can be used only when there are no strct tmng requrements for the processng of measurement data due to the ncreased delay resultng from the bundlng process. IV. COMPARISON TO RELATED WORK Unsynchronzed Clocs: The proposed tme synchronzaton scheme s smlar to the post-facto synchronzaton scheme descrbed n [17] n that sensor nodes clocs are unsynchronzed and that the tme of the occurrence of an event s frst recorded wth respect to the local cloc. The dfference s that tme translaton occurs at sensor nodes just after recevng synchronzaton messages from a thrd-party node (e.g., a head node n the post-facto synchronzaton, whle, n the proposed scheme, the tme of the occurrence s transmtted to the head node wthout watng for a synchronzaton message and the tme translaton s done at the head node. As for the recever cloc sew compensaton, the use of networ tme protocol (NTP [11] was suggested n [17]; because local-cloc resoluton and sew are mnmzed by the controlfeedbac desgn (e.g., based on phase-loced loop (PLL, however, the NTP s not proper for low-power/complexty sensor nodes. In ths regard the cloc sew s handled by a oneway, low-complexty asynchronous SCFR scheme descrbed n [9] for the proposed scheme. Reverse Two-Way Message Exchanges: The ReversePTP proposed n [18] for tme dstrbuton n software defned networs (SDNs s also smlar to our wor because the relatonshp and the drecton of cloc dstrbuton between one master (.e., the head node and many slaves (.e., sensor nodes are reversed n both schemes. Ther contexts (.e., SDNs vs WSNs, however, are qute dfferent from each other. As a consequence, the energy effcency s not the man focus of the ReversePTP. For example, the slaves n the ReversePTP dstrbute ther tmes to the master by perodcally sendng Sync messages. In the proposed scheme, t s the head node thastrbutes ts tme (.e., cloc frequency; f there s no data measurement, a sensor node does not transmt any message bac to the head node for tme synchronzaton. Passve Lstenng: The dea of passve lstenng to tmestamped messages from other nodes s also dscussed n [19] and [7]. The recever-only synchronzaton (ROS scheme proposed n [19] enables a subset of sensor nodes to acheve synchronzaton by overhearng the conventonal two-way synchronzaton message exchanges of a par of sensor nodes,

9 9 called super nodes; because of the addtonal nformaton from the overheard two-way messages exchanges, the nodes can estmate both cloc offset and cloc sew wthout any message transmsson. The applcaton of ROS scheme, however, s lmted to the sensor nodes located n the Regon of Parwse Synchronzaton where the messages from both super nodes can reach. Also, two super nodes performng two-way message exchanges are needed for ROS scheme. The jont estmaton algorthm based on prewhtenng of observaton models and least squares (LS n [7] can be mplemented ether n a centralzed or dstrbuted ways, but the dstrbuted approach ncurs addtonal broadcast messages to dstrbute tmestamps. In the proposed scheme, because the cloc sew and offset estmaton procedures are separated from each other, a dstrbuted recursve algorthm can be appled for the cloc sew estmaton based on passve lstenng wthout ncurrng addtonal message broadcastng. V. CONCLUSIONS In ths paper we have proposed an energy-effcent tme synchronzaton scheme for asymmetrc WSNs, whch s based on the asynchronous SCFR for one-way cloc sew estmaton/compensaton at sensor nodes and reverse two-way message exchanges for cloc offset estmaton/translaton at the head node n order to mnmze the number of message transmssons at sensor nodes. Tang notce of the asymmetry between the head node and sensor nodes n terms of processng power and suppled energy, we move most of tme synchronzaton operatons to the head node n the proposed scheme and thereby reduce the complexty and power consumpton of sensor nodes for tme synchronzaton. As for the one-way cloc sew estmaton/compensaton,.e., the only major operaton to be done at sensor nodes, because the low complexty of the algorthm s crtcal for power effcency, we use the CR estmator proposed n [9] for asynchronous SCFR, whch s an unbased estmator attanng the lower bound for a Gaussan delay as stated n Proposton. The comparatve analyss presented n Secton III-A shows that the CR estmator can provde relatvely good performance wth the lowest complexty among the estmators consdered rrespectve of the delay dstrbuton. The smulaton results n Secton III-B n fact verfes that the proposed tme synchronzaton scheme wth the CR as the cloc sew estmator provdes the best performance n terms of the accuracy of measurement tme estmaton and, as an ndrect measure of energy consumpton, the number of message transmssons and receptons at a sensor node. In addton, f there are no strct tmng requrements for the processng of measuremenata, we can further reduce the number of message transmssons by bundlng measuremenata n a Report/Response message wthout sgnfcantly affectng the tme synchronzaton performance, whch s demonstrated by the smulaton results n Secton III-C. Note that n ths paper we have focused on the essental aspects of the proposed tme synchronzaton scheme,.e., the one-way cloc sew estmaton/compensaton and the reverse two-way message exchanges combned wth measuremenata report messages, whch leaves room for further extensons of the proposed scheme; examples nclude the operatons wth multple head nodes for redundancy and possbly localzaton as well [7] and the applcaton of adaptve synchronzaton nterval []. Of possble extensons, we have already shown how the proposed scheme can be extended to a herarchcal structure for networ-wde, mult-hop synchronzaton through ether smple pacet-relayng or tme-translatng gateway nodes. The tradeoff between the complexty of gateway nodes and the performance of tme synchronzaton for the two approaches s an mportant topc for further research and wll be addressed n a followup wor. APPENDIX A PROOF OF PROPOSITION 1 We frst rephrase (9 as follows: t a, ( = R ( + θ + d + n(, (17 whee d s a fxed porton of delay, whch s assumed to be nown, and n( s a random porton of delay followng a zeromean whte Gaussan dstrbuton wth varance σ. From (17, we can construct a lelhood functon and a log-lelhood functon for the collectve observaton t a, {t a, (j} and the parameters θ and R as follows: L(θ, R ; t a, = ( 1 exp (t a,(j (R (j + θ + d πσ σ, and j=1 ln L(θ, R ; t a, = ln ( πσ 1 σ (18 {t a, (j (R (j + θ + d}. (19 Then the frservatves of the log lelhood functon wth respect to θ and R are gven by ln L(θ, R ; t a, = θ 1 σ {t a, (j (R (j + θ + d}, (0 ln L(θ, R ; t a, = R 1 σ (j {t a, (j (R (j + θ + d}. (1 From (0 and (1, therefore, we obtan the jont one-way E of cloc offset ( and sew ( as follows: ˆθ ˆR ( = ˆθ ˆR ˆR ( = t a, ( d, ( { (ˆθ } t a, (j ( + d t d (j. (3

10 10 After some manpulatons, we can obtan closed-form expressons for ˆθ ( and ˆR ( gven n (10 and (11, respectvely. To obtan the CRLBs, we need the second dervatves of the log lelhood functon wth respect to θ and R, whch are gven by ln L(θ, R ; t a, = σ, (4 ln L(θ, R ; t a, R ln L(θ, R ; t a, θ R θ = d(j σ = σ t d, (5 = d(j σ = σ. (6 Note that there s no need to tae expectatons because there are no terms related wth t a, (j. The Fsher nformaton matrx n ths case becomes I(θ, R = 1 σ. (7 Because the CRLBs are the dagonal elements of the nverse of the Fsher nformaton matrx [16, p. 40], we obtan I 1 σ (θ, R = {t d } t d. (8 1 From (8, therefore, we can obtan the CRLBs of cloc offset and sew for the Gaussan delay model gven n (1 and (13,.e., σ Var (ˆθ t d ( {t d }, ˆR Var ˆR ( σ t d {t d }. ˆθ Note that the jont Es of cloc offset ( and sew ( are effcent estmators: Frst, from (11, we obtan [ ] E ˆR ( = E [ ] t a, td E [ ] t a, t d. (9 From (17, we have E [ t a, ] = R t d + (θ + d, (30 E [ ] t a, = R + (θ + d. (31 [ ] Insertng (30 and (31 nto (9, we obtan E ˆR ( =R, whch shows that ˆR ( s an unbased estmator. Secondly, from (11, we also obtan Var ˆR ( = Var t a, td Var ta, {t d }. (3 From (17, we have Var ( t a, = (j Var (n( = σ t d, (33 Var ( t σ a, = = σ. (34 Insertng (33 and (34 nto (3, we obtan σ Var ˆR ( = {t d }, whch s the CRLB n (13. Smlar procedures can be taen for ( to show that t s also an effcent estmator. ˆθ APPENDIX B PROOF OF PROPOSITION CR ˆR It s straghtforward to see that [ ] E ˆRCR ( = E [ R + 1 d( ( ( s unbased because ] = R. (35 As for the varance, we frst construct from (14 a lelhood functon for the observaton t a, ( and the parameter R : 1 ( t L(R ; t a, ( = exp a, ( R ( ( 4πσ 4σ. (36 Then the frst and the second dervatves of the log lelhood functon are gven by ln L(R ; t a, ( R = t a, ( R ( σ (, (37 ln L(R ; t a, ( R = ( σ. (38 From (37, we can see that ˆRCR ( s n fact the estmator maxmzng the log-lelhood functon. Also from (38, applyng the same steps to derve a CRLB, we can obtan the lower bound of cloc sew as follows Var ˆRCR ( σ (. (39 Because the varance of Var ˆRCR ( = Var CR ˆR ( s gven by ( R + d( ( = σ (, (40 ˆR CR ( attans the lower bound gven by (39.

11 11 REFERENCES [1] Y.-C. Wu, Q. Chaudhar, and E. Serpedn, Cloc synchronzaton of wreless sensor networs, IEEE Sgnal Process. Mag., vol. 8, no. 1, pp , 011. [] M. Ahlaq and T. R. Sheltam, RTSP: An accurate and energy-effcent protocol for cloc synchronzaton n WSNs, IEEE Trans. Instrum. Meas., vol. 6, no. 3, pp , Mar [3] D. Mac, A. Ageev, and A. Somov, Power consumpton reducton n wreless sensor networs through optmal synchronzaton, n Proc. IMTC 09, May 009, pp [4] J. Elson, L. Grod, and D. Estrn, Fne-graned networ tme synchronzaton usng reference broadcasts, n Proc. of 5th Symp. Operatng System Desgn and Implementaton, Dec. 00, pp [5] T. Trump, Maxmum lelhood trend estmaton n exponental nose, IEEE Trans. Sgnal Process., vol. 49, no. 9, pp , Sep [6] R. T. Rajan and A.-J. van der Veen, Jont rangng and cloc synchronzaton for a wreless networ, n Proc. of CAMSAP 011, Dec. 011, pp [7] S. P. Chepur, R. T. Rajan, G. Leus, and A.-J. van der Veen, Jont cloc synchronzaton and rangng: Asymmetrcal tme-stampng and passve lstenng, IEEE Sgnal Process. Lett., vol. 0, no. 1, pp , Jan [8] S. B. Moon, P. Selly, and D. Towsley, Estmaton and removal of cloc sew from networ delay measurements, n Proc IEEE INFOCOM, New Yor, NY, Mar [9] K. S. Km, Asynchronous source cloc frequency recovery through aperodc pacet streams, IEEE Commun. Lett., vol. 17, no. 7, pp , Jul [10] A. Manwarng, J. P. R. Szewczy, D. Culler, and J. Anderson, Wreless sensor networs for habtat montorng, n Proc. WSNA 0. New Yor, NY, USA: ACM, 00, pp [11] D. L. Mlls, Internet tme synchronzaton: The networ tme protocol, IEEE Trans. Commun., vol. 39, no. 10, pp , Oct [1] S. Ganerwal, R. Kumar, and M. B. Srvastava, Tmng-sync protocol for sensor networs, n Proc. SenSys 03. New Yor, NY, USA: ACM, 003, pp [13] K. S. Km, Comments on "eee 1588 cloc synchronzaton usng dual slave clocs n a slave", IEEE Commun. Lett., vol. 18, no. 6, pp , Jun [14] K.-L. Noh, Q. M. Chaudhar, E. Serpedn, and B. W. Suter, Novel cloc phase offset and sew estmaton usng two-way tmng message exchanges for wreless sensor networs, IEEE Trans. Commun., vol. 55, no. 4, pp , Apr [15] A. Guchhat and R. M. Karth, Jont mnmum varance unbased and maxmum lelhood estmaton of cloc offset and sew n oneway pacet transmsson, n Proc. IEEE 81st Vehcular Technology Conference (VTC015-Sprng, May 015, pp [16] S. M. Kay, Fundamentals of statstcal sgnal processng: Estmaton theory. Upper Saddle Rver, NJ, USA: Prentce-Hall, Inc., 1993, vol. 1. [17] J. Elson and D. Estrn, Tme synchronzaton for wreless sensor networs, n Proc. IPDPS 001, San Francsco, CA, USA, Apr [18] T. Mzrah and Y. Moses, Usng reverseptp to dstrbute tme n software defned networs, n Proc. of ISPCS 014, Sep. 014, pp [19] N. Kyoung-lae and E. Serpedn, Parwse broadcast cloc synchronzaton for wreless sensor networs, n Proc. WoWMoM 007, Jun. 007, pp. 1 6.