HONGWEI_corrected_Layout 6 A signal-strength model for gravity-field spectrum components in a local space Liu Hongwei1,*, Wang Zhaokui2, Zhang Yulin1,2 1 Institute of Aerospace Science and Engineering, National University of Defence Technology, Changsha, China 2 School of Aerospace, Tsinghua University, Beijing, China ANNALS OF GEOPHYSICS, 55, 5, 2012; doi: 10.4401/ag-5459 ABSTRACT In the Earth gravity-field model, potential coefficients are projections of the Earth gravitational potential on the corresponding spherical harmonics by degree and order. These coefficients reflect the spectral composition of the Earth gravitational potential in the whole space, but they cannot reflect the spectral composition in a local area. In this study, the Earth gravitational potential was projected on spherical harmonics in a local area using a window function introduced to facilitate the mathematical expressions. Then a model for the spectral structure of the Earth gravitational potential in the local area was established. The model can reflect the signal strength of any gravity-field spectrum component in the local area, which adds information for the description of the gravity field of the Earth, and which has great significance for Earth-sciences research and satellite gravity measurements. Taking the third-degree coefficients, for example, the signal-strength distributions of potential coefficients were computed. The data show that the signal-strength distributions of third- degree coefficients on the sphere surface are independent of the longitude. Among the third-degree coefficients, the zero-order, first-order and second- order coefficients are stronger near the two poles of the Earth, while the third-order coefficient is stronger near the equator. 1. Introduction The Earth gravity field is one of basic physical fields of the Earth, and it reflects information relating to matter dis- tribution inside the Earth and affects physical events on the Earth and in the neighboring space. Therefore, research into the Earth gravity field is always a hot issue in geodesy [Jin- sheng 2001, 2002]. The gravity field outside the Earth is caused by matter inside the Earth. As well as the shape of the Earth affecting the Earth gravity potential, this is also affected by the asymmetry of the matter density distribution inside the Earth [Jian 1994]. Researches into the Earth shape and density are two basic missions in Earth gravity research [Xiaoguang et al. 2000]. Due to anomalies in the shape of the Earth and to the asymmetry of the matter density distribution inside the Earth, the Earth gravity potential cannot be expressed as a simple an- alytical form. The Earth gravity field is usually approximated by a spherical harmonics series, and the set of all the gravita- tional potential coefficients in an infinite series expansion of spherical harmonics represents the Earth gravity field model. The idea that the Earth gravity potential can be ex- pressed as an infinite series of spherical harmonics has been widely used in geoscience research, and it is a basic method in Earth sciences and gravity satellite research. The literature contains conclusions about the relationships between dis- turbing potential coefficients and matter density anomalies inside the Earth [Jian 1994]. It was believed that lower degree potential coefficients reflect deep Earth matter density anomalies, while higher degree potential coefficients reflect shallow Earth matter density anomalies. By using the GRACE time-varying gravity-field model, global water stor- age changes were analyzed in Xuhua et al. [2006], Antarctic ice-sheet mass changes were studied in Guangbin et al. [2009], and application of a gravity-field model in tectonic detection was analyzed in Leyang et al. [2007]. At the same time, the spherical harmonics infinite series expansion of the Earth gravitational potential has been widely used in satellite gravity measurements. Satellites can measure the Earth grav- ity field, such as CHAMP, GRACE and GOCE, which were designed for recovering gravity potential coefficients in spe- cial bands. Among these three satellites, CHAMP has meas- ured the Earth gravity field using orbit perturbation theory, which has helped to detect lower degree potential coeffi- cients [Reigber et al. 2003]. GRACE measured the Earth grav- ity field by satellite-to-satellite tracking in a low Earth orbit, which helped in the detection of the lower- and middle-de- gree potential coefficients [Reigber et al. 2005, Zhao et al. 2011]. GOCE has used gravity gradient measurements that have helped to detect middle- and higher-degree potential coefficients [Muzi and Allasio 2003]. In the Earth gravity field model, gravitational potential coefficients are projections of gravitational potential on the Article history Received October 28, 2011; accepted June 11, 2012. Subject classification: Earth gravitational potential, Local area, Spectrum structure, Spherical cap harmonic analysis, Signal strength. 1017 RESEARCH ARTICLES corresponding spherical harmonics by degree and order in the whole space, and they reflect the spectral composition of the gravitational potential in the whole space. However, gravitational potential coefficients cannot reflect the spectral composition of gravitational potential in a local space, which is a shortcoming of conventional Earth gravity-field models. Although spherical cap harmonic analysis (SCHA) can give the spectral composition of gravitational potential on a spherical cap [De Santis et al. 1999], because of less gravity potential coefficients, the spectral structure of the spherical cap area reflected by a SCHA model is limited. In the present study, based on the results in Slepian and Pollak [1961], the Earth gravitational potential was projected onto spherical harmonic functions in a local area, and by appropriate math- ematical expressions, the signal-strength model for the grav- ity potential spectral component in a local space was established. This model will be helpful for Earth-science re- search and gravity satellite measurements. 2. A signal-strength model for gravity-field spectrum components in a local space We establish a spherical coordinate system with the cen- ter of the Earth as the origin, and we take the z-axis and the Earth rotation axis in the same direction, as is shown in Fig- ure 1. The spherical harmonic expansion of Earth's gravita- tional potential is: (1) where the definition of the spherical coordinates (t,i,m) is shown in Figure 1, and G is the gravitational constant, M is the Earth mass, a is the semi-major axis of the Earth ellipsoid, are the gravitational potential coeff icients, and is the fully normalized associated Legendre poly- nomial. The relationship between and the associ- ated Legendre polynomial, Pnk (cosi), is: (2) The nonspherical Earth perturbation potential function can be given from Equation (1): (3) By the orthogonality of spherical harmonics functions, the gravitational potential coefficients can be de- rived as: (4) (5) where X, the whole spherical surface, is the integral area. X can be divided into n parts, (6) where, Xi (i=1,2,...,n) are the local areas on the spherical sur- face, as shown in Figure 2. Xi can be described by a certain longitude and latitude range, and its center is of interest. In Equations (4) and (5), if the integral area X is changed into Xi, then projections of the nonspherical Earth pertur- bation potential R(t,i,m) on n-degree and k-monics in the local area Xi can be obtained: (7) (8) Then, (9) It is known that reflect projections of the nonspherical Earth perturbation potential on n-degree and k-order spherical harmonics {Pnk (cos i)cos km, Pnk (cos i)sin km} in the local areas Xi, which represent the signal strengths of n-degree and k-order spherical harmonic components in Xi, while reflect the sum of these projections in the LIU HONGWEI ET AL. 1018 , , , V cos sin cosGM a C k S k P1 nn n k2 0 t i m t t m m i = + + 3 = = nk nk nk ^ ^ ^ h h h; E` j// ! , , ! , . cos cos P n k k n n k P 2 0 1 2 1 2 k k nk i d i d + = = = + - nk k 0! ^ ^ ^ ^ ^ h h h h h ' , , sin .cos cos R C S k PGM a k k n n n 02 t i m t t m m i = + 3 == nknk nk ^ ^ ^ h h h` j// 2 ! ! , , sin cos d d ,kcos C GM a n n k n k R P 2 1 2 1n k nk m m r t t d t i m i i i = + + - X nk ^ ^ ^ ^ h h h h ` j ## 2 ! ! , , sin cos d d ,kcos C GM a n n k n k R P 2 1 2 1 ,nk i n k nk i m m r t t d t i m i i i = + + - X ^ ^ ^ ^ h h h h ` j ## 2 ! ! , , sin d d ,sinkcos S GM a n k n n k R P 2 1 2 1n nk m m r t t t i m i i i = + + - X nk ^ ^ ^ ^ ^ h h h h h ` j ## 2 ! ! , , sin d d .sinkcos S GM a n k n n k R P 2 1 2 1 ,nk i n nk i m m r t t t i m i i i = + + - X ^ ^ ^ ^ ^ h h h h h ` j ## ... ,n1 2X X X X= " , , .C C S Snk i n nk i n 1 1 = = = = , ,nk i nk i/ / Figure 1. Spherical coordinate system with the Earth center as its origin. C , Snk nk" , C , Snk nk" , C , Snk nk" , cosP ink^ h cosP ink^ h C , Snk, nk,i i" , 1019 whole area X, and represent the total signal strength of n-degree and k-order spherical harmonic components in the whole area. For purposes of satellite gravity field measurements, if the re- covery frequency bands involved contain n-degree potential co- efficient, sampling data at the places where the signal of the n-degree spectral component is stronger helps to increase the signal-to-noise ratio of n-degree spectrum components and in- crease the recovery accuracy of n-degree potential coefficients. The signal strength of n-degree and k-order spectrum component in the local areas Xi can be computed by Equa- tions (7) and (8), which is approximately equivalent to the computing of the data with perturbation potentials in Equa- tions (4) and (5), multiplied by a window function. The win- dow function takes the maximum in Xi, and the farther the distance away from Xi it is, the smaller the window function will be. Here, on the sphere surface of fixed radius t0, the signal strength of n-degree and k-order spectrum compo- nents in the local area X0 will be computed, where the cen- ter of X0 is (t0,i0,m0). The unit vector from the sphere center to the point (t0,i0,m0) is: (10) The unit vector from the sphere center to any point, (t0,i,m), on the sphere surface is: (11) Assuming that the angle between two unit vectors in the above is }, then: (12) Thus, (13) On the sphere surface of radius t0, making a window function with the point (i0,m0) as its center, (14) As } is the spherical distance from (t0,i,m) to (t0,i0,m0), the window function satisfies 00. The greater p is, the narrower H is, and the smaller p is, the wider H is. Multiplying Equations (4) and (5) by the window func- tion in Equation (14) and substituting Equation (3) gets the sine and cosine signal strength of n-degree and k-order spec- trum components in the local area X0 where the center is (t0,i0,m0), (15) (16) Then the signal strength of n-degree and k-order spec- trum components in the local area X0 is: (17) 3. Comparison among spherical harmonic analysis, spher- ical cap harmonic analysis, and translated origin spherical cap harmonic analysis From spherical harmonics expansion, it is known that the gravitational potential can be seen as the superposition of an infinite number of spherical harmonic waves, and po- tential coefficients are just the vibration amplitudes of these spherical harmonic waves, which reflect the energy of cor- responding spherical harmonic waves. By projecting gravi- tational potential on spherical harmonic functions in the global area, i.e., using spherical harmonic analysis (SHA), a conventional gravity model can be obtained, and one poten- tial coefficient reflects the total energy of the corresponding spherical harmonic in the global area. For the projection of gravitational potential on spherical harmonic functions in a local area, i.e., the method investigated here, we get the en- ergy of the spherical harmonics in the local area. This is the relationship between SHA and our method. Spherical cap harmonic analysis (SCHA) is a method for the gravitational potential on the spherical cap, which was first proposed by Haines, and the effectiveness of this method has been shown by geomagnetic field modeling [Haines 1985]. SCHA was first applied in geodesy by Jiancheng, with an SCHA model for a local gravity field given GRAVITY FIELD SPECTRUM IN LOCAL AREA .sin cos sin sin cose 0 0 0 0 0i m i m i=0 ^ h .sin cos sin sin cose i m i m i=^ h .cos sin sin cos cos cose e 0 00$} i i m m i i= = - +0 ^ h arc cos cos cos .sin sin cos0 0 0} i i m m i i= - +^ h6 @ , .H e,0 0 2 i m =i m p}-^ h ) 2 ! ! , sin cos d dk cos sin cos cos C a n n k n k a H C k S k P P 2 1 2 1 , n k m l l l lm lm lm nk 02 0 0 m m r t d t i m m m i i i i = + + - + 3 i m X == nk ^ ^ ^ ^ ^ ^ h h h h h h ` `j j// ## ) 2 ! ! , sin d d .k cos sin cos cos sin S a n k n n k a H C k S k P P 2 1 2 1 , n m l l l lm lm lm nk 02 0 0 m m r t t i m m m i i i i = + + - + 3 i m X == nk ^ ^ ^ ^ ^ ^ ^ h h h h h h h ` `j j// ## ) ) .E C Snk 2 2= +nk nk^ ^h h Figure 2. Sketch map of the local region Xi on the spherical surface. [Jiancheng et al. 1995, De Santis and Torta 1997]. The basic functions in SCHA are composed of noninteger degree as- sociated Legendre polynomials and trigonometric functions, with orthogonality and completeness, and they can approx- imate the Earth gravitational potential function on the spher- ical cap. The nonspherical perturbation potential on the spherical cap can be expressed as: (18) where Pnkm (cos i) is the noninteger degree associated Le- gendre polynomial. From a mathematical point of view, SCHA is the method that the potential on a spherical cap is projected on spherical cap harmonic functions. Compared to the SHA model, when the same spatial resolution for a local gravity field is achieved, fewer gravitational potential coeffi- cients are required in the SCHA model [Fuqing and Jinhai 2000], which means that the SCHA model contains less spec- trum components than SHA model. At the same time, with SCHA, the nonorthogonality of odd basis functions (k-m = odd) and even basis functions (k-m = even) can lead to com- putation instability of potential coefficients and deterioration in the SCHA model at the edge of the cap when the data dis- tribution is not uniform over the region [De Santis 1991]. For the computation instability of potential coefficients in SCHA, De Santis [1991] proposed translated origin spher- ical cap harmonic analysis (TOSCA), which was mathemati- cally constructed by moving the new origin in the direction of the cap region and applying SCHA in the new reference sys- tem. TOSCA can greatly improve the detail in the central re- gions and can smooth the observations at the edge of the cap. In addition, for the computation of more complex Legendre functions with non integer harmonic degrees in SCHA, De Santis proposed adjusted SHA (ASHA), in which the co-lati- tude interval was adjusted to a hemisphere, and then con- ventional spherical harmonics could be used [De Santis 1992]. Just like SCHA and TOSCA, the method in the present study can be used for analysis of a potential spectrum in a local area. However, in our method, the gravitational poten- tial in a local area is directly projected on the spherical harmonic functions, not on the spherical cap harmonic functions, which compared to SCHA, can obtain more potential coefficients and consequently a more comprehensive spectrum structure for the local area, although the computation cost will increase. To fa- cilitate the mathematical expressions, the product of gravita- tional potential and an exponential function called the window function, which decreases with distance from the center of the local area involved, is projected on the spherical harmonic func- tions in the global region, which indirectly obtains the projec- tion of the gravitational potential in the local region; i.e., Equations (15) and (16). It is in the global region for the projec- tion of gravitational potential that the observation data scope increases and the detail in the central regions is greatly im- proved, which is similar to De Santis 1991, and helps to reduce errors. In De Santis [1991], the details in the central regions were greatly improved by using more observation data under a new coordinate system that was obtained by moving the origin of the old coordinates upwards from the center of the Earth. The relationships between SHA, SCHA, TOSCA and ASHA and the method in the present study are shown in Figure 3. LIU HONGWEI ET AL. 1020 Figure 3. The relationships between SHA, SCHA, TOSCA, ASHA and the method in the present study. , , ,cos sin cos R mGM a C m S P m k k n km km n m 02 k k t i m t t m m i = + 3 X == ^ ^ ^ h h h` j// 1021 4 Numerical computations Take the first 10 degrees in the Earth Gravitational Model 1996 (EGM96) as a simulation of the Earth gravity field. The attenuation factor p determines the size of local area selected, and a right attenuation factor can be obtained by computa- tional results. According to Equations (15), (16) and (17), in the case of p = 0~10, the computed signal strength of the third- degree potential coefficient along the longitude m = r/4 at the altitude of 200 km is shown in Figures 4-7. From Figures 4-7, it is known that: (1) When attenuation factor p = 0, the signal strength of each spherical harmonic component does not change with the co-latitudes. The reason for this is that p = 0 means the local area is just the entire spherical surface, and their projections on spherical harmonics are certainly the same, making the signal strength of each spherical harmonic com- ponent equal in any area of the spherical surface. (2) With the increase in the attenuation factor, local fea- tures gradually appear in the signal strength of each spheri- cal harmonic component. With different attenuation factors, the sizes of selected local areas are different, and conse- quently there will be some differences among the signal strengths of each spherical harmonic component, although their tendencies with the co-latitudes are consistent. GRAVITY FIELD SPECTRUM IN LOCAL AREA Figure 4. The signal strength of third-degree and zero-order spherical harmonics at different co-latitudes (longitude, m = r/4). Figure 5. The signal strength of third-degree and first-order spherical harmonics at different co-latitudes (longitude, m = r/4). According to Figures 4-7, when p = 5, the change ten- dency for the signal strength of each spherical harmonic component with co-latitudes can be shown clearly. Take the attenuation factor of the window function as 5, as shown in Figure 8, and analyze the distribution law for the signal strength of the third-degree spherical harmonic component along the whole sphere surface. From Equations (15), (16) and (17), the signal strength dis- tributions of the third-degree potential coefficients on the sphere surface at the altitude of 200 km are shown in Figures 9 to 12. Figures 9 to 12 show that the signal strength distributions of the third-degree potential coefficients on the sphere surface are independent of the longitude. The zero-order, first-order and second-order components are stronger near the two poles of the Earth, and weaker near the equator, while the third-order is weaker near the two poles and stronger near the equator. From Equation (3), it is known that when the longitude m is unchanged, the associated Legendre polynomial, Pnk (cosi), reflects the change tendency of the n-degree and k-order spherical harmonic with the co-latitude. This represents the signal strength change tendency of n-degree spherical har- monics in the conventional gravity model, such as the EGM96 LIU HONGWEI ET AL. 1022 Figure 6. The signal strength of third-degree and second-order spherical harmonics at different co-latitudes (longitude, m = r/4). Figure 7. The signal strength of third-degree and third-order spherical harmonics at different co-latitudes (longitude, m = r/4). 1023 model. Figure 13 shows the changes of P30 (cosi), P31 (cosi), P32 (cosi) and P33 (cosi) with the co-latitudes, and Figure 14 shows the signal strengths of third-degree spherical harmonic components along the longitude m= r/4, with attenuation fac- tor p= 5. From the comparison of Figures 13 and 14, it is known that for the change tendency of first-order and third-order spherical harmonic components with the co-latitudes, the re- sults given by the EGM96 model and the method in the present study are the same, while for the zero-order and third-order, the results from the EGM96 model and the present method are dif- ferent. This shows that the method in the present study can give a more sophisticated structure for the gravitational potential, and increases the information describing the Earth gravity field, which cannot be seen from the conventional gravity model. 5. Conclusions Here, by projecting the Earth gravitational potential on the spherical harmonics in a local area, the signal strength model for each spectral component in the Earth gravity field model in a local area was established. Compared to the tradi- tional Earth gravity-field model, i.e., SHA model, this can pro- vide a more sophisticated gravity field spectrum structure in GRAVITY FIELD SPECTRUM IN LOCAL AREA Figure 8. The window function with an attenuation factor p = 5 on the Earth surface. Figure 9. The distribution of on the sphere surface.C S2 2+) ) 30 30 the local area and increase the information that describes the Earth gravity field, which cannot be seen from the conven- tional gravity model. Compared to SCHA, the present model can reflect more spectral components in the local area, which helps to design a more reasonable satellite orbit and sampling area for the given gravity field measurement band, and thus we can strengthen the gravity satellite sample data at the places where the signal of the potential coefficients is stronger, and improves the accuracy of the gravity field measurement. In the end, the signal strength distributions for third-de- gree coefficients were computed. These results show that the signal strength distributions of third-degree coefficients are independent of the longitude, and the zero-order, first-order and second-order coefficients have stronger signals nearer the two poles of the Earth, while the third-order coefficient has a stronger signal near the equator. LIU HONGWEI ET AL. 1024 Figure 10. The distribution of on the sphere surface. Figure 11. The distribution of on the sphere surface. C S2 2+) ) 31 31 C S2 2+) ) 32 32 1025 Acknowledgements. This study is supported by the National Natural Science Foundation of China (Grant No. 11002076) and National Defence Pre-Research (Grant No. 51320010201). The author would like to thank Zou Xiancai, Associate Professor in the School of Surveying and Mapping, Wuhan University, China, for useful discussions. The author gratefully acknowledges helpful suggestions and comments from Prof. Angelo De Santis, Istituto Nazionale di Geofisica e Vulcanologia (INGV), Rome, Italy. References De Santis, A. (1991). Translated origin spherical cap harmonic analysis, Geophysical Journal International, 106, 253-263 De Santis, A. (1992). Conventional spherical harmonic analy- sis for regional modeling of the geomagnetic field, Geo- physical Research Letters, 19 (10), 1065-1067. De Santis, A., and J.M. Torta (1997). Spherical cap harmonic analysis: a comment on its proper use for local gravity field representation, Journal of Geodesy, 71, 526-532. De Santis, A., J.M. Torta and F.J. Lowes (1999). Spherical cap GRAVITY FIELD SPECTRUM IN LOCAL AREA Figure 12. The distribution of on the sphere surface. Figure 13. The associated Legendre polynomials P30(cosi), P31(cosi), P32(cosi) and P33(cosi). C S2 2+) ) 33 33 harmonics revisited and their relationship to ordinary spher- ical harmonics, Phys. Chem. Earth (A), 24 (11-12), 935-941. Fuqing, P., and Y. Jinhai (2000). The characters and compu- tation of Legendre function with non-integral degree in SCHA, Acta Geodaetica et Cartographica Sinica, 29 (3), 204-208 (in Chinese). Guangbin, Z., L. Jiancheng, W. Hanjiang, W. Zhengtao (2009). Investigation on the mass variations of ice sheet in Antarctic with GRACE time-variable gravity models, Ge- omatics and Information Science of Wuhan University, 34 (10), 1185-1189 (in Chinese). Haines, G.V. (1985). Spherical cap harmonic analysis, J Geo- phys Res, 90, 2583-2591. Jian, F. (1994). Investigation of anomaly density by using satellite gravity data, Progress in Geophysics, 3, 60-65 (in Chinese). Jiancheng, L., C. Dingbo and N. Jinsheng (1995). Spherical cap harmonic expansion for local gravity field represen- tation, Manuscr. Geod., 20, 265-277. Jinsheng, N. (2001). Following the developments of the world, devoting to the study on the Earth gravity field, Geomatics and Information Science of Wuhan Univer- sity, 26 (6), 471-474 (in Chinese). Jinsheng, N. (2002). The satellite gravity gravity surveying technology and research of Earth's gravity field, Journal of geodesy and geodynamics, 22 (1), 1-5 (in Chinese). Leyang, W., Z. Jianjun and Z. Ziqiang (2007). Applicable study of EIGEN-GRACE02S gravity model on detection of Earth's crust tectonics, Journal of geodesy and geody- namics, 27, 31-34 (in Chinese). Muzi, D., and A. Allasio (2003). GOCE: the first core Earth explorer of ESA's Earth observation programme, Acta Astronautica, 54, 167-175. Reigber, Ch., P. Schwintzer, K.H. Neumayer, F. Barthelmes, R. König, Ch. Förste, G. Balmino, R. Biancale, J.M. Lemoine, S. Loyer, S. Bruinsma, F. Perosanz and T. Fa- yard (2003). The CHAMP-only Earth gravity field model EIGEN-2, Advances in Space Research, 31 (8), 1883-1888. Reigber, Ch., R. Schmidt, F. Flechtner, R. König, U. Meyer, K.-H. Neumayer, P. Schwintzer and S.Y. Zhu (2005). An Earth gravity field model complete to degree and order 150 from GRACE: EIGEN-GRACE02S, Journal of Geo- dynamics, 39, 1-10. Slepian, D., and H.O. Pollak (1961). Prolate spheroidal wave functions, Fourier analysis and uncertainty — I, The Bell System Technical Journal, 40, 43-64. Xiaoguang, H., X. Houze and L. Dajie (2000). Earth's density flat rate and the hypothesis of latitudinal normal density, Science in China (Series D), 30 (4), 436-441 (in Chinese). Xuhua, Z., W. Bin, P. Bibo, X. Houze (2006). Detection of global water storage variation using GRACE, Chinese J . Geophys., 49 (6), 1644-1650 (in Chinese). Zhao, Q., J. Guo, Z. Hu, Ch. Shi, J. Liu, H. Cai and X. Liu (2011). GRACE gravity field modeling with an investigation on correlation between nuisance parameters and gravity field coefficients, Advances in Space Research, 47, 1833-1850. *Corresponding author: Liu Hongwei, Institute of Aerospace Science and Engineering, National University of Defence Technology, Changsha, China; email: liuhw05@163.com. © 2012 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. LIU HONGWEI ET AL. 1026 Figure 14. The signal strength of third-degree spherical harmonics at different co-latitudes (longitude, m = r/4). << /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