The Relationship between Global Solar Radiation and Sunshine Durations in Cameroon ()

R. Mbiaké^{1*}, A. Beya Wakata^{2}, E. Mfoumou^{3}, E. Ndjeuna^{4}, L. Fotso^{1}, E. Tiekwe^{1}, J. R. Kaze Djamen^{5}, C. Bobda^{6}

^{1}University of Douala, Faculty of Sciences, Douala, Cameroon.

^{2}University of Yaoundé, Faculty of Sciences, Yaoundé, Cameroon.

^{3}Nova Scotia Community College, Division of Applied Research, Springhill, Canada.

^{4}Ecole Normale Supérieure d’Enseignement Technique (ENSET), Douala, Cameroon.

^{5}Centre de Physique Atomique, Moléculaire Optique et Quantique (CEPAMOQ), Douala, Cameroon.

^{6}Public Health of the San Diego University, San Diego, USA.

**DOI: **10.4236/ojap.2018.72006
PDF
HTML XML
1,000
Downloads
2,798
Views
Citations

Based on the well-known modified Angstrom formula on the relationship between the sunshine duration and the global solar radiation, this paper aimed to estimate the value of the constant a and b in Cameroon. Only five cities (Maroua, Garoua, NGaoundéré, Yaoundé and Douala) had the both available in-situ data recorded during the period of eleven years (1996-2006) beside which four others cities (Dschang, Koundja, Yoko and Manfé) had only the in-situ sunshine duration available data recorded during the period of twenty years (1986-2006). The 9 cities were grouped in 3 different climate regions. Based on the data of the 5 first cities belonging the 3 regions, the follow constant values a1 = -0.05, a2 = -0.02, a3 = -0.14 and b1 = 0.94, b2 = 0.74, b3 = 1.12 were obtained. The Root Mean Square Error (RMSE) Mean Bias Error (MBE) and correlation coefficient (r) were also determined. Then we used these values to estimate the global solar radiation for the other four remain cities. The constants a and b obtained values are in accordance with those of the West Africa region which Cameroon belongs to. So they can be employed in estimating global solar radiation of location in Cameroon paying attention only to the geographical location information.

Share and Cite:

Mbiaké, R. , Wakata, A. , Mfoumou, E. , Ndjeuna, E. , Fotso, L. , Tiekwe, E. , Djamen, J. and Bobda, C. (2018) The Relationship between Global Solar Radiation and Sunshine Durations in Cameroon. *Open Journal of Air Pollution*, **7**, 107-119. doi: 10.4236/ojap.2018.72006.

1. Introduction

Energy is the motive force behind the sustained technology development, especially for developing countries, where the renewable energy is capable of supplying the necessary energy for their rapid development. The awareness of the limited availability of non-renewable resources and their associated environmental problems is making it imperative that the world shift emphasizes the renewable energy resources. The environment consequences of harnessing these non-renewable energy sources are assuming alarming proportions.

So, the knowledge of global solar radiation and the sunshine duration received during the average day of each month are prerequisite in any solar energy exploitation for the optimal design and the prediction of the system performance. Obviously the best way of knowing the amount of global solar radiation at a given site is to install Pyrhéliomètre as possible at many locations in the region and look after them day to day.

A global study of the world distribution of the solar radiation has been carried out by Lof et al. [1] and Gueymard et al. [2] , while there have been many attempts to find common models applicable anywhere in the world [3] . Beside these studies, the Angstrom-Prescott relationship has been examined in many countries throughout the world, e.g. Australia [4] , Canada [5] , Caribbean countries as Guatemala [6] , the West Indies [7] , South Asian countries such as Sri Lanka [8] , Pakistan [9] in Africa continent Nigeria [10] , Sudan [11] , and South East Asia as Hong Kong [12] .

Nevertheless, one can notice that comparing to Northern countries, where there is a several hundreds of ground meteorological stations directly and indirectly measuring solar radiation, and continuous irradiance values deriving directly from the meteorological geostationary satellites (e.g.; METEOSAT), the countries in development face with the insufficient of these meteorology stations. Definitively in these regions, modeling is the better tool for estimating incoming global solar radiation from the sunshine duration, in a day of location where the measurements are not available [13] .

To reduce this gap in the African continent, this work aims to improve the global solar radiation in Cameroon by it estimation for the four cities (Manfé, Koundja, Yoko and Dschang) where only daily sunshine records are available.

2. Methodology

2.1. Theoretical Aspect

Theoretically, the global solar radiation that reaches at the earth’s surface is a result of complex interactions between the solar radiations at the top of the earth’s atmosphere evaluated by the solar constant G_{sc} = 1369 ± 6 W/m^{2} and the climatic conditions of the location. This evaluation does take also into account the rotation of the earth about its own axis and by its elliptical orbit about the sun and the solar declination (δ). Excluding the meteorological and climatic conditions, the daily extraterrestrial solar radiation H_{0} on a horizontal surface for a day in a month is well approximated by the given expression [14] :

${H}_{0}={G}_{sc}\times {\left(\frac{{R}_{0}}{R}\right)}^{2}\times \frac{24\times 3600}{\pi}\times \left[\mathrm{cos}\left(L\right)\mathrm{cos}\left(\delta \right)\mathrm{cos}\left(\omega \right)+\mathrm{sin}\left(L\right)\mathrm{sin}\left(\delta \right)\right]$ (1)

L is the latitude and $\omega $ is the solar hour angle. The solar declination $\left(\delta \right)$ is expressed by:

$\delta =23.45\mathrm{sin}\left[0.986\times \left(J+284\right)\right]$ (2)

Equation (1) overestimates the global irradiation as it does not take into account the interaction with atmosphere’s components. Indeed, the intensity of solar radiation through the atmosphere is attenuated by various atmospheric constituents, namely gases (air, molecules, ozone, CO_{2}, OH, etc.), liquid and solid particles (aerosols PM_{2.5}, PM_{10} and clouds). This attenuation is mainly due to clouds, and the path length through atmosphere is also critical.

Facing with the insufficiency and the quality of the in-situ measurement global solar radiation data, and the overestimation of the global irradiation by Equation (1), let us used in this paper the general empirical approach formula pioneered by angstrom in 1942 [15] , and later completed by Prescott 1940 [16] . This formula relates average daily global radiation with average daily sunshine hours:

$\stackrel{\xaf}{H}={\stackrel{\xaf}{H}}_{0}\left(a+b\frac{\stackrel{\xaf}{d}}{{\stackrel{\xaf}{d}}_{0}}\right)$ (3)

where a, b are the constants to be determined. $\stackrel{\xaf}{d}$ is the monthly mean daily bright sunshine duration obtained by the Campbell-Stokes recorder, and ${\stackrel{\xaf}{d}}_{0}$ is the highest monthly mean daily sunshine duration that can be calculated using the followed expression:

${\stackrel{\xaf}{d}}_{0}=\frac{2}{15}{\mathrm{cos}}^{-1}\left(\mathrm{tan}\delta \times \mathrm{tan}L\right)$ (4)

2.2. Solar Radiation Data in Cameroon

Despite the lowest availability of the radiation in-situ instruments at meteorological stations in Cameroon, we used in this paper, the sunshine duration of nine (9) stations and five (5) global solar radiations of them. These data were furnished by the Cameroon National Meteorological Office (DNM).

Based on the Trewarth’s climatic classification [17] , the 9 stations were grouped by region. Owing to its geographical position stretching from latitude 4˚N to 14˚N, Cameroon belongs to three following types of climate. Douala, Yaoundé, Manfé and Dschang cities belonging to latitude 3˚50'N and 5˚45'N make up the first region with a humid subtropical climate. The second region stretching from latitude 5˚50'N to 7˚30'N with a tropical wet climates contains the cities of Yoko, Koundja and N’Gaoundéré, while a warm and dry climate make up the third region with Garoua and Maroua. The 9 stations shared out in the 3 regions are presented in Table 1.

Among the 9 meteorological stations where these data were recorded, only 5 had the both types of data: the hour’s solar length in a day (d) and the global solar radiation (H) (Table 2).

Table 1. Geographical location of cities and records length.

Table 2. Stations with global solar radiation and sunshine duration availables.

For the 9 stations, the monthly average daily values over 10 to 20 years of measured data for the sunshine length hours were computed. Figure 1(a); Figure 2(a); Figure 3(a) are the curve of the measured values of the monthly average daily sunshine duration of the 3 climatic regions. And Figure 1(b), Figure 2(b), and Figure 3(b) are the mean value for each of them.

All the obtained plots brought out the similar trend of the monthly average daily sunshine. However, when we draw the 3 region curves on the same graph (Figure 4), it appears curiously that they have the same horizontal trend while their vertical positions is related to the latitudinal coordinate and the type of climate from the wet to the dry one.

Taking into account the above climatic classification, the only 5 stations that had in-situ recorded global solar radiation were grouped by region according to Table 2.

After the sunshine duration, we draw as one can see on Figure 5 the trend of the monthly average daily global solar radiation for the 5 stations that have these data. Once again, the 5 form curves have the same annual trend, while the graph

(a) (b)

Figure 1. (a) Monthly average daily sunshine duration for the 4 stations of region 1; (b) Mean value of the monthly average daily sunshine duration for region 1.

(a) (b)

Figure 2. (a) Monthly average daily sunshine duration for the 3 stations of region 2; (b) Mean value of the monthly average daily sunshine duration for region 2.

(a) (b)

Figure 3. (a) Monthly average daily sunshine duration for the 2 stations of region 3; (b) Mean value of the monthly average daily sunshine duration for region 3.

position depends to the region. Based on these curves trend and the way they gathered together in Figure 5, we draw in Figure 6 the blue curve to express the dry climate global solar radiation H_{d} and in red the humid subtropical one H_{h}.

2.3. Local Values of Angstrom-Prescott Constants

To estimate the unknown global solar radiation values of the 4 other stations, we computed the regression coefficients a and b using the following equations:

Figure 4. Monthly average daily sunshine.

Figure 5. Monthly average daily global solar radiation for the 5 stations.

Figure 6. Monthly average daily global solar radiation for dry and humid climate.

$a=\frac{\left[\left({\displaystyle \sum \frac{\stackrel{\xaf}{H}}{{\stackrel{\xaf}{H}}_{0}}}\right)\left({\displaystyle \sum {\left(\frac{\stackrel{\xaf}{d}}{{\stackrel{\xaf}{d}}_{0}}\right)}^{2}}\right)\right]-\left[\left({\displaystyle \sum \frac{\stackrel{\xaf}{d}}{{\stackrel{\xaf}{d}}_{0}}}\right)\left({\displaystyle \sum \left(\frac{\stackrel{\xaf}{d}}{{\stackrel{\xaf}{d}}_{0}}\right)\left(\frac{\stackrel{\xaf}{H}}{{\stackrel{\xaf}{H}}_{0}}\right)}\right)\right]}{\left[M{\displaystyle \sum {\left(\frac{\stackrel{\xaf}{d}}{{\stackrel{\xaf}{d}}_{0}}\right)}^{2}-{\left({\displaystyle \sum \frac{\stackrel{\xaf}{d}}{{\stackrel{\xaf}{d}}_{0}}}\right)}^{2}}\right]}$ (5)

$b=\frac{M\left({\displaystyle \sum \left(\frac{\stackrel{\xaf}{d}}{{\stackrel{\xaf}{d}}_{0}}\right)\left(\frac{\stackrel{\xaf}{H}}{{\stackrel{\xaf}{H}}_{0}}\right)}\right)-\left[\left({\displaystyle \sum \frac{\stackrel{\xaf}{d}}{{\stackrel{\xaf}{d}}_{0}}}\right){\displaystyle \sum \left(\frac{\stackrel{\xaf}{H}}{{\stackrel{\xaf}{H}}_{0}}\right)}\right]}{\left[M{\displaystyle \sum {\left(\frac{\stackrel{\xaf}{d}}{{\stackrel{\xaf}{d}}_{0}}\right)}^{2}-{\left({\displaystyle \sum \frac{\stackrel{\xaf}{d}}{{\stackrel{\xaf}{d}}_{0}}}\right)}^{2}}\right]}$ (6)

The obtained values were then used in Equation (3) to estimate the values of the global solar radiation ${\stackrel{\xaf}{H}}_{est}$ which were compared to its corresponding measured values ${\stackrel{\xaf}{H}}_{mes}$ at each of the 5 stations of Table 2. The deviation between the estimated and the measured values were summarized by calculating the following statistical parameters: Mean Bias Error (MBE), Root Mean Square Error (RMSE) and correlation coefficient r.

$MBE\left(\%\right)=100\left(\frac{1}{{\stackrel{\xaf}{H}}_{M}}\right)\left({\displaystyle \underset{i}{\sum}\frac{{E}_{i}}{M}}\right)$ $RMSE\left(\%\right)=100\left(\frac{1}{{\stackrel{\xaf}{H}}_{M}}\right)\sqrt{\left({\displaystyle \underset{i}{\sum}\frac{{E}_{i}^{2}}{M}}\right)}$ (7)

where ${E}_{i}={\stackrel{\xaf}{H}}_{est}-{\stackrel{\xaf}{H}}_{mes}$ with i = 1, M was the total number of the observations points.

$r=\frac{{\displaystyle \sum \left({\stackrel{\xaf}{H}}_{est}-{\stackrel{\xaf}{H}}_{E}\right)\left({\stackrel{\xaf}{H}}_{mes}-{\stackrel{\xaf}{H}}_{M}\right)}}{\sqrt{\left({\displaystyle \sum {\left({\stackrel{\xaf}{H}}_{est}-{\stackrel{\xaf}{H}}_{E}\right)}^{2}{\displaystyle \sum {\left({\stackrel{\xaf}{H}}_{mes}-{\stackrel{\xaf}{H}}_{M}\right)}^{2}}}\right)}}$ (8)

Angstrom constants and statistical parameters are compiled in Table 3.

Noticing that the both solar parameters were almost the same at a close latitudinal coordinate we determined the mean values of each region (Table 4).

To validate our model, we tested it by using Rietveld and Turton models [3] [18] that could be considered most likely to be appropriate.

Let’s recall that the Turton’s model analyzes on a long term monthly average sunshine duration and global solar radiation data from 25 stations especially in humid tropical countries around the world in the both hemispheres. The formulation of his model is:

$\frac{\stackrel{\xaf}{H}}{{\stackrel{\xaf}{H}}_{0}}=0.30+0.40\left(\frac{\stackrel{\xaf}{d}}{{\stackrel{\xaf}{d}}_{0}}\right)$ (9)

Table 3. The results of the angstrom and Prescott model applied for Cameroon.

Table 4. The results of the Angstrom and Prescott model applied for Cameroon.

While Rietveld’s model applied to 42 stations in different countries around the world that were selected to collect the same data, expecting to derive a unified correlation applicable anywhere in the world. The follow relation was established:

$\frac{\stackrel{\xaf}{H}}{{\stackrel{\xaf}{H}}_{0}}=0.18+0.62\left(\frac{\stackrel{\xaf}{d}}{{\stackrel{\xaf}{d}}_{0}}\right)$ (10)

Based on this two Equations (9) and (10) for the test, we used the length day ${\stackrel{\xaf}{d}}_{0}$ and the sunshine duration averages $\stackrel{\xaf}{d}$ to estimate the value of the global solar radiation for the 3 regions.

Afterwards we compared the outcome estimated values of global solar radiation with the local values of the regions (Table 5).

Table 7 contains the must values of a, b constants calculated around the world. These values let the constants of the region 1 to be the best to estimate the global solar radiation of the four remain cities namely Manfé, Koundja, Yoko and Dschang.

2.4. Estimation of the Unknown Solar Radiation

Amongst all these a, b constants in Table 7, it stand out that the best one to estimate the unknown global solar radiations of the Manfé, Koundja, Yoko and Dschang cities are those of the region 1 of our country model with the closeness of their geographical coordinates.

The outcome computed values are compiled in Table 8 and Figure 7. We have the graph of the estimated global unknown global solar radiation and the figure shows once again the trend coherence of the vertical position that it is linked to the geographic position.

3. Results and Discussion

Despite the fact that Driesse and Thervenard [19] pointed out the error prone in the measurement of the sunshine duration due for instance, the threshold that depends on the humidity of the recording card Iqbal [20] , or the human factors in the evaluation of the track burned by the Campbell-Stoke device introduces variability in sunshine recordings Baumgartner [21] , Painter [22] . The both limits can lead to the over or under estimation the sunshine duration at different times.

To go around this variability of expressing the sunshine duration in term of the average daily beam radiation H_{b}, the model suggested by Suehrcke (2000) [23] :

Table 5. Comparison between the 3 models: country, Turton’s and Rietveld’s models, based on the measured and estimated solar radiations.

$S=\frac{{\stackrel{\xaf}{H}}_{b}}{{\stackrel{\xaf}{H}}_{b,clear}}$ (11)

showed no better representation of the sunshine global radiation and let Angstrom and Prescott model be easier to use, after the analysis of the large data base of the World Data radiation Centre (WDRC).

So, as other authors who used widely the Angstrom-Prescott relation in many research area as, T. Dunne & L. B. Leopold [24] , J. M. Kowal et al. [25] and G. Guyot [26] in hydrological and agro meteorology applications, C. Augustine et al. [10] , and E. O. Falayi et al. [27] in renewable Energy to predict the solar radiation, we did the same calculating the regression coefficients a and b for the three selected regions summarized in Table 5.

To make sure of the quality of our model, we applied Turton’s, Rietveld’s and our model to the three select regions, and the obtained results were compared to one year in-situ measured solar radiation data, they are compiled in Table 6. Our model gives better RMSE, MBE and correlation coefficient r values compared with those of the Turton’s, Rietveld’s models (Table 7).

In Table 8, our model with its parameters between the others computed for different areas in the world show that the found values of a, b are clearly situated in the space defined by Davies in the case of West Africa [28] .

The regression coefficients a, b obtained from our model are clearly acceptable if compared with Turton’s and Rietveld’s models (Table 7), and be used successfully to estimate the monthly average daily global solar radiation of the 4 cities that have no solar radiation in-situ measurement data.

Figure 7. Monthly estimated daily global solar radiation for the 4 remained cities.

Table 6. Comparison of the regression a and b coefficients.

Table 7. Some published Angstrom constants and those computed in this work.

4. Conclusions

As others researchers who applied the Angstrom type regression model for estimating global solar radiation, we used this model to obtain, for the first time, the

Table 8. The estimated values of the solar radiation the unknown solar cities radiation.

local values of a and b for Cameroon three regions that will lead a better prediction of the solar potential energy at a given site.

It can be point out that a, b coefficient of the region 1, 2 and 3 on Table 6 let the estimated values of global solar radiation close to the measured one. Moreover, this study shows that countries without enough in-situ meteorological equipment’s to broadcast solar data can easily use Angstrom model to accurate the value of the predicted global solar radiation.

The better the predicted global solar are, the better evaluated potential energy for good exploitation will be.

Conflicts of Interest

The authors declare no conflicts of interest.

[1] | Lof, G.O.G., Duffie, J.A. and Smith, C.O. (1966) World Distribution of Solar Radiation. Solar Energy, 10, 27-37. https://doi.org/10.1016/0038-092X(66)90069-7 |

[2] | Gueymard, C.A. and Thevenard, D. (2009) Monthly Average Clear-Sky Broadband Irradiance Database for Worldwide Solar Heat Gain and Building Cooling Load Calculations. Solar Energy, 83, 1998-2018. |

[3] |
Rietveld, M.R. (1978) A New Method for Estimating the Regression Coefficient in the Formula Relating Solar Radiation to Sunshine. Agricultural Meteorology, 19, 243-252. https://doi.org/10.1016/0002-1571(78)90014-6 |

[4] |
Hutchinson, M.F., Booth, T.H., McMahon, J.P. and Nix, H.A. (1984) Estimating Monthly Mean Value of Daily Total Solar Radiation for Australia. Solar Energy, 32, 277-290. https://doi.org/10.1016/S0038-092X(84)80045-6 |

[5] | Haluoani, N., Nguyen, C.T. and Vo-Ngoc, D. (1993) Calculation of Monthly Average Global Solar Radiation, on Horizontal Surfaces Using Daily Hours of Bright Sunshine. Solar Energy, 50, 247-258. https://doi.org/10.1016/0038-092X(93)90018-J |

[6] |
Orozco, E.B. (1987) Gautemalan Solar Map. Solar and Wind Technology, 4, 381-388. https://doi.org/10.1016/0741-983X(87)90068-3 |

[7] |
Hawas, M.M. and Muner, T. (1983) Correlation between Global Radiation and Sunshine Data for India. Solar Energy, 30, 289-290. https://doi.org/10.1016/0038-092X(83)90158-5 |

[8] | Samuel, T.D.M.A. (1991) Estimation of Global Radiation for Sri-Lanka. Solar Energy, 47, 333-337. https://doi.org/10.1016/0038-092X(91)90026-S |

[9] | Raja, I.A. and Twidell, J.W. (1989) Distribution of Global Insolation over Pakistan. Solar Energy, 43, 355-357. https://doi.org/10.1016/0038-092X(89)90106-0 |

[10] | Augustine, C. and Nnabuchi, M.N. (2009) Relationship between Global Solar Relation Sunshine Hours for Calabar, Port Harcourt and Enugu, Nigeria. International Journal of Physic Sciences, 4, 182-188. |

[11] | Khogali, A. (1983) Solar Radiation over Sudan—Comparison of Measured and Predicted Data. Solar Energy, 31, 45-53. https://doi.org/10.1016/0038-092X(83)90032-4 |

[12] | Leung, C.T. (1980) The Fluctuation of Solar Irradiance in Hong Kong. Solar Energy, 25, 485-494. https://doi.org/10.1016/0038-092X(80)90080-8 |

[13] |
Rasmussen, M.S. (1998) Developing Simple, Operational, Consistent NDVI-Vegetation Models by Applying Environmental and Climatic Information: Part I. Assessment of Net Primary Production. International Journal of Remote Sensing, 19, 97-117. https://doi.org/10.1080/014311698216459 |

[14] | Duffie, J.A. and Beckman, W.A. (1994) Solar Engineering of Thermal Processes. 2nd Edition, John Wiley & Son, New York. |

[15] | Angström, A. (1924) Solar Terrestrial Radiation. Quarterly Journal of the Royal Meteorological Society, 50, 121-125. |

[16] | Precott, J.A. (1940) Evaporation from a Water Surface in Relation to Solar Radiation. Transactions of the Royal Society of South Australia, 64, 114-125. |

[17] | Trewartha, G.T. (1961) The Earth’s Problem Climates. The University of Wisconsin, Madison. |

[18] |
Turton, S.M. (1987) The Relationship between Total Irradiation and Sunshine Duration in the Humid Tropics. Solar Energy, 38, 353-354. https://doi.org/10.1016/0038-092X(87)90007-7 |

[19] | Driesse, A. and Thevenard, D. (2002) A Test of Suebrcke’s Sunshine-Radiation Relationship Using a Global Data Set. Solar Energy, 72, 167-175. |

[20] | Iqbal, M. (1983) An Introduction to Solar Radiation. Academic Press, Toronto. |

[21] | Baumgartner, T. (1979) Die Schwellenintensitat des Sonnenscheinautog-Raphen Campbell-Stokes an wolkenlosen Tagen, Arbeitsberichte der Schweizerischen Meteorologischen Zentralanstalt. Zürich. |

[22] | Painter, H.E. (1981) The Performance of a Campbell-Stokes Sunshine Recorder Compared with a Simultaneous Record of Normal Incidence Irradiance. Meteorological Magazine, 110, 102-109. |

[23] |
Suehrcke, H., Bowden. R.S. and Hollnds, K.G.T. (2013) Relationship between sunshine Duration and Solar Radiation. Solar Energy, 92, 160-171. https://doi.org/10.1016/j.solener.2013.02.026 |

[24] | Dunne, T. and Leopold L.B. (1978) Water in Environmental Planning. W.H. Freeman and Company, New York, 818. |

[25] | Kowal, J.M. and Kassam, A.H. (1978) Agricultural Ecology of Savanna: A Study of West Africa. Clarendon Press, Oxford, 403. |

[26] | Guyot, G. (1998) Physics of the Environment and Climate. John Wiley & Son, Chichester, 632. |

[27] | Falayi, E.O., Adepitan, J.O. and Rabiu, A.B. (2008) Empirical Models for Correlation of Global Solar Radiation with Meteorological Data for Iseyin, Nigeria. International Journal of Physical Sciences, 3, 210-216. |

[28] |
Davies, J.A. (1966) The Assessment of Evapotranspiration for Nigeria. Geografiska Annaler: Series A, Physical Geography, 48, 136-156. https://doi.org/10.1080/04353676.1966.11879735 |

Journals Menu

Contact us

customer@scirp.org | |

+86 18163351462(WhatsApp) | |

1655362766 | |

Paper Publishing WeChat |

Copyright © 2023 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.