PMD monitoring in traffic-carrying optical systems and its statistical analysis

Differential group delay (DGD) experienced by the optical signal in in-service terrestrial optical fiber systems has been monitored for the first time without the requirement of looping-back, in which the live traffic carried in the fiber was used as the probing signal. The relationship between the measured DGD using this technique and the actual fiber PMD parameter is formulated and verified by field experiments. ©2008 Optical Society of America OCIS codes: (060. 2330) Fiber optics communications; (260.5430) Polarization. References and links 1. B. L. Heffner, “Automated Measurement of Polarization Mode Dispersion Using Jones Matrix Eigenanalysis,” IEEE Photon. Technol. Lett. 4, 1066–1069 (1992). 2. R. M. Jopson, L. E. Nelson, and H. Kogelnik, “Measurement of Second-Order Polarization-Mode Dispersion Vectors in Optical Fibers,” IEEE Photon. Technol. Lett. 11, 1153–1155 (1999). 3. P. Williams, “PMD measurement techniques and how to avoid the pitfalls,” J. Opt. Fiber Commun. Rep. 1, 84-105 (2004). 4. M. Boroditsky, M. Brodsky, N. J. Frigo, P. Magill, and J. Evankow, “Estimation of eye penalty and PMD from frequency-resolved in-situ SOP measurements,” Proc. 17th Annual Meeting of the IEEE Lasers and Electro-Optics Society, (Piscataway, 2004), pp.8889. 5. S. X. Wang, A. M. Weiner, M. Boroditsky, and M. Brodsky, “Monitoring PMD-induced penalty and other system performance metrics via a high-speed spectral polarimeter,” IEEE Photon. Technol. Lett. 18, 17531755 (2006). 6. B. Fu and R. Hui, “Fiber chromatic dispersion and polarization-mode dispersion monitoring using coherent detection,” IEEE Photon. Technol. Lett. 17, 1561–1563 (2005). 7. R. Hui, R. Saunders, B. Heffner, D. Richards, B. Fu, and P. Adany, “Non-blocking PMD monitoring in live optical systems,” Electron. Lett. 43, 53–54 (2007). 8. M. Karlsson, J. Brentel, and P. A. Andrekson, “Long-Term Measurement of PMD and Polarization Drift in Installed Fibers,” J. Lightwave Technol. 18, 941–951 (2000). 9. H. Kogelnik, L. E. Nelson, and R. M. Jopson, “Polarization-mode dispersion,” in Optical Fiber Telecommunications IVB, I. P. Kaminov and T. Li, Eds. (Academic, New York, 2002). 10. G. Bosco, B. E. Olsson, and D. J. Blumenthal, “Pulsewidth distortion monitoring in a 40-Gb/s optical system affected by PMD,” IEEE Photon. Technol. Lett. 14, 307-309 (2002). 11. M. Karlsson and H. Sunnerud, “PMD impact on optical systems: Singleand multichannel effects,” J. Opt. Fiber Commun. Rep. 1, 123-140 (2004). 12. Y. Li and A. Yariv, “Solutions to the dynamical equation of polarization-mode dispersion and polarizationdependent losses,” J. Opt. Soc. Am. B. 17, 1821-1827( 2000). 13. J. P. Elbers, C. Glingener, M. Duser, and E.Voges, “Modelling of polarisation mode dispersion in singlemode fibres,” Electron. Lett. 33, 1894–1894 (1997). 14. C. Antonelli, A. Mecozzi, K. Cornick, M. Boroditsky, and M. Brodsky, “PMD-induced penalty statistics in fiber links,” IEEE Photon. Technol. Lett. 17, 1013-1015 (2005). 15. C. Antonelli, A. Mecozzi, M. Brodsky, and M. Boroditsky, “A Simple Analytical Model for PMD Temporal Evolution,” in Optical Fiber Communication Conference and Exposition and The National Fiber Optic Engineers Conference, Technical Digest (CD) (Optical Society of America, 2006), paper OWJ4. #95305 $15.00 USD Received 21 Apr 2008; revised 7 Aug 2008; accepted 22 Aug 2008; published 26 Aug 2008 (C) 2008 OSA 1 September 2008 / Vol. 16, No. 18 / OPTICS EXPRESS 14057


Introduction
In high speed optical fiber communication systems, polarization mode dispersion (PMD) is one of the most important factors of performance degradation.Traditionally the PMD parameter of a fiber can be measured by a number of techniques, such as Jones Matrix Eigenanalysis, Poincare Sphere Analysis and Mueller Matrix method [1][2][3].Fig. 1 shows the Poincare sphere representation of signal polarization vector.With the frequency change of the optical signal which propagates through an optical fiber, the output state of polarization (SOP) rotates on the Poincare sphere around the principle state of polarization (PSP) vector Ω .For the polarization states ) (ω S and of two the frequency components selected from the optical signal shown in Fig. 1, if the separation between their azimuth angles is Δθ, the DGD of the fiber can be found as , where ω Δ is the frequency difference between these two components.Obviously, Δθ has to be small enough so that this linearization is valid for the measurement of the 1 st -order PMD.In practice, in order to measure Δθ, several different SOP settings of the input optical signal have to be used to complete a Jones matrix or a Mueller matrix.In addition, both the Jones matrix and the Mueller matrix techniques require the synchronization between the PSP settings of the input optical signal and the polarimeter measurement at the output side.As the consequence, these traditional PMD measurement techniques require the accesses to both ends of the fiber, which prevents their application from monitoring in-service optical systems since the source and the receiver of live optical networks are at distance and usually are not accessible at the same time.However, there is clearly a need for a more practical approach that supports a network provider's planning and route design process for possible capacity upgrading and the system characterization has to be done without disrupting customer traffic.
Recently, several techniques were proposed for in-situ evaluation of PMD utilizing the optical signal carried in the fiber as the probe signal.A heterodyne polarimeter with an RF spectrum analyzer was used to estimate the PMD-induced system penalty by measuring the state of polarization 'string' length in a nonintrusive way [4].A similar technique with a higher measurement speed was proposed in [5], which uses direct detection with a high resolution optical spectrum analyzer consisting of an InGaAs line-scan camera and a virtually imaged phase array.We have proposed a simplified method to directly measure the DGD in traffic-carrying optical links using coherent detection and RF signal processing [6].For all these non-intrusive PMD monitoring techniques, the SOP of the input optical signal is not adjustable and the measurement of Δθ in Fig. 1 is therefore not feasible.In fact, the core angle α Δ shown in Fig. 1 is usually measured in these in-service monitoring techniques, because it only depends on the relative polarization walk-off between two frequency components within the optical spectrum of the probe.However, since it is Δθ instead of Δα which represents the DGD between the fast and the slow axis of the fiber, it is important to find the relationship between them.Δα is generally smaller than Δθ and ω α τ Δ Δ = p represents the actual DGD seen by the probing signal.In order to correctly interpret the results obtained by the in-service PMD monitoring technique, it is important to rigorously examine the relationship between the PMD parameter of the fiber and the DGD measured by the technique demonstrated in references [6] and [7].

Theoretical analysis
Figure 1 indicates that Δα is related to Δθ by, where, β represents the angle between point A and the PMD vector Ω .When θ Δ is small enough, which can be ensured by choosing appropriate frequency difference ω Δ , Eq.( 1) can be simplified to, In Stokes space, the well-known PSP model indicates that a long fiber can be regarded as a wave plate with the time retardation equals to the modulus of the PMD vector in the fiber, while the principle axis of the wave plate is aligned with the slow axis of the PMD vector.Thus, the angle between the input polarization state of the signal in S and the fiber PMD- vector is also equal to β and therefore, In a Cartesian coordinator, the PMD vector can be decomposed into three orthogonal components, , where x a , y a and z a are unit vectors, and thus 2 3 . When each of the three orthogonal components 1 Ω , 2 Ω and 3 Ω follows an independent Gaussian distribution with zero mean and the same standard deviation q, the statistics of PMD vector will exhibit a Maxwellian distribution [8], In general, the fiber PMD parameter is regarded as its mean DGD which is related to the parameter q by, π τ where, the mean DGD f τ is the average value of the Maxwellian distribution shown in Eq. (4).In practice, using the live traffic carried in the fiber as the probing signal is critical for the in-service monitoring of live optical systems.Since the SOP of the input optical signal is determined by the laser in the transmitter, it is relatively stable.Without losing generality, one can arbitrarily assume that the SOP of the input optical signal is in S = (1, 0, 0), then The combination of Eqs. ( 2 sin Note that, a Maxwellian distribution is referred to as a Chi distribution with 3 degrees of freedom because it is related to three independent components 1 Ω , 2 Ω and 3 Ω .In our case, Eq. ( 7) indicates that τ p is only related to two of the three independent orthogonal components, and therefore, it should follow a Chi distribution with 2 degrees of freedom, which is also known as Rayleigh distribution and its probability density function can be expressed as [9], The mean value of this distribution is 2 / π τ From Eq. ( 5) and ( 9), the relationship between f τ and p τ can be easily found as, π τ τ It is worth noting that quantity τ p in eq.( 7) is a "partial" DGD, which is in fact the projection of the fiber PMD vector perpendicular to the signal SOP direction.However, the mean and the statistic distribution of the actual fiber DGD can be derived from the measured τ p , and it is sufficient for most of the practical applications since it is directly related to system eyeclosure penalty [10,11].
In practical fiber-optic systems, polarization-depend loss (PDL) may exist in addition to PMD.When PDL is taken into accounted, the output polarization state S will vary with optical frequency as [12], where, Λ is the differential attenuation slope (DAS) vector which is related to PDL vector Γ .Λ can be decomposed into three orthogonal and independent random Gaussian components with the same standard deviation q′ when PDL is small enough: Where, ϕ is the angle between the vectors S × Ω and S S × × Λ ) ( .Equation (12) indicates that the distribution of ω α Δ Δ still follows Rayleigh statistics and its mean value is, Where, ( ) 2 / q q L ′ = As an example, with a 2dB PDL, the value of L will be 0.021 and the difference between the mean values of τ p with and without PDL is only 1.03%.Therefore one can generally conclude that the impact of PDL on PMD measurement is negligible when system PDL is less than 2 dB.

Experimental setup and results
We have assembled an experimental setup for in-service PMD monitoring using coherent detection [6] as schematically shown in Fig. 2, where a small portion of the optical signal is tapped from the transmission link for the measurement., where P L and P s are the optical powers of the signal and the LO [7].After the heterodyne IF spectrum is amplified, two RF filters with 1GHz bandwidth are used to select two different frequency components of the signal and their central frequency difference is 10 GHz.The measurement is relatively independent of modulation format of the optical signal since the signal average power is used for measurement.The two frequency components selected by the RF bandpass filters can be any part within the modulated signal spectrum.By measuring the differential polarization walk-off between the two frequency components, the first-order DGD experienced by the optical signal can be evaluated.The accuracy of the PMD measurement in the laboratory environment was verified by using a PMD emulator and setting 90 = β [6].The smallest DGD that can be measured by the current setup is about 0.3 ps which was verified in a system without DGD.This measurement error is believed to be mainly caused by electrical circuit noise.In a previous field trial, we have also demonstrated that the system Q margin was inversely proportional to the instantaneous DGD measured by this technique [7].
In the current measurement apparatus, since polarization scrambling is used for LO, the variation of signal SOP at the fiber output cannot be monitored.If the signal SOP needs to be measured, one can programmatically switch the SOP of the LO between three orthogonal polarization states on the Poincare sphere and performing Stokes parameter analysis of the detected IF signal.
In order to verify the statistical distribution predicted by our analysis, we have recently carried out a number of field trials in various long-distance terrestrial fiber-optic systems carrying DWDM traffics at 10Gb/s data rate with non-return-to-zero (NRZ) modulation.Fig. 3 shows the results of DGD measurements at Sprint's Kansas City switch site and -20dBm of signal optical power was tapped to perform the measurement.Figure 3 [13] where the statistics of PMD-induced system impairments was numerically simulated.Since the eye-closure penalty in the receiver depends on the alignment between the SOP of the optical signal and the PSP of the optical fiber [14,15], a Rayleigh distribution was expected.From in-service system monitoring point of view, our coherent detection technique evaluates τ p which has a mean value of p τ .Eq.( 10) can be used to convert this result to the more traditionally defined mean DGD of the fiber, f τ .The PDL of this system was estimated to be approximately 0.997dB through another measurement, therefore the impact of PDL in the PMD measurement is only about 0.26% and is negligible.To the best of our knowledge, this is the first PMD measurement reported in commercial DWDM systems carrying live traffic and without the requirement of loopingback.

Conclusion
In conclusion, PMD monitoring in traffic-carrying DWDM optical fiber systems is reported for the first time without the requirement of looping-back.The simple relationship between the partial DGD measured with the coherent detection technique and the actual PMD parameter of the fiber is theoretically derived and verified, which allows the accurate evaluation of the PMD parameter in installed fiber systems without disturbing the commercial traffic.The measured partial DGD statistics fits well with a Rayleigh distribution as predicted by the theory.

Fig. 1 .
Fig. 1.Poincare sphere representation of polarization vectors and output SOP rotation with optical frequency change

( 14 )
Under the small PDL assumption, the relationship between PMD, PDL and DAS vectors is, #95305 -$15.00USD Received 21 Apr 2008; revised 7 Aug 2008; accepted 22 Aug 2008; published 26 Aug 2008 (C) 2008 OSA 1 September 2008 / Vol.16, No. 18 / OPTICS EXPRESS 14060 A tunable laser is used as a local oscillator (LO) for coherent heterodyne detection and transmission channel selection.A polarization controller is placed at the output of the local oscillator to randomly scramble the SOP of the local oscillator.When the SOP of the LO is aligned with the received optical signal at the frequency component ω , the relative angular walk-off Δα between ) of LO.Under this condition the angle can obtained through the measured IF intensity ) (

Fig. 2 .
Fig. 2. Block diagram of the coherent PMD monitor (a) shows the result of 268-hour continuous measurement of partial DGD, τ p as the function of time for a fiber link between Kansas City and Chicago which is approximately 900 km, while Fig.3(b) shows the statistical distribution of τ p which is composed of approximately 480,000 data points.The correlation time of this link is about 0.5 h as shown in the inset of Fig. 3(b).Our 268 hours of monitoring is equivalent to 536 uncorrelated samples, which is reasonably sufficient to reconstruct a statistic distribution.The solid line in Fig. 3(b) is a Rayleigh distribution which fits well to the measured partial DGD, while as a comparison the dotted line in the same figure shows a Maxwellian distribution which is obviously not a good fit.It is noticed that Fig. 3(b) looks very similar to Fig. 1 in

Fig. 3 .
Fig. 3. (a) Normalized partial DGD versus time measured over a 900km link; (b) normalized statistic distribution of (a).Solid lines in (b): Rayleigh distribution, dotted lines: Maxwellian distribution with the same mean value.Inset in (b) is the autocorrelation function.Time (h)