152: Belarus - temperature trends STABLE before 1980

When it comes to analysing the temperature data of Belarus the biggest problem is the low quantity of data. There is no data before 1880 and only three stations have data before 1950, two of which are long stations with over 1200 months of data (for a full list of stations see here). On a more positive note, there are seventeen medium stations with over 480 months of data and these are quite evenly distributed across the country. That means it should be possible to construct a reliable measure of the overall mean temperature change for Belarus, at least for the last sixty years. What this data then shows is that the climate of Belarus appears to have been fairly stable over the hundred years prior to 1980, but then in 1988, like much of Europe, the temperature suddenly increased by about 1.1°C (see Fig. 152.1 below).

Fig. 152.1: The mean temperature change for Belarus since 1880 relative to the 1981-2010 monthly averages. The best fit is applied to the monthly mean data from 1881 to 1980 and has a statistically insignificant positive gradient of +0.23 ± 0.23 °C per century. After 1980 there is an abrupt warming of 1.1°C.

In order to quantify the changes to the climate of Belarus the temperature anomalies for all stations with over 480 months of data before 2014 were determined and averaged. This was done using the usual method as outlined in Post 47 and involved first calculating the temperature anomaly each month for each of the nineteen valid stations relative to its own monthly reference temperature (MRT). Then those anomalies were averaged to determine the mean temperature anomaly (MTA) for the whole country for each month. The MRTs for each station in Belarus were calculated using the same 30-year period, namely from 1981 to 2010. The resulting MTA for Belarus is shown as a time series in Fig. 152.1 above and clearly shows that temperatures were rising slowly (at about 0.23°C per century) for about 100 years up until 1980 but that this rise was only comparable to the uncertainty in the trend and so is not significant. After 1980, however, the MTA suddenly increases by about 1.1°C. Such behaviour is seen in the MTA for many other European countries and for Europe as a whole (see Post 44).

The total number of stations included in the MTA in Fig. 152.1 each month is shown in Fig. 152.2 below. The peak in the frequency after 1970 indicates why the 1981-2010 interval was determined to be the most appropriate to use for calculating the MRTs in this case.

Fig. 152.2: The number of station records included each month in the mean temperature anomaly (MTA) trend for Belarus in Fig. 152.1.

The locations of the nineteen stations used to calculate the MTA in Fig. 152.1 are indicated on the map in Fig. 152.3 below. These stations appear to be distributed very evenly across the country with no significant clusters. This means that a simple average of their temperature anomalies should be just as accurate as any of the gridding or homogenization processes that are used by the main climate science groups in their analyses.

Fig. 152.3: The (approximate) locations of the 19 longest weather station records in Belarus. Those stations denoted with squares are long stations with over 1200 months of data, while diamonds denote medium stations with more than 480 months of data.

If we next consider the change in temperature based on Berkeley Earth (BE) adjusted data we get the MTA data in Fig. 152.4 below. This again was determined by averaging the anomalies each month from the nineteen longest stations and also suggests that the climate was warming very slightly before 1980 and then more rapidly thereafter. In this case, however, the post-1980 warming is more continuous and gradual in nature than that seen in the raw data in Fig. 152.1.

Fig. 152.4: Temperature trends for Belarus based on Berkeley Earth adjusted data. The best fit linear trend line (in red) is for the 12-month moving average data over the period 1881-1980 and has a positive gradient of +0.25 ± 0.08°C/century.

If we compare the curves in Fig. 152.4 with those from the published Berkeley Earth (BE) version for Belarus shown in Fig. 152.5 below, we see that there is excellent agreement between the two sets of data at least as far back as 1880. This indicates that the simple averaging of adjusted anomalies used to generate the BE MTA in Fig. 152.4 is as effective and accurate as the more complex gridding method used by Berkeley Earth in Fig. 152.5. In which case simple averaging should be just as effective and accurate in generating the MTA using raw unadjusted data in Fig. 152.1. What is more difficult to explain is how Berkeley Earth was able to determine the mean temperature of Belarus as far back as 1750 when there appears to be no reliable data before 1880.

Fig. 152.5: The temperature trend for Belarus since 1750 according to Berkeley Earth.

But if we next compare the adjusted data in Fig. 152.4 with the raw data shown in Fig. 152.1 we see that there is excellent agreement between these two sets of data as well (see Fig. 152.6 below).

Fig. 152.6: A comparison of the 5-year mean temperature change for Belarus since 1880 between the original raw data from Fig. 152.1 (in blue) and the Berkeley Earth adjusted data from Fig. 152.4 (in red).

The small differences between the MTA from the raw data in Fig. 152.1 and that from the BE adjusted data in Fig. 152.4  are mainly due to the data processing procedures used by Berkeley Earth. These include homogenization, gridding, Kriging and most significantly breakpoint adjustments. These lead to changes to the original temperature data, the magnitude of these adjustments being the difference in the MTA values seen in Fig. 152.1 and Fig. 152.4. The magnitudes of these adjustments are shown graphically in Fig. 152.7 below.

Fig. 152.7: The contribution of Berkeley Earth (BE) adjustments to the anomaly data in Fig. 152.4 after smoothing with a 12-month moving average. The blue curve represents the total BE adjustments including those from homogenization. The linear best fit (red line) to these adjustments for the period 1911-2010 has a negative gradient of -0.242 ± 0.016 °C per century. The orange curve shows the contribution just from breakpoint adjustments.

The blue curve in Fig. 152.7 is the difference in MTA values between adjusted (Fig. 152.4) and unadjusted data (Fig. 152.1), while the orange curve is the contribution to those adjustments arising solely from breakpoint adjustments. Both are relatively small with the former adding a slight cooling to the data since 1911 of about 0.24°C. The large offset between the blue curve and the orange curve in Fig. 152.7 is due to the different MRT intervals used for calculating the anomalies in Fig. 152.1 (1981-2010) and in Fig. 152.4 (1961-1990).

Fig. 152.8: A comparison of the 5-year mean temperature change for Belarus (blue curve) with that of the Baltic States in Post 51 (red curve).

As mentioned at the start of this post, the main weakness of the data for Belarus is the lack of good data before 1950 which obviously raises questions over its accuracy and reliabilty. One way to test the accuracy is to compare the Belarus data with that from neighbouring states to see what the level of similarity is. This has been done in Fig. 152.8 above where the comparator set is data from the Baltic States in Post 51. As can clearly be seen, the level of agreement between the two datasets is very good, particularly after 1950. However, even before 1950 there is good agreement even though the Belarus temperature trend is based on data from three stations at most. This suggests that the MTA for Belarus in Fig. 152.1 is reliable and likely to reflect the true climate of Belarus as far back as 1880.

Summary

The raw temperature data for Belarus clearly shows that the climate was stable up until 1980. Any warming in the trend over this period was significantly less than the natural variation in the mean temperature anomaly (MTA) that was seen for all timescales up to ten years in duration (see Fig. 152.1 and Fig. 152.4).

After 1980 the climate has warmed sharply by about 1.1°C. This behaviour is similar to patterns seen across Europe. The reason for this abrupt temperature increase is still unknown as it does not correlate with increases in carbon dioxide levels in the atmosphere.

There appears to be a very good level of correlation between the MTA trend for Belarus and that for the Baltic States reported in Post 51. This allows each to in effect corroborate the other and therefore strengthen the validity of each.



Acronyms

BE = Berkeley Earth.

MRT = monthly reference temperature (see Post 47).

MTA = mean temperature anomaly.

Long station = a station with over 1200 months (100 years) of data before 2014.

Medium station = a station with over 480 months (40 years) of data before 2014.

List of all stations in Belarus with links to their raw data files.

150: Malta - temperature trends STABLE before 1980

The difficulty in determining the extent of climate change in Malta is the lack of data. According to Berkeley Earth (BE) there are only two stations with over 480 months of data and one of those (A. M. D.) has virtually no data after 1934. The other station is Luqa (BE ID:156721) which is situated near the main airport and actually has over 1800 months of data, so it could be a good indicator of the climate of Malta as a whole. Its monthly mean temperature anomaly (MTA) relative to the period 1981-2010 is shown in Fig. 150.1 below.

 

Fig. 150.1: The mean temperature change for Luqa since 1840 relative to the 1981-2010 monthly averages. The best fit is applied to the monthly mean data from 1881 to 1980 and has a slight positive gradient of +0.05 ± 0.11 °C per century.

 

The data in Fig. 150.1 indicates that there was virtually no climate change in Malta before 1980. Then after 1980 the local temperature appears to have risen by about 0.6°C, although this depends on how one interprets the data. If one uses the 5-year average (yellow line) as a guide the temperature appears to increase by over 1°C from 1880 to 2010. However, average temperatures in 2010 are only about 0.6°C above the 1881-1980 best fit line (in red) which is virtually horizontal. So this would indicate that temperatures in 2010 are only about 0.6°C above average.

 

Fig. 150.2: Temperature trends for Luqa based on Berkeley Earth adjusted data. The best fit linear trend line (in red) is for the period 1886-2010 and has a positive gradient of +0.56 ± 0.03°C/century.

 

A similar temperature trend for Luqa is seen in its Berkeley Earth (BE) adjusted data (see Fig. 150.2 above) which appears to have zero adjustments after 1870. However, the actual published BE trends for Malta shown in Fig. 150.3 below appear to exhibit more warming over the period from 1890 to 2000 (~1.6°C) than is seen in the Luqa data in Fig. 150.2 (~1.4°C). It is also interesting that the temperature trend in Fig. 150.3 extends back to 1750 even though there appears to be no temperature data for Malta before 1850 (for a list of stations in Malta see here).

 

 Fig. 150.3: The temperature trend for Malta since 1750 according to Berkeley Earth.

 

If we compare the raw data for Luqa with the BE adjusted version we find that there is virtually no difference between the two (see Fig. 150.4 below). It is actually very unusual for such a long dataset not to undergo at least one breakpoint adjustment so this result is quite surprising.

 

Fig. 150.4: The 5-year average of the monthly mean temperature change for Luqa since 1820 based on the original raw data from Fig. 150.1 (in blue) and the Berkeley Earth adjusted data from Fig. 150.2 (in red).

 

As there is no data from Malta to compare the Luqa data against we could instead compare it with temperature data from the nearby Italian islands of Lampedusa (BE ID: 155869) and Pantelleria (BE ID: 175525 ). The locations of these islands relative to Malta and Sicily is shown in the map in Fig. 150.5 below. Lampedusa is about 150 km from both Malta and Pantelleria.

Fig. 150.5: The (approximate) locations of the weather stations in Malta (Luqa), Lampedusa and Pantelleria. Those stations denoted with squares are long stations with over 1200 months of data, while diamonds denote medium stations with more than 480 months of data.

 

Unfortunately the temperature data from Lampedusa and Pantelleria only start around 1960 but both sets show temperature rises after 1980 that are similar to that seen in the Luqa data. There is also good correlation in the high frequency (less than 12 months) fluctuations, but the medium timescale fluctuations with peak widths of more than 12 months are not well correlated.

 

Fig. 150.6: A comparison of the 5-year average temperature change for the islands of Malta, Lampedusa and Pantelleria since 1960.

 

The other data that we could compare Luqa with is that of Italy. This is shown in Fig. 150.7 below and here the medium timescale fluctuations in the two datasets correlate much better.

 

Fig. 150.7: A comparison of the 5-year average temperature change for Malta and Italy.

 

 

Summary

The temperature data for Malta clearly shows that the climate warmed by over 0.6°C after 1980 (see Fig. 150.1).

Before 1980 there is a large amount of natural variation in the temperature data but no overall upward trend.

The temperature trend for Malta can only be determined using one set of station data: Luqa (BE ID:156721). This makes it very unreliable. Comparing it with the neighbouring islands of Lampedusa and Pantelleria only confirms the temperature trend after 1960 (see Fig. 150.6).

There does appear to be a good correlation between the temperature trend of Luqa in Malta and that of Italy in Fig. 148.2 of Post 148 extending back as far as 1890 (see Fig. 150.7). This would suggest that the temperature trends of Malta and Italy have been very similar over at least the last 150 years.



Acronyms

BE = Berkeley Earth.

MRT = monthly reference temperature (see Post 47).

MTA = mean temperature anomaly.

Long station = a station with over 1200 months (100 years) of data before 2014.

Medium station = a station with over 480 months (40 years) of data before 2014.

List of all stations in Malta with links to their raw data files.

148: Italy - temperature trends STABLE before 1980

Whereas France only has four long stations with over 1200 months of data, Italy has at least twelve long stations, most with data stretching back over two hundred years. This means it is possible to determine  with a high degree of certainty the extent of climate change in Italy as far back as 1820. What this climate data shows is that the climate of Italy was stable for almost two hundred years up until 1980. Then over the last forty years it has warmed by about 1°C.

In addition to its twelve long stations, Italy also has 81 medium stations with over 480 months of data (for a full list see here). The locations of these 93 stations are shown on the map in Fig. 148.1 below. While the stations are generally spread evenly, there is a higher concentration of long stations in the north of the country compared to the south, and some clustering around Milan, Venice and Rome.

Fig. 148.1: The (approximate) locations of the 93 longest weather station records in Italy. Those stations with a high warming trend between 1911 and 2010 are marked in red while those with a cooling or stable trend are marked in blue. Those denoted with squares are long stations with over 1200 months of data, while diamonds denote medium stations with more than 480 months of data.

In order to quantify the changes to the climate of Italy the temperature anomalies for all stations with over 480 months of data before 2014 were determined and averaged. This was done using the usual method as outlined in Post 47 and involved first calculating the temperature anomaly each month for each station relative to its monthly reference temperature (MRT), and then averaging those anomalies to determine the mean temperature anomaly (MTA) for the whole country for each month. The MRTs for each station in Italy were calculated using the same 30-year period, namely from 1961 to 1990. 

The resulting MTA is shown as a time series in Fig. 148.2 below and clearly shows that temperatures were stable for over two hundred years up until 1980. Then they appear to increase rapidly by about 1.0°C over thirty years. That said, the change is comparable in size and speed to natural variations seen in earlier parts of the trend such as in 1940. On the other hand the MTA after 1980 is based on data from many more stations (up to ninety) and so is likely to be more accurate.

Fig. 148.2: The mean temperature change for Italy since 1740 relative to the 1961-1990 monthly averages. The best fit is applied to the monthly mean data from 1881 to 1980 and has a slight positive gradient of +0.03 ± 0.12 °C per century.

The total number of stations included in the MTA in Fig. 148.2 each month is shown in Fig. 148.3 below. The peak in the frequency around 1970 suggests that the 1961-1990 interval was indeed the most appropriate to use for the MRTs. It also indicates that data from about ten stations were used to calculate the MTA for almost every month back to 1820. As fifteen stations appears to the minimum number needed to provide an accurate MTA, this suggests that the trend in Fig. 148.2 is reliable at least as far back as 1820.

Fig. 148.3: The number of station records included each month in the mean temperature anomaly (MTA) trend for Italy in Fig. 148.2.

If we next consider the change in temperature based on Berkeley Earth (BE) adjusted data we get the MTA data in Fig. 148.4 below. This again was determined by averaging each month the anomalies from the 93 longest stations but here the picture is slightly different from that depicted in Fig. 148.2. Instead of a stable MTA we see a positive trend of over 0.5°C per century suggesting that the climate was warming before 1980. This clearly contradicts the raw data in Fig. 148.2.

Fig. 148.4: Temperature trends for Italy based on Berkeley Earth adjusted data. The best fit linear trend line (in red) is for the period 1881-1980 and has a positive gradient of +0.53 ± 0.04°C/century.

But if we next compare the curves in Fig. 148.4 with those from the published Berkeley Earth (BE) version for Italy shown in Fig. 148.5 below, we see that there is excellent agreement between the two sets of data at least as far back as 1770. This indicates that the simple averaging of adjusted anomalies used to generate the BE MTA in Fig. 148.4 is as effective and accurate as the more complex gridding method used by Berkeley Earth in Fig. 148.5. In which case simple averaging should be just as effective and accurate in generating the MTA using raw unadjusted data in Fig. 148.2.

Fig. 148.5: The temperature trend for Italy since 1750 according to Berkeley Earth.

The differences between the MTA in Fig. 148.2 and the BE versions using adjusted data in Fig. 148.4  are instead mainly due to the data processing procedures used by Berkeley Earth. These include homogenization, gridding, Kriging and most significantly breakpoint adjustments. These lead to changes to the original temperature data, the magnitude of these adjustments being the difference in the MTA values seen in Fig. 148.2 and Fig. 148.4. The magnitudes of these adjustments are shown graphically in Fig. 148.6 below. The blue curve is the difference in MTA values between adjusted (Fig. 148.4) and unadjusted data (Fig. 148.2), while the orange curve is the contribution to those adjustments arising solely from breakpoint adjustments. Between 1880 and 1980 both are considerable with the former leading to an additional warming since 1880 of over 0.5°C.

Fig. 148.6: The contribution of Berkeley Earth (BE) adjustments to the anomaly data in Fig. 148.4 after smoothing with a 12-month moving average. The blue curve represents the total BE adjustments including those from homogenization. The linear best fit (red line) to these adjustments for the period 1881-1980 has a positive gradient of +0.510 ± 0.009 °C per century. The orange curve shows the contribution just from breakpoint adjustments.

The overall impact of the BE adjustments can be seen more clearly if we compare the 5-year averages for the raw data (from Fig. 148.2) and the BE adjusted data (from Fig. 148.4). This comparison is shown in Fig. 148.7 below. It clearly shows that the trend based on adjusted data (red curve) exhibits considerably more warming since 1840 but slightly less since 1990.

Fig. 148.7: The 5-year mean temperature change for Italy since 1740 based on the original raw data from Fig. 148.2 (in blue) and the Berkeley Earth adjusted data from Fig. 148.4 (in red).

Summary

The raw unadjusted temperature data for Italy clearly shows that the climate was stable from 1780 to 1980 (see Fig. 148.2).

In contrast, the BE adjusted data claims that the climate first cooled and then warmed by 0.5°C from 1880 to 1980 (see Fig. 148.4).

After 1980 the climate has clearly warmed by about 1°C but this is still only of the same magnitude as the natural variations in the climate, so it may be too early to state definitively how much of the warming is permanent and how much more is still to come.



Acronyms

BE = Berkeley Earth.

MRT = monthly reference temperature (see Post 47).

MTA = mean temperature anomaly.

Long station = a station with over 1200 months (100 years) of data before 2014.

Medium station = a station with over 480 months (40 years) of data before 2014.

List of all stations in Italy with links to their raw data files.

147: France - temperature trends STABLE before 1980

What is surprising about France is the poor quality of its temperature data. Its data is worse than for Spain, Italy, Germany, the UK, Norway, Sweden and Finland. And yet France is the de-facto home of metrology, the country that gave us SI units.

The longest temperature record in France is for Bourges in the middle of the country. It is one of only four long stations with over 1200 months of data in France. There are a further 91 medium stations in France with over 480 months of data (for a full list see here). The locations of all these 95 stations are shown on the map in Fig. 147.1 below. In this analysis I have also included two stations in the Channel Islands; a long station in Guernsey and a medium station in Jersey. The reason for this is that they are much closer to France than England and so it is more reasonable for their data to be combined with data from France than with that from the UK. It also results in another long station being included in the analysis for France thereby improving the reliability of its long-term trend.

Fig. 147.1: The (approximate) locations of the 97 longest weather station records in France and the Channel Islands. Those stations with a high warming trend are marked in red while those with a cooling or stable trend are marked in blue. Those denoted with squares are long stations with over 1200 months of data, while diamonds denote medium stations with more than 480 months of data.

 

In order to quantify the changes to the climate of France the temperature anomalies for all stations with over 480 months of data before 2014 were determined and averaged. This was done using the usual method as outlined in Post 47 and involved first calculating the temperature anomaly each month for each station relative to its monthly reference temperature (MRT), and then averaging those anomalies to determine the mean temperature anomaly (MTA) for the whole country for each month. The MRTs for each station in France were calculated using the same 30-year period, namely from 1961 to 1990. The resulting MTA is shown as a time series in Fig. 147.2 below and clearly shows that temperatures were fairly stable for over 120 years up until 1980. Then they appear to increase suddenly by over 0.8°C.

Fig. 147.2: The mean temperature change for France since 1820 relative to the 1961-1990 monthly averages. The best fit is applied to the monthly mean data from 1861 to 1980 and has a slight positive gradient of +0.22 ± 0.11 °C per century.

The total number of stations included in the MTA in Fig. 147.2 each month is shown in Fig. 147.3 below. The peak in the frequency around 1990 suggests that the 1961-1990 interval was an appropriate one to use for the MRTs as it enabled all but five of the 97 datasets to be included in the MTA. All five of these medium station datasets had no data after 1900 and at least four exhibited a negative temperature trend. This suggests that in the 19th century the climate of France was cooling not warming.

Fig. 147.3 also indicates that data from less than ten stations were used to calculate the MTA for almost every month before 1945. As fifteen stations appears to the threshold number needed to provide an accurate MTA, this suggests that the trend in Fig. 146.2 is reliable only as far back as 1950. In fact most of the MTA trend before 1930 is dependent on data from only six stations. Of these, two are based in large cities (Paris and Marseille) and appear to exhibit strong linear warming trends (over 1.7°C since 1900) consistent with an urban heat island effect, and three have fragmented data (Bourges, Guernsey, and Montpellier). So only one (Chateauroux) is of any real quality and it has warmed by only about 0.3°C from 1900 to 2013.

Fig. 147.3: The number of station records included each month in the mean temperature anomaly (MTA) trend for France in Fig. 147.2.

If we next consider the change in temperature based on Berkeley Earth (BE) adjusted data we get the MTA data in Fig. 147.4 below. This again was determined by averaging each month the anomalies from the 97 longest stations and suggests that the climate was slowly warming before 1980 by about 0.3°C since 1860. While this is slightly more than the 0.2°C suggested in Fig. 147.2, the difference is not that significant.

Fig. 146.4: Temperature trends for France based on Berkeley Earth adjusted data. The best fit linear trend line (in red) is for the period 1861-1980 and has a positive gradient of +0.31 ± 0.04°C/century.

Next, if we compare the curves in Fig. 147.4 with those from the published Berkeley Earth (BE) version for France shown in Fig. 147.5 below, we see that there is excellent agreement between the two sets of data at least as far back as 1825. This indicates that the simple averaging of adjusted anomalies used to generate the BE MTA in Fig. 147.4 is as effective and accurate as the more complex gridding method used by Berkeley Earth in Fig. 147.5. In which case simple averaging should be just as effective and accurate in generating the MTA using raw unadjusted data in Fig. 147.2.

Fig. 147.5: The temperature trend for France since 1750 according to Berkeley Earth.

Any differences between the MTA in Fig. 147.2 and the BE versions using adjusted data in Fig. 147.4  are mainly due to the data processing procedures used by Berkeley Earth. These include homogenization, gridding, Kriging and most significantly breakpoint adjustments. These lead to changes to the original temperature data, the magnitude of these adjustments being the difference in the MTA values seen in Fig. 147.2 and Fig. 147.4. The magnitudes of these adjustments are shown graphically in Fig. 147.6 below. The blue curve is the difference in MTA values between adjusted (Fig. 147.4) and unadjusted data (Fig. 147.2), while the orange curve is the contribution to those adjustments arising solely from breakpoint adjustments. Both are minimal with the former leading to an additional warming since 1860 of between 0.1°C and 0.2°C.

Fig. 147.6: The contribution of Berkeley Earth (BE) adjustments to the anomaly data in Fig. 147.4 after smoothing with a 12-month moving average. The blue curve represents the total BE adjustments including those from homogenization. The linear best fit (red line) to these adjustments for the period 1861-1980 has a positive gradient of +0.079 ± 0.010 °C per century. The orange curve shows the contribution just from breakpoint adjustments.

The overall impact of the BE adjustments can be seen more clearly if we compare the 5-year averages for the raw data (from Fig. 147.2) and the BE adjusted data (from Fig. 147.4). This comparison is shown in Fig. 147.7 below. It clearly shows that the trend based on adjusted data (red curve) exhibits more warming before 1900 but very little extra afterwards.

Fig. 147.7: The 5-year mean temperature change for France since 1780 based on the original raw data from Fig. 147.2 (in blue) and the Berkeley Earth adjusted data from Fig. 147.4 (in red).

Summary

The raw unadjusted temperature data for France clearly shows that the climate warmed by less than 0.2°C from 1880 to 1980 (see Fig. 147.2)

In contrast, the BE adjusted data claims that the climate warmed by 0.3°C over the same period (see Fig. 147.4).

After 1980 the climate has clearly warmed by almost 1°C (see Fig. 147.7).

The temperature data before 1980 presented here clearly disagree with that of Spain in Post 146. However, as the MTA for Spain over the 100 years before 1950 is based on data from between twelve and fifty different stations compared to only about five for France, that would suggest that the Spain data is the more accurate.



Acronyms

BE = Berkeley Earth.

MRT = monthly reference temperature (see Post 47).

MTA = mean temperature anomaly.

Long station = a station with over 1200 months (100 years) of data before 2014.

Medium station = a station with over 480 months (40 years) of data before 2014.

List of all stations in France with links to their raw data files.

146: Spain - temperature trends COOLING before 1980

The analysis of temperature data for Portugal in the previous post appeared to indicate that the local climate had warmed continuously since 1870 by over 1°C in total. The caveat to this was the poor data quantity before 1940 with only two stations of significance contributing data to the regional trend, and one of them, Lisbon, clearly displayed characteristics in its data that were suggestive of the urban heat island effect. 

In this post I will look at the corresponding temperature data for Portugal's neighbour, Spain, to see if the temperature trends seen in Post 145 are repeated, as would be expected of neighbouring territories. The results will in fact show that they are not, and that the trends for Portugal before 1940 are probably wrong. In fact the data for Spain indicates that the climate cooled over the one hundred years before 1980 and has only recently begun to warm.

Spain has many more weather stations compared to Portugal, which given the difference in size is not surprising. Despite this, the station densities of the two countries are roughly the same. Spain has fourteen long stations with over 1200 months of data before 2014 and another 69 medium stations with over 480 months of data. The locations of these stations are indicated on the map in Fig. 146.1 below. They include two stations in Gibraltar and four in the Balearic Islands. They are distributed fairly evenly across the country but there are significant clusters around Madrid and south of Seville. 

Fig. 146.1: The (approximate) locations of the 83 longest weather station records in Spain. Those stations with a high warming trend between 1911 and 2010 are marked in red while those with a cooling or stable trend are marked in blue. Those denoted with squares are long stations with over 1200 months of data, while diamonds denote medium stations with more than 480 months of data.

In order to quantify the changes to the climate of Spain the temperature anomalies for all stations with over 480 months of data before 2014 were determined and averaged. This was done using the usual method as outlined in Post 47 and involved first calculating the temperature anomaly each month for each station relative to its monthly reference temperature (MRT), and then averaging those anomalies to determine the mean temperature anomaly (MTA) for the whole country for each month. The MRTs for each station in Spain were calculated using the same 30-year period, namely from 1961 to 1990. The resulting MTA is shown as a time series in Fig. 146.2 below and clearly shows that temperatures were decreasing for over 110 years up until 1980. After which they appear to increase suddenly by about 0.8°C.

Fig. 146.2: The mean temperature change for Spain since 1780 relative to the 1961-1990 monthly averages. The best fit is applied to the monthly mean data from 1871 to 1980 and has a slight negative gradient of -0.19 ± 0.10 °C per century.

The total number of stations included in the MTA in Fig. 146.2 each month is shown in Fig. 146.3 below. The peak in the frequency around 1970 suggests that the 1961-1990 interval was indeed the most appropriate to use for the MRTs. It also indicates that data from at least ten stations were used to calculate the MTA for almost every month back to 1865. As fifteen stations appears to the minimum number needed to provide an accurate MTA, this suggests that the trend in Fig. 146.2 is reliable at least as far back as 1865.

Fig. 146.3: The number of station records included each month in the mean temperature anomaly (MTA) trend for Spain in Fig. 146.2.

If we next consider the change in temperature based on Berkeley Earth (BE) adjusted data we get the MTA data in Fig. 146.4 below. This again was determined by averaging each month the anomalies from the 83 longest stations and suggests that the climate was warming before 1980. This clearly contradicts the raw data in Fig. 146.2.

Fig. 146.4: Temperature trends for Spain based on Berkeley Earth adjusted data. The best fit linear trend line (in red) is for the period 1871-1980 and has a positive gradient of +0.32 ± 0.04°C/century.

But if we next compare the curves in Fig. 146.4 with those from the published Berkeley Earth (BE) version for Spain shown in Fig. 146.5 below, we see that there is excellent agreement between the two sets of data at least as far back as 1865. This indicates that the simple averaging of adjusted anomalies used to generate the BE MTA in Fig. 146.4 is as effective and accurate as the more complex gridding method used by Berkeley Earth in Fig. 146.5. In which case simple averaging should be just as effective and accurate in generating the MTA using raw unadjusted data in Fig. 146.2.

Fig. 146.5: The temperature trend for Spain since 1750 according to Berkeley Earth.

The differences between the MTA in Fig. 146.2 and the BE versions using adjusted data in Fig. 146.4  are instead mainly due to the data processing procedures used by Berkeley Earth. These include homogenization, gridding, Kriging and most significantly breakpoint adjustments. These lead to changes to the original temperature data, the magnitude of these adjustments being the difference in the MTA values seen in Fig. 146.2 and Fig. 146.4. The magnitudes of these adjustments are shown graphically in Fig. 146.6 below. The blue curve is the difference in MTA values between adjusted (Fig. 146.4) and unadjusted data (Fig. 146.2), while the orange curve is the contribution to those adjustments arising solely from breakpoint adjustments. Both are considerable with the former leading to an additional warming since 1865 of over 0.5°C.

Fig. 146.6: The contribution of Berkeley Earth (BE) adjustments to the anomaly data in Fig. 146.4 after smoothing with a 12-month moving average. The blue curve represents the total BE adjustments including those from homogenization. The linear best fit (red line) to these adjustments for the period 1871-1980 has a positive gradient of +0.506 ± 0.013 °C per century. The orange curve shows the contribution just from breakpoint adjustments.

The overall impact of the BE adjustments can be seen more clearly if we compare the 5-year averages for the raw data (from Fig. 146.2) and the BE adjusted data (from Fig. 146.4). This comparison is shown in Fig. 146.7 below. It clearly shows that the trend based on adjusted data (red curve) exhibits considerably more warming since 1870.

Fig. 146.7: The 5-year mean temperature change for Spain since 1820 based on the original raw data from Fig. 146.2 (in blue) and the Berkeley Earth adjusted data from Fig. 146.4 (in red).

Summary

The raw unadjusted temperature data for Spain clearly shows that the climate cooled by about 0.2°C from 1870 to 1980 (see Fig. 146.2)

In contrast, the BE adjusted data claims that the climate warmed by 0.3°C over the same period (see Fig. 146.4).

After 1980 the climate has clearly warmed by about 0.8°C (see Fig. 146.7).

The results presented here clearly disagree with those for Portugal in Post 145. However, as the MTA for Portugal before 1960 is based on data from less than five stations compared to two or three times as many stations for Spain, that would suggest that the Spain data is the more accurate.

Acronyms

BE = Berkeley Earth.

MRT = monthly reference temperature (see Post 47).

MTA = mean temperature anomaly.

Long station = a station with over 1200 months (100 years) of data before 2014.

Medium station = a station with over 480 months (40 years) of data before 2014.

List of all stations in Spain with links to their raw data files.

145: Portugal - temperature trends WARMING

The biggest problem in quantifying the extent of climate change in mainland Portugal is the lack of data before 1960. There are only two long stations with over 1200 months of data before 2014, the most significant of which is located in the capital, Lisbon (see here). This station appears to have recorded a temperature increase of over 2°C since 1850, yet the only other station with a comparable length of data in Coimbra shows a more modest increase with significant natural variability (see here) while data from Porto appears to suggest temperatures in the 19th century were warmer than today (see here). As Lisbon is very likely subject to some degree of urban heat island effect due to its size, this makes any accurate determination of the climate changes in Portugal before 1960 problematic.

In addition to the two long stations at Lisbon and Coimbra, there are also fifteen medium stations with over 480 months of data (for a full list of stations see here). The locations of all seventeen of these stations are shown on the map in Fig. 145.1 below. This suggests that they are fairly evenly distributed across the country with no part of the country being more than 70 km from a weather station. This is consistent with the station density of seventeen in 92,212 km² (the area of Portugal), or one in every 5424 km².

Fig. 145.1: The (approximate) locations of the 17 longest weather station records in Portugal. Those stations with a high warming trend are marked in red while those with a cooling or stable trend are marked in blue. Those denoted with squares are long stations with over 1200 months of data, while diamonds denote medium stations with more than 480 months of data.

In order to quantify the changes to the climate of Portugal the temperature anomalies for all stations with over 480 months of data before 2014 were determined and averaged. This was done using the usual method as outlined in Post 47 and involved first calculating the temperature anomaly each month for each station relative to its monthly reference temperature (MRT), and then averaging those anomalies to determine the mean temperature anomaly (MTA) for the whole country for each month. The MRTs for Portugal were calculated using the 30-year period from 1961 to 1990. The resulting MTA is shown as a time series in Fig. 145.2 below.

Fig. 145.2: The mean temperature change for Portugal since 1850 relative to the 1961-1990 monthly averages. The best fit is applied to the monthly mean data from 1871 to 1980 and has a positive gradient of +0.96 ± 0.08 °C per century.

The MTA in Fig. 145.2 clearly shows temperatures rising continuously from 1850 to 2013. However, as has already been pointed out, much of the increase before 1940 is due to the influence of the station in Lisbon. This is demonstrated by the graph below in Fig. 145.3 which shows only two or three stations contributing to the MTA before 1940. After 1960 the MTA is much more reliable as it is the result of averaging data from up to fifteen stations each month.

Fig. 145.3: The number of station records included each month in the mean temperature anomaly (MTA) trend for Portugal in Fig. 145.2.

If we next consider the change in temperature based on Berkeley Earth (BE) adjusted data we get the MTA data in Fig. 145.4 below. This again was determined by averaging each month the anomalies from the seventeen longest stations and suggests that the climate warmed slowly by 0.4°C in the one hundred years prior to 1980 and then warmed by another 0.7°C in the next thirty years.

Fig. 145.4: Temperature trends for Portugal based on Berkeley Earth adjusted data. The best fit linear trend line (in red) is for the period 1871-1980 and has a positive gradient of +0.41 ± 0.04°C/century.

Comparing the curves in Fig. 145.4 with those in the published Berkeley Earth (BE) version for Portugal shown in Fig. 145.5 below indicates that there is good agreement between the two sets of data. This demonstrates that the simple averaging of anomalies used to generate the BE MTA in Fig. 145.4 is as effective and accurate as the more complex gridding method used by Berkeley Earth in Fig. 145.5. In which case simple averaging should be just as effective and accurate in generating the MTA using raw unadjusted data in Fig. 145.2 even though the geographical distribution of stations is not completely homogeneous, as was shown in Fig. 145.1.

Fig. 145.5: The temperature trend for Portugal since 1750 according to Berkeley Earth.

While the temperature trends for the raw unadjusted data in Fig. 145.2 and the BE adjusted data in Fig. 145.4 look very similar, there are some significant differences. These can be seen more clearly by comparing the 5-year average of each dataset as shown in Fig. 145.6 below. This shows how the adjustments in Fig. 145.4 have altered the shape of the trend between 1900 and 2010 but not the overall total temperature change. The impact is to reduce the warming between 1900 and 1980 and to increase it thereafter, so making it look more like the classical 'hockey stick'.

Fig. 145.6: The 5-year mean temperature change for Portugal since 1850 based on the original raw data from Fig. 145.2 (in blue) and the Berkeley Earth adjusted data from Fig. 145.4 (in red).

The adjustments themselves can be calculated by subtracting the MTA values of the raw data in Fig. 145.2 from the adjusted values used in Fig. 145.4. The result is shown in Fig. 145.7 below.

Fig. 145.7: The contribution of Berkeley Earth (BE) adjustments to the anomaly data in Fig. 145.4 after smoothing with a 12-month moving average. The blue curve represents the total BE adjustments including those from homogenization. The linear best fit (red line) to these adjustments for the period 1871-1980 has a negative gradient of -0.238 ± 0.017 °C per century. The orange curve shows the contribution just from breakpoint adjustments.

The blue curve in Fig. 145.7 is the difference in MTA values between the adjusted data (Fig. 145.4) and the unadjusted data (Fig. 145.2), while the orange curve is the contribution to those adjustments arising solely from breakpoint adjustments. Additional contributions to the blue curve come from other adjustments based on techniques such as homogenization, gridding, Kriging and also any difference in MRT interval. Overall these adjustments appear to reduce the warming between 1900 and 1980 by about 0.3°C and then increase it by a similar amount from 1980 to 2010.

Summary

Both the raw data and the BE adjusted data appear to show that the climate of Portugal has warmed by about 1°C since 1890 (see Fig. 145.6). Most of this warming has occurred since 1980.

The effect of BE adjustments is to modify the shape of the trend from 1900 onwards rather than to increase or decrease the amount of overall warming (see Fig. 145.7).

While the trend since 1960 is incontrovertible as it is based on data from about fifteen different stations each month (see Fig. 145.3), the trend for earlier data is less so as it is based on only three sets of station data that profoundly contradict each other (Lisbon, Coimbra and Porto). However, data from Portugal's neighbour Spain in the next post may resolve this issue.

Acronyms

BE = Berkeley Earth.

MRT = monthly reference temperature (see Post 47).

MTA = mean temperature anomaly.

Long station = a station with over 1200 months (100 years) of data before 2014.

Medium station = a station with over 480 months (40 years) of data before 2014.

List of all stations in Portugal with links to their raw data files.

143: England - temperature trends WARMING

It is probably not surprising that England has more weather stations of note than Scotland. After all it has about ten times the population and almost twice the area. Yet the difference is not as great as one might imagine. For while Scotland has nine long stations with over 1200 months of data before 2014, England has only a slight advantage with ten stations. For medium stations with over 480 months of data the difference is greater with England having 55 compared to 13 in Scotland. There is, however, more clustering of stations in England as the map in Fig. 143.1 below shows.

Fig. 143.1: The (approximate) locations of the 65 longest weather station records in England. Those stations with a high warming trend between 1911 and 2010 are marked in red while those with a cooling or stable trend are marked in blue. Those denoted with squares are long stations with over 1200 months of data, while diamonds denote medium stations with more than 480 months of data.

In order to quantify the changes to the climate of England the temperature anomalies for all stations with over 480 months of data before 2014 were determined and averaged. This was done using the usual method as outlined in Post 47 and involved first calculating the temperature anomaly each month for each station relative to its monthly reference temperature (MRT), and then averaging those anomalies to determine the mean temperature anomaly (MTA) for the whole country for each month. The MRTs for England were calculated using the same 30-year period as for the UK in Post 141, namely from 1956-1985. 

The resulting MTA is shown as a time series in Fig. 143.2 below and clearly shows that temperatures rose slightly over 150 years up until 1975 before increasing more rapidly thereafter. In this respect the MTA data for England more resembles that of Great Britain (see Fig. 141.2 in Post 141) than it does that of Scotland (see Fig. 142.2 in Post 142) or Ireland (see Fig. 140.2 in Post 140).

Fig. 143.2: The mean temperature change for England since 1760 relative to the 1956-1985 monthly averages. The best fit is applied to the monthly mean data from 1826 to 1975 and has a positive gradient of +0.35 ± 0.08 °C per century.

The temperature trend for England was calculated using the usual method as outlined in Post 47 and involved first calculating the temperature anomaly each month for each station relative to its monthly reference temperature (MRT), and then averaging those anomalies to determine the mean temperature anomaly (MTA) for the whole country for each month. The graph in Fig. 143.3 below indicates how many stations were available each month in order to contribute to that month's MTA.

The MRTs for England were calculated using the same 30-year period as for the UK in Post 141, namely from 1956-1985. The resulting MTA is shown as a time series in Fig. 143.2 above and clearly shows that temperatures were slowly increasing for over 150 years up until 1975. Then at some point in the 1980s (probably in 1988) the mean temperature appears to increase abruptly by about 1°C. This is a phenomenon that has been seen in many other temperature trends across Europe.

Fig. 143.3: The number of station records included each month in the mean temperature anomaly (MTA) trend for England in Fig. 143.2.

If we next consider the change in temperature based on Berkeley Earth (BE) adjusted data we get the MTA data in Fig. 143.4 below. This again was determined by averaging each month the anomalies from the 65 longest stations and also suggests that the climate was warming slowly before 1980 but then warmed more strongly by over 1°C thereafter.

Fig. 143.4: Temperature trends for England based on Berkeley Earth adjusted data. The best fit linear trend line (in red) is for the period 1826-1975 and has a positive gradient of +0.27 ± 0.03°C/century.

The difference between the MTA based on raw unadjusted data (from Fig. 143.2) and the MTA based on BE adjusted data (from Fig. 143.4) is shown in Fig. 143.5 below. The blue curve in Fig. 143.5 is the difference in MTA values between the adjusted data (Fig. 143.4) and the unadjusted data (Fig. 143.2) and represents the total of all the data adjustments made including those from homogenization, gridding, Kriging and most significantly breakpoint adjustments. The orange curve is the contribution to those adjustments arising solely from breakpoint adjustments.

Fig. 143.5: The contribution of Berkeley Earth (BE) adjustments to the anomaly data in Fig. 143.4 after smoothing with a 12-month moving average. The blue curve represents the total BE adjustments including those from homogenization. The linear best fit (red line) to these adjustments for the period 1871-2010 has a small positive gradient of +0.003 ± 0.003 °C per century. The orange curve shows the contribution just from breakpoint adjustments.

The overall impact of any adjustments can perhaps be seen more clearly if we compare the 5-year averages for the raw unadjusted data and the BE adjusted data as is shown in Fig. 143.6 below. This shows that the two datasets agree almost perfectly from 1870 onwards while before 1870 the adjusted data implies the climate is more stable. It should be noted though that the MTA trends before 1870 are based on data from less than five stations and this drops to less than two before 1850 (see Fig. 143.3). This is one reason why the MTA exhibits more natural variability over the earlier period before 1870.

Fig. 143.6: The 5-year mean temperature change for England since 1760 based on the original raw data from Fig. 143.2 (in blue) and the Berkeley Earth adjusted data from Fig. 143.4 (in red).

Summary

What the raw data for England shows is that the climate was warming slowly for most of the period up to 1975 (see Fig. 143.2). This is similar to the trend seen previously for Great Britain (see Fig. 141.2 in Post 141), but is different from both Scotland (see Fig. 142.2 in Post 142) and Ireland (see Fig. 140.2 in Post 140) where little warming was seen in this period. This suggests that England is the dominant country in determining the overall climate of the UK but is also the outlier. But why?

The obvious answer is that England has a much greater population density and so experiences much more urban or surface heating from human activities (see Post 14, Post 29, Post 127 and Post 134). As I pointed out in Section (iv) of Post 127, the energy consumption of Greater London is sufficient to raise the local temperature by over 4°C.


Acronyms

BE = Berkeley Earth.

MRT = monthly reference temperature (see Post 47).

MTA = mean temperature anomaly.

Long station = a station with over 1200 months (100 years) of data before 2014.

Medium station = a station with over 480 months (40 years) of data before 2014.