SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS ON OCTOBER 15TH 2014

1

Dual Polarized Modulation and Reception for Next
Generation Mobile Satellite Communications

arXiv:1803.07314v1 [eess.SP] 20 Mar 2018

Pol Henarejos, Student, IEEE,, Ana Perez-Neira, Senior, IEEE

Abstract—This paper presents the novel application of Polarized Modulation (PMod) for increasing the throughput in mobile
satellite transmissions. One of the major drawbacks in mobile
satellite communications is the fact that the power budget is
often restrictive, making unaffordable to improve the spectral
efficiency without an increment of transmitted power. By using
dual polarized antennas in the transmitter and receiver, the PMod
technique achieves an improvement in throughput of up to 100%
with respect to existing deployments, with an increase of less than
1 dB at low Eb /N0 regime. Additionally, the proposed scheme
implies minimum hardware modifications with respect to the
existing dual polarized systems and does not require additional
channel state information at the transmitter; thus it can be
used in current deployments. Demodulation (i.e. detection and
decoding) alternatives, with different processing complexity and
performance, are studied. The results are validated in a typical
mobile interactive scenario, the newest version of TS 102 744
standard (Broadband Global Area Network (BGAN)), which aims
to provide interactive mobile satellite communications.
Index Terms—Satellite Communications, Polarized Modulation, Dual Polarized Antennas, Interactive Services

I. I NTRODUCTION

M

ULTIPLE-input multiple-output (MIMO) schemes
were introduced as a promising way to notably increase the spectral efficiency (SE) using multiples antennas at
transmission and/or reception [1]–[3]. In contrast to terrestrial
communications, where the transmitter can obtain Channel
State Information at the Transmitter (CSIT), in satellite links
it is impossible to maintain CSIT updated due to the long
round trip time. By the time the satellite receives the feedback
parameters, the channel varies and the CSIT becomes outdated.
Among the different approaches that do not use CSIT,
Vertical Bell Laboratories Layered Space-Time (VBLAST)
scheme and successive improvements present a simple way
to increase the achievable rate with a relative increase of the
processing complexity [4]–[7] in the absence of CSIT. However, VBLAST introduces interference between the streams
since all signals are transmitted through all antennas without
any interference pre-cancellation. In consequence, the signals
must be transmitted with higher amplitude to obtain the same
error rates compared with the single stream case.
In contrast to VBLAST, Spatial Modulation (SM) appeared
recently to increase the SE [8]–[10]. In SM, the bits of
P. Henarejos is with the Communications Systems Division of Catalonian Technological Telecommunications Center, Barcelona, e-mail:
pol.henarejos@cttc.es.
A. Perez-Neira is with the Signal and Communication Theory Department
of Universitat Politcnica de Catalunya and with the Catalonian Technological
Telecommunications Center, e-mail: ana.isabel.perez@upc.edu
Manuscript submitted on October 15th 2014

information are split in blocks; some are coded as antenna
indices and the rest are transmitted through the antennas
that are selected by the indices. By estimating the antenna
indices, the receiver can recover the bits used for this purpose.
Nevertheless, this approach is very sensitive to the channel
variations and requires an accurate channel estimation as well
as spatially uncorrelated channels [11]–[13].
In satellite scenarios, due to the Line of Sight (LOS), the
spatial components become correlated at the receiver side
although the transmitting antennas may be separated at half
wavelength. Hence, in the absence of scatterers, the receiver
only discovers a single transmission path and the sensitivity of
the terminal is not enough to distinguish the different spatial
signatures and detect the antenna indices. Because of this,
SM does not seem suitable as it does not provide sufficient
diversity in satellite scenarios.
On the contrary, the polarization channel provides more
diversity and may be used for this kind of schemes [14].
Although dual polarized antennas were used for broadcasting,
where subscribers only tuned a single polarization, recent studies unveil that dual-polarized MIMO channel is richer in terms
of diversity [15]. Additionally, the use of dual polarized antennas is increasingly motivated by the new possibilities arising,
together with the newest standards including dual polarized
MIMO, such as Digital Video Broadcasting-Next Generation
broadcasting system to Hand-held (DVB-NGH [16]). Finally,
research projects such as [17], [18] reported that the throughput can be increased as in a conventional MIMO system
if more antennas, and the consequent radio frequency (RF)
chains [19]–[21], are added in order to multiplex polarizations.
The price to be paid is that the complexity of the satellite payload increases since interference among polarizations appears.
For instance, extending the VBLAST strategy to dual polarized
schemes requires higher transmit power to maintain the same
Quality of Service (QoS) in point-to-point clients [17].
The primary motivation behind the present work is to apply
the SM concept to dual polarized communications in the
challenging mobile satellite channels. Hereafter, this scheme
is referred to as polarized modulation (PMod). Indeed, it
is not polarization multiplexing since only one polarization
is activated at a time and therefore precludes the presence
of interference. Although this work has been triggered as
an attempt to apply a simple diversity technique as SM to
the satellite scenario, the paper could also be viewed as the
extension to satellite communications of the PMod idea that
has been reported previously in optical communications [22].
However, from the authors’ knowledge there is no literature
describing in detail PMod demodulation scheme (detection

SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS ON OCTOBER 15TH 2014

and decoding) for optical communications, being polarization
multiplexing far more common.
The scheme to be proposed contains the following contributions:
• The proposed PMod technique exploits the polarization
diversity in satellite scenarios, where the spatial diversity
is highly penalized.
• PMod does not require CSIT and increases the throughput
maintaining the robustness based on the polarization
diversity.
• Usually, satellite systems operate with dual polarized
antennas and thus the proposed scheme does not require
additional antennas. The minimum requirement is to use
a dual polarized feeder.
• The success of this scheme lies not only on the simplicity
of the transmission technique, but also on the receiver
design, which is also one main contribution of the present
work, together with the performance evaluation. Note that
the information is conveyed not only in the transmitted
bit stream, but also in the polarization.
• Finally, as we demonstrate in a maritime mobile satellite
L-band scenario, the result is an increase of the overall
performance in terms of throughput, whereas guaranteeing a minimum QoS and requiring a minimum increase
in power usage. The best performance is obtained for low
order modulations, where the proposed method achieves
a gain of 2 when it is compared with a basic system
without PMod.
II. S YSTEM M ODEL
Let us consider a MIMO system where transmitter and
receiver are equipped with a single antenna with dual polarization, and a Rician frequency flat channel. Each symbol
contains b + 1 bits of information, where b bits are modulated
within the constellation S and the remaining bit is used for
polarization selection. This remaining bit is denoted as c and
the modulated bits as s ∈ C. We would like to remark that the
information is conveyed through the symbol s as well as bit
c. The channels across the polarizations 1 and 2 are denoted
h11 ∈ C and h22 ∈ C, respectively, and their respective crosschannels as h21 ∈ C and h12 ∈ C. All channel coefficients hij
follow a Rice statistical distribution with (K,σh ) parameters
with a pair-wise correlation ρij = [0, 1]. The received signals
for polarizations 1 and 2 are denoted as y1 ∈ C and y2 ∈ C,
respectively, and ωi ∈ C follows Additive White Gaussian
Noise (AWGN), ω ∼ CN (0, N0 I2 ).
Depending on the value of c ∈ {0, 1}, the s symbols are
conveyed using one polarization or the other. Hence, we can
formulate the system model as follows:

 




y1
h11 h21
1−c
ω1
=
s+
(1)
c
y2
h12 h22
ω2
or in a more compact way as:
= [h1 h2 ] cs + ω

Since this scheme adds an additional bit to the transmission
by keeping the same power budget, the throughput gain of
PMod with respect to Single-input Single-output (SISO) case
is
1
b+1
=1+ .
(4)
G=
b
b
For higher order modulations, (4) is asymptotically 1 and
thus the proposed PMod scheme increases the gain for low
order modulations. For instance, the gain is 2 for Binary PhaseShift Keying (BPSK) modulation or 1.5 for Quadrature PhaseShift Keying (QPSK) modulation. As low order modulations
are used in low signal to noise ratio (SNR) regime, it is clear
that PMod increases significantly the throughput gain G in
low SNR systems. This is exactly the scenario for mobile
satellite communications where shadowing, fading and power
limitations cause low SNR.
III. D EMODULATION S CHEMES
The implementation of the receiver derives into several
approaches depending on the scenario constraints. Since PMod
transmits a single stream, we aim to extract this stream to be
processed into a SISO decoder. This scheme offers two main
advantages:
•
•

x̂ = arg max p (y|x, H)

(5)

x∈M

where x = cs and M is the symbol alphabet of x. Hereafter,
we assume that the noise contribution is Gaussian. Note that
there is no restriction with the characteristics of the channel
matrix. Thus, we do not take any assumption on the statistical
independence of H.
Hence, (5) can be simplified as
2

x̂ = arg min ky − Hxk .

(6)

x∈M

Based on the expression proposed in (6) four different
demodulators are proposed:
•

(3)

•

where hi is the channel corresponding to the ith polarization.

•

= Hcs + ω,

Reduces the complexity drastically since the signal processing is one dimensional.
Can be combined with existing SISO decoders, maintaining the compatibility with the current standards.

The reception scheme is illustrated in Fig. 1.
With this scheme, PMod−1 implements one of the four
demodulation schemes that are introduced in this section,
estimates the bit c and manipulates the received signal y to
produce the signal r, which is capable to be used by a SISO
decoder.
The optimal demodulation scheme is based on the Maximum a Posteriori (MAP) criteria, which is equivalent to
the Maximum-Likelihood (ML) criteria for the case where
the transmitted symbols are equiprobable. Thus, the optimal
receiver can be derived from the expression

(2)

y

2

•

First approach: zero forcer.
Second approach: per-symbol detection.
Third approach: per-hard-bit detection.
Fourth approach: per-soft-bit detection.

SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS ON OCTOBER 15TH 2014

A. First Approach: Zero Forcer
In this approach a zero forcing pre-filter is applied before
detecting the information in x. We note that the zero forcer is
2
df
the solution of dx
H = 0, where f (x) = ky − Hxk . This is
equivalent to apply the filter
−1 H
W = HH H
H .
(7)
Therefore we obtain the signal after the processing as
follows:
−1 H
z = Wy = cs + HH H
H ω.
(8)

In the case where the signal is being transmitted through
polarization 1, i.e. c = 0, z1 contains the signal plus noise
and z2 only receives the noise; and the reciprocal case for
c = 1, z1 contains only the noise and z2 conveys the signal
plus the noise. Therefore, to decide on c we propose a power
detector. It is denoted as:


(9)
ĉ = arg max |zi |2 − 1.
i

Once the receiver knows the used polarization path, it is able
to decode the symbol s based on the signal r = zĉ+1 .
This solution presents a simple implementation but it is not
a sufficient statistic to decode the whole problem (symbol s
and bit c). Since it computes the envelope of the vector, the
information conveyed through the phases is lost and therefore
it is not sufficient.
Nevertheless, this solution presents a major disadvantage.
HH H could be badly conditioned and thus it may produce an
excessive noise enhancement, doing impossible the demodulation.
B. Second Approach: Per-Symbol Detection
As we stated in the previous section, applying the W filter
may introduce important distortions. Since the solution have
to lie in the subset M, the solution of the first approach may
not be the optimal. Thus, the optimal approach to solve (6) is
performing an exhaustive search over the subset M.
In the particular case of PMod, the transmitted vector x can
be restricted to the subset x ∈ M, where the first vector [ 0s ] of
the set defines the transmission using the first polarization and
the second vector [ 0s ] of the set, using the second polarization.
Hence, the decision rule for demodulating the bit c is based
on ĉ = 0 if x̂ = [ 0s ] and ĉ = 1 if x̂ = [ 0s ]. In the same way,
the signal r can be written as r = xĉ+1 .
This scheme, however, presents a notable increase of computational complexity. The exhaustive search requires to find
the solution amongseveral
 possibilities. The complexity of the
2
initial search is O 2b but due to the restriction of the set

aforementioned, the complexity is reduced to O 2b+1 , which
is equivalent to the complexity of the SISO case where b + 1
bits are conveyed.
Furthermore, the previous demodulation schemes introduce
hard decisions that induce non-linearities, such as sign() or
abs() functions. In the presence of coded information, as it can
be seen in [23], soft decoding outperforms the previous ML
implementation. In the following sections we describe schemes
that introduce soft information.

3

C. Third Approach: Likelihood Ratio with Hard Decision
Usually, to deal with channel impairments, the transmitted
bits are coded. The channel decoder computes the metrics
based on the likelihood of the received signal and is able to
estimate the uncoded bits.
The third approach is based on this philosophy and bit c is
estimated based on the likelihood ratio. If the likelihood ratio
is defined as


2
P
2 s̃k
exp − ky−h
2
σw
P (c = 1|y)
P2
s̃∈S

 , (10)
=
= P
Λ (y) =
2
1 s̃k
P1
P (c = 0|y)
exp − ky−h
2
σ
s̃∈S

w

the decision rule for estimating c depends only on the sign
of log (Λ (y)). In the case where likelihood ratio is greater
than 1, it means that it is more probable that c = 1 and viceversa. Assuming that only b bits are coded, an estimator of
the uncoded c can be stated as:
1 + sign (log (Λ (y)))
.
(11)
ĉ =
2
Once the receiver obtains the estimation of c, it knows which
polarization is being used and thus it can recover the symbol
s using the signal r = yĉ+1 .
Although this scheme uses soft information in the decoding
of symbol s, the decision of bit c is still hard. Thus, if this bit
is also coded, the result is suboptimal. In the next section we
describe how to obtain a soft version of bit c.

D. Fourth Approach: Likelihood Ratio with Soft Decision
The three approaches described above perform hard decision
for the estimation of bit c. However, they can introduce errors
if the system conveys coded information as it was mentioned.
The soft version of bit c corresponds to the log-likelihood,
exactly as the bits b. That is ĉ = log (Λ (y)).
After that, the bit c is soft and can be passed to the soft
decoder. However, there is the problem of which polarization
to choose for decoding. In the previous schemes, since the bit
c is hard, it is possible to process the received signal on the
polarization indicated by c. In the present scheme it is not
possible to decide which polarization conveys information.
To solve this issue we compute the average received signal:
r = P1 y1 + P2 y2 .

(12)

Using the likelihood ratio Λ (y) computed as in (10), and
using
P2 = P (c = 1|y) = 1 − P (c = 0|y) = 1 − P1 ,

(13)

we can rewrite
P2 = P (c = 1|y) =

Λ (y)
.
1 + Λ (y)

(14)

Therefore, the receiver can recover the signal by weighting
the received signals from both polarizations by P1 = 1 − P2
and P2 , respectively. If we assume that the bit c is transmitted
with equal probability, the averaged received signal takes the
following form:
1
(y1 + y2 Λ (y)) . (15)
r = (1 − P2 ) y1 + P2 y2 =
1 + Λ (y)

SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS ON OCTOBER 15TH 2014

Finally, the combined signal r is passed to the decoder in order
to obtain the b bits.
IV. N UMERICAL R ESULTS FOR U NCODED BER
In this section we analyse the results of the proposed
schemes. To compare them, we deploy a system conveying
QPSK symbols in addition to the switching bit c. For this
purpose, we only examine the uncoded bit error rate (BER).
The channel model used corresponds to the Rician maritime
mobile channel model described in the experiment V in [24]
with a correlation factor of ρij . All parameters are summarized
in Table I.
In all results the following labels are used:
1) Reference denotes the reference scenario, i.e. the scenario where single polarization is used.
2) VBLAST is the polarization multiplexing VBLAST coding scheme.
3) PMod ZF is the first approach described in Section III-A.
4) PMod ML is the second approach described Section III-B.
5) PMod HD is the third approach described in Section III-C.
6) PMod SD is the fourth approach described in Section III-D.
7) OSTBC corresponds to the Orthogonal Space Time
Block Codes applied to polarization instead of spatial
diversity [25].
In Fig. 2, we compare the BER of the four PMod schemes.
The four curves are labelled in the same order that they have
been introduced (from the first to the fourth approach).
As we stated in Section III-B, the ML solution provides the
lowest error rate, immediately followed by the fourth solution.
As expected, in the absence of channel coding, the ML
receiver becomes the optimal solution. Although we remark

that ML uses a reduced search space of order O 2b+1 , the
computational complexity is sensibly higher with respect to
the other solutions.
Next to the curve of ML is the pure soft scheme (the fourth).
If we examine the magnified area, we observe the gap between
the ML solution with the pure soft is tight. Hence, we conclude
that the fourth demodulation scheme stays very close to the
optimal solution.
Finally, the third approach, which does not use the conditional mean of the signal, performs close to PMod SD, whereas
the first approach PMod ZF obtains the highest BER.
Hereinafter, we choose the fourth approach, PMod SD, to
compare it with other schemes different to PMod. The reason
is twofold:
• First, it performs a near-optimal ML solution, with a
small gap of 0.05 dB of Eb /N0 for a fixed BER of 10−6 .
• Second, the computational complexity is less expensive
than in ML.
Fig. 3 compares the PMod SD solution with the conventional OSTBC, VBLAST and reference scenario using a QPSK
constellation for all schemes. Note that even though we use
the same constellation for all schemes, the total SE is different
for each scheme. Thus, although we are comparing different

4

schemes with different SE, the most remarkable point is the
fact that the PMod is bounded by OSTBC (lower SE) and
VBLAST (higher SE) and therefore PMod achieves a tradeoff between OSTBC and VBLAST in terms of BER and SE.
In all these schemes, 2 bits/channel use are conveyed. As
expected, OSTBC obtains lowest BER, followed by PMod and
VBLAST. However, OSTBC does not allow to increase the
granularity of the adaptive bit rate. In other words, there is
no choice to transmit 3 bits/channel use. The next step is to
transmit a 16QAM with OSTBC, which is 4 bits/channel use.
Newest standards, such as DVB-S2X [26], aim to include new
modulation schemes to refine the rate adaptation curve.
A. Equal SE Analysis
In contrast to the previous section, where the comparison is
performed maintaining the same constellation, in this section
we analyse the performance of PMod compared with the other
schemes but constraint to the same SE. To do this, we use the
following transmission schemes:
• PMod SD with BPSK constellation.
• VBLAST with BPSK constellation.
• OSTBC with QPSK constellation.
• Reference with QPSK constellation.
In all schemes, 2 bits per channel use are conveyed. Fig. 4
describes the curves of the different throughputs and it is clear
that all curves tend to the same throughput for high SNR.
Fig. 5 depicts the BER for the different techniques. In this
case, OSTBC obtains the lowest BER, followed by PMod SD,
Reference and VBLAST, respectively. As expected, OSTBC
exploits the full diversity of the channel and is closely followed
by PMod. However, one of the advantages of PMod in front of
OSTBC is the ability to increase the granularity of the throughput adaptation. Whereas OSTBC increases the throughput by
a powers of two, PMod can increase the throughput by small
fractions, as it is seen in (4).
V. R ESULTS IN A R EALISTIC S YSTEM C ONTEXT
In this section we describe the implementation of the PMod
solution inside the Broadband Global Area Network (BGAN)
standard. In more detail, we deploy the downlink of the
Next Generation Satellite Communications standard, currently
being redacted at the European Telecommunications Standards
Institute (ETSI) committee (more detail at [27]). This part of
the standard defines the scrambling, turbo coding and mapping
stages, among other procedures. In order to offer flexibility in
terms of data rate, several bearers and subbearers are detailed.
They are different profiles with many combinations of coding
rate and constellations. Focusing in the downlink part, the
symbol rate is 33.6 ksps and the frame length is 80 ms,
where the blocks of coded symbols are not interleaved. In
order to simplify the model, QPSK bearers will be used in all
simulations.
A. Next Generation Satellite Communications Simulation
Framework
We simulate a L-band geostationary satellite with 7 beams
(one desired beam and six interfering beams) and dual polarization. Since the beams are not perfectly orthogonal, we

SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS ON OCTOBER 15TH 2014

consider six adjacent beams at the same frequency subband
as interferences, as well as the cross polarization couplings.
All these values are summarized in Table II and are obtained
via realistic multibeam antenna pattern during the project
Next Generation Waveforms for Improved Spectral Efficiency
(NGW), whose results are summarized in [17]. In more detail,
Fig. 6 illustrates the beam pattern, where the working beam is
marked with a red circumference and the interfering beams as
yellow circumferences. It is important to remark that not all
beams induce the same levels of interference. Depending on
the position of the satellite and the geometry of the reflectors,
the power of interferences varies between beams. In more
detail, Fig. 7 and Fig. 8 illustrate the co-polar and crosspolar coverage for the forward link with contours at 3 dB
(red lines) and 4.5 dB (blue lines). One of the relevant aspect
is the asymmetry of the co-polar and cross-polar gains in each
beam. From these figures, it is clear that gains are different
for each beam spot. Finally, Fig. 9 shows the block diagram
used for the simulations described hereafter.
In Fig. 9 the identified blocks are:
• Forward Error Correction (FEC) Encoder: encodes the bit
stream using the specifications of [27].
• P M od: groups the bits in blocks of size b + 1, modulates
the symbols s and uses the c bits to select the polarization
for each symbol.
• Framing: encapsulates the symbols of each polarization
in a frame defined in [27]. It inserts pilots for channel
estimation, a preamble for synchronization and a header
for modulation-code identification.
• Interference matrix Bi : models the cross polarization
by a factor defined in Table II. B0 corresponds to the
cross-polarized matrix of intended data and B1 , . . . , B6
correspond to the cross-polarized matrices of interfering
beams.
√
P . It is
• P : the signal is amplified by a factor of
important to remark that this is possible due to the fact
that, for each symbol, only a single polarization is active
and thus all power budget P is available, whereas in the
case of VBLAST and OSTBC this factor is divided by 2.
• Li , i = 0, 1: equivalent path-loss for each polarization.
• Hi , i = 0, 1: convolves the signal using the Rician fast
fading channel model.
• Noise: adds the AWGN.
−1
• P M od
: implements one of the schemes.
• FEC Decoder: performs the inverse operation of FEC
Encoder. It implements a Turbo Coder with Systematic
Recursive Convolutional Codes (SRCC).
We consider the Rician maritime mobile channel model
described in the experiment V in [24] and the parameters
described in Table I.
The aim is to evaluate the basic transmission and reception
concepts and schemes; thus, in this work it is assumed perfect
synchronization at the receiver side as well as perfect channel
estimation. Prior to detection of symbol s, one of the four
approaches is performed in order to estimate the bit c and
filter the received signal.
We remark that this scenario includes nongaussian interference. Thus, as we described the PMod solution under this

5

assumption, we need to cope with the interference to minimize
it. To achieve that, the receiver implements a MMSE linear
filter. This configuration mitigates the interferences from the
other beams as well as the other polarization for the detection
of symbol s.
One important aspect is the Faraday Rotation (FR), which
appears at L-band. This effect is caused by the free electrons
in the ionosphere and causes a rotation of the polarization.
Since it changes the polarization, FR may be critical in order
to estimate the bit c. Fortunately, this effect can be reduced
using a circular polarization or performing an estimation and
assuming that the FR remains invariant during the time slot.
An estimation of FR is described in [28] and it can be
applied using the pilot symbols used by the channel estimation.
Nevertheless, for the simulations, we assume that this effect
is corrected.
Finally, in the next stage, the demodulated soft bits are
passed to the turbo decoder and scrambled to obtain the
information bits. In contrast to the previous section, since we
consider interferences in this scenario, we use the signal to
interference plus noise ratio (SINR) in the x-axis rather than
SNR.
B. Comparing PMod Solutions
We compare the four proposed demodulation schemes. In
contrast to Fig. 10, although the ML solution is the optimal in
absence of channel coding, this is not the case in the presence
of coded information. Certainly, the PMod SD scheme produces the lowest BER, followed by PMod HD. Both schemes
use soft bits and, thus, their performance is better than the
hard solutions (PMod ZF and PMod ML).
In order to compare the proposed schemes with the existing
ones, we compare the performance in terms of throughput,
which corresponds to the average rate of successful information delivery and is defined as
T = R(1 − BLER)G.

(16)

This is equivalent to the bitrate (R) of the particular bearer
weighted by the probability of no error in the whole block (1−
BLER) and the throughput gain (G), defined in (4). During all
simulations, a fixed modulation-code is simulated with coding
rate of 0.625 (R = 40 kbps without frame overhead). BLER
is obtained by simulations and corresponds to the number of
erroneous blocks divided by the total number of blocks.
Fig. 11 describes the throughput achieved using the four
schemes. We observe that the four curves are grouped in the
soft and hard receivers. In contrast to the previous section
where the gap between the solutions is tight, in this case the
gap increases notably, making clearer the performance of the
PMod SD/HD.
C. Comparing PMod SD with Other Solutions
In this section we compare the performance of the PMod
SD with OSTBC and VBLAST in the same interference scenario. With the following comparisons we examine different
strategies to increase the throughput.

SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS ON OCTOBER 15TH 2014

Fig. 12 illustrates the coded BER for the different schemes.
As with the uncoded BER (see Fig. 3)), in this case, PMod
SD lies between OSTBC and VBLAST. One important aspect
is that improves the BER of the Reference scenario. This is
positive since PMod increases the SE but also the error rate.
Finally, Fig. 13 illustrates the throughput achieved by each
scheme. The interesting part of this figure is the adaptation of
the rate. For very low Eb /N0 the most effective scheme is
OSTBC. From 3.5 dB, the PMod SD increases the throughput
by a factor of 1.5, followed by VBLAST from 5.5 dB.
This motivates the use of PMod in Adaptive Modulation and
Coding Schemes (AMC).

D. XPD Analysis
In addition to prior comparisons, we also include a crosspolarization discrimination (XPD) analysis for the PMod. The
results are extremely encouraging and reveal that the PMod
scheme is robust in front of cross-polarization impairments.
The reason is twofold:
•

•

For high XPD values, only one polarization carries the
data symbol whereas the other only contains noise. In
this case, the system will decode the symbols s and the
switching bits c correctly.
For low XPD values, both polarizations carry the same
symbol but only one polarization is decoded. In this case,
the probability of error of decoding bit c increases as the
XPD decreases but the probability of error of decoding
the symbol s remains the same. This is motivated by the
fact that, even in the case where the c bit is erroneous
and the decoded polarization is the wrong one, it also
contains the symbol s and thus, is able to decode the s
symbols as if it was decoded from the other polarization.

To analyse the XPD of the PMod technique, the XPD is
defined as follows:


|yc |
,
(17)
XP D = 20 log
|y1−c |
where yc is the signal received at the polarization where the
symbol is transmitted and y1−c is the other one.
Fig. 14 compares the throughput of the four proposed
schemes (PMod ZF, PMod ML, PMod HD and PMod SD) for
different values to the XPD with the reference (Reference).
Note that PMod HD and PMod SD are overlapped in the
fig., although the PMod SD has slightly higher robustness.
Particularly, only for these simulation a fixed SNR of 20 dB
was set, whilst the other parameters remain the same as in
previous figures. As aforementioned, the PMod technique is
robust in front of XPD as it can exploit the fact that the 2/3
of the bits are transmitted through the both polarizations.

6

channel norm. Indeed, the power of the error ξ is defined as
follows:
o
n
2
E h − h̄
2
o
n
(18)
|ξ| =
2
E |h|
where E {x} is the expectation value of the variable x and h̄
is the estimated coefficient. The results can be examined in
Fig. 15.
The PMod scheme becomes more robust in front of the
reference scheme (Reference). Particularly, three of the four
schemes (PMod ZF, PMod HD and PMod SD) offer the
same tolerance, but with the difference that PMod ZF in
addition is able to decode the bit c correctly. This means that
although the scheme may be inaccurate, it is always capable to
decode the bit c. This motivates hierarchical modulations. For
example, using the PMod ZF we could establish a hierarchical
BPSK+QPSK and always succeed on decoding the BPSK
scheme at least.
Finally, it is worth mentioning that the PMod technique
presents a good trade-off between robust techniques but with
less throughput, such as OSTBC, and more throughput available techniques but more power consuming such as VBLAST,
as Fig. 15 depicted.
VI. C ONCLUSIONS
This paper introduces a novel application to mobile satellite
communications of the entitled Polarized Modulation, which
is based on dual polarized antennas. The work shows that with
dual-polarized modulation the throughput can be increased
by a factor of 1 + b−1 in the absence of CSIT in low
Eb /N0 regime and that the transmission results robust to
cross-polarization and imperfect channel estimation. Performance depends on the implemented receiver, that is why in
this paper different alternatives are proposed that trade-off
computational complexity vs. performance. One of the demodulation schemes is based on probabilities, which involves soft
detections and to the authors knowledge it is novel in the
context of either spatial or polarized modulation. Finally, the
proposed techniques have been thoroughly tested and validated
using a maritime mobile satellite scenario and the newest
implementation of the novel ETSI’s standard TS 102 744 [27],
known as BGAN, which is used for interactive mobile satellite
communications. It validates the PMod scheme and demonstrates the enhancement of throughput and the robustness.
Further work is to extend the results and receiver architectures
to more than two polarizations and investigate the PMod in
aeronautical and urban channels. PMod exploits the diversity
of the channel and therefore, whereas the polarization channel
has diversity, the PMod will work as expected. Additionally,
although the union bound for Rayleigh channel is provided,
an interesting action is to study the impact of the averaged
probability of error in Rician channels as well as the mutual
information and capacity analysis.

E. Imperfect Channel Estimation
In this section we analyse the impact of an imperfect
channel estimation. To study this behaviour we introduced an
error ξ into the channel estimation which is normalized by the

ACKNOWLEDGEMENT
The simulator tool was developed in the context of the
European Space Agency (ESA) contract for the project Next

SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS ON OCTOBER 15TH 2014

Generation Waveform for Improved Spectral Efficiency and
for the Satellite Network of Experts (SatNEx IV). The real
measurements used as the interference gains as well as the
coverage patterns were provided by ESA.
This work has received funding from the Spanish Ministry
of Economy and Competitiveness (Ministerio de Economia
y Competitividad) under project TEC2014-59255-C3-1-R and
from the Catalan Government (2014SGR1567).
We would also like to especially thank P. D. Arapoglou
for his valuable comments and appreciations. We also kindly
thank Prof. Miguel Angel Lagunas for his fruitful discussions.
R EFERENCES
[1] Schwarz, A. Knopp, B. Lankl, D. Ogermann, and C. A. Hofmann,
“Optimum-capacity MIMO satellite broadcast system: Conceptual design for los channels,” in Advanced Satellite Mobile Systems, 2008.
ASMS 2008. 4th, 2008, pp. 66–71.
[2] R. T. Schwarz, A. Knopp, D. Ogermann, C. A. Hofmann, and B. Lankl,
“Optimum-capacity MIMO satellite link for fixed and mobile services,”
in Smart Antennas, 2008. WSA 2008. International ITG Workshop on,
2008, pp. 209–216.
[3] P. Arapoglou, P. Burzigotti, M. Bertinelli, A. Bolea Alamanac, and
R. De Gaudenzi, “To MIMO or not to MIMO in mobile satellite
broadcasting systems,” IEEE Trans. Wireless Commun., vol. 10, no. 9,
pp. 2807–2811, 2011.
[4] P. W. Wolniansky, G. J. Foschini, G. D. Golden, and R. Valenzuela,
“V-BLAST: an architecture for realizing very high data rates over the
rich-scattering wireless channel,” in Signals, Systems, and Electronics,
1998. ISSSE 98. 1998 URSI International Symposium on, 1998, pp. 295–
300.
[5] G. D. Golden, C. J. Foschini, R. Valenzuela, and P. W. Wolniansky,
“Detection algorithm and initial laboratory results using V-BLAST
space-time communication architecture,” Electron. Lett, vol. 35, no. 1,
pp. 14–16, 1999.
[6] C. Shen, H. Zhuang, L. Dai, and S. Zhou, “Detection algorithm
improving V-BLAST performance over error propagation,” Electron.
Lett, vol. 39, no. 13, pp. 1007–1008, 2003.
[7] Y. Xue and T. Kaiser, “Adaptive detection of V-BLAST in slowly
time-varying channels,” in Signal Processing Advances in Wireless
Communications, 2003. SPAWC 2003. 4th IEEE Workshop on, 2003,
pp. 477–481.
[8] R. Y. Mesleh, H. Haas, S. Sinanovic, C. W. Ahn, and S. Yun, “Spatial
modulation,” IEEE Trans. Veh. Technol., vol. 57, no. 4, pp. 2228–2241,
2008.
[9] Y. A. Chau and S.-H. Yu, “Space modulation on wireless fading
channels,” in Vehicular Technology Conference, 2001. VTC 2001 Fall.
IEEE VTS 54th, vol. 3, 2001, pp. 1668–1671.
[10] R. Mesleh, H. Haas, C. W. Ahn, and S. Yun, “Spatial modulation
- a new low complexity spectral efficiency enhancing technique,” in
Communications and Networking in China, 2006. ChinaCom ’06. First
International Conference on, 2006, pp. 1–5.
[11] M. Di Renzo and H. Haas, “Bit error probability of sm-MIMO over
generalized fading channels,” IEEE Trans. Veh. Technol., vol. 61, no. 3,
pp. 1124–1144, 2012.
[12] M. Di Renzo, H. Haas, and P. M. Grant, “Spatial modulation for
multiple-antenna wireless systems: a survey,” IEEE Commun. Mag.,
vol. 49, no. 12, pp. 182–191, 2011.
[13] R. Mesleh, S. Ganesan, and H. Haas, “Impact of channel imperfections
on spatial modulation OFDM,” in Personal, Indoor and Mobile Radio
Communications, 2007. PIMRC 2007. IEEE 18th International Symposium on, 2007, pp. 1–5.
[14] P. Horvath and I. Frigyes, “Investigation of the polarization properties of
satellite channels with multiple antennas,” in Antennas and Propagation,
2006. EuCAP 2006. First European Conference on, 2006, pp. 1–4.
[15] K. P. Liolis, J. Gomez-Vilardebo, E. Casini, C. Peel, and A. I. PerezNeira, “Statistical modeling of dual-polarized MIMO land mobile satellite channels,” IEEE Trans. Commun., vol. 58, no. 11, pp. 3077–3083,
2010.
[16] Digital Video Broadcasting (DVB); Next Generation broadcasting system
to Handheld, physical layer specification (DVB-NGH), DVB Std.

7

[17] P. Henarejos, M. A. Vazquez, G. Cocco, and A. Pérez-Neira, “Forward
link interference mitigation in mobile interactive satellite systems,” in
in Proceedings of AIAA International Communications Satellite Systems
Conference (ICSSC), 14-17 October 2013.
[18] G. Gallinaro, E. Tirrò, F. D. Cecca, M. Migliorelli, N. Gatti, and
S. Cioni, “Next generation interactive s-band mobile systems: challenges
and solutions,” International Journal of Satellite Communications and
Networking, vol. 32, no. 4, pp. 247–262, 2014.
[19] H. Vasseur, “Degradation of availability performance in dual-polarized
satellite communications systems,” IEEE Trans. Commun., vol. 48, no. 3,
pp. 465–472, 2000.
[20] Arapoglou, P. Burzigotti, A. B. Alamanac, and R. de Gaudenzi, “Capacity potential of mobile satellite broadcasting systems employing dual
polarization per beam,” in Advanced satellite multimedia systems conference (asma) and the 11th signal processing for space communications
workshop (spsc), 2010 5th, 2010, pp. 213–220.
[21] Zorba, M. Realp, M. Lagunas, and A. I. Perez-Neira, “Dual polarization
for MIMO processing in multibeam satellite systems,” in Signal Processing for Space Communications, 2008. SPSC 2008. 10th International
Workshop on, 2008, pp. 1–7.
[22] M. Karlsson and E. Agrell, “Which is the most powerefficient modulation format in optical links?” Opt. Express,
vol. 17, no. 13, pp. 10 814–10 819, Jun 2009. [Online]. Available:
http://www.opticsexpress.org/abstract.cfm?URI=oe-17-13-10814
[23] P. Fertl, J. Jalden, and G. Matz, “Performance assessment of MIMObicm demodulators based on mutual information,” IEEE Trans. Signal
Processing, vol. 60, no. 3, pp. 1366–1382, 2012. [Online]. Available:
http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=6093979
[24] M. Sellathurai, P. Guinand, and J. Lodge, “Space-time coding in mobile
satellite communications using dual-polarized channels,” IEEE Trans.
Veh. Technol., vol. 55, no. 1, pp. 188–199, 2006.
[25] A. Perez-Neira, C. Ibars, J. Serra, A. del Coso, J. Gomez, and M. Caus,
“MIMO applicability to satellite networks,” in Signal Processing for
Space Communications, 2008. SPSC 2008. 10th International Workshop
on, 2008, pp. 1–9.
[26] Digital Video Broadcasting (DVB); Second generation framing structure,
channel coding and modulation systems for Broadcasting, Interactive
Services, News Gathering and other broadband satellite applications;
Part 2: DVB-S2 Extensions (DVB-S2X), DVB Std.
[27] TS 102 744-2-1: Satellite Component of UMTS (S-UMTS); Family SL
satellite radio interface; Part 2: Physical Layer Specifications; Sub-part
1: Physical Layer Interface, ETSI Std.
[28] F. J. Meyer and J. B. Nicoll, “Prediction, detection, and correction of
faraday rotation in full-polarimetric l-band sar data,” Geoscience and
Remote Sensing, IEEE Transactions on, vol. 46, no. 10, pp. 3076–3086,
2008.

SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS ON OCTOBER 15TH 2014

TABLE I
S CENARIO M AIN PARAMETERS
Maritime
Rician flat fading
10
2 Hz
Jakes
ρij = 0.5
35786 km
187.05 dB
200 kHz
−12.5 dB/K
L (1.59 GHz)
0.625
40 kbps

10 0
10 -1
10 -2
10 -3

BER

Profile
Channel model
Rician K factor
Doppler shift
Doppler spectrum
Stream correlation
Path distance
Path loss
Bandwidth
Terminal G/T
Carrier band
Code rate
Bitrate

8

10 -4
10 -5
10 -6
10

10 -8

TABLE II
D ATA C OUPLING P OLARIZATION M ATRIX + I NTERFERENCE M ATRICES

Data

0
1

Interference

2
3
4
5
6

Interference

40.8
−11.6

3.7
−12.3
8.7
 −13
3.6
 −6.7
13.4
 −8.9
8.9
−4.7

11.6
−3.7

matrix (dB)

−11.6
40.8

−12.3
3.7 
−13
8.7 
−6.7
3.6 
−8.9
13.4 
−4.7
8.9 
−3.7
11.6

-2

0

2

4

6

8

10

12

14

SNR [dB]

Fig. 2. Comparison of the uncoded BER of the four proposed PMod
techniques conveying a QPSK constellation.

10 0
10 -1
10 -2

BER

Index

PMod ZF
PMod ML
PMod HD
PMod SD

-7

10 -3
10 -4
10 -5
Reference
OSTBC
VBLAST
PMod SD

10 -6
10 -7

0

2

4

6

8

10

12

14

16

SNR [dB]

Fig. 1. Reception scheme. PMod−1 applies one of the following demodulation schemes to estimate the bit c and prepare the signal r to be processed
by a common SISO decoder.

Pol Henarejos was worn in Barcelona, Catalonia.
He received the M.Sc. degree in telecommunication engineering from UPC in May 2009 and the
European Master of Research on Information and
Communication Technologies (MERIT) in 2012.
He joined CTTC in January 2010 as a research engineer. In 2010, he participated in European projects
implementing real receivers using the filterbank multicarrier approach. After that, he was involved in
industrial projects based on implementations of the
physical layer aspects of LTE in prototypes. That
research permitted the deployment of the new standard for interactive services
in satellite communications and testing innovative techniques.
His interests comprise the implementations of the physical layer of radio
communications in real devices to theoretical studies in resource management
in multiuser scenarios.

Fig. 3. Comparison of the uncoded BER of the PMod SD with other existing
techniques conveying a QPSK constellation.

Ana I. Pérez-Neira is full professor at UPC (Technical University of Catalonia) in the Signal Theory
and Communication department.
She has been the leader of 20 projects and has
participated in over 50 (10 for ESA). She is author
of 50 journal papers (20 related with Satcom) and
more than 200 conference papers (20 invited). She
is co-author of 4 books and 5 patents (1 on satcom). Since 2008 she is member of EURASIP BoD
(European Signal Processing Association) and since
2010 of IEEE SPTM (Signal Processing Theory
and Methods). She has been guest editor in 5 special issues and currently
she is editor of IEEE Transactions on Signal Processing and of Eurasip
Signal Processing and Advances in Signal Processing. She has been the
general chairman of IWCLD09, EUSIPCO11, EW14 and IWSCS14. She has
participated in the organization of ESA conference 1996, SAM04 and she
is the general chair of ASMS16. She has been in the board of directors of
ETSETB (Telecom Barcelona) from 2000-03 and Vicepresident for Research
at UPC (2010-13). She created UPC Doctoral School (2011). Currently, she is
Scientific Coordinator at CTTC (Centre Tecnolgic de Telecomunicacions de
Catalunya). She is the coordinator of the Network of Excellence on satellite
communications, financed by the European Space Agency: SatnexIV.
Her research topics are in: multi-antenna signal processing for satellite
communications and wireless and in physical layer scheduling for multicarrier

SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS ON OCTOBER 15TH 2014

9

68
66

Throughput [kbps]

64
62
60
58
56
PMod SD
VBLAST
OSTBC
Reference

54
52
0

2

4

6

8

10

12

14

16

18

20

SNR [dB]

Fig. 4. Comparison of the throughput of the PMod SD with other existing
techniques constraint to the same SE.

10 0

BER

10 -1

10

Fig. 6. Considered beam pattern to perform realistic simulations. Working
beam is marked with a red circumference and interfering beams as yellow
circumferences.

-2

10 -3
PMod SD
VBLAST
OSTBC
Reference

10 -4

0

2

4

6

8

10

12

14

16

18

20

SNR [dB]

Fig. 5. Comparison of the Uncoded BER of the PMod SD with other existing
techniques constraint to the same SE.

systems.

Fig. 7. Co-polar coverage for the forward link with contours at 3 dB (red
lines) and 4.5 dB (blue lines).

SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS ON OCTOBER 15TH 2014

10

60
PMod ZF
PMod ML
PMod HD
PMod SD

Throughput [kbps]

50

40

30

20

10

0
0

1

2

3

4

5

6

7

8

SINR [dB]

Fig. 11. Comparison of the throughput of the four proposed PMod techniques
conveying a QPSK constellation.
Fig. 8. Cross-polar coverage for the forward link with contours at 3 dB (red
lines) and 4.5 dB (blue lines).
10 0
Reference
OSTBC
VBLAST
PMod SD

10 -1

BER

10 -2
10 -3
10 -4
10 -5
10 -6
10 -7

Fig. 9. Block diagram of the simulation framework.

0

1

2

3

4

5

6

7

8

SINR [dB]

Fig. 12. Comparison of the coded BER of the PMod SD with other existing
techniques conveying a QPSK constellation.
10 0

80
10 -1

Reference
OSTBC
VBLAST
PMod SD

70
60

Throughput [kbps]

BER

10 -2

10 -3

10 -4

PMod ZF
PMod ML
PMod HD
PMod SD

10 -5

10 -6

0

1

2

50
40
30
20
10

3

4

5

6

7

8

SINR [dB]

0
0

1

2

3

4

5

6

7

8

SINR [dB]

Fig. 10. Comparison of the coded BER of the four proposed PMod techniques
conveying a QPSK constellation.

Fig. 13. Comparison of the throughput of the PMod SD with other existing
techniques conveying a QPSK constellation.

SUBMITTED TO IEEE TRANSACTIONS ON COMMUNICATIONS ON OCTOBER 15TH 2014

60

Throughput [kbps]

50

40

30

20
PMod ZF
PMod HD
PMod SD
Reference
PMod ML

10

0
0

1

2

3

4

5

6

7

8

XPD [dB]

Fig. 14. Comparison of the throughput with respect of XPD of the different
techniques conveying a QPSK constellation.

80
Reference
VBLAST
PMod ZF
PMod HD
PMod SD
OSTBC
PMod ML

70

Throughput [kbps]

60
50
40
30
20
10
0
0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

|ξ|2

Fig. 15. Impact of the imperfect channel estimation in the different techniques
conveying a QPSK constellation.

11

