Ionospheric trends in Southern Hemisphere stations due to the increasing greenhouse gases concentration ANNALS OF GEOPHYSICS, 58, 3, 2015, A0332; doi:10.4401/ag-6707 A0332 Ionospheric trends in Southern Hemisphere stations due to the increasing greenhouse gases concentration José V. Venchiarutti1, Blas F. de Haro1,*, Ana G. Elías1,2 1 Universidad Nacional de Tucumán, Departamento de Física, Tucumán, Argentina 2 Consejo Nacional de Investigaciones Científicas y Técnicas, CONICET, Argentina ABSTRACT The lower and middle atmosphere present long-term variations in tem- perature and other parameters linked to anthropogenic sources, such as the increase in greenhouse gases concentration since the start of the in- dustrial era. Some examples are the well known temperature increase in the troposphere and stratosphere cooling. Upper atmosphere parameters also present long-term variations. While in the case of the middle and lower atmosphere it is debated whether the origin of the trends is prima- rily anthropogenic or solar, in the upper atmosphere other sources are also able to induce long-term changes, such as long-term variations in geo- magnetic activity and secular variations of the Earth’s magnetic field. In this paper, trends of the F2 layer critical frequency, foF2, measured at three Southern Hemisphere stations (Brisbane, Canberra and Christchurch) are analyzed to determine the importance of increasing greenhouse gases concentration effect. According to our results for the pe- riod covering solar cycles 18 to 21 (period 1944-1986), it is more important than natural forcings. Update estimates including solar cycle 23 are pre- sented although the difficulty due to two of the three stations present big data gaps during cycle 22 and traditional ionospheric filtering is no longer entirely reliable for cycle 23. The aim of this study is to contribute both to an active area of aeronomy as is the study of trends in the upper atmos- phere, and to the understanding of climate change. 1. Introduction Since the beginning of the 1990s there has been a great interest in the study of long-term trends in the upper atmosphere and several papers have been pub- lished since then on this topic (see Lastovicka et al. [2012] for a very good review). Some of these studies associate trends in the ionosphere with cooling in the thermos- phere due to increased concentration of greenhouse gases [Roble and Dickinson 1989, Rishbeth 1990]. If the concentration of CO2 is doubling, according to theo- retical models a cooling of 30-40 K in the thermosphere should be expected together with a decrease in air den- sity, a decrease of 15-20 km of the peak height of the F2 layer, hmF2, and a decrease in the critical frequency of the F2 layer, foF2, less than 0.5 MHz globally. However the overall pattern of hmF2 and foF2 trends estimated with experimental data of several worldwide stations is very complex and would not completely agree with the greenhouse gases hypothesis. Another possible cause of the ionospheric trends would be the long-term variation of geomagnetic ac- tivity [Mikhailov and Marin 2000] and the secular vari- ations of the Earth’s intrinsic magnetic field [Foppiano et al. 1999, Elias 2009, Cnossen 2014]. In this paper we estimate foF2 trends measured at three Southern Hemisphere ionospheric stations, which have records since solar cycle 18, and analyze them in terms of the cooling expected from the anthropogenic effects. Being able to determine the causes of the ob- served ionospheric trends, would contribute to the study of variability in the ionosphere for predictive purposes and the problems of global change directly linked to in- creased concentration of greenhouse gases. 2. Data analysis Monthly median values of foF2 from Damboldt and Suessmann database [2012] (ftp://ftp.ips.gov.au/ wdc-data/iondata/medians/Damboldt/) were used at 12 LT. Three Southern Hemisphere stations which are listed in Table 1 were selected which had long enough records and few data gaps. Canberra is the only that presents an almost unin- terrupted record for the period 1941-2009. Christchurch presents a first almost complete period 1941 to1986, and then returns in 1993-2009, while in the case of Bris- bane the almost complete data periods are 1943-1986 and 1998-2009. Due to the big data gap in solar cycle 22, trends are estimated first for the period 1944-1986 which covers four complete solar cycles, 18 to 21. Then Article history Received December 3, 2014; accepted May 22, 2015. Subject classification: Ionospheric trend, EUV proxies, Forecasts, Dynamics. they are updated to include solar cycle 23 with special care in the filtering process required in ionospheric trend estimation. 2.1. Trends during period covering solar cycles 18 to 21 In order to minimize spurious effects introduced by the filtering process here applied [Lastovicka et al. 2012, Elias 2014], trends should be assessed for the same period in all the cases analyzed and include complete solar cycles. In our case this period corresponds to 1944- 1986, that is solar cycles 18 to 21. The steps followed to estimate foF2 trend are: 1. Seasonal variation was filtered out assessing 12- month running mean of each data series analyzed. 2. Periods with Rz> 160 were eliminated due to saturation effect occurring in the F2 region. 3. Solar cycle 23 was not considered due to lack of a close agreement between Rz (or F10.7) and solar EUV [Emmert et al. 2010, Chen et al. 2011, Tapping and Valdes 2011]. Its inclusion is analyzed in the following section. 3. Solar activity effect was filtered out estimating the residuals from the regression between the experimen- tal foF2 values and the solar activity index Rz, that is: foF2res =foF2exp − (a + b Rz) (1) where foF2exp is the experimental foF2 value, Rz is the sunspot number used as a proxy for the solar EUV ra- diation, and a and b are the regression constants as- sessed with least squares. 4. The linear trend a (in MHz/year) is calculated from foF2res = a t + b (2) Figure 1 shows foF2res together with the fitted lin- ear trend. Table 1 shows the values of the trends ob- tained. Trend errors in all the cases are around 8-11% of trend values, that is ~6-9×10-4 Mhz/year. All trend val- ues are significantly different from zero at a 95% confi- VENCHIARUTTI ET AL. 2 Figure 1. foF2 residual for (a) Brisbaen (27.5°S, 152.9°E), (b) Can- berra (35.3°S, 149.0°E) and (c) Christchurch (43.6°S, 172.8°E), after filtering seasonal and solar activity effects using Rz as solar EUV proxy (black dots) and fitted linear trend (solid line). Station Geographic latitude Geographic longitude foF2 trend [MHz/year] Period cycles 18-21 Period cycles 18-23 Rz as proxy F10.7 as proxy Rz as proxy F10.7 as proxy Brisbaen 27.5°S 152.9°E -0.0086 -0.0120 +0.0004* -0.0102 Canberra 35.3°S 149.0°E -0.0074 -0.0107 -0.0008* -0.0089 Christchurch 43.6°S 172.8°E -0.0073 -0.0108 +0.0015 -0.0071 Table 1. Geographic coordinates of the Southern Hemisphere ionospheric stations here analyzed and linear foF2 trend in MHz/year cal- culated for four different conditions: filtering seasonal and solar activity effects using Rz (4th and 6th columns) and F10.7 (5th and 7th columns), and for periods covering solar cycles 18 to 21 (4th and 5th columns) and 18 to 23 (6th and 7th columns) (values with * mean not significantly different from 0). a b c 3 dence level. The significance was assessed calculating the statistic t as a/v, where (3) v is the standard deviation of trend and r the correlation coefficient of the regression between foF2res and time. Rz was chosen as solar EUV proxy instead of F10.7, which is the one recommended [Lastovicka et al. 2012, Elias 2014], since F10.7 is not available since the begin- ning of cycle 18. In addition, since we are not including here cycle 23, the difference between filtering solar ac- tivity using Rz or F10.7 should not be large. As a test, we linearly extrapolated F10.7 values to 1944, using a linear regression between Rz and F10.7 for the period 1947-1986. Trends, listed in Table 1, were assessed fol- lowing the same steps as in Rz case. It can be noticed that although the difference is ~40-50% for the three stations, they are all significant decreasing trends. 2.2. Extending period to cycle 23 The three stations here analyzed have available data records until year 2009. However two of them, Christchurch and Brisbane, present a big data gap dur- ing solar cycle 22, which cannot be easily interpolated for trend detection purposes. As a first approach, with- out interpolation, trends for the period covering cycles 18 to 23 were estimated. Special care should be taken in foF2 filtering process when solar cycle 23 is included [Elias et al. 2014]. The period 2007-2009 that corresponds to minimum activity level, is characterized by lower EUV solar radiation than previous solar cycles [Solomon et al. 2013]. In ad- dition it is different than that deduced from traditional solar EUV proxies as Rz and F10.7. Emmert et al. [2010] suggested that the long-term relationship between EUV irradiance and F10.7 has changed markedly since around 2006 with EUV levels decreasing more than ex- pected from F10.7 proxy. Rz instead, underestimates EUV during the maximum and declining phase of the cycle [Lukianova and Mursula 2011]. In addition, Rz and F10.7, which were used interchangeably as EUV proxy, present a significant change in their relationship since solar cycle 23 [Tapping and Valdes 2011]. Taking all this into account, trends were esti- mated for each station considering Rz and F0.7 to fil- ter solar activity effect. Trend values are listed in Table 1. Figure 2 shows foF2res together with the fitted linear trend, considering both Rz and F10.7 separately as solar activity filters. It can be clearly noticed that the inclusion of solar cycle 23 marks a significant differ- ence not only trend values, but also in the sign, be- tween both EUV proxies. 3. Geomagnetic activity effect Danilov and Mikhailov [1999], Mikhailov and Marin [2000] and Danilov [2002] with a new approach, obtained a dependence of the trend magnitude on ge- omagnetic latitude (with more negative values at higher latitudes) indicating that F2-layer trends might be re- lated to long-term changes in geomagnetic activity and t n 2 r r 1 r n 2 1 12 2v a a a = = - - = - - IONOSPHERIC TRENDS IN SOUTHERN HEMISPHERE Figure 2. foF2 residual for (a) Brisbaen (27.5°S, 152.9°E), (b) Can- berra (35.3°S, 149.0°E) and (c) Christchurch (43.6°S, 172.8°E), after filtering seasonal and solar activity effects using Rz (black dots) and F10.7 (red dots) as solar EUV proxy and fitted linear trend (corre- sponding solid lines). a b c F2-layer storm mechanisms. They show those periods with negative and positive foF2 trends exist, which cor- respond to the periods of long-term increasing or de- creasing geomagnetic activity. However, more recently, according to Danilov [2012], the assumptions made at the end of 1990’s that the observed trends in ionospheric parameters are a manifestation of secular changes in ge- omagnetic activity could be proved only for the data be- fore 1980-1990. In our case, the three stations here analyzed are mid-latitude stations, so the geomagnetic activity effect should not be large. To prove this, aa index was added to Equation (1) to filter foF2, that is foF2res =foF2exp − (a + b Rz + c aa) (4) This index is a global geomagnetic activity index which measures the disturbance level of the Earth’s magnetic field. The index, in nT units, is available from http://www.wdcb.ru/stp/data/geomagni.ind/aa/ and from http://isgi.cetp.ipsl.fr/source/indices/. Following the same steps described in Section 2.2, trends were assessed considering period 1944-1986, and 1944-2009, using Rz and F10.7 as solar EUV proxies sep- arately. Trend values are shown in Table 2. As can be noticed, in all the cases trend values are almost the same as those seen in Table 1. This can be expected since the correlation coefficient between foF2 and Rz or F10.7, is already greater than 0.97 meaning that Rz (or F10.7) explain more than 94% of foF2 variance in all the cases. Including aa into the linear regression model, only increases the correlation by ~0.005, that is it adds only 1% to the explained variance. Taking into account the result of Danilov [2012] that there is a different foF2 negative trend after ~1980- 1990, and of Lastovicka et al. [2012] that observes a change of dependency of F2 trends after 2000, we made a rough analysis in the case of Canberra, which presents complete data set covering the period 1944- 2009. We estimated the aa coefficient, c in Equation (4), considering F10.7 as EUV proxy, for the period 1944- 1989 and 1990-2009 separately. In the first case it results −0.15±0.03 and in the latter −0.04±0.03, while b results 0.395±0.003 and 0.382±0.003. That is, after 1990 the co- efficient of aa is almost null. This could be evidence of a decline in importance of aa, however when we look at trends estimated with and without aa in foF2 filter- ing, we cannot arrive to a definitive conclusion through this type of analysis. Anyway as it has been observed by Bremer [1992, 1998] and Mikhailov and Marin [2000] the inclusion of a geomagnetic index in the regression in fact does not remove the geomagnetic dependence. 4. Earth’s magnetic field effect The possibility of ionospheric trends induced by the Earth’s magnetic field secular variations was first suggested by Foppiano et al. [1999], and followed then by other papers [Elias 2009, Cnossen 2014]. Earth’s magnetic field, generated in the Earth’s core, presents long term variations in the field’s strength and orienta- tion. A simple mechanism through which trends in the Earth’s magnetic field would affect the ionosphere is through changes in the dip angle (I). The sin(I)cos(I) fac- tor, associated to effects of neutral winds on hmF2 [Rishbeth 1998] will also change. The horizontal ther- mospheric wind U drives ions and electrons, up during the night and down during the day, along the geomag- netic field lines at speed U cos(I). The vertical compo- nent U sin(I) cos(I) raises the F2-peak during night time (when U blows from Pole to Equator) and lowers it dur- ing daytime (when U blows from Equator to Pole), in- creasing or decreasing the peak electron density. An increase in the sin(I)cos(I) factor would produce an ad- ditional lowering of the F-region with a decrease in foF2, during daytime, and an additional raise of the re- gion with an increase in foF2 during the night. A de- crease in the sin(I) cos(I) factor would produce the opposite effect. The region of strongest variations of foF2, lies be- tween 10°N and 30°S in latitude and between 20°E and 80°W in longitude, which is also the region of strongest changes in I and sin(I) cos(I) factor [Elias 2009, Cnossen VENCHIARUTTI ET AL. 4 Station foF2 trend [MHz/year] Period cycles 18-21 Period cycles 18-23 Rz + aa F10.7 + aa Rz + aa F10.7 + aa Brisbaen -0.0083 -0.0118 +0.0001* -0.0101 Canberra -0.0068 -0.0105 -0.0007* -0.0089 Christchurch -0.0068 -0.0103 +0.0016 -0.0070 Table 2. foF2 trend in MHz/year calculated for four different conditions: filtering seasonal and solar activity effects using Rz and aa (4th and 6th columns) and F10.7 and aa (5th and 7th columns), and for periods covering solar cycles 18 to 21 (4th and 5th columns) and 18 to 23 (6th and 7th columns) (values with * mean not significantly different from 0). 5 2014]. At the mid-latitude zone, corresponding to the locations of the ionospheric stations here analyzed, the expected trends due to changes in I, should not be sta- tistically significant different from zero. A qualitative analysis can be made through the sin(I) cos(I) factor. Figure 3 presents the time variation of this factor calculated from the International Geo- magnetic Reference Field (IGRF) model (http://ccmc. gsfc.nasa.gov/modelweb/models/igrf_vitmo.php), and for the three stations here analyzed, it can be no- ticed an increase for the period 1944-1986. This implies that although it is very small, an in- creasing trend in daytime foF2, is contrary to the val- ues obtained from measured foF2 data. 5. Conclusions All stations present a decreasing trend in foF2 dur- ing the period 1944-1986, covering complete solar cycles from cycle 18 to 21, consistent with the consequences ex- pected in the ionosphere F2 region due to the increased concentration of greenhouse gases. For a doubling of CO2 concentration Qian et al. [2009] obtain a maximum decrease of 18% to 20% in foF2. Making a linear extrapolation to the actual increase in CO2 during the period analyzed in this paper, and as- suming an average foF2 value of 10 MHz, a foF2 decrease of −0.005 MHz/year should be expected. According to the values here obtained when Rz is used as solar EUV proxy (~0.008 MHz/year), the greenhouse effect would account for 60% of these trends. When F10.7 is used as EUV proxy, extrapolating to cover the complete period 1944-1986, this account diminish to 50%. With regard to natural forcings, in the case of long- term variation in geomagnetic activity, our results hardly vary when considering the aa index as another linear term in Equation (1). Probably in a different kind of analysis, and considering the result of Danilov [2012] this could be more noticeable. In the case of secular variations of the intrinsic magnetic field of the Earth, its ability to induce trends depends on the geographic location. The values of the expected trends are very small (even lower than those expected by anthropogenic effect) except in the region between 10°N and 40°S and 60°W and 10°E [Elias 2009, Cnossen 2014] . The stations considered are located outside this region. In addition a rough qualitative estimation of the expected effect does not even agree in the sign of the expected trends. When using the complete data sets that go until solar cycle 23, the most reliable trends are those ob- tained using F10.7 as EUV proxy [Lastovicka et al. 2012, Elias 2014]. Trend values are close to those correspon- ding to period 1944-1986, but less negative, very likely due to a lack of agreement between F10.7 and the solar EUV radiation which F10.7 is supposed to represent. The importance of understanding and measuring trends in the upper atmosphere lies in: 1) the importance of understanding and predict- ing the behavior of the whole atmosphere, which is es- sential to human life. 2) the effect of the decreasing density on space de- bris. There are approximately 22,500 objects larger than 10 cm in Earth orbit, and 95% of them is classified as “space debris”. These objects represent a significant risk to spacecraft operations due to possible collisions. Re- cent studies [Lewis et al. 2011] suggest that the decrease in the density of the upper atmosphere makes the space debris remains in orbit longer than expected. The results obtained here are intended to be a con- tribution to both topics. References Bremer, J. (1992). Ionospheric trends in mid-latitudes as a possible indicator of the atmospheric greenhouse effect, J. Atmos. Terr. Phys., 54, 1505-1511. Bremer, J. (1998). Trends in the ionospheric E and F re- gions over Europe, Ann. Geophysicae, 16, 986-996. Chen, Y., L. Liu and W. Wan (2011). Does the F10.7 index correctly describe solar EUV flux during the deep solar minimum of 2007-2009?, J. Geophys. Res., 116; doi:10.1029/2010JA016301. Cnossen, I. (2014). The importance of geomagnetic field changes versus rising CO2 levels for long-term change in the upper atmosphere, J. Space Weather Space Clim., 4, A18. Damboldt, T., and P. Suessmann (2012). Consolidated Database of Worldwide Measured Monthly Medians of Ionospheric Characteristics foF2 and M(3000)F2, INAG Bulletin on the Web, INAG-73; www.ips.gov. au/IPSHosted/INAG/web-73/2012/damboldt_con solidated_database.pdf. IONOSPHERIC TRENDS IN SOUTHERN HEMISPHERE Figure 3. Standardized sin(I) cos(I) factor estimated with International Geomagnetic Reference Field (IGRF) model (http://ccmc.gsfc.nasa. gov/modelweb/models/igrf_vitmo.php) for Brisbaen (black solid line), Canberra (red solid line) and Christchurch (dashed black line). Danilov, A.D., and A.V. Mikhailov (1999). Long-term trends in the parameters of the F2- region: A new approach, Geom. Aeron., 39, 473-479. Danilov, AD. (2002). Overview of the trends in the ionospheric E and F2 regions, Phys. Chem. Earth, 27, 579-588. Danilov, A.D. (2012). Long-Term Trends in the Upper Atmosphere and Ionosphere (a Review), Geom. Aeron., 52, 271-291. Elias, A.G. (2009). Trends in the F2 ionospheric layer due to long-term variations in the Earth’s magnetic field, J. Atmos. Solar Terr. Phys., 71, 1602-1609. Elias, A.G. (2014). Filtering ionosphere parameters to detect trends linked to anthropogenic effects, Earth, Planets and Space, 66, 113. Elias, A.G., B.F. de Haro Barbas, S. Kiyoto and J.R Souza (2014). Effect of solar cycle 23 in foF2 trend estima- tion, Earth, Planets and Space, 66, 111. Emmert, J.T., J.L. Lean and J.M. Picone (2010). Record- low thermospheric density during the 2008 solar minimum, Geophys. Res. Lett., 37; doi:10.1029/2010 GL043671. Foppiano, A.J., L. Cid and V. Jara (1999). Ionospheric long-term trends for South American mid-latitudes, J. Atmos. Solar Terr. Phys., 61, 717-723. Lastovicka, J., S.C. Solomon and L. Qian (2012). Trends in the Neutral and Ionized Upper Atmosphere, Space Sci. Rev., 168, 113-145. Lewis, H.G., A. Saunders, G. Swinerd and R.J. Newland (2011). Effect of thermospheric contraction on re- mediation of the near Earth space debris environ- ment, J. Geophys. Res., 116. Lukianova, R., and K. Murusla (2011). Changed rela- tion between sunspot numbers, solar UV/EUV ra- diation and TSI during the declining phase of solar cycle 23, J. Atmos. Solar Terr. Phys., 73, 235-240. Mikhailov, A.V., and D. Marin (2000). Geomagnetic con- trol of the foF2 long-term trends, Ann. Geophysi- cae, 18, 653-665. Qian, L., A.G. Burns, S.C. Solomon and R.G. Roble (2009). The effect of carbon dioxide cooling on trends in the F2-layer ionosphere, J. Atmos. Solar Terr. Phys., 71, 1592-1601. Rishbeth, H. (1990). A greenhouse effect in the iono- sphere?, Planet. Space Sci., 38, 945-948. Rishbeth H. (1998). How the thermospheric circulation affects the ionospheric F2 layer, J. Atmos. Solar Terr. Phys., 60, 1385-1402. Roble, R.G., and R.E. Dickinson (1989). How will changes in carbon dioxide and methane modify the mean structure of the mesosphere and thermos- phere?, Geophys. Res. Lett., 16, 1441-1444. Solomon, S.C., L. Qian and A.G. Burns (2013). The anomalous ionosphere between solar cycles 23 and 24, J. Geophys. Res., 118, 6524-6535. Tapping, K.F., and J.J. Valdes (2011). Did the Sun change its behavior during the decline of cycle 23 and into cycle 24?, Sol. Phys., 272, 337-350. Corresponding author: Blas F. de Haro, Universidad Nacional de Tucumán, Departamento de Física, Tucumán, Argentina; email: blasdeharo2000@yahoo.com.ar. © 2015 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. VENCHIARUTTI ET AL. 6 << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles false /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Warning /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.1000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams true /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments false /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile (None) /AlwaysEmbed [ true /AndaleMono /Apple-Chancery /Arial-Black /Arial-BoldItalicMT /Arial-BoldMT /Arial-ItalicMT /ArialMT /CapitalsRegular /Charcoal /Chicago /ComicSansMS /ComicSansMS-Bold /Courier /Courier-Bold /CourierNewPS-BoldItalicMT /CourierNewPS-BoldMT /CourierNewPS-ItalicMT /CourierNewPSMT /GadgetRegular /Geneva /Georgia /Georgia-Bold /Georgia-BoldItalic /Georgia-Italic /Helvetica /Helvetica-Bold /HelveticaInserat-Roman /HoeflerText-Black /HoeflerText-BlackItalic /HoeflerText-Italic /HoeflerText-Ornaments /HoeflerText-Regular /Impact /Monaco /NewYork /Palatino-Bold /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman /SandRegular /Skia-Regular /Symbol /TechnoRegular /TextileRegular /Times-Bold /Times-BoldItalic /Times-Italic /Times-Roman /TimesNewRomanPS-BoldItalicMT /TimesNewRomanPS-BoldMT /TimesNewRomanPS-ItalicMT /TimesNewRomanPSMT /Trebuchet-BoldItalic /TrebuchetMS /TrebuchetMS-Bold /TrebuchetMS-Italic /Verdana /Verdana-Bold /Verdana-BoldItalic /Verdana-Italic /Webdings ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 150 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.10000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 150 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.10000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.08250 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName (http://www.color.org) /PDFXTrapped /Unknown /CreateJDFFile false /SyntheticBoldness 1.000000 /Description << /ENU (Use these settings to create PDF documents with higher image resolution for high quality pre-press printing. The PDF documents can be opened with Acrobat and Reader 5.0 and later. These settings require font embedding.) /JPN /FRA /DEU /PTB /DAN /NLD /ESP /SUO /NOR /SVE /KOR /CHS /CHT /ITA >> >> setdistillerparams << /HWResolution [2400 2400] /PageSize [595.000 842.000] >> setpagedevice