Covariance-Based Estimation for Clustered Sensor Networks Subject to Random Deception Attacks

In this paper, a cluster-based approach is used to address the distributed fusion estimation problem (filtering and fixed-point smoothing) for discrete-time stochastic signals in the presence of random deception attacks. At each sampling time, measured outputs of the signal are provided by a networked system, whose sensors are grouped into clusters. Each cluster is connected to a local processor which gathers the measured outputs of its sensors and, in turn, the local processors of all clusters are connected with a global fusion center. The proposed cluster-based fusion estimation structure involves two stages. First, every single sensor in a cluster transmits its observations to the corresponding local processor, where least-squares local estimators are designed by an innovation approach. During this transmission, deception attacks to the sensor measurements may be randomly launched by an adversary, with known probabilities of success that may be different at each sensor. In the second stage, the local estimators are sent to the fusion center, where they are combined to generate the proposed fusion estimators. The covariance-based design of the distributed fusion filtering and fixed-point smoothing algorithms does not require full knowledge of the signal evolution model, but only the first and second order moments of the processes involved in the observation model. Simulations are provided to illustrate the theoretical results and analyze the effect of the attack success probability on the estimation performance.


Introduction
Nowadays, communication networks are widely used to cope with signal estimation problems in engineering, economy, health or security, among others, since they generally provide more robust and precise estimators of the target signal than a single sensor. Conventional centralized and distributed fusion estimation architectures have been widely studied under a state-space approach (see e.g., [1][2][3][4] and references therein). In addition, assuming that the evolution model of the signal is not fully known and only covariance information is available, centralized and distributed fusion estimation algorithms have been proposed for sensor networks affected by different network-induced random uncertainties (see e.g., [5][6][7]). A comprehensive survey of recent developments in estimation and fusion for networked systems with randomly occurring phenomena can be found in [8]. The key theories and methodologies of distributed multisensor data fusion are comprehensively reviewed in [9] and a survey of distributed fusion estimation algorithms with applications in networked systems, including an interesting analysis of some network-induced uncertainties, is provided in [10].
These conventional centralized and distributed fusion architectures require connecting all sensor nodes to a central processor, which sometimes can involve a serious communication burden, especially when the number of sensors is large. Additionally, in the distributed fusion estimation architecture, equipping each single sensor with a local estimator may be sometimes unaffordable. To overcome these shortcomings, a usual practice in engineering consists of arranging the sensors in clusters and selecting one node per cluster (called cluster-head node) acting as a local processor. Each cluster head first collects measurements from the regular sensor nodes of its cluster to generate a local estimator and then all local estimators are collected in the central processor, where the distributed fusion estimator is generated. This structure, called hierarchical cluster-based estimation structure, provides an efficient way to analyze big data in a great variety of application fields [11]. Clustering is also widely used, for example, in underwater sensor networks established for military purposes, such as anti-submarine warfare, communications, positioning and guidance [12]. An essential issue when dealing with complex networked systems consisting of multiple clusters or individual subsystems is coordination control and the design of appropriate consensus protocols; relevant results concerning these problems in different kinds of networked systems are found in [13][14][15][16]. A state-of-the-art and comprehensive survey on clustering approaches can be found in [17] and a detailed outline of some modern energy-efficient clustering approaches to improve the lifetime of wireless sensor networks can be seen in [18]. Recently, a comprehensive review with comparisons and classifications of different optimized clustering approaches according to different metrics has been carried out in [19].
A remarkable disadvantage of using sensor networks for signal estimation is that their reliability can be compromised by possible cyber-attacks from adversaries. For this reason, cyber-security of networked systems is drawing considerable research attention. The most typical kinds of attacks are the denial-of-service attacks and the deception attacks [20]. While the first ones strike at data availability by obstructing the flow of information through the network, the second ones violate data integrity by injecting false information that modifies the real data packets. Such false information can include a wrong sensor measurement or control input, an incorrect time-stamp or a wrong identity of the sending device. The fusion estimation problem for stochastic signals from measured data coming from a sensor network subject to deception attacks is currently an important focus of research (see e.g., [20][21][22][23][24][25][26] and references therein). In [20], new insecurity conditions for the state estimation problem under false data injection attacks are proposed. In [21], the variance-constrained distributed filtering problem is studied for time-varying systems subject to multiplicative deception attacks with bounded attack noises. A distributed recursive filtering algorithm is designed in [22] for a class of discrete time-delayed systems subject to both uniform quantization and intermittent deception attacks. In [23], the security-guaranteed filtering problem is addressed for a class of nonlinear discrete time-delayed systems with both stochastic sensor saturations and deception attacks. The centralized security-guaranteed filtering problem for linear discrete time-invariant stochastic systems with multi-rate sensor nodes is addressed in [24], when deception attacks are launched during the transmission of information from the sensors to the centralized filter. In [25], the cluster-based covariance intersection fusion estimation problem under stochastic deception attacks is investigated. An integrated analysis of event-triggered fault detection and fault estimation is proposed in [26] for a class of discrete-time stochastic systems subject to unknown disturbances and stochastic deception attacks.
Research motivation.
In this paper, distributed fusion filtering and fixed-point smoothing algorithms are designed for clustered sensor networks subject to stochastic linear deception attacks. This study is motivated by the following challenges: (i) To find out fusion estimation structures for multi-sensor networked systems that reduce transmission burdens (with respect to both centralized and distributed fusion) and the cost of embedding a local processor in each sensor (with respect to distributed fusion).
(ii) To design recursive and easily implementable estimation algorithms with high estimation accuracy, fusing unreliable multi-sensor data, corrupted by possible deception attacks, under general assumptions on the target signal that do not require the full knowledge of the system state-space model.
Paper contributions. In light of our motivations, the main contributions of this paper are highlighted as follows: (i) The fusion estimation problem in multi-sensor networked systems when the sensors are grouped into clusters is investigated assuming, on the one hand, that the measurements are subject to stochastic deception attacks and, on the other, that the signal evolution model is not necessarily known, but only information about its first and second order statistical properties (covariance information) is available. This covariance-based estimation approach is more general than the conventional one based on the full knowledge of the state-space model of the system since, in that case, the mean and autocovariance function of the signal process can be calculated. Hence, this estimation approach provides a comprehensive framework to deal with a great variety of stochastic signals.
(ii) A two-stage fusion estimation algorithm is designed for both filtering and fixed-point smoothing problems. In the first stage, each local processor collects measurements (subject to random deception attacks) from its cluster to generate local least-squares linear estimators; in this stage, the innovation approach is used to design the local estimation algorithms, which are recursive and computationally simple. In the second stage, in order to improve the estimation performance, all these local estimators are transmitted to the fusion center, where fusion estimators are obtained by matrix-weighted linear combinations of the local estimators; the weight matrices are computed by minimizing the mean squared estimation error, for which the cross-correlation matrices between any two local estimators need to be previously calculated.
In contrast with previous papers concerning the fusion estimation problem under a covariance-based approach, this is the first time that a cluster structure of the sensor network is used and stochastic deception attacks are considered. In comparison with the existing literature about the estimation problem in clustered sensor networks, either subject to random deception attacks or not, the main difference is the kind of information required for the derivation of the algorithms: full knowledge of the state-space model in the existing literature and only covariance information (which covers the conventional formulation based on state-space model and also more general situations) in the current paper.
Paper structure. The measurement model framework is presented in Section 2, where the equations of the clustering sensor measurements and the stochastic deception attacks, together with the hypotheses under which the distributed estimation problem will be addressed, are specified. In Sections 3 and 4 the first and second stages of the proposed cluster-based fusion estimation structure are developed. The local least-squares linear filtering and fixed-point smoothing estimators, as well as formulas for the error covariance matrices, are obtained by recursive algorithms in Section 3. The proposed distributed fusion estimators are proposed in Section 4 by a matrix-weighted linear combination of the local estimators using the mean squared error as optimality criterion; also, to measure the estimation accuracy, recursive formulas for the error covariance matrices are derived. The performance of the proposed estimation algorithms is analyzed in a numerical simulation study, carried out in Section 5. Finally, some concluding remarks are provided in Section 6.
Notation. Unless otherwise stated, the notation used in the paper is fairly standard. R n denotes the n-dimensional Euclidean space and R m×n the set of all m × n real matrices. When the dimensions of vectors or matrices are not specified, they are assumed to be compatible with algebraic operations. 1 n denotes the n-dimensional column vector with all ones, while the identity and zero matrices are denoted by I and 0, respectively. For any function Γ h,k , depending on h and k, and Ξ (rs) , depending on r and s, we will write Γ k = Γ k,k and Ξ (r) = Ξ (rr) , respectively, for simplicity. The Kronecker product of the matrices A ∈ R m×n and B ∈ R m ×n , which is a matrix in R mm ×nn , is denoted by A ⊗ B. The Hadamard product of the matrices C, D ∈ R m×n , which is also a matrix in R m×n , is denoted by C • D and it is defined by (C • D) ij = C ij D ij . δ k,s is the Kronecker delta function, i.e., δ k,s = 1 when k = s and zero otherwise. The autocovariance function of a second-order process {α k } k≥1 is defined as stands for the mathematical expectation operator.

Estimation Problem Formulation and Measurement Model
In this paper, we consider the distributed fusion estimation problem of a discrete-time stochastic signal from a set of measurements provided by multiple sensors. A cluster-based estimation structure will be adopted: the sensors are grouped into different clusters, each of them connected with a local processor which, in turn, is connected with a global fusion center (FC). The distributed fusion estimation process operates in two stages. In the first one, the measured outputs provided by the sensors of each cluster are sent to the corresponding local processor, where local least-squares (LS) linear estimators are obtained; it is assumed that, during the transmission to the local processors, the clustering sensor measurements are subject to stochastic linear deception attacks. In the second stage, the local estimators received from all the local processors are gathered and fused in the FC, where the proposed distributed signal estimators are generated by a matrix-weighted linear combination of the local estimators using the mean squared error as optimality criterion.
The aforementioned estimation problem will be addressed without requiring full knowledge of the evolution model generating the signal process. Instead, it will be assumed that the signal mean function is zero and its covariance function is factorizable (covariance-based estimation approach). More precisely, the following assumption is required ( [7]):

(A1)
The n x -dimensional signal {x k } k≥1 is a zero-mean second-order process and its autocovariance function is expressed in a separable form; namely, E The signal process of the most common signal evolution models (e.g., the signal of linear systems, or that of uncertain systems with a sum of multiple multiplicative noise terms) meets this assumption (A1) and, hence, the covariance-based estimation approach provides a comprehensive framework to cope with different signal evolution models, thus overcoming the necessity of deriving specific algorithms for each situation.

Clustering Sensor Measurements and Stochastic Deception Attacks
Consider a sensor network and assume that the sensor nodes are grouped into L clusters to measure the stochastic signal of interest. Specifically, assume that each cluster r = 1, . . . , L is made up of m r sensors that provide measurements of the signal according to the following model: where C (r(i)) k are known matrices and z (r(i)) k ∈ R n z is the signal measured output from the i-th sensor of the r-th cluster at time k, which is transmitted to the r-th local processor to obtain the local LS linear signal estimators. The following assumption is required on the measurement noises, v (r(i)) k k≥1 , i = 1, . . . , m r , r = 1, . . . , L: (A2) The measurement noises of different clusters are independent and, for each r = 1, . . . , L, the noises v (r(i)) k k≥1 , i = 1, . . . , m r , are zero-mean second-order white processes with known covariance matrices For each cluster r = 1, . . . , L, the transmissions of the measured outputs z (r(i)) k , i = 1, . . . , m r , to the r-th local processor are affected by random linear deception attacks and the deceptive signal injected by the attackers, ξ (r(i)) k , is described by: This signal involves two parts: the first one neutralizes the true information and the second one is the blurred information (noise) added by the attackers. These noises, w (r(i)) k k≥1 , i = 1, . . . , m r , r = 1, . . . , L, are assumed to satisfy the following requirement: (A3) The attack noises of different clusters are independent and, for each r = 1, . . . , L, the noises w (r(i)) k k≥1 , i = 1, . . . , m r , are zero-mean second-order white processes with known covariance matrices

Measurements Received by the Local Processors
Usually, in practice, the attacks may randomly succeed or not. So, taking into account this random nature of the attacks, for every r = 1, . . . , L, the measurements received by the r-th local processor are modelled by introducing different sequences of Bernoulli random variables, λ For each r = 1, . . . , L, and i = 1, . . . , m r , the value λ (r(i)) k = 1 models a successful attack into the i-th communication channel from the r-th cluster, meaning that only noise w (r(i)) k arrives to the r-th local processor; conversely, the value λ (r(i)) k = 0 models a failed attack, which means that the real measured output z (r(i)) k is received by the r-th local processor. Taking these considerations into account, the following model for y (r(i)) k , the measurements received by the r-th local processor, is considered: T are known and their entries are easily calculated taking into account that Finally, the following independence hypothesis is also assumed: (A5) For r = 1, . . . , L and i = 1, . . . , m r , the signal process {x k } k≥1 and the processes v and λ (r(i)) k k≥1 are mutually independent.

First Stage: Local LS Linear Estimators
To start with, as indicated previously, our aim is to calculate in every local processor, r = 1, . . . , L, LS linear estimators of the signal based on the measurements received from all the sensors of the r-th cluster. Therefore, to estimate the signal x k in the local processor r at time k + N, we consider all the measurements received from all the sensors i = 1, . . . , m r , of the r-th cluster, up to time k + N; that is, the measurement set y , made up of all the measurements received by the r-th processor at each sampling time h, the problem at hand is formulated as that of determining local LS linear estimators of the signal x k based on the vectors y

Stacked Model for the Measurements Received by the Local Processors
Taking into account Equation (4) for the measurements received by the r-th local processor, the following model for the above-defined vectors y (r) k is clearly deduced: where The following statistical properties of the processes involved in model (6), which will be used to address the LS linear estimation problem, are easily inferred from the model assumptions (A1)-(A5) stated in Section 2: , r = 1, . . . , L, are independent zero-mean noise processes with , r = 1, . . . , L, are independent sequences of independent random matrices with known means, Λ Moreover, the Hadamard product properties guarantee that, for any random matrix G independent of Λ (r) where K λ (r) k is calculated from (5).

Recursive Local LS Linear Filtering and Fixed-Point Smoothing Algorithms
This subsection is devoted to the design of recursive algorithms, at each local processor r = 1, . . . , L, for the LS linear filtering and fixed-point smoothing estimators based on the measurements received from all the sensors of the r-th cluster. In other words, for the r-th cluster, r = 1, . . . , L, the aim in this subsection is the design of algorithms to obtain the local LS linear estimators, x (r) k/k+N , N ≥ 0, of the signal x k based on the vectors y (r) h , h ≤ k + N, given by (6); specifically, a recursive algorithm for the local LS filter x (r) k/k , k ≥ 1, and a recursive algorithm for the local LS smoother x (r) k/k+N , for fixed k ≥ 1 and N = 1, 2, . . ., will be derived.
For this purpose, we will use the measurement innovations rather than the raw measurements (innovation approach), where the innovation at time k is defined as µ k/k−1 , being y h , h ≤ k − 1). Taking into account the properties of the innovation process, the following expression for the LS linear estimators is derived: Taking orthogonal projections in (6) and using properties (P2)-(P4), we get that the one-stage observation predictor is y h/h−1 and, consequently, the innovation is given by Furthermore, from property (P5), it is clear that the innovation covariance matrix, Π Using now (6) for y Then, by defining we obtain that the estimators x L , and, by substitution in (9), we obtain that µ Finally, since (11) guarantees that . Bearing in mind the above results, the following filtering algorithm is derived.

Recursive Local LS Linear Filtering Algorithm.
For the r-th cluster, r = 1, . . . , L, the local filtering estimators, x (r) k/k , and the error covariance matrices, Σ (r) The vectors o and the innovation covariance matrices, Π (r) where the matrices Σ y (r) k are given in (7). The general expression (8) for the LS linear estimators is also the starting point to derive the following covariance-based recursive fixed-point smoothing algorithm. The derivation is omitted, since this algorithm can easily be deduced by an analogous reasoning to that used in Theorem 2 of [6].

Recursive Local LS Linear Fixed-point Smoothing Algorithm.
For the r-th cluster, r = 1, . . . , L, starting from the filter, x (r) k/k , the local LS linear fixed-point smoothers, x (r) k/k+N , N ≥ 1, are calculated as where the matrices X k,k+N , obey the following recursive formula Starting from the error covariance matrix of the filter Σ (r) k/k , the fixed-point smoothing error covariance matrices,

Second Stage: Distributed Signal Estimators
As we have already mentioned, once the local LS linear estimators, x (r) k/k+N , r = 1, . . . , L, have been obtained, they are sent to the FC and our goal is to fuse these local estimators to obtain distributed estimators, x k/k+N , N ≥ 0, as matrix-weighted linear combinations that minimize the mean squared estimation error.
Defining the stacked vectors by comprising all the local estimators, X k/k+N = x and applying the LS criterion, the proposed distributed fusion estimators are given by: where E X k/k+N X T k/k+N = E x . Hence, the derivation of the distributed estimators in (16) requires to obtain the cross-covariance matrices between the local ones Since the initial condition of the recursive Formula (15) for the local smoothing estimators is the local filter, the cross-covariance matrices, Σ x (rs) k/k+N , N ≥ 1, between the local smoothers will be recursively obtained by starting from the cross-covariance matrices, Σ x (rs) k/k , between the local filters.
Next, for t = r, s, we use (6) for y (t) k and (14) for y (t) k/k−1 ; then, using again that, from the OPL, k−1 , the following expression is clear: where the matrices Σ y (rs) k are given in (7).

Cross-Covariance Matrices between Local Smoothing Estimators
k+N , it is clear that the cross-covariance matrices between any two local smoothing estimators x • To obtain the expectations Φ Next, we use that, from the OPL, the estimator is orthogonal to the estimation error, to write k+N−1 , the following expression for Φ (rs) k,k+N is immediately derived: • Finally, the expectations Ψ k , k ≥ 1; r, s = 1, . . . , L.

Distributed Filtering and Fixed-Point Smoothing Estimators
The distributed estimators, x k/k+N , are calculated from (16), while a formula for the error covariance matrices, Σ k/k+N ≡ E (x k − x k/k+N ) (x k − x k/k+N ) T , is easily derived from assumption (A1) and (16). More specifically, we can state the following results: T be the vectors constituted by the local estimators calculated from the recursive algorithms in Section 3.2; then, the distributed filtering and smoothing estimators are given by k/k+N , r, s = 1, . . . , L, are obtained by (18)) for N ≥ 1, with initial condition given in (17). The error covariance matrices of the distributed estimators are computed by

Numerical Simulation Study
A numerical simulation example is presented to illustrate the application of the algorithms designed in the current paper to estimate a two-dimensional signal. Specifically, such algorithms have been implemented by using MATLAB software and both local and distributed filtering and fixed-point smoothing error variances of the signal components have been calculated, for 100 iterations, in order to examine the estimation accuracy. The effect of the stochastic deception attacks on the estimation performance has also been analyzed.
Consider a two-dimensional signal, x k , whose evolution is described by the following model: The multiplicative noise, ε k k≥1 , and the additive noise, α k k≥1 , are standard white gaussian scalar noises. The initial signal x 0 is a gaussian two-dimensional random vector with zero mean and covariance matrix E[x 0 x T 0 ] = I. The noise sequences and initial signal vector are assumed to be mutually independent; then, it is easy to see that the signal covariance function is given by This example shows that assumption (A1) on the signal autocovariance function is fulfilled by uncertain systems with state-dependent multiplicative noise.
Suppose that a sensor network, comprising 12 sensors, is deployed to measure the stochastic signal x k . These sensors are grouped into three clusters and the number of sensors in each cluster is m 1 = 3, m 2 = 4 and m 3 = 5, respectively. Scalar measurements are provided by the 12 sensors, according to model (1), with the following time-invariant observation matrices: k } k≥1 , r = 1, 2, 3, are independent zero-mean Gaussian white process with variance 10; so the measurement noises in each cluster r = 1, 2, 3, are correlated with R r(ij) k = 10(v (r) ) 2 , k ≥ 1; i, j = 1, . . . , m r . Assume that the transmissions to the local processors are subject to linear deception attacks and the data injected by the attackers are described by (2), where the attack noises are defined as k , for all i = 1, . . . , m r , with w (1) = 0.1, w (2) = 0.25, w (3) = 1, and {ζ (r) k } k≥1 , r = 1, 2, 3, independent standard Gaussian white processes. Clearly, in each cluster r = 1, 2, 3, the attack noises are correlated with S r(ij) k = (w (r) ) 2 , k ≥ 1; i, j = 1, . . . , m r . According to the theoretical study, for r = 1, 2, 3, the measurements received at the r-th local processor are modelled by (3); that is, the attacks are considered to take place randomly and the statuses of the attacks are described by independent sequences of independent Bernoulli random variables with known probabilities P(λ (r(i)) k First, considering λ (r(i)) k = 0.1i, for r = 1, 2, 3, and i = 1, . . . , m r , the performance of the local filtering estimators has been compared with that of the local fixed-point smoothing estimators; for r = 1, 2, 3, these local estimators are computed in each r-th processor using only the measurements received from the sensors of the corresponding r-th cluster. For the first signal component, Figure 1 displays the error variances of the local filtering estimators, Σ (r) k/k 11 , and the local fixed-point smoothing estimators Σ (r) k/k+N 11 , for N = 1, 3, 5, 7, are displayed. For the three processors, Figure 1 shows that, as expected, the error variances corresponding to the local smoothers are less than those of the local filters. It is also observed that the accuracy of the smoothers at each fixed-point, k, becomes better as the number of available observations, k + N, increases, although this improvement is practically imperceptible for N > 3 in processor 1, N > 5 in processor 2 and N > 7 in processor 3. Note that the best estimation accuracy is obtained in processor 1 and the worst one in processor 3; therefore, we can conclude that, when the filtering performance is better, the improvement of the smoothers is less significant and, actually, it becomes practically negligible for lower values of N. Similar results and conclusions are obtained for the second signal component.
Next, considering again λ (r(i)) k = 0.1i, for r = 1, 2, 3, and i = 1, . . . , m r , the performance of the distributed filtering and fixed-point smoothing estimators have been compared with that of the local filtering estimators of the three processors. Figure 2 displays, for both first and second signal components (a = 1, 2), the error variances of the local filtering estimators Σ (r) k/k aa , r = 1, 2, 3, and the distributed filtering and smoothing estimators Σ k/k+N aa , N = 0, 1, 3. On the one hand, this figure shows that the error variances of the distributed fusion filtering estimators are smaller than those of every local filter; consequently, the distributed fusion filtering estimators outperform all the local ones, agreeing with what theoretically is expected, since all the information of the three clusters is available to the distributed filters, while the local ones are based on the information of a single cluster. On the other hand, agreeing with the comments made about Figure 1, from Figure 2 it is also concluded that the error variances of the distributed fusion smoothers are less than those of the filters and, at each fixed-point, k, the performance of the smoothers is better as the number of available observations, k + N, increases (note that, in this case, the improvement is practically imperceptible for N > 1).    Finally, assuming the same success probability of attack during the transmissions of the measured outputs of all the sensors in the three clusters, λ (r(i)) k = λ, r = 1, 2, 3, i = 1, . . . , m r , we analyze how this probability affects the error variances of the distributed filtering estimators. Taking into account that these variances stabilize after a sufficiently large number of iterations, only the results of the last iteration (k = 100) are shown in Table 1. Specifically, the distributed filtering error variances of the first and second signal components, Σ 100/100 aa , a = 1, 2, together with the corresponding percent variation rates are shown in Table 1, when λ varies from 0.1 to 0.9. From this table, we conclude that, for both signal components, as λ increases, the distributed filtering error variances become higher, meaning that, as expected, the lower the probability of successful attacks is, the better estimators are obtained. In addition, for both signal components, it is observed that the deterioration of the estimators is more significant (the percent variation rate of the error variance is larger) for high attack probabilities (λ ≥ 0.7).

Concluding Remark
This paper contributes to the literature on the distributed fusion filtering and fixed-point smoothing problems in clustering sensor networks subject to random deception attacks. Actually, the sensor nodes of the network are grouped into some clusters and each cluster is connected to a local processor, which gathers the measured outputs of all the sensors in the cluster. During this transmission, the measured outputs sent by the sensors to the local processor may be deceptively modified by malicious attackers, who launch random deception attacks with known probabilities of success that may be different at each sensor. Aggregating the information received at every local processor, and using an innovation approach, recursive and easily implementable local filtering and fixed-point smoothing algorithms have been designed without requiring full knowledge of the signal evolution model, but only the first and second order moments of the processes involved in the measurement equations. Once the local estimators are available, they are transmitted to the global fusion center where they are fused to obtain the optimal matrix-weighted fusion filter and fixed-point smoother, under the minimum mean-squared error criterion. The simulation results have illustrated the effectiveness of the proposed algorithms, showing that the distributed fusion estimators outperform the local ones. In addition, the influence of the probability of successful attacks on the estimation performance has been analyzed, concluding that the estimation error variance is larger (and, consequently, worse estimations are obtained) for higher attack probabilities of success. Some new research directions stemming from this work would be: -Considering other different attack strategies (not necessarily linear ones) to cover more general scenarios in practice. Another interesting generalization would be studying the case where the Bernoulli random variables modelling the success or failure of the attacks are not necessarily independent, but they are correlated at consecutive sampling times or they obey a Markovian dependence structure, for example. - Addressing the estimation problem in clustering sensor networks with a topology represented by a directed or undirected graph, thus generalizing the current study by allowing the exchange of information between sensors and between clusters.
-Investigating other interesting issues related with clustering sensor networks (cluster consensus, sensor mobility, overlap of clusters, cluster position awareness, efficient and energy, uniform clustering or stability of the cluster, among others).
Author Contributions: All the authors contributed equally to this work. R.C.-Á., A.H.-C. and J.L.-P. provided original ideas for the proposed model and collaborated in the derivation of the estimation algorithms; they participated equally in the design and analysis of the simulation results; and the paper was also written and reviewed cooperatively.

Conflicts of Interest:
The authors declare no conflict of interest.