ACTA UNIVERSITATIS CAROLINAE
AUC GEOGRAPHICA 55 2/2020
CHARLES UNIVERSITY • KAROLINUM PRESS AUC Geographica is licensed under a Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
© Charles University, 2020 ISSN 0300-5402 (Print) ISSN 2336-1980 (Online) Original Article 137 Symbolic and social boundaries of the integration of Russian immigrants in Czechia Ekaterina Ignatyeva*
Charles University, Faculty of Science, Department of Social Geography and Regional Development, Centre for Urban and Regional Research, Czechia * Corresponding author: [email protected]
ABSTRACT This article focuses on the integration process of Russian immigrants into Czech society. The integration of immigrants into Czech society is a key topic in the public debate as well as a political issue. Ukrainians, Slovaks, Vietnamese, and Russians are the most numerous groups within the half-million migrant population. Czechia is therefore predominantly attractive to non-EU immigrants. Representing highly educated and financially well-secured migrants who come as entire families, the Russians are distinct from other Eastern European immigrants. However, various factors hinder their integration. The article discusses the factors that shape symbolic and social boundaries in this integration process: (1) the development of Czech-Russian relationships that have been influ- enced by dramatic past events, (2) the representation of Russians in Czech media, (3) their specific socio-economic status, and (4) Czech immigration and integration policies. Negative experience, socio-economic inequalities, strict implementation of immigration policies towards third-country immigrants, and an unfavourable media discourse affect the attitudes of the majority toward the Russians and limit meaningful encounters.
KEYWORDS symbolic and social boundaries; immigration and integration policies; media discourse; Russian immigrants; Czechia
Received: 24 April 2019 Accepted: 20 May 2020 Published online: 3 July 2020
Ignatyeva, E. (2020): Symbolic and social boundaries of the integration of Russian immigrants in Czechia. AUC Geographica 55(2), 137–148 https://doi.org/10.14712/23361980.2020.10 © 2020 The Author. This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0). 138 Ekaterina Ignatyeva
1. Introduction countries, Czechia is the country with the most neg- ative attitudes towards immigrants. Concerning the With the rapid growth of immigrants in Czechia, the situation inside Czechia, the assessment of majority issue of their integration becomes important at the attitudes indicates that 64 percent of Czechs consid- institutional level through the implementation of er immigrants a problem (Public Opinion Research immigration policy, as well as at the social level due Centre 2017). Interestingly, this indicator has never to its impact on public opinion in society. During the since 2003. With regard to Russians, the Czech major- way in different areas of life within a host country. itydropped manifests below more fifty antipathy percent overthan thesympathy past ten towards years, Meanwhile,integration process,members immigrants of the majority try to may find not their hinder, own them, and, moreover, the attitude of Czechs towards but rather, they may aid and support societal new- the Russians has grown worse in recent years (Public comers. Hence, integration is considered a two-way Opinion Research Centre 2017). as well as a reciprocal process (Bosswick, Heckmann Pettigrew (1998) argues that not only social con- 2006), where both immigrants and the majority may tacts, but also cultural, economic, political factors, participate. On both sides, participation is necessarily - conditioned by mutual interaction that mostly mani- ence majority attitudes towards immigrants. Similar- fests in coexistence at the local level, where encoun- ly,media several discourse, different and factors dramatic that contribute social events to making influ ters between foreigners and members of the majority - society come about most often. sians, and therefore, slow down or block integration Prague increasingly attracts foreigners and belongs processes.boundaries In influence this article, Czech I discuss attitudes factors towards (historical, Rus to regions with their highest concentration (Czech discursive, socio-economic, and institutional) that capital of Czechia into a cosmopolitan city. The most of Russian immigrants into Czech society. numerousStatistical groupsOffice 2019). of up to Immigrants half a million transform migrants the in haveThe a significantaim of the articleimpact ison to the explore integration the context process of whole Czechia include Ukrainians (23%), Slovaks the integration process through the concept of sym- (21%), Vietnamese (11%), and Russians (7%)1. In bolic and social boundaries (Lamont, Molnár 2002) this article, I focus on Russian immigrants living in with an emphasis on (1) the historical development Czechia. Russians are particularly distinct from other of Czech-Russian relationships, which has undergone foreigners from Eastern Europe, especially in terms of both positive and negative events in the past; (2) rep- their composition. They represent a group of highly resentation of Russians in Czech media discourse; educated, economically well-resourced immigrants of (3) the socio-economic status of Russians that dis- working age who migrate as whole families (Drbohlav tinguish them from ‘others’; and (4) modes of Czech et al. 2010; Drbohlav, Janská 2004). immigration and integration policies that disadvan- According to the intergroup contact hypothesis tage the legal status of Russians as third-country (Pettigrew 1998; Allport 1954), in countries with - larger numbers of immigrants, there are better oppor- cantly contributes to understanding how symbolic tunities for encountering and social interaction that andimmigrants social boundaries over EU immigrants. between Russian This article immigrants signifi improve the majority attitudes towards immigrants through prejudice reduction. Even though today the everyday negotiation by limiting encounters that hin- number of immigrants in Czechia is 14 times higher derand the the integration Czech majority of immigrants. are created The andresults influence will be than in 1989 and is continuously growing, members useful predominantly for integration policymakers in of the Czech majority society, particularly of the older the setting, development, and improvement of inte- generation, still have negative attitudes towards for- gration activities for foreigners. eigners. In contrast, the younger generation comes First, the article discusses the theoretical frame- into contact with foreigners more so; and most of work, which is built around the concept of symbolic them have foreigners as work colleagues or school- and social boundaries (Lamont, Molnár 2002). After mates – a common circumstance today. Some of the reasons for this impact on the attitudes of the older the methodological point of view, in regard to research majority generation are (1) a lack of experience with onthis, the I briefly development introduce of Czech-Russianthe research methods. relationships, From international migration in Czechia, and thus, locals I study the available literature on this subject. Next, have not become used to the presence of ‘others’ yet, an analysis of the Czech press serves as a methodi- and (2) relatively recent dramatic events in Czech cal approach to investigating the discourse of media history, such as the Soviet military intervention and coverage concerning Russian immigrants. A descrip- occupation of Czechoslovakia, that are remembered by elders to this day. The results of the 2017 Euro- immigrants follows. Finally, there is an attempt to barometer survey show that, in a comparison of EU explaintion of the the specific institutional socio-economic status of Russians status of throughRussian an analysis of the development of Czech immigration 1 According to data from the Directorate of the Alien Police Service from 31 December 2018. are in the following section, and in conclusion, there and integration policies. The article’s main findings Integration of Russian immigrants in Czechia 139 is a summary of the research results and an explana- Caviedes 2015; Estrada et al. 2016) which limit mutu- tion of how the discussed factors impact the attitudes of the majority towards Russians and, consequently, Boundaries may become exaggerated by frames how they hinder the integration process of Russian basedal encounters on stereotyping and influence or which immigrant emphasise incorporation. the ‘oth- immigrants into Czech society. erness’ of the immigrant (Estrada et al. 2016). In
may lead to distinguishing society as ‘us’ and ‘them’. 2. Theoretical and methodological Asparticular, a result, definingthere is no who favourable ‘we’ are environmentand who ‘they’ where are frameworks fruitful social interactions can arise. Importantly, the immigrant voice is generally almost always missing 2.1 Symbolic and social boundaries in media coverage as well. Immigrants belong to the ‘others’ as well as poor, handicapped people, sexual In this article, I use the concept of symbolic and social and ethnic minorities, who exclude from the ‘normal boundaries by Lamont and Molnár (2002) to describe values of Western culture’ (Spivak 1999). Simultane- barriers in the process of Russian integration into ously immigrants belong to the discourse powerless Czech society as well as to evaluate the incorpora- salient social group, that less quoted or that have less tion process through the historically developed atti- to say (van Dijk 1988). As van Dijk (1988) asserts, people who are neglected in the press are people socio-economic status of Russian immigrants, and who are neglected in social life. According to van Dijk thetudes implementation of the majority, of mediaimmigration discourse, policies. the specificBy cat- (1988: 140) ‘there is not much difference between egorising objects, people, practices, time, and space, the free press of the Western countries or the more people draw symbolic boundaries that help individ- controlled press in most communist and many Third uals and social groups identify each other and, at the World countries’. Thus, by ignoring the opinion of same time, to determine their belonging within a par- ticular group during the negotiation of everyday life and, on the other hand, provides readers a biased (Lamont, Molnár 2002). Thus, symbolic boundaries viewimmigrants, of immigrants. the media deprives them of justification, contribute to the emergence and maintenance of a How immigrants are described or labelled in media dichotomy of ‘us’ and ‘them’ or the in-group and the coverage creates public opinion and policy outcomes out-group. (Sohoni, Sohoni 2014; Bleich et al. 2015; Estrada et Symbolic boundaries, however, could transform al. 2016); in other words, media discourse may have into social boundaries, which are broadly accepted an impact on immigration policies. In particular, rep- in society (Lamont, Molnár 2002). Unequal access resentations of restrictive immigration legislation to resources, their distribution, and different social in media contribute to the maintenance of symbolic opportunities shape social boundaries in a society, boundaries, especially when readers rely principally on the media discourse in place of the original legisla- (Lamont, Molnár 2002: 168). Social boundaries are tive resources or actual law texts (Estrada et al. 2016). oftenwhich institutionalised are ‘objectified (Bail forms 2008; of social Heizmann differences’ 2016) Nevertheless, I will discuss immigration legislation as and policy-oriented (Neumann, Moy 2018), determin- a separate dimension of boundary-making below. ing who in a society can or cannot access resources According to Estrada et al. (2016), the government - continues to be a key player in drawing boundaries nities (Lamont, Molnár 2002), and consolidate social among groups of immigrants and natives through inequalities.(material and Even non-material), though symbolic define and social social opportu bound- immigration law-making. Social boundaries are man- aries are closely connected and ‘should be viewed ifested by restrictive immigration policies that ‘draw equally as real’ (Lamont, Molnár 2002: 169), symbolic a sharper line between “us” and “them”, highlighting or brightening boundaries’ (Heizmann 2016: 1793). for the existence of social boundaries’ (Lamont, Mol- Immigration legislation provides immigrants with boundariesnár 2002: 169). are a ‘necessary but insufficient condition unequal rights in comparison to those that citizens Symbolic boundaries are formed discursively based enjoy. Immigrants may have different rights even within one country. For instance, in the European by representations from the outside. Today, the depic- context, there is a distinction between immigrants tionon our of immigrantssubjective perceptions in media discourse and are also plays influenced an active from the EU and non-EU countries and third countries role in shaping and reinforcing symbolic boundaries. (Heizmann, Böhnke 2018). This distinction between Media might represent immigrants in different ways – immigrants and natives results in the creation of positive, negative, or neutral. However, negative por- social boundaries. However, for third-country immi- trayals mostly contribute to the drawing of symbolic grants living in the EU, there is an extra layer of social boundaries. Frames, where immigrants are linked with crime or terrorism, represent them as a threat to as non-EU immigrants. Contrarily, integration policies the receiving society, creating negative attitudes in the endeavourboundaries tothat reduce define inequalities them as both and as foreignersachieve equi and- majority towards newcomers (Sohoni, Sohoni 2014; librium between immigrants and non-immigrants. 140 Ekaterina Ignatyeva
Here it is necessary to mention the importance of implementation of Czech immigration and integration the time aspect, respectively the time of arrival and policies that shape social boundaries by distinguish- ing among immigrants from EU and non-EU countries difference in the time when migrant arrived in a new and reinforces them at the institutional level. In the country,duration preciselytime that itinfluence depends boundary-making. on how many migrants It is a this country has? How more restrictive immigration policy towards newcomers? How many experienc- next2.2 Research section, Idesign briefly discuss the research methods. es locals have with living in a multicultural society and which attitudes they have towards foreigners? The article offers the results of three-part research. etc. Another important aspect is the duration time of Firstly, a historical development overview of the foreigner’s stay. First, as time goes on, it comes Czech-Russian relationships provides a summary of about a mutual habit and adaptation, on both part, previous detailed academic research and studies built of foreigners and locals. Second, and more important upon the memories of immigrant descendants, analy- is, with increasing a duration time in a new country, ses of documents from historical archives, as well as a the migrant acquires more and more rights in society, discussion of the current situation. and at the same time, has to overcome less and less - boundaries. covered through an analysis of Russian immigrant This article examines symbolic boundaries found coverageSecondly, in the the online representation version of Czech findings newspapers. were dis in the Czech media coverage of Russian immigrants. The media analysis focused on four of the most read The perception of historical experiences with Rus- Czech national daily newspapers, such as Lidové novi- sians as well as the current conditions of Czech-Rus- ny, Deník, Blesk, and Reflex, whose articles have been - under observation throughout an eight-year period tion of Russians. Moreover, Czechia, as a post-socialist (2011–2018). The important reason for choosing country,sian relations is in transition influence andthe Czechis still media’sdealing withrepresenta its her- these newspapers was the simple logical claim that - the most read sources have a wider target audience tions of Czech society, especially human mentality and which could be affected by reading its content. behaviouritage. Večerník patterns (2002) that points are not to the completed social transforma yet (Sýko- Newspapers articles were searched using the key- - words ‘Russian’ and ‘Russians’ on newspapers’ web- self with the shaping of symbolic boundaries between sites which archives were available online. The key- Russianra, Bouzarovski speaking 2012). immigrants Klvaňová and (2018) Czechs concerns through her an words ‘migrant’ and ‘immigrant’ were not using for examination of collective memory – in particular, the the purpose to reduce or eliminate the useless inci- cultural trauma caused by the Soviet military inter- dence of found articles that deal with immigrants in vention and occupation of Czechoslovakia in 1968. general. The selection of articles had to concern infor- She found that Czech society perceives immigrants mation about Russian immigrants as a main require- from the former Soviet Union as ‘colonizers’, devolves ment. Each found article was read and its information responsibility for past acts of occupational violence value detected. Articles that did not directly deal with to current immigrants, marks them as ‘others’, and Russian immigrants living in Czechia were excluded. keeps them at a distance due to the former dominance For instance, articles about tourism, sport, culture, and international politics. If the article contained trauma in Czech society contributes to the stigmatisa- relevant information, it was included in the database tionof the of Soviet newcomers regime through (Klvaňová past 2018). negative Thus, experienc cultural- of articles in the form of a simple Excel spreadsheet. Articles were categorised by topic in order to ascer- Russia’s current geopolitical behaviour (for example, tain what is written and spoken of in the Czech press itses andviolation draws of symbolic international boundaries law in the(Klvaňová annexation 2018). of in relation to Russian immigrants. The total number of found articles revealed whether the Russians are the majority attitudes towards its citizens. popular in the Czech media discourse or not. Subse- anotherSocial state’s boundaries territory) are primarily also significantly represented influences by an immigrant’s socio-economic status. Most Russians for each article: geographic level (national, regional, come to Czechia with a high-level socio-economic local),quently, type the followingof narrative characteristics with the numbers were identified of each status and try to keep or enhance it upon arrival. This type, and the possible participation of Russians. points to the differences between Russian immigrants Finally, the article offers (1) a discussion of the and other Russian-speaking foreigners from the East. socio-economic status of Russian immigrants based On the other hand, Russians belong to third-coun- try immigrants; thus, unlike immigrants from EU with emphasis on three features – income, educa- - tion,on 2011 and occupation census data – and from (2) Czech an assessment Statistical of Office their ditions, access to the labour market, and limited legal status via an evaluation of Czech immigration countries, they have much more difficult entry con and integration policies, which was based on an over- access to public health insurance, the right to vote, view of their general development from 1990 and any etc.).possibilities On that to account, use other it isstate necessary benefits to (for examine instance, the
relevant changes. The analysis focuses specifically Integration of Russian immigrants in Czechia 141 on immigrants from the third countries (including state organised an unprecedented humanitarian Russians) and their institutional position in society action for Russians (the so-called ‘Russian Action’) as a result of implemented immigration policies. To that set up and secured not only the basic needs of compare with EU immigrants, I discuss the limits and the immigrants but also provided temporary asylum disadvantages of Czech immigration policies towards and the opportunity to study and work. The respon- immigrants from non-EU countries that contribute to sive Czechoslovak policy emphasised the cultural the creation of social boundaries in the integration enrichment of Europe due to the presence of Russian process. refugees who carried with them the traditional Rus- Nevertheless, there are some limitations of the sian culture. The relationships between Czechs and Russians were mostly friendly at the beginning when development overview is limited by the lack of num- all these actions appeared as a temporary situation. bersstudy. of The academic first part studies which about deals Czech-Russian with the historical rela- However, when Russians realised the way back was tionships. The Russians are not much explored immi- grants’ groups in Czechia for instance in contrast to to their new home while the Czechs had to come to Ukrainians or Vietnamese minorities. The analysis termsdefinitely with closed, the permanent they had presenceto accustom of new themselves citizens of media discourse about Russian immigrants based and learn to live together. At this point, Russians start- on the research of newspapers only – it is the second ed to build their own community, establishing soci- limitation of this research. It may be better to include eties, organisations, and institutions, some of which the other media sources – TV, radio, and social net- still function today (the Slavic Library in Prague, for works on the Internet (e.g. Facebook, Twitter, Insta- instance). Sládek (1999) notes a disadvantage in the gram, etc.). On the other hand, this would lead to the existence of these societies: the hermetic closeness of overshoot an extent of the article and would make it the Russians to the host society represented a major possible to create another independent article, hence barrier to the process of their integration. Keeping the media discourse of this article is represented only their own Russian culture and traditions without an by analysis of the press. effort to assimilate to that of Czechs led to the shap- ing of symbolic boundaries by Russians themselves towards the receiving society. 3. Results Top among positive Czech-Russian relationships was a victory in World War II, and especially the lib- 3.1 Heritage of the past: Development eration of Czechoslovakia by the Soviet Army; Czech of Czech-Russian relationships people appreciated and were grateful towards Rus- sian as well as Soviet soldiers. After the end of the Sec- - ond World War in 1945, there was a certain euphoria tudes of Czechs towards Russian immigrants. Since theDifferent establishment historical eventsof the influenceCzechoslovak the current and Czech atti of the Soviet Union helped the Communist Party gain Republics, there have been three waves of Russian politicalin Czech-Russian power in relations. Czechoslovakia The significant and establish influence a immigration. Sládek (2010) gave names to these communist totality in 1948. Although according to waves according to the periods in which they took - place: the First Republic wave (1918–1948), the cessful and could not enter deep inside into the Czech socialist wave (1948–1989), and the post-socialist Večerník (2002) the communist regime was not suc wave (1989–present). The individual waves differ Czech population, and these negative experiences are from each other by volume and structure of immi- stillsociety, passed it significantly on from generation influenced to the generation. behaviour During of the grants, their motivation to move, and the attitudes of socialism, Czechs and Russians were connected by the Czech receiving society. Moreover, the migration a lot of common things, for instance, economic and history of Russians is characterised itself by alternat- cultural relations (Sládek 2010). Eventually, previous ing voluntary and involuntary migration periods. open and fruitful mutual everyday interaction and Czech-Russian relationships arose in 1918, during collaboration have changed into pragmatic economic - - dents, professors, scientists, and wealthier intelligen- the first Czechoslovak Republic, when Russian stu Sovietcooperation power. (Kratochvíl The subordination et al. 2006), and and the subsequent increasing Russia for political reasons (Sládek 2010) after the dependencely, it has moved on the towards Soviet theUnion negative led to influencethe gradual of Bolsheviktsia (Kopřivová coup. 2001) Later, were they forced were tojoined flee from by Russian Tsarist destruction and backwardness of the Czechoslovak soldiers who did not want to return after the First economy, which once was one of the most advanced World War, and therefore, stayed in Czechoslovakia. in post-war Europe. The largest concentration of Russian students and Unfortunately, in 1968, the Prague Spring was professors was in the capital of Czechoslovakia where followed by the tragic Soviet occupation of Warsaw most universities were located. Hence the reason Pact troops, which ruled Czechoslovakia under orders 1920s Prague was nicknamed the ‘Russian Oxford’ from Moscow from that point on. Gratitude to the (Sládek 2010). The newly established Czechoslovak Russians for the liberation of Czechoslovakia in the 142 Ekaterina Ignatyeva
Second World War was replaced by hatred towards any that had Russian roots due to the invasion. The the representations of Russians in the Czech media negative experiences of the communist totality peri- willThe beinfluence discussed of thesein the negativefollowing past section. experiences on od form the foundation of the symbolic boundaries towards Russians and everything of Russian genesis. 3.2 Media coverage of Russian immigrants - in the Czech press ma of communism in the majority society tends to According to Klvaňová (2018), this collective trau The media discourse analysis found 68 articles: 30 in responsibility onto contemporary Russian and/or Lidové noviny, 19 in Deník, 11 in Blesk, and 8 in Reflex. Russian-speakingreflect these events immigrants. into the present and to devolve Considering this fact, we can claim that Russian immi- After the events of 1989 took place, a third Russian grants are not the dominant object of Czech newspa- migration wave began which continues to this day per coverage. In terms of geographical scale (Fig. 1), (Sládek 2010). At present, negative attitudes persist 37% of the articles reported on Russians in the Czech within Czech-Russian relations. And there are rea- national context. The regional level was represented sons for that. First, today’s Russian immigrants inher- in 56% of the articles; most of them dealt with Rus- it a ‘collective guilt’ for the 1968 occupation. A sec- sian immigrants in Prague and Karlovy Vary, which is towards foreigners in general within Czech society. in these two cities. Written about less frequently were Thisond reasonnegative reflects reaction the triggers typical axenophobic connection attitudesbetween Russiansexplained in by other the significant Czech cities, concentration such as Brno, of Russians Hradec Králové, and Kunovice. In terms of the local level, Russians as agents of Putin or the Kremlin. only 7% of articles in the case of Prague focused on Russians and the mafia, espionage, and perception of threat of contemporary Russian imperialism, which - is intenselyKlvaňová perceived (2018) alsoand monitored points to in the Central potential and ice.Russian In a comparison immigrants betweenliving in thenational districts and oflocal Bubeneč, levels, Eastern Europe. The Russian occupation of Crimea in Nové Butovice, Zličín, Stodůlky, Letňany, and Vršov 2014 and the subsequent war in the Donbas region claiming that media discourse on immigrants is gen- evokes Czech memories of the 1968 invasion when erallythe same much findings more were nationalised discovered and by mostly Lawlor neglects (2015), every Czechoslovak was considered a victim and at the local context. However, in my analysis, the region- the same time, every Russian had been perceived al level emerged and was represented more often (56%) than others. deep-seeded post-communist collective trauma, contemporaryas a perpetrator Russia’s (Klvaňová geopolitical 2018). As behaviour a result ofin the the participation of Russians in the media debate. The questionThe next is howfinding often in the Russians discourse are analyses given the concerns oppor- attitudes of the Czech majority towards Russian immi- tunity and space to express themselves in the Czech grantsinternational due to anarena equation is largely of Russians reflected with by Russia. negative In press? The media discourse analysis found only 16% August, Czech people annually remember the tragic of all articles include the opinions of Russian entre- events of 1968 – the 50th anniversary of the Soviet preneurs, students, and journalists living in Czechia. invasion of Czechoslovakia passed in 2018, stirring a Most of these articles are informative narratives, great response in Czech public discourse. In the fol- where Russians descript why Czechia is attractive to lowing section, the issues of the portrayal of Russians them. For instance, the high achiever Russian entre- in Czech media will be discussed in more detail. Not- preneur presents: withstanding, it seems the more time passes, the more the events of 1968 have been made into a represen- Czechia is attractive to the Russians because there is a tative reminder of the post-communist trauma, thus close mentality … Western European societies are very sustaining the symbolic boundaries which are then reinforced and transferred onto the next generation. communism. New traditions and relationships are now beingmuch built based in onEastern tradition Europe. and Here, are notRussian influenced integrates by common. Today, Russian students, intellectuals, and into society faster. In Western Europe, no matter how wealthyThe first entrepreneurs and last waves choose of migration Czechia as have a migration much in much money Russian has, he will always be a foreigner - there. (Lidové noviny, 23.12.2018) lic period. This tradition, however, today contains a destination, similarly to the first Czechoslovak Repub 7% eventssignificant of the critical communism point, where period. positive The experiencescommunist 37% National regimefrom the and past the have Soviet been occupation influenced of byCzechoslovakia the negative 56% Regional - Local slovak nation, which is still one of the primary causes havebehind left the a significant existence trail of symbolic in the history boundaries of the Czecho in the process of negotiation between Czechs and Russians. Fig. 1 Percentage of all articles by geographical levels. Integration of Russian immigrants in Czechia 143
Media which offers the majority an immigrant’s (2015: 861–862) noted that ‘media outlets, especially - print media in Europe, are often associated with par- ation and improvement of an individual’s perception ticular political viewpoints’. As a result, an equation asperspective well as public directly opinion has aboutsignificant immigrants value inin general.the cre takes place between Russian immigrants living in Cze- Thus, an immigrant voice in media discourse can be chia and the political force in their motherland. These used as a tool in the elimination of symbolic and/or representations of Russians as spies create a partic- social boundaries as well as in the prevention of their ular perception within the majority society, which is formation. accompanied by feelings of suspicion and mistrust Found articles about Russians covered various towards the whole Russian immigrant population. narratives (Tab. 1), the most often discussed is their This is yet another example of how boundaries can be cohesion with the majority (21%), entrepreneurship created in the negotiation between the majority and (16%), activity on the real estate market (15%), and immigrants. It is also worth mentioning that almost criminality (15%). Articles dealing with everyday no attention was paid to the themes such as the Rus- negotiations between Russians with Czechs never - tion (1%), discrimination (1%), and emigration from - Russiasian financial (1%). crisis (3%), debts (1%), church restitu providedtions. Nonetheless, information it is about remarkable conflicts that or tensions,articles con but- rathercerning they the representedeveryday life non-conflictual of Russians in Czechia social interac sever- Tab. 1 Percentage of all articles about Russian immigrants al times (29%) made mention of the 1968 occupation. by narrative type. For instance, a resident from Carlsbad narrates about Narrative % of articles cohabitation Russians with locals, and about some- Cohesion with the majority, everyday life1 21 times provocative behaviour from the Russian side: Entrepreneurship 16
The Russians still claim that Carlsbad is beautiful and Real estate market activity 15 that they are very well here. And gradually the local Criminality 15 people from Carlsbad are getting used to them. Never- Stereotypes 12 theless, the Russians can dial locals almost reliably. For Participation and election preferences 8 of Russians abroad August 21st or when the Russians defeat the Czechs Espionage 6 example, when they (Russians) indulge in fireworks on in hockey. The locals here do not forgive that. Russian financial crisis 3 10.12.2011) (Deník, Debts 1 Church restitution 1 This demonstrates that media discourse maintains and reinforces symbolic boundaries which were ini- Discrimination 1 tially created by negative past experiences. In articles Emigration 1 about the housing or real estate market, Russians Source: own research were usually described as owners of expensive lux- 1 29% of these narrative type articles mentioned the 1968 occupation. ury apartments – mostly in Karlovy Vary but Prague as well. This title of the article is a very good example that captures the nature of fondness for luxury prop- collectivises Russians when covering other Rus- erty: ‘Russians love Czechia. Castles are cheaper here’ sian-speakingAnother interesting immigrants finding from former is that Soviet Czech Union print (Lidové noviny, 15.11.2011) Czech media discourse countries. Therefore, the narratives in these articles largely accuses Russian immigrants of raising prop- lead to a misrepresentation of reality, shaping media erty prices in such a way that others cannot afford bias towards Russians. Only three articles of this sort them. Less attention was paid to topics dealing with were discovered in the research, but it is noteworthy the presence of Russians in Czechia generally and that two of them deal with criminality. For instance, which related to outdated stereotypes from time to the title in Blesk newspaper introduces the crime as time (12%). committed by Russians: ‘Two Russians raided a mon- Only 6% of found articles connected the Russians ey truck: They neutralized the drivers with tear gas’, with espionage. This happened thanks to informa- but the content of the article tells us that they were tion within reports of the Security Information Ser- not actual Russians and instead were possibly Rus- vice of Czechia (BIS), which mostly have a political sian-speaking foreigners or even people who speak a context and are built on fears of Russia’s geopolitical language similar to Russian: ‘According to witnesses, behaviour in the international arena (for more infor- one of them was nervously and loudly telephoning in mation about the securitisation of European media Russian or similar language’ (Blesk, 22.6.2015). discourse see, for example, Caviedes 2015). 8% of the The question is how many people read only titles articles dealt with the Russian presidential elections, with this distorted reality, which therefore create a or rather they portrayed the participation and elec- negative public perception of Russians and shape the tion preferences of Russians living abroad. Bleich et al. majority attitudes towards them? Bleich et al. (2015) 144 Ekaterina Ignatyeva note that articles about immigrant individuals with a from other foreigners but also from the Czech major- criminal or economic threat context can lead to repre- ity. For instance, a survey provided by Schebelle et al. sentation of the whole group of immigrants as deeply (2015) found that Russians had the highest monthly problematic for society. income and lowest debt in comparison to Ukrainian
3.3 Socio-economic status of Russian immigrants in Czechia averageand Vietnamese and median immigrants gross monthly in Czechia. salary Vavrečkováof Russians inand Czechia Dobiášová exceeded (2015) the discovered average and that median in 2013 gross the - monthly salary of domestic inhabitants. ries is the socio-economic status (SES) of immigrants, whichA significant can be factorseen as in a thesource shaping of various of social inequalities. bounda status discussed above rank them among the most Russian immigrants living in Czechia are perceived The specific features of Russian socio-economic them from others, including the Czech majority. This 2010) as an elite group of immigrants. They differ socio-economicself-sufficient immigrants division of in ‘us’ Czechia, and ‘them’ distinguishing leads to a fromby some the immigrantsscholars (e.g. of Drbohlavother former et al. Soviet 2010; countries Janičko deepening of inequalities in society and, therefore, to in the following ways. First of all, there is a high lev- the creation of social boundaries to which Russians el of education among Russians. Most of them (43%) contribute themselves. carry a university degree, which exceeds the Czech national average by more than three times (Czech 3.4 Czech immigration policies that foreigners from Western Europe and the United The history of Czech immigration policy is thirty StatesStatistical also Officecome with2011). a high However, education it should level. beRussians added also care about the education of their children and pay close attention to it. Believing in the European years old. Some scholars (e.g. Barša, Baršová 2005;- education system, which in their opinion is better eralDrbohlav and restrictive et al. 2010; approaches Kušniráková, cyclically Čižinský alternated. 2011) and cheaper than Russian, they send their children distinguish five historical periods during which lib when Czechia did not regulate or limit entry to its ter- share of Russian students attending Czech univer- The first period took place between 1990 and 1996 sitiesto study from abroad. 4.9% This in 2007 is confirmed to 12.9% by in an 2019 increase (Czech of 222) approach towards all foreigners in Czech immi- grationritory. This policy ‘liberal enabled tolerance’ free entry (Barša, to the Baršová country 2005: but, high level of education provides better opportunities as Drbohlav et al. (2009: 46) note, ‘without a legal way inStatistical obtaining Office well-paying 2019). Russians jobs and aregeneral confident well-being. that a for permanent residence or naturalization, except for - marriage with a Czech citizen’. sians in Czechia is their high economic activity and Between 1996 and 1999, due to the deterioration typeThe of nextoccupation. specific Mostfeature of themthat characterises run a business Rus or of the socio-economic situation in Czechia as well as do a highly-skilled job in a position that corresponds increasing numbers of illegal foreign workers, Czech to their education level. As 2011 census data shows, immigration policy turned to a restrictive approach 16.2% of Russians are frequently employed in whole- through a tightening of the rules. At the same time, sale or retail; 13.1% in the real estate sector; almost Czechia became an EU candidate country and there- fore sought to adapt its entry requirements accord- technical activities; 8% are in manufacturing, 7.5% ingly. As a result, an amendment to Act No. 326/1999 in9% information of Russians and carry communication out qualified, technologies; scientific, and Coll., on the Residence of Foreign Nationals in the 6.4% carry out administrative and support activities. Czech Republic, came into force, which complicated This differs among Russians; for instance, Ukrainians, the lives of immigrants via the implementation of a who are mostly employed in Czechia as construction visa requirement before entering Czechia (Drbohlav workers or in manufactories, do lower-skilled jobs in comparison to the jobs they performed in their home for a permanent residence permit were permissible country (Drbohlav, Janská 2004). Furthermore, it is afteret al. 2009;ten years Kušniráková, of continuous Čižinský stay 2011). in CzechiaApplications and typical of Russians to create an immigrant economy, - which focuses on their compatriots or other Rus- ment/entrepreneurship, and humanitarian cases. In sian-speaking foreigners, allowing them to remain only for the purpose of family reunification, employ- relatively independent and, at the same time, limiting alised (Bail 2008; Heizmann2016) and politicised their interaction with the Czech majority. In this case, (Neumann,essence, this Moy was the2018) first social step in boundaries creating institution towards Russians initiate the shaping of social boundaries newcomers through the regulation of their entrance themselves through the creation of their own small and residence. world with strong inner ties that help them to sepa- The third period – 2000 to 2004 according to rate and close off from others, including Czechs. Several studies show that Russian immigrants in Czechia are distinguished by their incomes not just partialBarša and liberalisation. Baršová (2005) In 2004, or 2006Czechia in joinedkeeping the with EU Kušniráková and Čižinský (2011) – was marked by a Integration of Russian immigrants in Czechia 145 and the most important change in the Czech immi- to be employed as a pedagogical/academic worker gration policy came into force: a division of all immi- at a Czech university; or (c) have been posted to Cze- grants into foreigners from EU and non-EU countries. chia for the provision of services by his or her foreign Free entry, movement, and access to the Czech labour employer based in some other EU state. market was given to EU citizens and their family mem- Two years later, in 2016, the next amendment to the Alien Act was implemented. Two types of permit - residence newly came into existence: (1) a short- cantbers change– essential to bring benefits disadvantages distinct fromleading immigrants to institu- term visa for seasonal workers and (2) a long-term tionallycoming fromrooted non-EU social countries. boundaries, It was especially the first towards signifi residence permit for investment purposes. However, third-country immigrants (include Russians). there were also restrictive changes. For instance, a further restrictive step was taken in terms of acquir- ing a permanent residence permit by children in cas- wasThe called next a neoliberal period, from immigration 2005 to policy. 2007 A (Barša, conse- Baršováquence of 2005) economic or 2008 growth, (Kušniráková, Czechia Čižinskýwas faced 2011), with underage children. Such an unhappy implementation a labour shortage, and the solution to this problem ofes immigrationof family reunification policy leads when to the those division eligible of family were was the implementation of a green card as a way to members and seriously impacts Russians, who often attract a cheap labour force quickly (Drbohlav et al. move with the whole family. Therefore, one of the 2010). The green card project existed from 2007 to - 2009. As Drbohlav et al. (2010) note, Russians as well cure Russians is to send their child to be educated in as foreigners from Vietnam, Moldavia, and Mongo- most common migration strategies for financially-se lia could not apply for a green card. Given the high In 2018, a recent amendment of the Alien Act socio-economic status of Russians, it was likely that broughtCzechia first some and liberal then move changes to them. related to students they would not have been interested in this type of visa, which primarily targeted those doing low-skilled remain in Czechia for nine months for the purpose of and poorly paid jobs. Despite this fact, this deprival of seekingand scientists employment who, after or practising finishing theirentrepreneurship studies, may the opportunity to apply was the next brick in the wall of disadvantages and the drawing of social boundaries of foreign nationals into which it invests consider- towards immigrants from Russia. In the middle of this able– finally funds Czechia (a foreigner has considered may study the in Czech human language capital period, in 2006, the length of a continuous stay in Cze- free of charge and after graduation they now have free access to the national labour market). Today, at et al. 2009). least in small steps, the permeability of social bound- chiaThe was latest shortened period fromof Czech ten immigrationto five years policy(Drbohlav dis- aries has begun to be relaxed, even though it is only cussed in academic literature started in 2008 and con- for select groups of foreigners. Additionally, every for- eigner has an obligation to complete an integration call this time a neo-restrictive period during which entrytinues and to thisstay day.requirements Kušniráková have and been Čižinský tightened (2011) and arrival in Czechia. entrance to Czechia for some nationals, such as Mon- and adaptation course during the first year after their- golians, Moldavians, Thais, Ukrainians, and Vietnam- ty is related to the division in voting rights of immi- ese, has been temporarily cancelled. grantsAnother who significantare living in type Czechia. of institutional In accordance inequali with In 2011, the next amendment brought new, stricter the Election Acts (No. 491/2001 Coll., No. 62/2003 application requirements. Every applicant must now Coll.), EU citizens with a permanent residence permit provide proof of secure accommodation, health insur- have the right to vote in municipal elections as well ance, and funds for their stay in the country. More- as the European Parliament elections. On the oth- over, personal attendance when applying, as well as er hand, citizens of third countries do not have any voting rights, except holders of Czech citizenship. process longer and more apprehensive than before. It means not all immigrants who live long-term in Finally,an interview the newly with implemented a police officer, permit makes card the with whole bio- metric data increases expenditures for immigrants. living conditions in the receiving country. The govern- So as to attract a high-skilled labour force to Cze- mentCzechia enables have anforeigners equal possibility to come, work, of influencing and live in their the chia, a blue card was later implemented and, in 2014, country while limiting opportunities and withholding an employee card for all types of labour (includ- the right to change and enhance them until citizen- ing low-skilled) replaced the previous green card. ship is obtained. According to the Ministry of Interior, there are two Even though Czechia does not have a self-standing modes of the employee card: (1) dual, which contains integration law, its integration strategy has recently residence and employment permits, and (2) non-dual, been intensively developing. The Ministry of the Inte- which offers a residence permit only – for foreign- - ers with free access to the Czech labour market who tion policy in 2000 when the Alien Act came in force. (a) have obtained secondary, tertiary, tertiary profes- Focusingrior drafted on theequal first opportunities version of an and immigrant non-discrimina integra- sional, or university education in Czechia; (b) wish tion, immigrant integration policy struggles towards 146 Ekaterina Ignatyeva similar rights for long-term residents as those means they belong to a powerless salient social group received by Czech citizens. In 2006, 2011, and 2016, (van Dijk 1988) or in general to the ‘others’ (Spivak there were fundamental updates in the integration policy which paid attention primarily to social inter- Russians as a general term, pertaining not only to for- action between the majority population and immi- eigners1999; Estrada from Russia et al. 2016). but also And other finally, Russian-speaking the media use grants by supporting good relations in everyday life immigrants; such a generalisation can lead to the dis- negotiation. Since 2010, the Czech government has tortion of reality. annually published an action plan which contains pri- In summary, negative experiences from the past feed the present majority perception through media of the action plan in the previous year. representations, which help identify and determine orities, goals and means, and reports on the fulfilment in-group and out-group members. These practices contribute to the ‘otherness’ (Spivak 1999; Estrada et 4. Conclusion people in society as ‘us’ and ‘them’ and draw symbol- This article has investigated the symbolic and social ical. boundaries2016) of Russians (Lamont, and Molnár enable 2002). the classification Linking with of boundaries that hinder the integration process of a crime as well as a constant return to the past and Russian immigrants in Czechia. This study contributes a reminder of the events of 1968 classify Russians in to the literature on boundary-making, and particu- Czech media as a threat. Thus, the power of media larly the case of Russian immigrants, in Czechia. The and strengthening of symbolic boundaries in major- itydiscourse attitudes has towards a significant Russians impact that might on the lead creation to the werefindings (1) of the this development study are summarised of Czech-Russian henceforth. relation The- reinforcement of negative attitudes (Sohoni, Sohoni significantships, (2) the factors depiction that create of Russians boundaries in Czech discussed media 2014; Caviedes 2015; Estrada et al. 2016), limit their discourse, (3) the socio-economic status of Russian encounters, and therefore, obstruct their fruitful immigrants, and (4) the implementation of immigra- integration. tion and integration policies. The higher education of Russians provides them The positive attitudes towards the presence of the opportunity to get high-paying employment, and Russians in Czechia as well as Czech-Russian collab- therefore, a better position on the social ladder. This oration were disturbed by the negative experiences combination of high education level and high econom- during the communist regime. In particular, the inva- ic activity ranks Russians as an elite immigrant group sion and military occupation of Czechoslovakia by the that differentiates them from others and the Czech Soviet Army in 1968 left a dramatic footprint in the majority as well, leading to the consolidation of social souls of the Czech people. Unfortunately, this collec- inequalities. It is a next example of the ‘otherness’ tive trauma persists in the minds of elderly Czechs to (Spivak 1999; Estrada et al. 2016), but in this case to - which Russians contribute themselves. Thus, they are temporary Russia’s geopolitical behaviour saturates the presentnegative day. majority As Klvaňová attitudes (2018) due to indicated, a fear of conhis- However, Russians, as third-country immigrants, tory repeating itself. Collective trauma turns to col- haveable toan create unequal social legal boundaries status in comparisonby their specificity. to for- lective guilt, for which current Russian immigrants eigners from EU countries. The results of the research are deemed responsible in Czechia. Memories of the further show that the restrictive implementation of negative historical events, as well as the majority per- Czech immigration policy towards Russians, as well as other third-country immigrants, regulates their depiction of Russians. entrance, limits their access to the labour market, and ceptionAccording of the to current the media status, discourse can influence analysis, the ‘Russian media immigrants’ are not a popular topic in the Czech press. of previous studies (e.g. Heizmann, Böhnke 2018), a However, the frequency is not so important in con- detaileddefines them investigation as ‘they’ orof ‘others’.Czech immigration As with the findingspolicies trast to the narrative context, which created an overall - impression on readers. In this regard, based on the leged EU immigrants in comparison with immigrants most common narratives in the media, we can com- fromconfirms third that countries there is who an emphasisare legally on disadvantaged. legally privi pile a typical image of a Russian immigrant in Czechia. In sum, on the one hand, Russian immigrants create Probably he will be an entrepreneur who operates in the real estate market and owns a large number of - luxury apartments or castles, he is maybe involved in tunitiessocial boundaries in society (Lamont, themselves Molnár through 2002). their On thespecific oth- crime, or even he is a Russian spy; he gets along with erself-sufficient hand, other statussources that of inequalitiesdefines their that social lead oppor to the the majority without any problems, even though he creation of social boundaries are the implementation is sometimes able to provoke them, for example on of immigration policies that are broadly accepted in the anniversary of the occupation of Czechoslovakia. society (Lamont, Molnár 2002). A further result indicated that Russians rarely receive Based on this study I suggest some concluding space in the Czech media to express their opinions, it remarks that would help to make the process of Integration of Russian immigrants in Czechia 147 integration easier for foreigners and fruitful for soci- and Migration Studies 41, 857–873, https://doi.org ety in whole. First, members of the Czech majority /10.1080/1369183X.2014.1002197. society should intensify social contact with foreign- Bosswick, W., Heckmann, F. (2006): Integration of migrants: ers living in Czechia and behave towards newcom- Contribution of local and regional authorities. European ers with understanding and tolerance, regardless of Foundation for the Improvement of Living and Working Conditions, Dublin. their country of origin. Current immigrants cannot be Caviedes, A. (2015): An Emerging ‘European’ News responsible for past mistakes made by their predeces- Portrayal of Immigration? Journal of Ethnic and sors. Moreover, in everyday negotiations and attitudes Migration Studies 41, 897–917, https://doi.org/10.1080 towards foreigners, it is incorrect to connect immi- /1369183X.2014.1002199. grants with political affairs taking place in their ori- gin countries as they have left their motherland and Republic. Prague. live abroad. Second, the media should pay attention Drbohlav,Czech Statistical D., Janská, Office E. (2004):(2019): CurrentForeigners Ukrainian in the Czech and to how they represent immigrants who live among us Russian Migration to the Czech Republic: Mutual and how that may impact their lives. At the same time, Similarities and Differences. In: Górny, A., Ruspini, P. foreigners should be offered more avenues in which (eds.): Migration in the New Europe, East-West Revisited. to express their opinions, and interest in them should Palgrave Macmillan, London, 49–64. be shown. In the process of building attitudes towards foreigners, readers should rely on their own experi- Drbohlav,on its way D., Lachmanová-Medová,from emigration to immigration L., Čermák, country. Z., Janská, IDEA E., ences, not a mediated perception of news served by WorkingČermáková, Papers. D., Dzúrová, D. (2009): The Czech Republic: quick changes in legislation, Czech immigration andthe media.integration Third, policies although might it is adapt difficult to allto imagineforeign- EuropeanDrbohlav, D.Comission et al. (2010): (2017): Migrace Special a (i)migranti Eurobarometer v Česku. 469 – ers staying on its territory and intending to remain IntegrationKdo jsme, odkud of Immigrants přicházíme, in thekam European jdeme? Praha, Union. SLON. here. /index.cfm/survey/getsurveydetail/instruments /special/surveyky/2169.https://ec.europa.eu/commfrontoffice/publicopinion Acknowledgements Heizmann, B. (2016): Symbolic boundaries, incorporation policies, and anti-immigrant attitudes: what drives This research was supported by Grant Agency of exclusionary policy preferences? Ethnic and Racial Studies 39(10), 1791–1811, http://dx.doi.org/10.1080 Charles University (grant project ‘Multilevel separa- /01419870.2015.1124128. tion of Russian immigrants in the post-socialist city’, Heizmann, B., Böhnke, P. (2018): Immigrant life satisfaction - in Europe: the role of social and symbolic boundaries. tre program UNCE/HUM/018. I would like to thank Journal of Ethnic and Migration Studies 45(7), thereg. two № 192418) anonymous and reviewersCharles University whose comments/sug Research Cen- 1027–1050, https://doi.org/10.1080/1369183X gestions helped improve this paper. .2018.1438252.
Bachelor’s thesis. Charles University in Prague, References Janíčko,Department M. (2010): of Sociology. Pražští Rusové (The Prague Russians).
Act No. 326/1999 Coll., on the Residence of Foreign boundaries and cultural trauma in immigration from Nationals in the Czech Republic. Klvaňová,the former R. (2018): Soviet ‘TheUnion Russians to the Czech are back’: Republic. Symbolic Act No. 491/2001 Coll., on Elections to Municipal Councils Ethnicities 0(0), 1–20. https://doi.org/10.1177 and on Amendments to Some Acts. /1468796817752740. Act No. 62/2003 Coll., on Elections to the European Parliament and on Amendments to Some Acts. Allport, G. W. (1954): The Nature of Prejudice. Reading, Kopřivová, A. (2001): Střediska ruského emigrantského MA: Addison-Wesley. života v Praze (1921–1952). Praha: Národní knihovna rhetoricalČR, Slovanská action knihovna. and the idea of otherness: The Czech Boundaries against Immigrants in Europe. American Kratochvíl,Republic P., and Cibulková, Russia. Communist P., Beneš, V. and(2006): Post-Communist Foreign policy, Bail,Sociological C. A. (2008): Review The Configuration 73, 37–59, https://doi.org of Symbolic Studies 39(4), 497–511, https://doi.org/10.1016 /10.1177/000312240807300103. /j.postcomstud.2006.09.002.
Baršová, A., Barša, P. (2005): Přistěhovalectví a liberální Kušniráková, T., Čižinský, P. (2011): Dvacet let české politologickýstát. Imigrační ústav, a integrační Brno, 205–255. politiky v USA, západní migrační politiky: Liberální, restriktivní, anebo ještě Evropě a Česku. Masarykova univerzita, Mezinárodní Lamont,jiná? Geografie M., Molnár, 116(4), V. (2002): 497–517, The Study https://doi.org of Boundaries in the/10.37040/geografie2011116040497. Social Sciences. The Annual Review of Sociology 28, Barvíková, H. (2005): Ruští emigranti v Československu 167–195, https://doi.org/10.1146/annurev.soc.28 Bleich,v letech E., Bloemraad, 1918–1938. I., In: de ExilGraauw, v Praze E. (2015): a Československu Migrants, .110601.141107. Minorities1918–1938. and Pražská the Media: edice, Information, Praha, 9–21. Representation, Lawlor, A. (2015): Local and National Accounts of and Participation in the Public Sphere. Journal of Ethnic Immigration Framing in a Cross-national Perspective. 148 Ekaterina Ignatyeva
Journal of Ethnic and Migration Studies 41, 918–941, The Sociological Quarterly 55(1), 49–71, https:// https://doi.org/10.1080/1369183X.2014.1001625. doi.org/10.1111/tsq.12039. Neumann, R., Moy, P. (2018): You’re (Not) Welcome: The Spivak, G. S. (1999): A Critique of Postcolonial Reason: Impact of Symbolic Boundaries, Intergroup Contact, Toward a History of the Vanishing Present. Cambridge: and Experiences with Discrimination on Immigration Harvard University Press. Attitudes. American Behavioral Scientist 62(4), 458–477, Sýkora, L., Bouzarovski, S. (2012): Multiple https://doi.org/10.1177/0002764218760370. transformations: Conceptualising the post-communist Pettigrew, T. F. (1998): Intergroup contact theory. Annual urban transition. Urban Studies 49(1), 43–60. Review of Psychology, 49(1), 65–85, https://doi.org van Dijk, T. A. (1988): News analysis. Case studies of /10.1146/annurev.psych.49.1.65. international and national news in the press. Lawrence Erlbaum Associates, Publishers, Hillsdale, New Jersey.
Sládek, Z. (1999): České prostředí a ruská emigrace (1918–1938). In: Běloševská, L.: Duchovní proudy ruské Vavrečková,Russian community J., Dobiášová, into K. majority (2015): Začleněnísociety). Výzkumný ruské a ukrajinské emigrace v Československé republice komunity do většinové společnosti (Incorporation of (1919–1939). Slovanský ústav AV ČR, Praha, 7–46. Sládek, K. (2010): Ruská diaspora v České republice: ústav práce a sociálních věcí, v.v.i., Univerzita Karlova Sohoni,Sociální, D., Sohoni, politická T. aW. religiózní P. (2014): variabilita Perceptions ruské of migrace, v Praze, Fakulta sociálních věd. ImmigrantPavel Mervart, Criminality: Červený CrimeKostelec. and Social Boundaries. Večerník, J. (2002): Výzkum společenské transformace a česká sociologie (Transition Studies and Czech sociology). Sociologický časopis 38(1–2), 55–77. Original Article 149
Flood susceptibility mapping in Erythropotamos river basin with the aid of Remote Sensing and GIS Christos Domakinis1,*, Antonios Mouratidis1, Kostas Voudouris1, Theodore Astaras1, Maria Chara Karypidou2
1 Aristotle University of Thessaloniki, Department of Environmental and Physical Geography, Greece 2 Aristotle University of Thessaloniki, Department of Meteorology and Climatology, Greece * Corresponding author: [email protected]
ABSTRACT Erythropotamos is a tributary of river Evros and during the last decade its drainage basin flooded many times, causing extensive damage on properties. In order to assess flood susceptibility in the aforementioned study area, the inundated areas of floods that occurred in 2010, 2017 and 2018 were initially delineated with the aid of SAR (Synthetic Aperture Radar) imagery by applying an established flood delineation methodology. Subsequently, flood susceptibility mapping was conducted for the study area by apply- ing the Analytical Hierarchy Process (AHP). Topographical, hydrological and meteorological factors were used and each one of them was classified into three (3) flood susceptibility categories (low, medium and high). The determination of the importance for each factor over the others, which is the main objective of this research, was decided according to the proportion of the 2010 inundated area, captured by ENVISAT/ASAR imagery, which intersected with each factor’s high susceptibility class. Finally, the resulting flood susceptibility map was validated according with the inundated areas of the 2017 and 2018 flood events, captured by SENTINEL – 1 A/B imagery, indicating that approximately 60% of both of these areas intersected with the map’s high susceptibility zone.
KEYWORDS GIS; Susceptibility mapping; Analytical Hierarchy Process (AHP); floods; remote sensing
Received: 15 May 2019 Accepted: 7 May 2020 Published online: 31 July 2020
Domakinis, C., Mouratidis, A., Voudouris, K., Astaras, T., Karypidou, M. C. (2020): Flood susceptibility mapping in Erythropotamos river basin with the aid of Remote Sensing and GIS. AUC Geographica 55(2), 149–164 https://doi.org/10.14712/23361980.2020.11 © 2020 The Authors. This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0). 150 Christos Domakinis, Antonios Mouratidis, Kostas Voudouris, Theodore Astaras, Maria Chara Karypidou
1. Introduction to become more effective. Inundation and suscepti- bility mapping are among the main procedures that Floods can potentially cause fatalities, displacement of people and damage to the environment, to severely its goals. compromise economic development and to under- floodRegarding hazard assessmentinundation follows,mapping, in SARorder systems to achieve are mine the economic activities of every community that particularly suitable, thanks to the synoptic view, the suffers the effects of these destructive environmental capability to operate in almost all-weather conditions hazards (Patrikaki et al. 2018; Zhong et al. 2018; Birk- and during both day-time and night-time, as well as holz et al. 2014; Mouratidis and Sarti 2013; Yésou et the sensitivity of the microwave radiation to water al. 2013; Astaras et al. 2011). During the last decade, (Pierdicca et al. 2013). Furthermore, various methods such phenomena have also plagued Greece, with their majority occurring in the eastern part of the region of water from SAR data. Change detection highlights the Thrace (Kazakis et al. 2015; Nikolaidou et al. 2015; temporalhave been changesused within in theland literature cover by to comparing delineate flood the Mouratidis 2011; Mouratidis et al. 2011). Most such cases are attributed to the river Evros, which is the Psomiadis 2016; Schlaffer et al. 2015). The differ- natural borderline between Greece and Turkey, and, enceflood between scene to thea previous images candry be image combined (Li et withal. 2018; oth- along with its tributaries, has burst its banks on sev- er image segmentation techniques to identify areas eral occasions during the aforementioned time peri- producing an unusually low backscatter response, od. Erythropotamos is one of Evros’ tributaries and when compared to the single image methodologies been observed within its river basin, there is a lack (Matgenimproving et al. the 2011). reliability of the flood delineation although, in many occasions, flood phenomena have - catchment. ing to the current literature, the contemporary of floodDuring hazard the last assessment few decades, studies advances referring in remote to this trendConcerning involves flood mostly susceptibility the creation mapping, of ensemble accord models, which are based on the combination of sensing and GIS have helped flood hazard assessment
Fig. 1 a) Location of the study area (drainage basin of Erythropotamos), b) Geological formations within the drainage basin of Erythropotamos river (CoG 1989, I.G.M.E. 2002), c) Spatial distribution of elevation within the catchment of Erythropotamos river (E.E.A. 2017) and d) Spatial distribution of slope angle values within the drainage basin of Erythropotamos river. Flood susceptibility mapping in Erythropotamos river basin 151 different data-driven (machine learning), statistical or Tab. 1 Distribution of elevation into categories according to Dikau’s multi-criteria methods. This approach aims in achiev- classification (Dikau 1989). ing higher accuracy of the delineated susceptibility Elevation Description Area (km2) Percent (%) <150 Lowland 429.6 26.543 mapping methodologies that employ a single meth- odzones, or model in comparison (Costache etwith al. 2020;the flood Kanani-Sadat susceptibility et al. 150–600 Hilly 969.1 59.876 2019; Wang et al. 2019; Khosravi et al. 2016). How- 600–900 Semi-mountainous 197.3 12.190 ever, there is a plethora of methodologies that can be 900> Mountainous 22.5 1.390 - tistical and data-driven approaches (Ettinger et al. 2016;used in Nandi flood et susceptibility al. 2016; Tehrany assessment, et al. 2015; such Tehrany as sta Tab. 2 Slope angle categorization within the study area according et al. 2014; Pulvirenti et al. 2011). Among them, the to Demek’s classification (Demek 1972). analytic hierarchy process (AHP) (Kazakis et al. 2015; Slope Angle Description Area (km2) Percent (%) Stefanidis and Stathis 2013) is considered as the most (°) widely used and, because of its simplicity, continues 0–2 Plain to slightly sloping 147.5 9.11 Seejata et al. 2018; Tang et al. 2018). Additionally, this 2–5 Gently inclined 390.3 24.11 methodologyto be popular haseven proved in recent many works times (Lyu that et it al.can 2018; han- 5–15 Strongly inclined 864.9 53.44 dle sparse or poor quality data and that it can operate 15–35 Steep 215.7 13.33
Wang et al. 2011). >35 Precipitous 0.1 0.01 efficientlyThe main in regionalaim of this studies research (Chen iset toal. 2015,introduce 2013; a methodology that deals with the subjectivity that involves the determination of the hierarchy of factors 2. Location of the study area
To this end, the proposed methodology employed the Erythropotamos is a tributary of Evros River, which resultsin flood of susceptibility SAR-based inundation mapping with mapping the use that of delin AHP.- is the longest river that runs solely in the interior
- of 1,618.5 km2. The largest part of its river basin eatedtors was the determined flood extent according of the 2010 to flood.the proportion Specifically, of belongsof the Balkans, to Greece and and its particularly catchment tocovers the geographic an extent the aforementionedhierarchy between inundated the flood areas susceptibility that intersected fac region of Thrace in Northern Greece, while the rest with each factor’s high susceptibility class. of its drainage basin belongs to Bulgaria (Figure 1a).
Fig. 2 Spatial distribution of land cover in the catchment of Erythropotamos river (Copernicus 2017). 152 Christos Domakinis, Antonios Mouratidis, Kostas Voudouris, Theodore Astaras, Maria Chara Karypidou
Fig. 3 General flowchart of the applied methodology.
Regarding administrative distribution within nation- - al borders, the Greek part of Erythropotamos’ river nicus 2020) data layer, is chronologically closer to the basin belongs to the Prefectures of Evros and Rhodopi. which, instead of the Corine Land Cover 2018 (Coper- The drainage basin of Erythropotamos River potamos River is dominated by forests and semi nat- belongs to both the Circum – Rhodope geotectonic uralgauged areas 2010 (Figure flood 2).event, Along the with catchment agricultural of Erythro areas, zone and the Rhodope massif. The geological for- mation that covers the largest part of the study area consists of orthogneissess and augen gneisses (Fig- these two land cover categories occupy approximately ure 1b). the 98% of the total extent of the study area (Table 3). Elevation in the drainage basin of Erythropotamos 3. Materials and methodology river ranges from 16 m to about 1,258 m above mean The materials that were used and the methodology area can be described as hilly according to Dikau’s that was followed in order to achieve the aims of this sea level (M.S.L.), and the largest part of the study Additionally, the spatial distribution of slope angle inundation mapping with the use of SAR images, while valuesclassification within (Dikauthe study 1989) area, (Figure indicates 1c and that Table most 1). of study can be divided into two parts. The first involves its terrain belongs to the strongly inclined category isthe given second in Figure part is 3. concerned with flood susceptibility angles (Demek 1972) (Figure 1d and Table 2). mapping with the use of AHP. The general flowchart (5°–15°)Finally, according based on toinformation Demek’s classification provided by theof slope data 3.1 Inundation mapping with SAR imagery
Eleven ENVISAT/ASAR and twenty seven SENTI- Tab.layer 3 Distributionof Corine Landof Land Cover within 2012 the (Copernicus river basin of 2017), - Erythropotamos according to Corine Land Cover 2012 (Copernicus 2017). NEL – 1 A/B images of VV (Vertical transmit – Verti cal receive) polarization were used to map the flood Land Cover Area (km2) Percent (%) - isticsextents appear of the on February Table 4. 2010, April 2017 and March Artificial surfaces 18.4 1.137 2018The flood aforementioned events. Their SARdetailed images product were character pre-pro- Agricultural areas 586.2 36.219 cessed with the aid of ESA’s SAR satellite image anal- Forest and semi natural areas 1007.9 62.274 ysis software SNAP (Sentinel Application Platform). Wetlands 0.1 0.006 - Water bodies 5.9 0.365 Initially, they were calibrated to σ° backscatter coef ficient values and despeckled using a 3 × 3 Gamma map filter. Regarding the co-registration step, the SAR Flood susceptibility mapping in Erythropotamos river basin 153
Tab. 4 Product information of ENVISAT/ASAR and SENTINEL – 1 A/B imagery.
Satellite ENVISAT SENTINEL – 1 andcategorizes non-water as flooded land cover only classes. areas that On arethe temporarilyother hand, Flood image: Flood images: 18/4/2017 NDFVIcovered was by water,used in excluding order to permanent detect shallow water water bodies in 16/2/2010 (1) & 26/3/2018 (2) low vegetation. Dates Reference images: Reference images: According to Cian et al. (2018) NDFI values that 5/8/2008 to 8/10/2016 to are greater than 0.7 and NDFVI values that are great- 27/4/2010 (10) 28/9/2018 (25) er than 0.75 can be used to delineate inundated areas Spatial 11.1 m × 11.1 m 8.8 m × 8.8 m in open land and in low vegetation respectively. How- Resolution - Pass Ascending Descending cessing according to the following criteria: Mode N/A IW ever, the resulting flooded areas require further pro Type N/A GRD 1. becauseFlooded theyareas can with be extent considered smaller as thanspurious the size (Cian of Level 1 1 et10 al. pixels 2018). in NDFI and NDFVI images were excluded Polarization VV VV Relative Orbit 14 109 < 0.015), which correspond to permanent water 2. Pixels with σο(mean) values less valuesthan 0.015 greater (σο(mean) than - images were co-registered with the use of EU-DEM, sistentlybodies, and decrease pixels with their σ ο(min)backscatter during the which is the Digital Surface Model (DSM) of European 0.03 (σο(min) > 0.03) that represent pixels that con Environment Agency (EEA) member and cooperating - - flood, indicating that something happened, but not- nated by the sensors. It is a hybrid product based on enoughdation maps to reach (Cian a σ etο(min) al. 2018). value typical of water pix countriesSRTM and that ASTER represents GDEM datathe first fused surface by a asweighted illumi Moreover,els, have to the be filteredadverse out weather from theconditions resulting during inun averaging approach (EEA 2017). Its horizontal spa- optical imagery and aerial vehicles from capturing the while its absolute and relative vertical accuracy are the 2010, 2017 and 2018 floods prevented satellite- 3.6tial mresolution and 5.3 m, is respectively1 arc second (Mouratidis (approximately et al. 2019).25 m), tion of their SAR-based inundation mapping results The Change Detection And Thresholding (CDAT) wasextents not offeasible. the corresponding floods and thus valida methodology by Cian et al. (2018), based on the work - 3.2 Flood susceptibility mapping with the use of AHP and with the aid of GIS and satellite imagery events.of Long This et al. procedure (2014), was involved applied the in calculation order to delinof the eate the inundated areas of the aforementioned flood The compilation of the susceptibility map can be achieved by conducting multi-criteria analysis (MCA), (NDFVI),Normalized which Difference are based Flood on the Index multi-temporal (NDFI) and stathe- which involves the selection of criteria whose weights tisticalNormalized analysis Difference of two setsFlood of inimages, low Vegetation one containing Index will be determined via the AHP. In this process, the only the images before or after the event, and anoth- selection of criteria is very important. A plethora of er one containing images both before or after the event and during the occurrence of the event. NDFI criteria has been used in previous research on flood susceptibility mapping (Hong et al. 2018; Lyu et al.
Tab. 5 Details regarding the data from which each factor was compiled.
Original Map scale or spatial Primary input data Source Derived map format resolution EU-DEM Raster 25 m × 25 m EEA Elevation EU-DEM Raster 25 m × 25 m EEA Slope Angle Corine Land Cover 2012 Vector Better than 100 m Copernicus Land Cover EU-DEM Raster 25 m × 25 m EEA Drainage Density EU-DEM Raster 25 m × 25 m EEA TWI 1) Geologic Map of SE 1) 1:200,000 Institute of Geology and 1) 1:200,000 Rhodope – Thrace Raster Mineral Exploration (IGME) of Greece Geology 2) 1:50,000 2) Geologic map of Bulgaria 2) 1:50,000 Committee of Geology (CoG) WorldClim Raster 825 m × 825 m Fick et al. 2017 Rainfall Distance from EU-DEM Raster 25 m × 25 m EEA Streams 154 Christos Domakinis, Antonios Mouratidis, Kostas Voudouris, Theodore Astaras, Maria Chara Karypidou
2018; Seejata et al. 2018; Tang et al. 2018; Xiao et 3.2.2 Hydrological indicators al. 2018; Zhao et al. 2018; Kazakis et al. 2015). Their main characteristics being that they should be con- 3.2.2.1 Topographic Wetness Index (TWI) - tion mechanism, they can be measured or quickly cal- the indices of soil erosion, since it is used to relate the nectedculated with for the the whole physical study process area andof the that flood they genera ought effectsThis index of runoff according with togeomorphometry. Miliaresis (2011) It belongs is used toin to have simple interpretability (Papaioannou et al. 2015). Beven and Kirkby (1979) equation: According to Xiao et al. (2018) and Zhao et al. order to assess soil moisture and it is defined by the = ln (2) tan( ) availability, three types of indicators were considered in(2018), the present the factor’s research, effect i.e.on thetopographical, flood hazard hydrolog and data- � � area draining through a certain point per unit con- indicators (Xiao et al. 2018) provide information of In equation (2), α stands for the local upslope ical and meteorological. Specifically, topographical High values of TWI indicate areas more susceptible to the impact of the terrain. In the current study they tour length and tan(β) is the local slope in radians. consistthe flow of or elevation, stagnating slope of the angle water and on drainage the ground densi due- ty. On the other hand, hydrological indicators (Xiao et to3.2.2.2 flooding. Distance from streams - The drainage network of the drainage basin of Eryth- tration and accumulation of the water and the river ropotamos river has been produced by the EU-DEM al. 2018) provide information of the intercept, infil data layer with the use of raster-processing routines (TWI), distance from streams, land cover and geolo- (Voudouris et al. 2007). Furthermore, the distance gy.network. Finally, They the meteorologicalconsist of Topographic indicators Wetness (Zhao Index et al. from the streams of the drainage network data lay- 2018) provide information on the spatial distribution er was compiled through the use of geoprocessing of precipitation in the study area and were represent- - ed by the annual total rainfall. ceptibility mapping, because areas that are closer to The input data, their original format, the source buffer routines. This factor is crucial to flood sus from which the input data originated and their map event. scale or spatial resolution for each factor are synopti- streams are more likely to be inundated during a flood cally presented in the following table (Table 5). 3.2.2.3 Geology The synoptic geologic map of SE Rhodope – Thrace 3.2.1 Topographical indicators - tion (IGME) of Greece, at a scale of 1:200,000 (I.G.M.E. 3.2.1.1 Elevation 2002),from the was Institute used in of order Geology to produce and Mineral the part Explora of the Elevation is considered as an important factor for data layer that belongs to Greece. Accordingly, the geologic map of Bulgaria from the Department of drainage basin areas with low elevation values. The Geophysical Prospecting and Geological Mapping datafloods, layer because of elevation flood-prone was derived areas from tend EU-DEM. to occupy of the Committee of Geology (CoG 1989), at a scale of 1:50,000, was used in order to produce the part of the 3.2.1.2 Slope Angle data layer that belongs to Bulgaria. The slope angle data layer was produced by the - EU-DEM data layer with the aid of raster-processing - routines. Slope angle is also an important factor when calGeology formations is considered favor surface a significant runoff. factorOn the in other determin hand, ing flood-prone areas, because impermeable geologi it comes to discerning flood-prone areas, because 3.2.2.4permeable Land geological Cover formations favor infiltration. surfaceareas in terrain. a river basin that occupy flat terrain surfaces tend to flood more easily than areas with more steep 2017) was used to determine the land cover classes 3.2.1.3 Drainage Density withinThe data the layer limits of Corine of the Landstudy Cover area. 2012 It is worth(Copernicus men- length per unit area, which can be calculated as shown inThe the drainage following density equation is defined(1) (Zhou as et the al. total2014): stream ittioning depicts that more Corine closely Land the Cover surface 2012 relief’s was chosen land cover over Corine Land Cover 2018 (Copernicus 2020), because 1 = (1) by the station on Didymoteicho’s bridge. Addition- ally,conditions regarding during the the catchment 2010 flood, of which Erythropotamos was gauged DD stands∑ for drainage density, while S represents River, the differences between the aforementioned S the area of the grid and Li represents the length of river i within the grid. Areas with high drainage den- 2 - lateddata layersland cover are insignificant classes with sinceManning’s they covern roughness a total extent of approximately 2 km . Vieux (2004) corre sity indicate high flood susceptibility. Flood susceptibility mapping in Erythropotamos river basin 155
annual precipitation sums were considered as more formula: coefficient1 (Table2 1 6), which participates in Manning’s = 3 2 (3) prone3.3 Analytical to flooding. Hierarchy Process (AHP)
In equation (3), V 3/s), 3.3.1 Determination of flood susceptibility classes n A is the “wetted” for each factor cross-sectional area stands (m2), rfor stands discharge/flow for the hydraulic (m In order to apply the AHP methodology, which was radius is Manning’s and S is roughnessthe slope of coefficient, hydraulic grade or the lin- introduced by Saaty (1980), the data layer of each ear head loss (m/m). Moreover, Manning’s n rough- which means that low Manning’s n values correspond factor was classified into three classes according to toness high coefficient discharge is inverselyvalues. In proportional that way, areas to discharge, suscepti- assignedhow prone a ratingeach one of three of these (3), classes while those was to that flooding. are of n values. mediumClasses that susceptibility are highly were susceptible assigned to aflooding rating of were two (2) and those of low susceptibility were assigned a 3.2.3ble to Meteorological floods can be related indicators to low Manning’s rating of one (1).
3.2.3.1 Rainfall 3.3.2 Determination of the hierarchy between The annual total rainfall layer was derived using raw the flood susceptibility factors with the aid data that were retrieved from the WorldClim data- of the results of SAR-based inundation mapping base (Fick et al. 2017). The raw data involve monthly Having to deal with the subjectivity that often accom- precipitation totals, which refer to the climatological panies this step of AHP, the importance of each factor - was determined according to the proportion of the byperiod 824 1970–2000m) grid (Fick and et are al. 2017).available The as totalan approxi annual of 6.84 km2 for both open water and shallow water precipitationmately 30 seconds layer by was 30 constructedseconds (approximately by summing 824 all ininundated low vegetation) areas of that the intersected2010 flood withevent each (total factor’s area 12 monthly precipitation totals with the aid of map high susceptibility class (Figure 4 and Table 7). This algebra. Subsequently, the aforementioned rainfall concept was based on the idea that a SAR image that
water is concentrated. Moreover, the factors or indi- oflayer 25 wasm, was converted derived. to The a point downscaling shapefile, of from the whichorigi- is taken during a flood indicates the areas where flood- the final rainfall data layer, with a spatial resolution to the layer that was eventually used in the current cators of flood susceptibility all coexist in these are nal WorldClim layer (824 × 824 m grid resolution) as and it is known how each factor influences floods. by employing the universal kriging spatial interpola- For example it is known that, regarding e.g. slope analysis (25 × 25 m grid resolution), was performed whereangle, flatthe areasmost tendfavourable to flood conditions more easily. for Thus,most facthe- were the total annual precipitation values obtained at areas where flood water is concentrating are those tion method (Li et al. 2014). The interpolated values- for most factors or indicators intersect. Subsequent- iary variables used were elevation, slope, aspect and ly,tors the coexist, more i.e.a high where susceptibility the high susceptibility class of a factor classes or distanceeach point from of the the original sea. The WorldClim elevation grid.data Theused auxil was indicator is encountered in inundated areas, the more
Monitoring Services data portal (EEA 2017) and is susceptibility. providedthe EU-DEM on a obtained 25 by 25 mfrom grid. the Slope COPERNICUS and aspect Landwere influential this factor or indicator is in terms of flood derived from the EU-DEM using the available ras- ter-processing routines. Distance from the sea was - Tab. 7 Proportion of the total inundated area of the 2010 flood tines, at a spatial resolution of 25 m. Areas with high event that intersects with each factor’s high susceptibility class. also computed, by applying proximity analysis rou Extent of inundated Percent Factor area (km2) ratio (%) Tab. 6 Manning’s n roughness coefficients for certain land cover Land Cover 0.15 2.19 types according to Vieux (2004). TWI 0.12 1.75 Land Cover Manning’s n coefficient Geology 0.39 5.70 Artificial surfaces 0.015 Distance from streams 3.44 50.29 Agricultural areas 0.035 Rainfall 0.01 0.15 Forest and semi natural areas 0.100 Slope Angle 5.59 81.73 Wetlands 0.700 Drainage Density 0.99 14.47 Water bodies 0.030 Elevation 6.66 97.37 156 Christos Domakinis, Antonios Mouratidis, Kostas Voudouris, Theodore Astaras, Maria Chara Karypidou
Fig. 4 The data layer of the inundated areas of the 2010 flood event has been superimposed upon the flood susceptibility classes of the factors’ data layers: a) Elevation, b) Slope angle, c) TWI, d) Distance from streams, e) Geology, f) Land cover, g) Rainfall and h) Drainage density. Flood susceptibility mapping in Erythropotamos river basin 157
Tab. 8 Pairwise comparison of the factors that affect flood susceptibility.
Distance Slope Drainage Land Elevation from Geology TWI Rainfall angle Density Cover streams Elevation 1 2 3 4 5 6 7 8 Slope angle 1/2 1 2 3 4 5 6 7 Distance from streams 1/3 1/2 1 2 3 4 5 6 Drainage Density 1/4 1/3 1/2 1 2 3 4 5 Geology 1/5 1/4 1/3 1/2 1 2 3 4 Land Cover 1/6 1/5 1/4 1/3 1/2 1 2 3 TWI 1/7 1/6 1/5 1/4 1/3 1/2 1 2 Rainfall 1/8 1/7 1/6 1/5 1/4 1/3 1/2 1 Total 2.718 4.593 7.5 11.28 16.08 21.83 28.5 36
Tab. 9 Calculation of the factor weights with the use of the arithmetic mean method.
Distance Slope Drainage Land Elevation from Geology TWI Rainfall Mean angle Density Cover streams Elevation 0.368 0.435 0.403 0.355 0.311 0.275 0.246 0.222 0.327 (32.7%) Slope angle 0.184 0.218 0.268 0.266 0.249 0.229 0.211 0.194 0.227 (22.7%) Distance 0.123 0.109 0.134 0.177 0.187 0.183 0.175 0.167 0.157 (15.7%) from streams Drainage 0.092 0.073 0.067 0.089 0.124 0.137 0.140 0.139 0.108 (10.8%) Density Geology 0.074 0.054 0.045 0.044 0.062 0.092 0.105 0.111 0.073 (7.32%) Land Cover 0.061 0.044 0.034 0.030 0.031 0.046 0.070 0.083 0.050 (5%) TWI 0.053 0.036 0.027 0.022 0.021 0.023 0.035 0.056 0.034 (3.4%) Rainfall 0.046 0.031 0.022 0.018 0.016 0.015 0.018 0.028 0.024 (2.4%)
- because the measurements from the gauging station The 2010 flood extent was chosen for that purpose, theresulting corresponding matrix table, column an arithmetic factor with value which of 8 it indi has beencates compared,that a row whilefactor an is mucharithmetic more value significant of 1 means than wason Didymoteicho’s taken on 16/2/2010, bridge confirmed Erythropotamos that during indeed the date and time that the ENVISAT/ASAR’s flood image completion of Table 8, the arithmetic mean method lower spatial resolution when compared with SENTI- hasthat been both applied factors to are its equally results significant. and the weights After thefor flooded. Additionally, ENVISAT/ASAR’s imagery has- each factor were calculated (Table 9). ing station went out of order in 2012, the only way To sum up, Table 10 presents synoptically the fac- NEL-1 A/B imagery. Since the aforementioned gaug events was to rely on statements from members of the Departmentto collect information of Civil Protection for the 2017of the and region 2018 of Evrosflood classtors, theand classes the weight of flood that susceptibility was calculated into for which each they fac- (C. Papapostolou, Department of Civil Protection of torwere via classified, the application the rating of AHP that methodologywas assigned (Kazakisfor each the region of Evros, personal communication, 2018). et al. 2015).
3.3.3 Pairwise comparison between the flood 3.3.4 Consistency ratio susceptibility factors and determination In order to check the consistency of the eigenvector of their weights with the use of the arithmetic mean method according to the following formula: The factors were paired with each other and follow- matrix of AHP, the consistency ratio was calculated ing that, each factor was given an arithmetic value = (4) agreement with Table 7, when compared to the other In mathematic formula (4), CR stands for consis- factor,between with 1 which and 8, it according formed the to pair its (Table significance, 8). In the in tency ratio, CI RI
stands for consistency index, and 158 Christos Domakinis, Antonios Mouratidis, Kostas Voudouris, Theodore Astaras, Maria Chara Karypidou
Tab. 10 Synoptic table presenting the factors, their flood the weigh for each factor and Xi are the rating values susceptibility classes, the rating that was assigned for each class and the weight for each factor that was assigned through AHP methodology. referred. forThe each resulting pixel according susceptibility to the map factor was to validated which it by is Factor Class Rating Weight calculating, with the aid of geoprocessing routines, >600 1 the proportion of the inundated areas of the April Elevation 150–600 2 0.327 (m) 0–150 3 with its high susceptibility areas. This procedure indi- cated2017 thatand 59%March and 2018 58% flood of the events inundated that intersectedareas of the >15 1 Slope angle 2–5 2 0.227 (°) - 0–2 3 ing2018 map. and 2017 floods respectively coincided spatially Distance >500 1 with the high flood susceptibility zones of the result from 200–500 2 0.157 streams (m) 0–200 3 4. Results
Drainage 2.15–9.68 1 density 9.68–17.21 2 0.108 - 2 (m/m ) 17.21–24.74 3 According to the results that were produced by flood Permeable formations 1 event,extent covermapping, a total the of inundated 6.84 km2 , areaswhile withinthe inundated Erytro Geology Semi – permeable formations 2 0.073 potamos’ drainage basin, regarding the 2010 flood 2 and Impermeable formations 3 areas20.60 ofkm the2 respectively. flood events The that proportions occurred on of April the inun2017- Forests and wetlands 1 and March 2018 cover a total extent of 18.23 km Land Agricultural areas and water bodies 2 0.050 Cover areas and as shallow water in low vegetation areas Artificial surfaces 3 aredated presented areas that in morewere detaildetected in Table in open-land 11. flooded Regarding susceptibility mapping, the resulting −0.22–6.06 1 - TWI 6.06–12.34 2 0.034 tained areas of high, medium, and low susceptibility 12.34–18.62 3 (Figuremap was 5). classified Moreover, into the three application categories, of AHP which method con- 579.34–623.25 1 ology produced the criteria weight for each indicator. Rainfall 623.25–667.18 2 0.024 According to these results elevation was considered (mm) as the most important indicator with the weight value 667.18–711.10 3 of 0.327, followed by the slope angle with the weight value of 0.227. Distance from stream and drainage density are respectively considered as the third and RI depends on the number fourth most important criteria, and their weight val- of factor that are used to perform AHP and in our case, ues are 0.157 and 0.108, respectively. The weights of stands for random index.RI= 1.41 (Saaty 1980), while RI the remaining indicators are below 0.1, which indi- can be calculated by the following equation: cates that they present less important than aforemen- for an 8 by 8 matrix, tioned four indicators. The criteria weight value of = max (5) geology, land cover, TWI and rainfall are 0.073, 0.05, 1 − 0.034 and 0.024, respectively (Tables 8 and 9). Finally, by superimposing the delineated areas of − - the April 2017 and March 2018 inundation mapping In equation (5), λmax is the = maximum 8.41 and n eigenvalue = 8, there- onto the resulting susceptibility map, with the aid of offore the CR comparison = 0.042. According matrix andto Saaty n is the(1980), number if CR ofis lessfac geoprocessing routines, calculations indicated that tors. In the current study, λmax
3.3.5than 0.1,Calculation then the weights’of flood consistencysusceptibilityand is affirmed. the largest part of the aforementioned estimated flood validation of the results Tab. 11 Flood extents of the inundated areas for February 2010, Finally, the data layers for each factor were added April 2017 and March 2018 flood events. together in accordance with the mathematical equa- NDFI based NDFVI based Total tion (6): Flood event inundated area inundated area inundated area (km2) (km2) (km2) = (6) 2010 February 6.49 0.35 6.84 =1 2017 April 17.52 0.71 18.23 � ∗ 2018 March 19.28 1.32 20.60 i is In equation (6), S is the value for each pixel of the final flood susceptibility map of the study area, w Flood susceptibility mapping in Erythropotamos river basin 159
Fig. 5 Flood susceptibility map upon which are superimposed: a) the data layer of the April 2017 inundated areas (both NDFI and NDFVI based), and b) the data layer of the March 2018 inundated areas (both NDFI and NDFVI based).
class, while high and medium susceptibility classes susceptibility of the resulting map (Figure 5). Regard- intersected with 19.73 km2 (95.8%) of the resulting extents coincided with the areas of high and medium2 of the total inundated area in both open-land and in inglow thevegetation, April 2017 10.6 flood, km2 from(58.17%) the total intersected of 18.23 withkm flood extent (Table 12). 5. Discussion and medium susceptibility classes intersected with 17.54the high km flood2 susceptibility class. Additionally, high Correspondingly, from the total of 20.60 km2 of the appeared that the inundated areas of 2010 are consid- March 2018 (96.22%) inundated of area,the resulting 12.22 km flood2 (59.33%) extent. erablyBy comparing smaller thethan flood the inundatedextents of theareas flood of 2018events, and it 2, were included within the high flood susceptibility 2017. The 2010 flood covered an extent of 6.84 km 160 Christos Domakinis, Antonios Mouratidis, Kostas Voudouris, Theodore Astaras, Maria Chara Karypidou
Tab. 12 Area extent and percentage of the part of the March such as TWI (Arabameri et al. 2019; Das 2018; Tang 2018 and April 2017 inundated areas, which intersect with high to medium classes of the susceptibility map. Das 2018; Mahmoud et al. 2018; Kazakis et al. 2015), Inundated area Percentage Flood susceptibility classes drainageet al. 2018), density flow (Arabameriaccumulation et (Vojtekal. 2019; et Souissi al. 2019; et (km2) (%) al. 2019; Vojtek et al. 2019; Das 2018; Mahmoud 2018 March flood et al. 2018; Seejata et al. 2018) and rainfall (Souis- High 12.22 59.33 si et al. 2019; Mahmoud et al. 2018; Seejata et al. High and medium 19.73 95.80 2018; Tang et al. 2018; Kazakis et al. 2015) appear in most works. On the contrary, curvature (Arabameri 2017 April flood et al. 2019; Das 2018), NDVI (Normalized Difference High 10.60 58.17 High and medium 17.54 96.22 (Curve Numbers) (Vojtek et al. 2019; Mahmoud et al. Vegetation Index) (Arabameri et al. 2019), runoff/CN 2019) and groundwater depth (Souissi et al. 2019) 2018), SPI (Stream Power Index) (Arabameri et al. in comparison with the inundated areas of the 2017 Moreover, the number of the factors that are used in appear rarely on flood susceptibility mapping works. km2 and 20.60 km2 respectively; however, as it was et al. 2018; Rahmati et al. 2015) with the most com- andmentioned 2018 flood earlier, events, validation which of covered these results areas ofwas 18.23 not monflood numbersusceptibility of factors mapping ranging varies from greatly seven (Mahmoud to nine. feasible due to unfavourable weather conditions. At This paper tries to employ the most common and - first place, these differences in flood extents might- included,important since factors its datathat layercan be was used indirectly in flood employed suscep thermore,indicate that it is theworth 2010 mentioning flood had that probably according a lowerto the tibility mapping. Thus, flow accumulation was not gaugingreturn period station than on Didymoteicho’sthe other two flood bridge, events. the 2010 Fur 3/s, utilizedby the TWIbecause, factor. according Likewise, to Miliaresis Stream Power(2006), Indexthese while, on the other hand, such gauges were not avail- are,(SPI) like and TWI, Sediment indices Transport of soil erosion Index with(STI) very were simi not- flood reached a peak discharge of 1,255.05 m similar results. Curvature is considered to have a inundationable for the mapping2017 and with 2018 the floods. use of Remote Sensing lar mathematical expressions and thus they produce andRegarding especially theSAR, uncertainties it has to be mentioned that exist that in floodsuch and therefore it was not included in the factors that techniques and methodologies suffer mostly from minor impact on the occurrence of a flood (Das 2018)-
- were utilized in the assessment of flood susceptibili rentlyspeckle there and fromis no undermethodology or over-detection that can overcome of flood ty. Additionally, Topographic position index (TPI) and- extents especially in urban and vegetated areas. Cur lyTopographic appear to be roughness more important index (TRI) than theare factorsrarely usedwith mapping is still considered appropriate for validation whichin flood they susceptibility are compared, mapping so they and were even too more omitted. rare these difficulties entirely. However, flood inundation- Finally, the curve numbers (CN) data layer was not ping. (Giustarini et al. 2015a; Giustarini et al. 2015b; feasible to be compiled since there were no available Schumannin cases of floodet al. 2015)susceptibility and flood hazard map maps depicting the spatial distribution of the hydro- logical soil groups. wide variety of works that utilize the AHP methodol- The determination of the importance between the ogyConcerning for its implementation. flood susceptibility The main mapping, differences there and is a similarities of these works with the present research susceptibility mapping can be achieved by applying focus on the following points: 1) the factors that are variousfactors thatprocedures. are used Many in AHP researches when conducting use sensitivity flood employed by the research, 2) the determination of the analysis in order to overcome the subjectivity of AHP importance between the factors that are used by the (Souissi et al. 2019; Mahmoud et al. 2018; Tang et al. 2018), while weight linear combination is also a pop- factors of the study area and 4) the validation of the ular approach (Vojtek et al. 2019; Kazakis et al. 2015). AHP procedure, 3) the dominant flood susceptibility There is also a great number of works that employ The number and type of factors that are used in resulting flood susceptibility map. (Das 2018; Seejata et al. 2018; Rahmati et al. 2015). susceptibility, with the use of AHP, depend heavily on However,expert opinion the current in dealing trend with involves the hierarchy the use of of factors train- dataorder availability to determine (Xiao the et spatial al. 2018; distribution Zhao et al. of 2018). flood ing algorithms over a part of the elements that will be However, it can be observed that certain factors such used for the validation of the resulting susceptibility as elevation, slope angle, land cover, lithology and map, which usually involves a database of historical distance from streams are used in the vast majori- ty of works due to being easily produced via Digital The present paper is introducing the use of the results Elevation Models (DEMs), geological maps and the points where floods occurred (Arabameri et al. 2019).
of SAR-based inundation mapping, of a confirmed various Corine Land Cover data layers. Other factors via gauges flood event, in the determination of the Flood susceptibility mapping in Erythropotamos river basin 161
- data layer of rainfall had to be downscaled in order to reach the spatial analysis of the EU-DEM data lay- importance of the factors that affect flood susceptibil er. However, the aforementioned spatial variations of ity in AHP. Specifically, the aforementioned hierarchy- edwas with determined each factor’s according highest to susceptibility the part of the class, extent thus of effect on the resulting susceptibility map since they overcomingthe inundated the area subjectivity of the 2010 of AHP. flood that intersect these data layers do not appear to have a significant- When it comes to the determination of the most thewere study ranked area among involved the least a considerably important flood large suscep drain- results of AHP in various researches indicate that agetibility basin, factors which, (Table in terms 10). Additionally,of size, allowed the the extent use of thereimportant is no factorfactor into floodhave clearsusceptibility dominance mapping, over other the small scale data, which are widely used in likewise factors. Many papers indicate slope angle as the most - important factor (Arabameri et al. 2019; Vojtek et al. 2015).cases according to the existing literature regard elevation (Souissi et al. 2019) or even rainfall (Seejata ing AHP flood susceptibility mapping (Kazakis et al. et2019), al. 2018) while have flow accumulationbeen determined (Kazakis as the et dominant al. 2015), 6. Conclusions flood susceptibility factors in certain regions and by The present research paper introduced the idea to susceptibilitycertain methodologies. factor, but Likewise it can be the observed present researchthat the resultsdetermined depend elevation heavily as on the both most the important procedure flood that SAR imagery in order to determine the importance is employed in the determination of the hierarchy of use the extent of a flood that has been captured by factors and the conditions that lie within the studied with the subjectivity that involves the determination region. ofbetween the hierarchy flood susceptibility of factors in AHP.factors The and larger thus the dealing part of the inundated area that intersects with the factor’s susceptibility mapping with the use of AHP, the vast high susceptibility zone, the more important the fac- majorityRegarding of works the validationinvolves the of compilationthe results of aflood his- tor is considered over the others. torical database that includes, in the form of points, According to the results of the applied methodolo- sites where according to eye-witnesses or Remote susceptibility factor in the catchment of Erythropota- 2019; Souissi et al. 2019; Vojtek et al. 2019; Mahmoud mos.gy, elevation However, was this found has to to be be the further most dominantascertained flood by etSensing al. 2018). techniques The present floods work occurred handles (Arabameri this matter et byal. utilizing the results of SAR-based inundation map- and by taking advantage of the current and prospec- considering more future flood events in the same area- over, the resulting susceptibility map appeared to be wereping fornot specific involved flood in the occurrences. determination To ofthis the end, impor the- intive consistency availability with of SENTINEL-1 the, April 2017 imagery and Marchdata. More 2018 tanceinundated between areas the of factors the 2017 in AHP,and 2018were floods,used in whichorder - - bilityflood classextents, of the since resulting the aforementioned map. inundated tibilityto provide zones the of proportions the resulting of map. their Furthermore, respective flood the areasIt appears coincided that mostly the suggested with the methodology, high flood suscepti regard- extents that intersected with the high flood suscep- - ed that the areas of high susceptibility are located on ceptibility factors, via the results of SAR-based inun- resulting flood susceptibility map (Figure 5), indicat dationing the mapping,determination in AHP of producedthe hierarchy some of interesting flood sus results. Nevertheless, more thorough testing of this the basineastern mouth. part of the study area, specifically in the proposed methodology is required, while it also firstMoreover, half of main the streamscores thatand wereappear achieved increased by thetoward val- remains to be seen if its application on other drainage idation of the susceptibility map were quite high. In basins shall indicate each time another factor as more
- withparticular, the high approximately susceptibility 60% zones of theof the inundated map. The areas per- tibilityprevalent are in unique flood susceptibility,for each catchment. thus maintaining the from the April 2017 and March 2018 floods intersect argument that the conditions that affect flood suscep the aforementioned inundated areas intersect with thecentage map’s rises high to to approximately moderate susceptibility 96% in the zones. case that References Finally, it is worth mentioning that the source data - Arabameri, A., Rezaei, K., Cerdà, A., Conoscenti, C., bility mapping, in terms of scale and spatial resolu- Kalantari, Z. (2019): A comparison of statistical methods tion,layers were of the quite factors consistent that were since used the in majority flood suscepti of them were derived from EU-DEM that has a spatial resolu- susceptibility in Northern Iran. Science of the Total tion of 25 m. The data layers of geology and Corine Environmentand multi-criteria 660, decision443–458, making https://doi.org/10.1016 to map flood hazard /j.scitotenv.2019.01.021.
Land Cover 2012 involved a smaller scale, while the 162 Christos Domakinis, Antonios Mouratidis, Kostas Voudouris, Theodore Astaras, Maria Chara Karypidou
vulnerability to hydro-geomorphic hazards: estimating Astaras, T., Oikonomidis, D., Mouratidis, A. (2011): Psifiakí damageArguedas, probability A., Macedo, from L., Llerena,qualitative N. M.vulnerability (2016): Building Chartografía kai Geografiká Systímata Pliroforión [Digital assessment using logistic regression. Journal of Beven,Cartography K., Kirkby, and M. GIS].(1979): Thessaloniki: A Physically Ekdóseis Based, Dísigma Hydrology 541, 563–581, https://doi.org/10.1016 Variable[Disigma Contributing Publishing Group]. Area Model of Basin Hydrology. /j.jhydrol.2015.04.017. Hydrological Science Bulletin 24, 43–69, https://doi Fick, S. E., Hijmans, R. J. (2017): Worldclim 2: New 1-km .org/10.1080/02626667909491834. spatial resolution climate surfaces for global land areas. Birkholz, S., Muro, M., Jeffrey, P., Smith, H. M. (2014): International Journal of Climatology 37(12), 4302–4315, https://doi.org/10.1002/joc.5086.
EnvironmentRethinking the 478, relationship 12–20, https://doi.org/10.1016 between flood risk Matgen, P. (2015a). Flood hazard mapping combining /j.scitotenv.2014.01.061.perception and flood management. Science of the Total Giustarini,hydrodynamic L., Chini, modeling M., Hostache, and multi R., Pappenberger, annual remote F., sensing data. Remote Sensing 7(10), 14200–14226, https://doi.org/10.3390/rs71014200. Chen, Y., Liu, R., Barrett, D., Gao, L., Mingwei, Z., Renchullo, climates.L., Emelyanova, Science I. of(2015): the Total A spatial Environment assessment 538, Hostache, R., Matgen, P., Verhoest, N.E.C. De Baets, B. framework for evaluating flood risk under extreme Giustarini,(2015b): L., Accounting Vernieuwe, for H., image Verwaeren, uncertainty J., Chini, in M.,SAR-based .2015.08.094. Chen,512–523, Y., Yu, J., http://dx.doi.org/10.1016/j.scitotenv Khan, S. (2013): The spatial framework for Observation and Geoinformation 34, 70–77, https://doi weight sensitivity analysis in AHP-based multi-criteria .org/10.1016/j.jag.2014.06.017.flood mapping. International Journal of Applied Earth decision making. Environmental Modelling and Software 48, 129–140, https://doi.org/10.1016/j.envsoft.2013 A-X., Chen, W., Kougias, I., Kazakis, N. (2018): Flood .06.010. Hong,susceptibility H., Panahi, assessmentM., Shirzadi, in A., Hengfeng Ma, T., Liu, area J., Zhu,coupling Cian, F., Marconcini, M., Ceccato, P. (2018). Normalized adaptive neuro-fuzzy inference system with genetic algorithm and differential evolution. Science of the Taking advantage of EO big data. Remote Sensing of total Environment 621, 1124–1141, https://doi.org EnvironmentDifference Flood 209, Index 712–730, for rapid https://doi.org/10.1016 flood mapping: /10.1016/j.scitotenv.2017.10.114. /j.envsoft.2013.06.010. CoG (Committee of Geology) (1989). Geologic map of (2002): Synoptic geologic map of SE Rhodope – Thrace, P.R. Bulgaria, Scale 1:50 000. I.G.M.E.Scale (Institute 1:200 000, of Xanthi. geology and mineral exploration) COPERNICUS, Europe’s eyes on earth – land monitoring Kanani-Sadat, Y., Arabsheibani, R., Karimipour, F., Nasseri, service. Available online: https://land.copernicus.eu /pan-european/corine-land-cover/clc-2012?tab assessment in data-scarce and ungauged regions based =download (accessed on 25/11/2017). onM. (2019):GIS-based A new hybrid approach multi criteria to flood decision-making susceptibility COPERNICUS, Europe’s eyes on earth – land monitoring method. Journal of Hydrology 572, 17–31, https:// service. Available online: https://land.copernicus.eu doi.org/10.1016/j.jhydrol. 2019.02.034. /pan-european/corine-land-cover/ clc2018 (accessed Kazakis, N., Kougias, I. Patsialis, T. (2015): Assessment on 17/4/2020).
applicationof flood hazard in Rhodope–Evros areas at a regional region, scale Greece. using anScience Costache,(2020): R., Flash-Flood Pham, Q. B., Susceptibility Sharifi, E., Linh, Assessment N. T. T., Abba, Using ofindex-based the Total Environment approach and 538, Analytical 555–563, Hierarchy https:// Process: S. I., Vojtek, M., Vojteková, J., Nhi, P. T. T., Khoi, D. N. doi.org/10.1016/j.scitotenv.2015.08.055. Supported by Remote Sensing and GIS Techniques. RemoteMulti-Criteria Sensing Decision 12(1), 106,Making https://doi.org/10.3390 and Machine Learning /rs12010106. Khosravi,its mapping K., Nohani, in Iran: E., a Maroufinia, comparison E., between Pourghasemi, frequency H. R. Das, S. (2019): Geographic information system and ratio(2016): and A weights-of-evidence GIS-based flood susceptibility bivariate statistical assessment and models with multi-criteria decision-making technique. Maharashtra, India. Arabian Journal of Geosciences 11, Natural Hazards 83, 947–978, https://doi.org/10.1007 576,AHP-based https://doi.org/10.1007/s12517-018-3933-4. flood hazard zonation of Vaitarna basin, /s11069-016-2357-2. Demek, J. (1972). Manual of detailed geomorphological mapping. Academia, Prague. in the environmental sciences: A review. Environmental Dikau, R. (1989): The application of a digital relief model to Li, ModellingJ., Heap, A. and(2014): Software Spatial 53, interpolation 173–189, https://doi.org methods applied landform analysis. In: J. F. Raper, ed. Three Dimensional /10.1016/j.envsoft.2013.12.008 Application in Geographical Information Systems. Taylor
European Environment Agency (EEA) (2017): Li, Sentinel-1Y., Martinis, SAR S., Plank, data. International S., Ludwig, R. Journal(2018): of An Applied automatic and Francis, London, 51–77. Earthchange Observation detection approach and Geoinformation for rapid flood 73, mapping 123–135, in EU-DEM. Published 7/12/2017 from: https:// https://doi.org/10.1016/j.jag.2018.05.023. www.eea.europa.eu/data-and-maps/data/Copernicus Land Monitoring Service – copernicus-land-monitoring-service-eu-dem. mapping for Namibia using change detection and Longthresholding S., Fatoyinbo, with T. E.,SAR. Policelli, Environmental F. (2014): Research Flood extent Thouret, J.-C., Manville, V., Negulescu, C., Zuccaro, Ettinger, S., Mounaud, L., Magill, C., Yao-Lafourcade, A.-F., /1748-9326/9/3/035002. Letters 9, 35002–35009, https://doi.org/10.1088 G., De Gregorio, D., Nardone, S., Uchuchoque, J.A.L., Flood susceptibility mapping in Erythropotamos river basin 163
approach: The case of Mouriki, Greece. Flood risk assessment in metro systems of mega-cities Geosciences 8(50), 1–13, https://doi.org/10.3390 Lyu,using H. M., a GIS-basedSun, W. J., Shen,modeling S. L., approach.Arulrajah, ScienceA. (2018): of the /geosciences8020050. Total Environment 626, 1012–1025, https://doi.org /10.1016/j.scitotenv.2018. 01.138. Mahmoud, S., Gan, T. Y. (2018): Multi-criteria approach Pierdicca,gained fromN., Pulvirenti, COSMO-SkyMed L., Chini, observations. M., Guerriero, Acta L. Candela, L. (2013): Observing floods from space: Experience of Middle East. Journal of Cleaner Production 196, /j.actaastro. 2012.10.034. 216–229,to develop https://doi.org/10.1016/j.jclepro.2018 flood susceptibility maps in arid regions Astronautica 84, 122–133, http://dx.doi.org/10.1016 .06.047. Psomiadis,https://doi.org/10.1117/12.2241055. E. (2016): Flash flood area mapping utilising SENTINEL-1 radar data. Proceedings of SPIE, 10005, Matgen, P., Hostache, R., Schumann, G. J. P., Pfister, L., learnedHoffmann, from L. Savenije,two case H.studies. H. G. (2011): Physics Towards and Chemistry an Pulvirenti,Synthetic L., Aperture Pierdicca, Radar N., Chini, (SAR) M., data Guerriero, using fuzzy L. (2011): logic. ofautomated the Earth SAR-based 36 (7–8), 241–252, flood monitoring https://doi.org/10 system: lessons NaturalAn algorithm Hazards for andoperational Earth System flood mappingSciences 11,from 529–540, .1016/j.pce.2010.12.009. https://doi.org/10.5194/nhess-11-529-2011. Miliaresis A. (2006). Eidikés Efarmogés sto AcGIS. Rahmati, O., Zeinivand, H., Besharat, M. (2015): Flood hazard zoning in Yasooj region, Iran, using GIS and multi-criteria decision analysis. Geomatics, Natural Mouratidis[Specialized A. (2011): applications Contribution in ArcGIS]. of Earth Ekdóseis Observation «ION» Hazards and Risk 7(3), 1000–1017, https://doi.org [ION Publishing Group], Athens. /10.1080/19475705.2015.1045043.
data and GIS to mapping and managing flood evens in McGraw-Hill, New York. Mouratidis,Greece. Final A., Ampatzidis, Report of “GEOGIS D. (2019): floods” European research Digital project Schlaffer,Saaty, T. L. S., (1980): Matgen, The P., AnalyticHollaus, M.Hierarchy Wagner, Process. W. (2015): Flood Elevationfunded by Model the John (EU-DEM) S. Latsis validationPublic Benefit against Foundation. detection from multi-temporal SAR data using harmonic analysis and change detection. International Journal comparison with SRTM and ASTER-GDEM in Central of Applied Earth Observation and Geoinformation 38, Macedoniaextensive Global (Greece). Navigation ISPRS, InternationalSatellite Systems Journal data of and Geo-Information 8(3), 108, https://doi.org/10.3390 Schumann, G. J. P., Moller, D. K. (2015): Microwave remote /ijgi8030108. 15–24, http://dx.doi.org/10.1016/j.jag.2014.12.001. the Earth 83–84, 84–95, https://doi.org/10.1016/j.pce F. Tsakiri-Strati, M. (2011): Flood studies in Greece using .2015.05.002.sensing in flood inundation. Physics and Chemistry of Mouratidis,Earth Observation A., Nikolaidou, data andM., Doxani, Geographical G., Lampiri, Information M., Sarti, Seejata, K., Yodying, A., Wongthadam, T., Mahavik, N., Systems. Geological Remote Sensing Group (GRSG) Annual General Meeting, 7–9 December 2011, ESA /ESRIN, Frascati, Italy. Basin,TantaneE, Sukhothai S. (2018): Province. Assessment Procedia of flood Engineering hazard areas212, Mouratidis, A., Sarti, F. (2013): Flash-Flood Monitoring 340–347,using Analytical https://doi.org/10.1016/j.proeng.2018 Hierarchy Process over the Lower Yom and Damage Assessment with SAR Data: Issues and .01.044. Future Challenges for Earth Observation from Space Sustained by Case Studies from the Balkans and Zghibi, A., Dlala, M. (2019): GIS-based MCDM – AHP Souissi, D., Zouhri, L., HammamI, S., Msaddek, M. H., eds. Earth Observation of Global Changes (EOGC). southeastern Tunisia. Geocarto International, https:// Eastern Europe. In: Krisp J., Meng L., Pail R., Stilla U., doi.org/10.1080/10106049.2019.1566405.modelling for flood susceptibility mapping of arid areas, Springer, Berlin, Heidelberg, 125–136, https://doi. org/10.1007/978-3-642-32714-8_8.Lecture Notes in Geoinformation and Cartography. hazard based on natural and anthropogenic factors Nandi, A., Mandal, A., Wilson, M., Smith, D. (2016). Flood Stefanidis,using analytic S., Stathis, hierarchy D. (2013): process Assessment (AHP). Natural of flood hazard mapping in Jamaica using principal component Hazards: Journal of the International Society for analysis and logistic regression. Environmental the Prevention and Mitigation of Natural Hazards, Earth Sciences 75, 465, https://doi.org/10.1007 Springer;International Society for the Prevention and /s12665-016-5323-0. Mitigation of Natural Hazards 68, 569–585, https:// doi.org/10.1007/s11069-013-0639-5. Tang, Z., Zhang, H., Yi, S., Xiao, Y. (2018): Assessment Nikolaidou,database M.,for Mouratidis,Greece supported A., Doxani, by Earth G., Oikonomidis, Observation dataD., Tsakiri-Strati, and GIS. Proceedings M., Sarti, F.of (2015): the 10th An International on-line flood probabilistic multi-criteria decision analysis. Journal of Congress of the Hellenic Geographical Society, 22–24 Hydrologyof flood susceptible 558, 144–158, areas https://doi.org/10.1016using spatially explicit, October 2014, Thessaloniki, Greece, 671–688. /j.jhydrol.2018.01.033. Tehrany, M. S., Pradhan, B., Jebur, M. N. (2014): Flood susceptibility mapping using a novel ensemble weights- Papaioannou,prone areas G., mapping. Vasiliades, Water L., Loukas, Resoures A. (2015):Management 29, of-evidence and support vector machine models in GIS. Multi-criteria analysis framework for potential flood Journal of Hydrology 512, 332–343, https://doi.org -0817-6. /10.1016/j.jhydrol.2014.03.008. Patrikaki,399–418, O., http://dx.doi.org/10.1007/s11269-014Kazakis, N., Kougias, I., Patsialis, T., Tehrany, M. S., Pradhan, B., Mansor, S., Ahmad, N. (2015): Theodosiou, N. Voudouris, K. (2018): Assessing Flood susceptibility assessment using GIS-based support vector machine model with different kernel types.
flood hazard at river basin scale with an index-based 164 Christos Domakinis, Antonios Mouratidis, Kostas Voudouris, Theodore Astaras, Maria Chara Karypidou
Catena 125, 91–101, https://doi.org/10.1016/j.catena method considering spatial heterogeneity of risk .2014.10.017. preference. Science of the Total Environment 599–600, 1034–1046, https://doi.org/10.1016/j.scitotenv using GIS. Kluwer Academic Publishers, Netherlands, .2017.04.218. Vieux,123–124. B. E. (2004). Distributed Hydrologic Modelling Yésou, H., Sarti, F., Tholey, N., Mouratidis, A., Clandillon, S., Voudouris, K., Mavrommatis, T. Antonakos, A. (2007): Huber, C., StudeR, M., De Fraipont, P. (2013): Addressing Hydrologic balance estimation using GIS in Korinthia prefecture, Greece. Advances in Science and Research 1, 1–8, https://doi.org/10.5194/asr-1-1-2007. emergency flood mapping and monitoring of inland 2013,water Edinburgh,bodies with UK, sentinel ESA SP-722. 1-2. Expectations and Mapping on a National Scale in Slovakia Using the Zhao,perspectives. G., Pang, B., Living Xu, Z., Planet Yue, J., Symposium, Tu, T. (2018): 9–13 Mapping September Vojtek,Analytical M., Vojteková, Hierarchy J. (2019): Process. Flood Water Susceptibility 11, 364, https:// doi.org/10.3390/w11020364. scale in China. Science of the Total Environment 615, 1133–1142,flood susceptibility https://doi.org/10.1016/j.scitotenv in mountainous areas on a national .2017.10.037. Wang,GIS Y.,Multi-Criteria Hong, H., Chen, Decision-Making W., Li, S., Pamucar, Method D., Gigovic,for Flood SusceptibilityL., Drobnjak, S., Mapping Bui, D. T., at Duan, Shangyou, H. (2019): China. A Remote Hybrid Crompton, D., Fitzgerald, G., Huang, C. (2018): The Sensing, 11, 62, https://doi.org/10.3390/rs11010062. Zhong,long-term S., Yang, physical L., Toloo, and S., psychological Wang, Z., Tong, health S., Sun, impacts X.,
Environment 626, 165–194, https://doi.org/10.1016 Wang, Y., Li, Z. Tang, Z. (2011): A GIS-based spatial /J.SCITOTENV.2018.01.041.of flooding: A systematic mapping. Science of the Total Resourcesmulti-criteria Management approach for25, flood 3465–3484, risk assessment https:// in Zhou, Z., Psomiadis, E., Dercas, N. (2014): Analysis of the doi.org/10.1007/s11269-011-9866-2.the Dongting Lake Region, Hunan, Central China. Water networks. Earth Surface Dynamics 2, 105–116, https:// assessment based on spatial ordered weighted averaging doi.org/10.5194/esurf-2-105-2014.drainage density of experimental and modelled tidal Xiao, Y., Yi, S., Tang, Z. (2018): Integrated flood hazard Original Article 165 A comparison of actual evapotranspiration estimates based on Remote Sensing approaches with a classical climate data driven method Soghra Tofigh, Dariush Rahimi*, Reza Zakerinejad
University of Isfahan, Faculty of Geographical and planning, Iran * Corresponding author: [email protected]
ABSTRACT The knowledge of actual evapotranspiration at farm level is a prerequisite for irrigation planning, farm management, to increase production and reduce water consumption. To accomplish this, comprehensive and accurate assessment methods should be applied. In order to evaluate accurately evapotranspiration processes we compared lysimeter evapotranspiration data with MODIS (Aqua and Terra satellites) and LANDSAT (SEBAL algorithm) satellite images as well as with the FAO Penman-Montith method. The findings indicate the low error rate, high correlation (1) and appropriateness of SEBAL in estimating actual evapotranspiration. The error values MAD, MSE and RMSE between lysimeter and the SEBAL algorithm were 0.59, 0.36 and 0.60 respectively. The second best performance was established for the FAO Penman-Montith method. The obtained error values MAD, MSE and RMSE between the lysimeter and FAO-Penman-Montith method are 0.91, 1.29 and 1.13, respectively.
KEYWORDS actual evapotranspiration; SEBAL algorithm; Landsat; MODIS; Penman-Montith; Wheat; lysimeter
Received: 4 October 2019 Accepted: 8 May 2020 Published online: 18 September 2020
Tofigh, S., Rahimi, D., Zakerinejad, R. (2020): A comparison of actual evapotranspiration estimates based on Remote Sensing approaches with a classical climate data driven method. AUC Geographica 55(2), 165–182 https://doi.org/10.14712/23361980.2020.12 © 2020 The Authors. This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0). 166 Soghra Tofigh, Dariush Rahimi, Reza Zakerinejad
1. Introduction pixels. The low value of MSE, MAE, MAD, and RMSE - Water need is one of the most important parame- ings mentioned above (e.g. Karbasi et al. 2016; Ghor- ters in crop cultivation and in terms of planning the baniobtained et al. in 2015; different Morshedi other et studies al. 2016; sustain Rezaei the Banaf find- sheh et al. 2014; Kamali and Nazari 2018; Atasever are one of the major challenges in dry and semi-ar- and Ozkan 2018; Fu et al. 2018; Rawat et al. 2017). idirrigation regions calendar. like Iran Waterwhich Sourceis due todeficit the lowestimations amount The results obtained by Wagle et al. (2017) on the of precipitation (248 mm), high temperatures (aver- - age temperature 18 °C, which is 3 °C higher than the els of (SEBAL), (METRIC), (SEBS), (S-SEBI (SSEBop)) global average) long dry season (in some areas up to foroperation evapotranspiration of five of surface of sorghum energy predictionbalance SEB, indicate mod 8 months), high evaporation, inappropriate cultiva- that the S-SEBI, SEBAL and SEBS outperform METRIC tion pattern and irregular irrigation methods (Aliza- and SSEBop models with higher accuracy. deh 2016; Zakerinejad and Masoudi 2019). In the In this study we estimate evapotranspiration current condition, the phenomenon of global warm- through: i) the SEBAL algorithm, ii) the FAO-Pen- ing and the occurrence of severe and continuous man-Montith method and iii) MODIS evapotranspira- tion products and compare the obtained results con (Zakerinejad and Maerker 2015). Consequently, due observed lysimeter data. The study area is located in todroughts the high and impact desertification of evapotranspiration aggravate the processes problem the Shahrekord plain that is characterized by a tem- in plant and water resources management we assess perate climate and wheat cropping as dominant agri- different models for the Iranian conditions. Especially, cultural production. we focus in our study on the main crop cultivation in Iran, which is namely wheat. Evapotranspiration is highly affecting the hydro- 2. Study area logical cycle and water balance equations. Meas- uring, calculating and estimating the evapotran- The Centre of Shahrekord plain is located at 32°29' to spiration volume are essential in water resource 32°38' N and 50°46' to 50°55' E at 2066 m above sea management. Different methods exist on the marked level (Fig. 1). The annual average precipitation of the such as: direct measurements (Lysimeter), multiple plain is 330 mm and the annual temperature average models like Priestley-Taylor, Jensen-Haise, Thornth- is of 12 °C. The test farm is located at Farrokhshahr Waite, Blaney-Criddle, FAO-Pennman Monteith, Har- Agricultural Meteorological Research Center (AMRC). greaves-Samani, Turc, Making and Ritchie (Allen et al. The farm is equipped with a drainage lysimeter with 1998; Zare Haghi et al. 2016). In addition, evapotran- a diameter of 3 m and area of 7.60 m2 and a cover spiration can be estimatedusing remotely sensed data consisting of 1200 wheat seeds. This farm is one of and respective modelling approaches such as SEBI, the experimental farms that estimats the actual evap- SEBAL, S-SEBI, SEBS, METRIC, S-TSEB and P-TSEB otranspiration data through the SEBAL algorithm and (Alizadeh et al. 2016). Moreover, the MODIS sensor we compare the results with that of the nearest wheat also measures evapotranspiration, that is represented farm (Fig. 1). in a 8-day composite dataset. The results of many studies in different countries like China, Poland, Slovakia, Iraq, and Brazil indi- 3. Method and materials cate that the SEBAL algorithm is suitable to estimate evapotranspiration even in areas with climate data The method applied in this study is comparative and shortage (Santos et al. 2017; Ndou et al. 2018; Li et illustrated in the following. al. 2013; Santos 2017; Jaber 2016; Jian 2015; Bezer- ra 2015; Sun 2011). The MODIS evapotranspiration 3.1 Method - tions of evapotranspiration over a wide area. Exten- According to the available databases, the observation- siveproduct research provides has been significant done in information this context onby e.g.varia Yu al data (lysimeter) is used a reference data. We esti- et al. (2019), Rasmussen et al. (2014), or Sun et al. mate the evapotranspiration processes following the (2012). SEBAL algorithm given by Eq. 1–20, and the FAO-Pen- In this context, the results obtained through the man-Monteith model reported by Eq. 21. To test the SEBAL algorithm where compared with experimen- accuracy of these models and select the optimal mod- tal methods like Hargrevi-Samani, Blaney-Criddle, el the RMSE Eq. 23, MES Eq. 24, MAD Eq. 25 and R Eq. FAO Penman-Monteith, Metric, SWAT and lysimeter 22 indexes are calculated. data in countries like Turkey, India and some cities in Iran (Zanjan, Mazandaran, Neyshabur). Particularly, 3.2 Materials the SEBAL algorithm has been compared with actu- al lysimeter data showing small errors, which main- The data bases consist of three data groups: i) mete- ly related to the determination of the cold and hot orological data collected from Farokhshahr station Actual evapotranspiration 167
Fig. 1 Geographical locationof the study area and sample farm.
2 which include average, minimum and maximum tem- Where λET ), Rn is the perature, average, minimum and maximum precipita- 2), G is the soil tion, relative humidity (RH), wind speed, sunny hours, is2 the latent heat flux (W/m ii) evapotranspiration data measured by the lysime- net radiation2). flux at the surface (W/m ter, iii) Landsat satellite images (2016–2017), and vi) heatThe flux surface (W/m ),energy and H isbudget the sensible equation heat isflux further to the the MODIS evapotranspiration product (2016–2017) explainedair (W/m in part 4 of this section. Table 1. n) represents the actual radiant energy available at the surface. It is The net radiation flux at the surface (R Tab. 1 Specification of the applied satellite images - (LANDSATandMODIS). tratedcomputed in the by surface subtracting radiation all outgoing balance equation:radiant fluxes Image type Imaging time(D/M/Y) Julian Day from all incoming radiant fluxes (Figure 2) as illus R + R )R (2) Landsat7 25/07/2017 205 n 0 Landsat7 22/05/2017 143 Where = (1 −R α)R iss� the L�incoming − RL� − (1 shortwave − ε L� radiation Landsat7 11/11/2016 314 2), α is the surface albedo (dimensionless), R S� 2), R is the 2 L� (W/m o is the 3.2.1 SEBAL algorithm issurface the incoming thermal long emissivity wave radiation (dimensionless) (W/m (WatersL� In the SEBAL model, ET is computed from satellite outgoinget al. 2002). long wave radiation (W/m ), and ε images and weather data using the surface energy balance. Since the satellite image provides informa- Surface Albedo (α): The albedo at the top of the tion for the overpass time only, SEBAL computes an atmosphere is compute as follows: is calculated for each pixel of the image as a “residual” αtoa ωλ × ρλ) (3) ofinstantaneous the surface energy ET flux budget for the equation: image time. The ET flux = Σ( is a weighting n Where ρλ is the reflectivity and ωλ λET = R − G − H (1) coefficient for each band compute as follows: 168 Soghra Tofigh, Dariush Rahimi, Reza Zakerinejad
Tab. 2 Component of atmospheric transmissivity.
(4) 2 Station τsw τ sw z (m) Where ESUN is elevation of the sun. Albedo is Shahrekord 0.79 0.62 2066
Tab. 3 Component of RS↓ Equation. thedefined incident as the light ratio emitted of the byelectromagnetic the sun. Surface radiation albedo 2 2 Date GCS (w/m ) cos θ dr τsw RS↓ (w/m ) isreflected computed from by the correcting surface of the the α soiltoa forand atmospheric the plant to transmissivity: 2016/11/11 1367 0.79 1.02 0.79 871.68 2017/05/22 1367 0.93 0.97 0.79 981.65 (5) 2017/07/25 1367 0.91 0.96 0.79 954.29
path_radiance range between 0.025 and 0.04 and for SEBAL we recommend a value of 0.03 basedValues on Bastiaanssen for α (1998). 2). Its value is computed as follows: sw includes the transmissivity of both direct solar radiation flux that actually reaches the earth’s surface beam radiation and diffuse (scattered) radiation to (W/mR GCS × cos θ × dr × τsw (7) theτ surface. We calculate τsw assuming clear sky and s� 2 relatively dry conditions using an elevation-based Where = Gsc ), cos relationship from FAO-56: above, dr is the is inverse the solar squared constant relative (1367 earth-sun W/m dis- –5 τsw × z (6) θ is the cosinesw isof thethe atmosphericsolar incidence transmissivity. angle as defined The value R Where = 0.75 z is × the 2 × elevation 10 above sea level (m). tance, and τ S� Incoming Shortwave Radiation (RS ): Incoming is computed for the days specified. shortwave radiation is the direct and diffuse solar � Sun Elevation = 37.84438276 (Metadata file) � θ = 90 − 37.84438276= 52.15561724 � cos θ = 0.79
Fig. 2 Surface Albedo of wheat farm (2016/11/11). Actual evapotranspiration 169
Fig. 3 Surface Albedo of wheat farm (2017/22/05).
The SAVI is an index that attempts to “subtract” the effects of background soil from NDVI so that impacts Sun Elevation = 68.00616546 (Metadata file) � of soil wetness are reduced in the index. It is comput- θ = 90 − 68.00616546= 21.99383 � cos θ = 0.93 ed as: Sun Elevation = 65.69948076 (Metadata file) � θ =Outgoing 90 − 65.69948076 Long wave = 24.300522 Radiation �(R cosL ): θThe = 0.91 outgo- SAVI L)(R1 2 1 + R2) (9)
� Where; = (1 L +is a constant − R )/(R for SAVI. If L is zero, SAVI ing long2). waveIt is computed radiation inis SEBALthe thermal through radiation the follow flux- becomes equal to NDVI. A value of 0.5 frequently emitteding steps: from the earth’s surface to the atmosphere appears in the literature for L. (W/m The LAI is the ratio of the total area of all leaves on 1. Computation of vegetation indices of Normal- a plant to the ground area represented by the plant. It ized Difference Vegetation Index (NDVI), Soil Adjust- is an indicator of biomass and canopy resistance. LAI ed Vegetation Index (SAVI), and Leaf Area Index (LAI) is computed for southern Idaho using the following - empirical equation: ties for the near-infrared band (5) () and the red band (4)The () to NDVI their is sum: the ratio of the differences in reflectivi (10)
1 2 1 + R2) (8) Where; SAVIID is the SAVI calculated from Equation (9). NDVIThe NDVI = (R is− Ra sensitive)/(R indicator of the amount and condition of green vegetation. Values for NDVI between 0 and 1 and water and cloud are usually less energy2. Computation radiated by of the Surface surface emissivity to the thermal (ε) energy thanrange zero. between −1 and +1. Green surfaces have a NDVI radiatedSurface by emissivity a blackbody (ε) at is the the same ratio temperature. of the thermal 170 Soghra Tofigh, Dariush Rahimi, Reza Zakerinejad
Fig. 4 NDVI of wheat farm (2016/11/11).
Fig. 5 NDVI of wheat farm (2017/22/05). Actual evapotranspiration 171
Fig. 6 LST of wheat farm (2016/11/11).
Fig. 7 LST of wheat farm (2017/22/05). 172 Soghra Tofigh, Dariush Rahimi, Reza Zakerinejad
3. Computation of corrected thermal radiance (Rc) Tab. 4 Constant of K for Equation 12. The corrected thermal radiance (Rc) is the actual Landsat K1 K2 radiance emitted from the surface. Landsat 8 band 10,11 666.09 1282.71 4. Computation of surface temperature (Ts) The surface temperature (Ts) is compute using the Tab. 5 Components of the RL↓ Equation for different days. following equation: 4 2 Date Tα (K) ε0 RL↓ (W/m ) 2016/11/11 284.09 0.74 275.30 (11) 2017/05/22 296.15 0.74 326.47 2017/07/25 295.70 0.74 323.30 Where TB is the Brightness Temperature, wave- length of emitted radiance (the value of λ for bands 4 RL ε × σ × Tα (14) h × × 10 - � ο 6stant the (6.62607015 Landsat 4,5,7 × 10is−2 11.45.), e is emissivity, c2 = = (1.380649 c/s =1.4388 × 10 −34m K, where h = Planck’s con× Net surface radiation (Rn) is calculated is comput- 8 10 −23 J s), s = Boltzmann constant ed using Equation (2). 0.004 × J/K), c = velocity of lightmin (2.998max –NDVI m/s).min. The value of e is obtained from this relation Soil Heat Flux (G): Brightness Pv +0.986, Temperature where Pv = is NDVI obtained – NDVI from/NDVI the fol- storage into the soil and vegetation due to conduction. lowing relation: n for agriculture Soil heat flux surfaces is the rateis between of heat 0.05–0.15. Estimates of G/R (12) (15)
Where K1 and K2 are constants for Landsat imag- Sensible Heat Flux (H): es, Lλ (Lλ LQcal + AL) spectral radiance where ML rate of heat loss to the air by convection and conduc- tion, due to a temperature difference. Sensible heat It is flux compute is the metadata, = Q Mcal is quantized and calibrated standard using the following equation for heat transport: productis band-specific pixel value multiplicative and AL rescaling factor from- caling factor from metadata.in this paper from Bright- (16) ness Temperature and wavelength is band-specific of emitted additive radi res- 3 ance recorded by the sensor (thermal band) is used. Where ρ ), cp - 5. Computation of Outgoing Long wave Radiation ence (T1 is2) airbetween density two (kg/m heights (z 1is and air zspecific2), and (RL heatrah is (1004 the aerodynamic J/kg/K), dT (K)resistance is the temperature to heat transport differ This is computed using the Stefan-Boltzmann 1 −is T the height just above the zero plane dis- equation:�) placement (d 0.67 × height of vegetation) for the (s/m).surface z or crop canopy and z2 is some distance above 4 RL ε × σ × T� (13) the zero plane≅ displacement, but below the height of the surface boundary layer. Based on experimental � ο = is the “broad-band” surface emissivity analysis, values of 0.1 meter for z1 and 2.0 meters for z2 are assigned. Temperature difference (dt) is given ο 2 4 (5.67Where × 10 ε ), and Ts is the surface temper- 1 – Tz2. The air temperature at each pixel is (dimensionless),ature (K). −8 σ is the Stefan-Boltzmann constant unknown, along with explicit values for Tz1and Tz2. W/m /K Therefore,as dT = Tz only the difference dT is utilized. SEBAL Choosing the “Hot” and “Cold” Pixels: The “cold” computes dT for each pixel by assuming a linear rela- pixel is selected as a wet, well-irrigated crop surface tionship between dT and Ts s, where b and having full ground cover by vegetation. The sur- TS is the land sur- face temperature and near-surface air temperature face temperature. a is obtained: dt =by b subtracting+ aT the dT are assumed similar at this pixel. The “hot” pixel is a are the correlation coefficients and - assumed zero. and(dt hot cold pixel pixels − dt are cold separated pixel) and according the LST to(LST vegetation hot pix selected as a dry, bare agricultural field where ET is andel − LSTtemperature cold pixel). of theUsing pixels Envi and software, the dt arefirst calculat the hot- Incoming Long wave Radiation (RL ): The incom- ed based on the difference of two hot and cold pixels: ing long wave radiation is the downward thermal � 2). It is com- puted using the Stefan-Boltzmann equation: a = (dt hot pixel − dt cold pixel) / (LST hot pixel − LST radiation flux from the atmosphere (W/m cold pixel). b is obtained by multiplying −a in LST hot pixel and dt hot pixel: b = (−a) × LST(hot) + dt(hot). Actual evapotranspiration 173
Fig. 8 Rn of wheat farm (2016/11/11).
Fig. 9 Rn of wheat farm (2017/22/05). 174 Soghra Tofigh, Dariush Rahimi, Reza Zakerinejad
Fig. 10 Soil Heat Flux of wheat farm (G) (2016/11/11).
Fig. 11 Soil Heat Flux (G) of wheat farm (2017/22/05). Actual evapotranspiration 175
By replacing the unknowns in the dt equation, the Tab. 6 Results obtained from compute of for different dates. temperature difference of z and z is obtained. 1 2 Date ETinst inst), 2017/07/25 0.15 and Reference ET Fraction (ETrF) computation Latent Latent Heat Flux (λET), Instantaneous ET (ET - 2017/05/22 0.20 face due to evapotranspiration. It can be computed for 2016/11/11 0.03 eachheat fluxpixel is using the rate the followingof latent heatEquation: loss from the sur
n The Reference ET Fraction (ETr λET = R − G − H ratio of the computed instantaneous ET (ETinst) for 2 ). each pixel to the reference ET (ETF)r) iscomputed defined as from the Where λET is an instantaneous value for the time weather data: of theAn instantaneoussatellite overpass value (W/m of ET in equivalent evapo- ration depth is computed as: ETrF (19)
ETinst (17) Daily= values of ET (ET24) are often more useful than instantaneous ET. Where = 3600ETinst 3600 is the time conversion from seconds to hours, ET24 ETrF × ETr–24 (20) and λ is the latent is theheat instantaneous of vaporization ET or (mm/hr), the heat Where = ETr-24 is the cumulative 24-hour ETr for is computed as: the day of the image. This is calculated by adding the absorbed when a kilogram of water evaporates (J/kg) hourly ETr values over the day of the image (Waters 6 S (18) et al. 2002).
λ = [2.501 − 0.00236(T − 273)] × 10
Fig. 12 Flow chart of the computational steps of SEBAL algorithm (Bezerra et al. 2015). 176 Soghra Tofigh, Dariush Rahimi, Reza Zakerinejad
3.2.2 FAO Penman-Monteith method Evapotranspiration is obtained through Eq. 17 (Na RMSE (23) tional Irrigation and Drainage Committee 2008): = MSE (24) ETo (21) = MAD (25) = Where: ETo In all =these error detection indexes is the modeled Rn is the pure radiation entering the surface of the data and is the observational data and N is the data is the 2reference evapotranspiration (mm/day), count. 2 plantT is (MJ/m the mean/day), daily air temperature at 2 meters (°C),G is the soil heat flux (MJ/m /day), 4. Results u2 is the average daily wind speed at 2 m height In total the Chaharmahal and Bakhtiari province is es is the saturation vapor pressure (kPa), covered by circa 248,000 hectares of agricultural land, (m/s),ea is the real vapor pressure (kPa), 79,854 hectares of these are located in the Shahre- es a is the lack of saturation vapor pressure kord plain, of which, 58,553 hectares are under irri- (kPa), gation cultivation and 21,301 hectares are under rain Δ − e °C), fed cultivation. In total (in Shahrekord city), there γ are 29,917 hectares under cultivation of crops and is the slope of the vapor pressure curve (kPa/ 21,759 hectares dedicated to vegetable and garden 3.2.3 is MODIS the constant evapotranspiration coefficient. product farming. Evapotranspiration product is an 8-day composite According to the available statistics, the average dataset produced at 500-meter (m) pixel resolution. productivity of agricultural water use of the prov- The algorithm used for the Evapotranspiration data 3 while it should be increased to product collection is based on the logic of the Pen- 3 in the Iranian 5-years program and to man-Monteith equation, which includes inputs of ince is 0.873 in thekg/m Iranian 10-years planning. The water daily meteorological reanalysis data along with Mod- 1.51resources kg/m of this province decrease by 46 million m3 erate Resolution Imaging Spectroradiometer (MODIS) 1.9in average, kg/m annually (Ministry of Agriculture 2016). remotely sensed data products such as vegetation Due to global warming, water demand and relat- property dynamics, albedo, and land cover. The down- ed tensions regarding water supply increased. Being load data set is the MODIS evapotranspiration prod- aware of the fact that the knowledge about the actual uct for Aqua Satellite MYD16A2 and for Terra Satellite evapotranspiration is essential in water supply and MOD16A2. The MOD16A2 and MYD16A2 layers pro- management in this research, we assess the produc- vide the following products: tivity and determine the appropriate pattern of pro- i) The composited Evapotranspiration (ET), portional water resources under the present day cli- ii) Latent Heat Flux (LE), matic and hydrological conditions. iii) Potential ET (PET), iv) Potential LE (PLE) along with a quality control 4.1 Lysimeter layer. The pixel values for the two Evapotranspiration Evapotranspiration obtained for the initial growth layers (ET and PET) are the sum of all eight days with- in the composite period and the pixel values for the increases with the advance of the growth period and two Latent Heat layers (LE and PLE) are the average theperiod highest is 1.2 evapotranspiration mm/day. The evapotranspiration (4.09) is related volume to the of all eight days within the composite period. active growth period of the plants. Although the high- est temperature is recorded at the end of the growth 3.2.4 Model evaluation period, the evapotranspiration volume of this stage is To assess model performance a comparison is run lower (Table 8). between observations obtained from applying the three models and the SEBAL algorithm. These indica- 4.2 FAO Penman-Monteith method tors include R (Eq. 18), RMSE (Eq. 19), MSE (Eq. 20) and MAD (Eq. 21), as follows: The FAO Penman-Monteith estimates, before applying the crop factor, the highest evapotranspiration volume
R (22) Because at this time of the growth period, the plant is(7.43 completely mm/day) ripe in andthe noharvest irrigation time takes (2017/07/25). place. The = largest volume of the plant consists of seed and chaff. Actual evapotranspiration 177
Consequently, to estimate the most accurate evapo- Tab. 7 Estimated Evapotranspiration without crop coefficient rate. transpiration volume, the volumes obtained from the Vegetable FAO Penman-Monteith Date FAO Penman-Monteith model are corrected through coefficient method (mm/day)
- 11/11/2016 0.40 1.60 otranspirationcrop coefficients volume as shown is obtained in Tables for 7 andactive 8. plantAfter 22/05/2017 1.19 3.72 growthapplying time. the The crop minimum coefficient, evapotranspiration the maximum evap vol- ume is recorded in the initial period of growth. 25/07/2017 0.50 2.30
Fig. 13 Evapotranspiration of initial stage of wheat growth. 178 Soghra Tofigh, Dariush Rahimi, Reza Zakerinejad
4.3 MODIS evapotranspiration product - The results of the study of variations of wheat evap- otranspirationand 1.09 mm/day per on day the forTerra Aqua satellite satellites (Table and 8). otranspiration on the sample plots through MODIS 0.58The final mm stageevapotranspiration of wheat growth per with day 0.59 for Terramm evap sat- images show that the maximum evapotranspira- ellites had the highest evapotranspiration after the middle wheat growth stage. The lowest evapotran- growth stage. The average evapotranspiration at spiration occurs in the initial stages of growth (Fig- tion of wheat in the field occurred in the middle ures 13–14). this stage is 1.13 mm/day based on Aqua satellite
Fig. 14 Evapotranspiration of mid stage of wheat growth. Actual evapotranspiration 179
4.4 SEBAL algorithm 4.5 Data evaluation
The obtained SEBAL parameters are illustrated in The results of the error comparison and the correla- - the NDVI, Rn and ET are expressed in hot and cold pix- lated in Table 10. We show that the SEBAL algorithm Tableels indicting 9. As observed one of the in figurescritical (4–5,functions 8–9 andof this 16–17), algo- togethertion coefficient with satelliteof the above-cited imagery andmethods the FAO are tabuPen- rithm. man-Monteith method, after applying the vegetation
Fig. 15 Evapotranspiration of end stage of wheat growth. 180 Soghra Tofigh, Dariush Rahimi, Reza Zakerinejad
Fig. 16 Instantaneous ET of wheat farm (2016/11/11).
Fig. 17 Instantaneous ET of wheat farm (2017/22/05). Actual evapotranspiration 181
Tab. 8 Estimated Evapotranspiration with crop coefficient rate Tab. 10 Comparison of Lysimeter with Estimation of (mm/day). Evapotranspiration Methods. MODIS evapotranspiration FAO MODIS evapotranspiration FAO Evaluation SEBAL product Penman- product Penman Date Lysimeter criteria algorithm Monteith -Monteith Aqua Aqua method Aqua Terra R 1 0.86 0.98 0.98 11/11/2016 1.20 0.64 0.15 0.15 MAD 0.59 0.91 2.14 2.16 22/05/2017 4.09 4.42 1.13 1.09 MSE 0.36 1.29 5.24 5.34 25/07/2017 3.01 1.15 0.59 0.58 RMSE 0.60 1.13 2.29 2.31
Tab. 9 Results obtained from the farm for SEBAL algorithm index.
Date NDVI Rn (W/m2) Ts (K) H (W/m2) G (W/m2) Evp (mm/day) 11/11/2016 0.27 515.19 297.42 428.00 67.80 0.72 22/05/2017 0.66 569.53 307.30 419.00 55.00 4.80 25/07/2017 0.41 567.58 295.80 361.23 70.78 3.60
- method with MSE equal to 1.29, an MAD of 0.91, and a RMSE of 1.13 compared to the lysimeter output. coefficient,obtaining evapotranspiration are considered the through best methods the installation for situa In general, the results of this study indicate that oftions the where Lysimeter. sufficient instruments are not provided in applying remote sensing and satellite images, allows Although the correlation between the evapo- estimating evapotranspiration volume in areas with transpiration measured through MODIS images and Lysimeter is high (0.98), the error between their dif- particularly, a high correlation between the SEBAL ferent outputs is also quite high. The error of evapo- algorithmdata deficits. and The the results lysimeter reported data. in Consequently, this study reveal we transpiration measured by the Aqua Satellite in com- suggest using data derived with the SEBAL algorithm parison to the Lysimeter data is less than for the Terra in areas with similar environmental conditions where satellite products. no data are available and as input information for a further assessment of hydrological process dynamics. We show that the results of this study can be applied 5. Conclusion in studies of water resources management and appro- priate irrigation management on farm level. Agriculture development and food security are threatened by a decrease in rainfall, rising tempera- tures, droughts, a decrease in water table level and an increase in evaporation. Evapotranspiration is an References effective parameter in providing water balance and Allen, R. G., Pereira, L. S., Raes, D., Smith, M. (1998): Crop food security because of its contribution in deter- evapotranspiration – Guidelines computing crop water mining plant water need. Accurate estimations of the requirements. FAO Irrigation and drainage paper No. 56, water needs and water supply for plants, especially FAO. for wheat, as a strategic product in Iran, are of major Alizadeh, A. (2011): Principles of Applied Hydrology. interest. Evapotranspiration estimations by experimen- tal methods and different algorithms as used in this 33rd ed. Imam Reza University, Mashhad. paper is an important step forward especially in data Atasever,Actual U.Evapotranspiration H., Ozkan, C. (2018): Mapping A New and SEBAL On-Site Approach scarce areas. However, a careful validation of the dif- Application.Modified with Journal Backtracking of the Indian Search Society Algorithm of Remote for ferent methods should be carried out using observed data e.g. from Lysimeter. The evapotranspiration esti- Bezerra,Sensing B. 46,G., Da 1213–1222, Silva, B. B., https://doi.org/10.1007 Dos Santos, C. A. C., Bezerra, mations obtained through the experimental models G./s12524-018-0816-9. R. C. (2015): Actual Evapotranspiration Estimation and the SEBAL algorithm revealed that the SEBAL algorithm have highest correlation and the least error Approaches. Advances in Remote Sensing 4, 234–247, with the observed Lysimeter data. The error values Using Remote Sensing: Comparison of SEBAL and SSEB of MSE is 0.36, of is MAD is 0.59, and of the RMSE is Fu, Q., Liu, W., Li, T., Liu, D., Cui, S. (2018): Study of the water 0.60 in terms of the SEBAL algorithm. Comparable savinghttps://doi.org/10.4236/ars.2015.43019. potential of an irrigation area based on a remote good values we found for the FAO penman Monteith sensing evapotranspiration model. Arabian Journal of 182 Soghra Tofigh, Dariush Rahimi, Reza Zakerinejad
Rasmussen, M. O., Sørensen, M. K., Wu, B., Yan, N., Qin, H., Sandholt, I. (2014): Regional-scale estimation Ghorbani,Geosciences A., Faramarzi, 183, 1–10, M., https://doi.org/10.1007 Karami, J., Gholami, N., of evapotranspiration for the North China Plain Sobhani,/s12517-018-3504-8. B. (2015): Comparative Evaluation of SEBAL using MODIS data and the triangle-approach. and Metric algorithms in Estimating Evapotranspiration International Journal of Applied Earth Observation and (Case Study: Malayer City). Journal of Spatial Planning 19, 184–153. (In Persian). Jaber, H. S., Mansor, S., Pradhan, B., Ahmad, H. (2016): Rawat,Geoinformation K. S., Bala, A., 31, Singh, 143–153, S. K., Palhttps://doi.org/10.1016 R. K. (2017): Evaluation of SEBAL Model for Evapotranspiration /j.jag.2014.03.017. mapping: A casestudy from Bhiwani District of Haryana, International Journal of Applied Engineering Research India.Quantification Agricultural of wheat Water crop Management evapotranspiration 187, 200–209, and 11,Mapping 3950–3955, in Iraq Using Remote Sensing and GIS. Jian-Ying, Y., Xu-Rong, M., Zhi-Guo, H., Chang-Rong, Y., Hui, Rezaei Banafsheh, M., Jahanbakhsh, S., Valizadeh Kamran, J, Feng-Hua, Z. H., Qin L. (2015): Water consumption in K.http://dx.doi.org/10.1016/j.agwat.2017.03.015. (2014): Estimation of Actual Evapotranspiration in summer maize and winter wheat cropping system based on SEBAL model in Huang-Huai-Hai Plain, China. Journal 241–262. Santos,Sefidrood C. A. G., Basin Silva, with R. M. Satellite. D., Silva, Image A. M., Processing Neto, R. M. 13(44), B (2017): Estimation of evapotranspiration for different Kamali,of Integrative M. I., Nazari, Agriculture R. (2018): 14(10), Determination 2065–2076, of maizehttps:// land covers in a Brazilian semi-arid region: A case waterdoi.org/10.1016/S2095-3119(14)60951-5. need using remote sensing data and SEBAL study of the Brígida River basin, Brazil. Journal of South algorithm. Agricultural Water Management 209,
Karbasi, M., Moghaddam, M., Nikbakht, J., Kaviani, A. Sun,American Zh., Gebremichael, Earth Sciences M., Ardö, 74, 54–66, J., Nickless http://dx.doi.org A., Caquet, (2016):197–205, Estimation https://doi.org/10.1016/j.agwat.2018.07.035. of actual evapotranspiration of the B.,/10.1016/j.jsames.2017.01.002. Merboldh, L., Kutschi W. (2012): Estimation of daily plant using SEBAL Algorithm (Case Study: Khoramdareh Region in Zanjan Province). Journal of Ecohydrology 3, 10.22059 evapotranspiration over Africa using MODIS/Terra and . Sun,SEVIRI/MSG Zh., Wei, B., data. Su, W., Atmospheric Shen W., Wang, Research Ch., You, 112, D., 35–44, Liu Li, 427–437.Zh., Liu, X., (In Ma, Persian). T., Kejia, https://doi.org/ D., Zhou, Q., Yao, B., Niu T. (2013): Zh.https://doi.org/10.1016/j.atmosres.2012.04.005. (2011): Evapotranspiration estimation based on Retrieval/IJE.2016.60031 of the surface evapotranspiration patterns in the SEBAL model in the Nansi Lake Wetland of China, the alpine grassland–wetland ecosystem applying SEBAL Mathematical and Computer Modelling 54, 1086–1092, model in the source region of the Yellow River, China. Wagle, P., Bhattarai, N., Gowda P. H., Kakani, V. G. (2017): https://doi.org/10.1016/j.mcm.2010.11.039. Meyabadi,Ecological A., ModellingAlizadeh, A., 270, Sanaynejad, 64–75, http://dx.doi.org S. H., Ghahreman, B., estimating daily evapotranspiration in high biomass Davari/10.1016/j.ecolmodel.2013.09.004. K. (2017): Evaluation of SEBAL Algorithm for sorghum.Performance ISPRS of fiveJournal surface of Photogrammetry energy balance models and Remote for Estimation of actual Evapotranspiration of Neyshabour –
Research in Agriculture 30(4), 525–541. (In Persian). Waters,Sensing R., Allen128, 192–203, R., Tasumi http://dx.doi.org/10.1016 M., Trezza R., Bastiaanssen W. Rokh Basin Using SWAT Model. Journal of Water. (2002):/j.isprsjprs.2017.03.022. S E B A L (Surface Energy Balance Algorithms Morshedi, A., Tabatabaee, S. H., Naderi, M. (2016): Validationhttps://doi.org/ of Estimation10.22092/JWRA.2017.109014 of Evapotranspiration from Implementation. Yu,for B., Land)Shang AdvancedS., Zhu, W., Training Gentine andP., Cheng User’s Y. Manual, (2019): Idaho Data. Water and Soil Journal (Science and Technology Mapping daily evapotranspiration over a large irrigation SEBAL and Hargreaves-Samani Models Using Lysimeter district from MODIS data using a novel hybrid dual- source coupling model. Agricultural and Forest MinistryAgricultural) of Agriculture 23, 367–376 (2016). (In Persian), https://doi.org National/10.22067/jsw.v30i2.39762. Irrigation and Drainage Committee (2008). Ndou, N. N., Palamuleni, L. G., Ramoelo, A. (2018): Modelling ZareMeteorology Haghi, D., Neyshabouri, 276–277, https://doi.org/10.1016 M. R., Majnooni-Heris, A., depth to groundwater level using SEBAL-based dry Jalilian,/j.agrformet.2019.06.011. Z. (2016): Planning, Water and Soil Science, season potential evapotranspiration in the upper Molopo 26(1–2), 1–30. River Catchment, South Africa. Egyptian Journal of Zakerinejad, R., Maerker, M. (2015): An integrated assessment of soil erosion dynamics with special emphasis on gully erosion in the Mazayjan basin, Pourmohammadi,Remote Sensing S., and Rezaanian, Space Science M. T., 21,Mokhtari, 237–248, M. H.,https:// southwestern Iran. Journal of Natural Hazards 79, 25–50, Rahimian,doi.org/10.1016/j.ejrs.2017.08.003. M. H. (2010): Determination and zoning of actual evapotranspiration by remote sensing technique Zakerinejad, R., Masoudi, M. (2019): Quantitative mapping and SEBAL algorithm (Case study: Menshad watershed https://doi.org/10.1007/s11069-015-1700-3. in Yazd province). Iranian Journal of Watershed model: a case study in the Mazayejan Plain, Southwest Management Science and Engineering 13, 21–30. (In of desertification risk using the modified MEDALUS Persian). Iran. AUC Geographica 54(2), 232–239, https://doi.org /10.14712/23361980.2019.20. Original Article 183 Are the Czech or Slovak regions “closer to Europe”? Pro-Europeanness from a subnational perspective Martin Plešivčák*
Comenius University in Bratislava, Faculty of Natural Sciences, Department of Economic and Social Geography, Demography and Territorial Development, Slovakia * Corresponding author: [email protected]
ABSTRACT Based on the 2003–2019 electoral data, this article evaluates the level of pro-Europeanness in Czechia and Slovakia at the regional and sub-regional levels during and after their EU accession period. The TOPSIS multi-criteria evaluation method and cluster analysis were used to quantify the pro-European levels and to create the subsequent categories of territorial units. The results show sup- port for the ideas of European integration primarily in large urban regions (Prague, Brno, Bratislava, Košice), territorial units with a higher concentration of ethnic minorities, larger scale agricultural activities (southwestern Slovakia), and a high degree of religiosity (northeastern Slovakia). The low level of pro-Europeanness was predominant in the less developed north-western Czechia and parts of Moravia. In Slovakia, the Eurosceptic regions were mostly located in the northwest, where the values of statism, egalitarianism and nationalism have a strong tradition. This approach can be used to identify areas of weak support for the EU project at a spatially disaggregated level in other EU countries.
KEYWORDS Pro-Europeanness; TOPSIS method; EU referendum; European Parliament elections; subnational level; Czechia; Slovakia
Received: 4 October 2019 Accepted: 16 July 2020 Published online: 30 September 2020
Plešivčák, M. (2020): Are the Czech or Slovak regions “closer to Europe”? Pro-Europeanness from a subnational perspective. AUC Geographica 55(2), 183–199 https://doi.org/10.14712/23361980.2020.13 © 2020 The Author. This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0). 184 Martin Plešivčák
1. Introduction real assistance from European structural and invest- ment funds. Indeed, if we want to avoid disintegration Czechia, Slovakia and eight other mainly post-social- processes within the EU and the threat of its gradual ist countries of Central and Eastern Europe, joined decomposition, just in those regions which are char- the European Union on May 1, 2004. In the periods acterized by the highest degree of Euroscepticism, the before and after accession, the moods in the two EU’s contribution to the future should be the most vis- countries in relation to the European integration pro- - ject differed based on time, location, and the politi- - cal and socio-economic conditions. Public support ticismible. This at theis prevented regional and by the sub-regional knowledge levels. of “problem for the country’s accession to the EU, as well as the atic” regions and the reasons that can cause Euroscep referendum turnout itself differed for both countries, 2. Theoretical background can be said for the period after May 1, 2004, in the depending on their intraregional specifics. The same for Eurosceptic parties in the European Parliament politics, in the case of European elections, is the sec- (EP)case ofelections. the territorial specificities of turnout and votes ond-orderThe dominant national concept election of understanding theory (Reif and European Schmitt We have seen a gradual shift of Euroscepticism 1980). European Parliament elections were character- toward the centre of European politics since the ized as national elections taking place simultaneously in all the member states of the European Community. both enlargement and elections in 2004. For the new member1990s. Nevertheless, states, there it was gained high stronger public support influence for after the of major (general) national elections and are domi- European Union project; however, it began to decline natedThese byare the less same important parties elections that focus held more in theon national shadow after accession. In addition to contextual factors such interests than European issues or on the position of the national parties to the EU (Hix and Marsh 2011). rhetoric of several political actors could have con- Nevertheless, the results of the national and Europe- as the financial and migration crises, changes in the an elections are still different. Second-order elections right-wing political parties learned from the past and do not lead to national government formation and insteadtributed of to proclaiming the rise of radicalEuroscepticism. views (including The extreme xeno- are therefore of less interest to voters, the media and phobia, racism and anti-Semitism) they moved on to a somewhat milder and more tolerant populism and decision than if the national elections were organized thus Euroscepticism was addressed to a wider elec- political actors. This leads voters to make a different- ate votes on the basis of what they think about the values or Euroscepticism include citizens on the one country’son the same economic day (Hix situation, and Marsh the 2011). government’s The elector per- hand,torate and (Goodwin their political 2011). representatives,The bearers of pro-European political par- formance or the topics that move domestic politics at
- et al. 2011). For this reason, the campaigning and tac- tifyties andthe degreespecific ofcandidates pro-Europeaness on the other. of the public in ticsthat ofmoment political (Carrubba parties in and second-order Timpone 2005; elections De Vries are CzechiaThe main and Slovakiaobjective at of thethis time contribution of accession is to to iden the motivated by national themes (Reif and Schmitt 1980; EU and afterwards at the hierarchically lower, region- Cabada 2010). Also, media coverage of these elections al and sub-regional (district) levels, to create a typ- is usually limited. Furthermore, immediately after the election and the end of the campaign, the European long-term pro-European orientation and to charac- Parliament returns to obscurity (Lodge 2010). In the ificationterize the of groups partial formed territorial on theunits basis according of its indica to the- is emphasized the fact that one of the main aspects of decision making method regarding the value distance theanalysis second-order of the first elections European is that elections there is(1979), less at therestake oftors. given For territorial this reason, unit we to apply the most TOPSIS positive multi-criteria and most negative value within the set of units under study. We the national parliament or the government, and the also try to outline possible factors related to pro-Eu- electorate(Reif and Schmitt is not 1980).highly Themotivated election to does participate not involve in ropeanness at the regional and district levels, dis- the election or to vote differently than they would if - - tion theory predicts that elections to the European Pinkcussing 2012; with Voda findings 2015; of Kostelecký previous studies et al. 2016; on elector Przy- Parliamentnational elections follow werethree held.main The formulas: second a orderlower elec rate al behaviour (e.g. Krivý et al. 1996; Madleňák 2012; of participation (lower turnout), a more positive out- of the sub-national units of these countries, which is come for small and new parties, and a loss of support particularlybyla 2019). Thus, important the aim in relationis to map to the running “Europeanity” an effec- for government parties stemming from the location of general elections in the national election cycle (Reif promotionaltive, place-specific activities, campaign and in focusing the context on EU of relevance regional of the theory in the context of post-socialist countries (cohesion)and benefits, policy, the also meaningful due to the dissemination proper direction of itsof hasand beenSchmitt addressed 1980; Hix by severaland Marsh authors, 2007). e.g. The Linek validity et al. Pro-Europeanness from a subnational perspective 185
Šaradín (2008), Havlík and Hoskovec (2009),
(2007), Charlotdomestic 1986). and European Nevertheless, issues, their which importance was confirmed varies CabadaIn order (2010), to understand Klíma and theOutlý perception (2010), Kovář of the (2013), nature fromin the election first direct to election elections and in from 1979 place (Blumer to place. 1983; In Kovářof the Europeanand Kovář Parliament(2014). elections by the political general, however, the importance of European themes actors themselves (parties, politicians and voters), gains momentum over time, but this is in contrast to the connection of these elections with the concept of the declining voter participation in the European elec- Europeanization must be recognized. Europeaniza- - opment, food security, environment, foreign policy, and content of politics in such a way that the political immigration,tions. Traditional and economicissues, such and as industrial research policyand devel may andtion economicis defined dimensionsas a process ofof thechanging European the directionCommu- - ert Schuman Foundation 2004). Nevertheless, the cur- shape the content and scope of national policies (Lan- rentbe considered pan-European as “more themes European” of the (Budgegiven period, 2001; Robsuch drechnity influence 2002). It the is clear structure that European of national issues politics, have and an increasing effect on political debates at the national crisis, Brexit, migration, reviving debates on the con- level, but the response of national party systems is ceptas, the of aeconomic two-speed and Europe, financial the futurecrisis, ofGreece’s EU regional debt policy, etc., as well as the hot topics of the domestic - political scene, or analyses of election results in key very limited, without a more significant formation of- European integration countries play a no less impor- (new) parties built on a “European basis”. In this con tant role in this context. politicaltext, it is worthactors mentioninginvolved in the the term existing “European social struccleav- tureage”, should around ideally which be the shaped pre-election (ideally, strugglereal awareness of the In general, five basic forms of party Euroscepticism - majorcan be and identified government (Kopecký parties and with Mudde a critical 2002; attitudeTaggart trastof the to dimension the more traditional of the “Europeanness” concept of cleavages concerning (Lip- and Szczerbiak 2002). The first category concerns these elections). The term European cleavage, in con the functioning and direction of the Community itself. a concept for and against a territorial integration pro- towards the definition of European policies as regards jectset and within Rokkan Europe 1967), (the canEU as in a principle centrally be oriented defined and as - bureaucratic superstate, with a common economic or tationThis moderate of the socially form concerns oriented the EU attitude program, of whilethe social the strongerdemocratic form parties is represented toward the by inadequate conservative implemen or neo- one hand, or as a concept of a more lenient bundle of liberal parties who call for less regulation and inter- states,financial for (monetary, example onfiscal a customs or budgetary) union policybasis onor the vention by EU institutions. It therefore has a more ide- existence of a common market, on the other). Accord- also concerns government parties and deals with the on three levels: general (for and against the EU as a Europeanological character. integration The projectsecond typein terms of Euroscepticism of a reasona- territorialing to Bartolini integrator (2007), = thisindependence term can be / understoodintegration ble (or necessary) degree of Brussels interventions in dimension), constitutive (the cultural level – who can policies implemented at the national level. It is there- access, division of competences between EU insti- fore a pragmatic debate on the division of competenc- tutions and member states themselves, community es and agendas between EU institutions and member decision-making mechanism) and isomorphic (ide- - ological issues – liberalism, protectionism, welfare ed by a radical opposition to mainstream and govern- state, immigration policy, civil rights, etc.). We could mentstates. parties, The third in typethe senseof Euroscepticism of protest and is represent enforcing assume that those who evaluate the process of Euro- one’s own, far-right or far-left-oriented ideology, in pean integration in their particular life as a disap- the sense of maintaining national interests and iden- pointment tend to ignore the European Parliament tity on the one hand, and protesting against excessive elections, or they use the protest vote against the liberalization of the public sector and economic rela- mainstream parties and support Eurosceptic political forces of an extreme right- or left-wing orientation. - damentallytions with negative restoring social the functioning consequences of the on Communi the other.- electorate and the political parties in the European tyThe as fourth a whole, category the policy of Euroscepticism of exiting or refusing is aimed entry at fun to ParliamentThe question elections remains place regarding an emphasis how muchon address of the- the EU. However, in principle, parties of this type do ing internal political issues and how much the current challenges of European integration or the content of scene, but, somewhat paradoxically, in some countries - theynot have enjoy a significant a more prominent political impactrole particularly on the domestic in the ever, evidence that the Eurosceptic and pro-European - Europeanparties place policies importance themselves on European matter. issues,There is,and how this scepticism is relatively marginal in terms of political representationEuropean Parliament and time elections. span, andThe islast often type connect of Euro- in terms of better electoral results. An election cam- paignapproach for thehas European ideological Parliament and practical features significance a mix of ed to specific persons or “single use” parties. These formations do not require a fundamental reform of 186 Martin Plešivčák relations within the EU but rather focus on subtle ĽS-HZDS’s topics such as transparency and accountability for the staff of EU institutions, effective use of EU resources, laterEU at electorate, the time (“at along the withlast minute”); SNS and laterbut especially salaries of MEPs, etc. ĽS-NSpolicies in office as well as the profile of the eng.then Peo and- ple’s Party – Our Slovakia) proved to be Eurosceptic to the largest (Ľudová extent strana when – Našecomparing Slovensko, all of the major 3. Methodology in the case of the European Union membership ref- In the empirical part of the paper, we use two types of political parties in Slovakia (Gyárfášová 2007). Even electoral geographic data: the referenda on the coun- waserendum also presentin 2003, in a lowthe pastturnout (in comparison(52%) demonstrably with the European Parliament elections held in 2004, 2009, currentconfirmed ĽS-NS the factat the that time Euroscepticism of the impending in Slovakia mem- 2014try’s accession and 2019. to In the the EU second, in 2003 we (“yes also to work EU”) with and the bership in a more latent form). We do not consider results of the Eurosceptic parties. We put them into it to be strictly correct and the only possible to refer this group based on the work of several authors and to the divisions of various authors in this matter, as electoral programs of the parties themselves. In the Šárovec (2018) describe the SaS as case of the Czech political situation, we used several a Eurosceptic party, while a large part of its elector- atee.g. supportsHynčica and the EU project. We also consider it rel- Šaradín 2004; Baun et al. 2006; Fiala et al. 2006; Linek studies for the inclusion of the partiesík (Bradová and 2008; Hloušek and Pšeja 2009; Hricová 2009; Havlík theatively 2004 subjective European to Parliamentperceive parties elections. such With as KDH these or 2010;et al. 2007;Kovář 2014). Havlík In 2008; the case Havl of andSlovakia, Vykoupilová we also examplesSMER-SD (Hendersonwe wanted 2008)to illustrate as Eurosceptic, the fact thateven thefor
- problematic and to some extent subjective matter, in usedgraphical relevant context literature of electoral for this support purpose for (Gyárfášová the parties whichclassification different of partiesauthors as can Eurosceptic work with is differentreally a very set andin European Velšic 2004; elections Gyárfášová 2007). We refer to geo of Eurosceptic parties in research, as a certain degree differentiation of the extreme right support (Mikuš (Plešivčák 2015) and the spatial- avoid. ter of the electorate was the decisive criterion for of For subjectivity each election of classification of the period under is not study possible (2004– to theand inclusionGurňák 2016; of a political Mikuš partyet al. among2016). theThe group charac of - Eurosceptic parties. Recessive parties were not taken ceptic parties: into account. Given the above division of parties into 2019), we identified the following parties as Euros types according to the degree of Euroscepticism and Czechia their relevance within the party system, we would like to mention the cases when the inclusion of parties to 2004 European Parliament Elections the Eurosceptic, or their exclusion, was a problematic strana, eng. Workers’ Party matter. We emphasize that in classifying the parties, strana Čech a Moravy, eng. Communist –Party DS (Dělnická of Bohe- we primarily took into account the nature of the par- mia and Moravia ), KSČM (Komunistickáeng. National Coalition eng. Independent), then leaders. Among other parties, in the case of the ), NARKOA (Národníeng. koalice, Republicans ty’s electorate rather than the official rhetoric of its of Miroslav Sládek),) NEZ (Nezávislí, Úsvit (2014) to be the Eurosceptic parties (more pre- RMS2009 (Republikáni European Miroslava Parliament Sládka, Elections ciselyCzechia, parties we also with considered a predominantly KSČM (2004–2019) Eurosceptic elec and- eng. National Party), torate), although some authors label them soft Euro- eng. Party – DS, ofKSČM, Free sceptic (Havlík and Kaniok 2006; Kaniok and Havlík CitizensLibertas.cz, NS (Národní strana, - SSO (Strana svobodnýchČeskoslovenska, občanů, eng. Association for the Republic), SPR-RSČ – Republican (Sdružení Party pro republikuof Czechoslovakia – Repub), 2016). The electorate of these two parties is indeed Suverenitalikánská strana (eng. Sovereignty) strana,quite different eng. Civic from Democratic another partyParty often). In the associated case of 2014 European Parliament Elections – Česká suv- with Eurospeticism, ODS (Občanská demokratická erenita (formerly Suverenita, eng. Czech Sovereignty), Eurosceptic Party Group, based on several arguments (Baunthis party, et al. we 2006), finally as decided this party to not can include rather itbe in per the- eng. Work- ers’DSSS/SPE Party (formerlyof Social Justice/NoDělnická strana, to Brussels Dělnická Dictate! strana), sociální spravedlnosti/Ne diktátuČeskoslovenska, Bruselu!, eng. moreceived educated as “pro-European voters living with in the reservations”. urban environment Given Communist Party of Czechoslovakia morethe ODS in favourelectorate, of the which EU project, is largely we characterizedhave chosen not by KSČM, KSČ (Komunistickáeng. strana No to Brussels – Nation- to include the party into the Eurosceptic group. In the al Democracy ), ND (Ne BruseluČech, case of Slovakia, there was a problem with ĽS-HZDS Moravy– Národní a Slzska,demokracie, eng. Republican Party of Bohemia, and SNS, parties that even expressed support for the Moravia and Silesia), RSČMS (Republikánská strana
), SSO, SZR-NE (Strana zdravého Pro-Europeanness from a subnational perspective 187 rozumu – Nechceme Euro, eng. Party of Common Sense we decided to work with the index of pro-European- – We Don’t Want the Euro), Úsvit ness. When constructing, in addition to supporting (eng. Dawn of Direct Democracy) 2019 European Parliament příméElections demokracie – ANS supporting Eurosceptic parties, as is commonly used. eng. Alliance of National Forc- EU accession by Referendum 2003, we considered es Českou republiku, eng. per cent minus the support for the Eurosceptic par- (AlianceAlternative národních for Czech sil, Republic), Česká suverenita, We could not automatically work “with the rest” (100 ), APAČI (AlternativaČesko pro (eng. Independents/Joyful group of parties would include much more hetero- Czechia), ČSNS/Patrioti ČR geneousties) as a political % for pro-Europeanness, entities (in relation as toa much the level wider of Svobodní/RadostnéČeské republiky, eng. Czech Nation- al Social Party/Patriots of (ČeskáCzech Republic strana národně (Konzervativnísociální/Patrioti alternativa, eng. Conservative Alter- EU support) than for parties defined as Eurosceptic. native eng. Moravians), KOALPrvní pro-EuropeannessThe settings of the vs.model votes calculation for Eurosceptic in the parties)TOPSIS republika (eng. First Republic - method technically solve this “discrepancy” (index of losti České), KSČM, republiky, Moravané eng. Independence ( Party of the of input variables with sensitivity to their orientation Czech Republic ), SNČR (Strana nezávis very easily and is based on evaluating the influence eng. Workers’ Party decrease in the value of index of pro-European and of Social Justice/National), DSSS/NF Front(Dělnická strana sociální- (increasing value of something “negative” means a- muraspravedlnosti/Národní (formerly Úsvit fronta, tation of the variable – for our research a positive ), SPD –eng. Tomio Freedom Oka orientationvice versa. Thein supporting researcher the sets country’s the desired accession orien and Direct Democracy přímé – Tomio demokracie, Okamura Svoboda a to the EU, and a negative orientation in supporting přímá demokracie – Tomio eng.Okamura, Reasonables - the Eurosceptic parties, which in both cases means ), SPR-RSČ,eng. increasing the value of the index of pro-European- NationalRozumní/ND Democracy (formerly) SZR, /Národ ness). ní demokracie, formerly Právo a Spravedlnost, Slovakia - ness thus were as follows: 2004 European Parliament Elections – KSS (Komu- The variables entering the index of pro-European eng. Communist Party – the more the better of Slovakia), ĽS-HZDS b)a) Votes forfor Eurosceptic accession to parties the 2003 in the European 2004 European Union nistická strana Slovenska,eng. People’s Party – Move- referendum (%) – the less the better ment for Democratic Slovakia (Ľudová strana – Hnutieľudová za c) Votes for Eurosceptic parties in the 2009 European demokratickéstrana, eng. Slovak Slovensko, People’s Party), SNS/PSNS (Slov- Parliament elections (%) – the less the better ), SĽS (Slovenská - d) Votes for Eurosceptic parties in the 2014 European na, eng. Slovak National Party/True Slovak National Parliament elections (%) – the less the better Partyenská) národná strana/Pravá Slovenská národná stra e) Votes for Eurosceptic parties in the 2019 European 2009 European Parliament Elections – KSS, Parliament elections (%) – the less the better ĽS-HZDS, SNS 2014 European Parliament Elections – KSNS WeParliament obtained elections data for (%) the state, regional and dis- eng. Christian trict levels from the databases of the Czech Statistical Slovak National Party), KSS, ĽS-NS (KresťanskáNaše Slovensko, slovenská eng. People’s národná Party strana, – Our Slovakia), (Ľudová strana – SlovakOffice and Republic the Statistical 2019). Office of the Slovak Republic Nation and Justice – Our Party Úsvit (eng. (Czech Statistical Office 2019; Statistical Office of the NaS-NS (Národ a Spravodlivosťeng. – Defiancenaša strana, – Labour eng. Preference by Similarity to Ideal Solution) to evaluate Party) ), SĽS, SNS, theWe position use the of TOPSISthe regions method and (Techniquedistricts under for Orderstudy Dawn),2019 Vzdor European – strana Parliament práce ( Elections – Kotleba – in mutual comparison based on the values of the set of ĽSNS (formerly Ľudová strana – Naše Slovensko, eng. Kotleba – People’s Party Our Slovakia), KSS/Vzdor – score for the index of pro-Europeanness to rank the mentionedindicators mentionedterritorial units. above. For This the methodneed of thegenerates empir- eng. Slovak People’s Part of Andrej Hlinka - ical part of the paper, we decided to use this method, strana práce, SĽS Andrejaeng. We HlinkuAre Family (formerly – Boris SĽS, Kollár eng.), which in relation to the objectives of the work can be ), SME RODI - stranaNA – Boris vlastencov, Kollár ( eng. Slovak National Unit – Patriot essary to use other methods, e.g. factor analysis). Giv- PartySNJ-sv), (formerlySNS KSNS, Slovenská národná jednota – enassessed that in asthis adequate part of the (for paper this reason,we decided it was to evaluatenot nec a set of variables indicating the degree of pro-Euro- We wanted to approach the issue from a positive peanness across the regions and districts of Czechia perspective, based on support for the EU project the multicriteria evaluation tools can be considered and Slovakia, the use of TOPSIS method as one of (Euro-optimistic), and not Euro-sceptically. Therefore, 188 Martin Plešivčák desirable. In addition, if we work with several territo- and the negative ideal solution rial units, in this case 20 at the regional and 149 at the district level, the use of this method is the right choice, Dj = (minj wij), j = 1, 2, … , k because in the case of a given territorial unit it takes into account the level of each input variable to ideal and to the least desirable value within the set of units solution by using the formula as follows: (i.e. with respect to the value of the most successful 4. To calculate the distance from the positive ideal and the least successful region or district). , i = 1, 2, … , p, Accelerators increasing the value of the pro-Euro- peanness index were the high values of indicator a and from∑ the negative ideal solution by using the ∑ (the higher the better) and the low values of indica- formula below: tors b-e (the lower the better). When calculating index (in scale from 0 to 1), each input indicator (a-e) was , i = 1, 2, … , p, ∑ (Hwang and Yoon 1981) is The Euclidean∑ distance measure was utilised to cal- consideredequally weighted, one of by the 1/5 most (0.2). classical multi-criteria culate the distance. decisionThe TOPSIS making method methods ideal solution by using the formula below: It constitutes a collection(Opricovic of shortcut and Tyeng methods 2004; 5. To calculate the relative distance from the negative designedShih et al. to 2007; minimize Manokaran the distance et al. 2011).from the ideal solu- , i = 1, 2, … , p
Variants are then arranged in descending order istion. then These the onemethods which use according an ideal variantto the selected as the object met- according to the ci values. ricsof aspiration. is the closest The toselected the ideal “best” option. compromise variant - It provides a complete ordering of all variants. ate groups of districts based on the pro-Europeanness Subsequently, we used the cluster method to cre matrix as well as the weight vector of individual cri- data and information has led to the need to devel- To resolve the problem, the multi-criteria decision opindex methods (Hastie to et clarify al. 2016). and Theclassify increasing them. In amount addition of method is to identify the variant that is closest to the positiveteria has ideal to be solution, determined. and farthestThe main from principle the negative of this ideal solution. numberto other ofclassification clusters, with methods, objects incluster one cluster analysis having has similarbegun to properties, be used. Thisand objectsmethod in produces different a clusters certain having as many different properties as possible. N matrixThe calculation procedure is as follows. objects denoted by indexes 1 < i < N, which have d fea- 1. To calculate the normalized multi-criteria decision turesThe indexed input for as 1cluster < j < d analysis is represented by R = (rij) to the N × d matrix: . These data are used to write using the formula:
rij = , i = 1, 2, … , p, j = 1, 2, … , k
∑ After this transformation, the columns in the matrix are vectors of unit size by Euclidean metrics.
matrix Line d-dimensional vector xi is a vector of the i-th 2. To calculate the weighted multi-criteria decision object, while element xij denotes the value of the j-th W = (wij) feature of the i-th object. - how the j-th column is multiplied by the appropri- al steps. 1. Selecting and extracting the features, 2. ate weight, as follows The cluster analysis is comprised of four gener- uating the results. wij = (vjrij) Selecting the algorithm, 3. Verifying accuracy, 4. Eval to conduct the clustering. As a result of clustering, basedThe on IBM the SPSS values Statistics of the index 22 programme of pro-Europeanness, was used
3. HToj =determine (maxi wij), the j = positive1, 2, … , kideal solution generated. five groups of districts with internal similarity were Pro-Europeanness from a subnational perspective 189
4. Analysis, results and findings 2) are the least inhabited areas. 4.1 Levels of Territorial Units and the Banská Bystrica Region (68 inhabitants/km 4.1.2 Districts 4.1.1 Regions At the regional level, we analysed 22 spatial units, 14 regions in Czechia and eight regions in Slovakia caseAt the of districtSlovakia, level, the municipalwe worked districts with 149 of Bratislava units, 77 - (5)in Czechia and Košice and (4)72 inwere Slovakia connected (Fig. to2, Tab.one district1). In the in the entire city in order to strengthen the compara- (Fig.different, 1). These as the units smallest also regionrepresent has the an NUTSarea of 3 onlylev tive value of the analysis with the other districts of 496.10el. Their km territorial2 and population size is markedly 2 (Central Bohemia (the city Region). of Prague In Slovakia, as a separate the smallest NUTS 3 ofthe 1,945.69 countries. km 2The, while largest the Levicedistrict district in Czechia in the is Nitra the regionregion) is while the Bratislava the largest Region has an (2,052.5 area of 11,014.97 km2), and km the districtRegion (1,551.1of Klatovy km in2 the Plzeň Region, with an area 2). In other hand, the Brno-mě 2) is terms of population, the city of Prague is the smallest ) is the largest in Slovakia. On the- regionlargest inis theterms Banská of area, Bystrica but it Region has largest (9,454.4 population km 2) issto the district smallest (230.22 in Slovakia. km theIn Czechia, smallest the in Czechiacity of Prague and the has Kysucké the largest Nové Mespop- Vary Region has the smallest population (data as of to district (173.7 km (1,301,135 inhabitants) in Czechia, while the Karlovy- kia, the differences between regions are also smaller ulation (1,301,135 inhabitants); on the contrary,- inDecember this indicator, 31, 2018, as the 295,686 region withinhabitants). highest number In Slova of the Jeseník district in the Olomouc Region has the inhabitants with permanent residence is the Prešov smallest population (38,330 inhabitants as of Decem Region (825,022), and the region with the lowest theber Medzilaborce31, 2018). The district most populousin the Prešov area region in Slovakia is the is the city of Bratislava (432,864 inhabitants), while- - estnumber population of inhabitants density is are the the Trnava city ofRegion Prague (563,591 (2,622 (2,622least populous inhabitants/km (11,8962) dominates inhabitants in Czechia, as of Decem while inhabitants/kmas of December2 31,) in 2018).Czechia The and areas the Bratislava with the Region great theber Prachatice31, 2018). Indistrict terms in of the population South Bohemian density, PragueRegion 2 2) is the least populated. In Slo- 2) 2) has the (321 inhabitants/km ) in Slovakia. On the contrary, (37 inhabitants/km the South Bohemian Region (63 inhabitants/km vakia, Bratislava (1,177 inhabitants/km
Fig. 1 Territorial composition of NUTS 3 regions in Czechia and Slovakia. Source: Czech Statistical Office, Statistical Office of the Slovak Republic (2019). 190 Martin Plešivčák
Fig. 2 Territorial composition of districts in Czechia and Slovakia. Source: Czech Statistical Office, Statistical Office of the Slovak Republic (2019).
Tab. 1 Order number of districts in Czechia and Slovakia. Order District Region Country Order Number District Region Country Number 26 Karviná Moravian-Silesian Region Czechia 1 Benešov Central Bohemian Region Czechia 27 Kladno Central Bohemian Region Czechia 2 Beroun Central Bohemian Region Czechia 28 Klatovy Plzeň Region Czechia 3 Blansko South Moravian Region Czechia 29 Kolín Central Bohemian Region Czechia 4 Brno-město South Moravian Region Czechia 30 Kroměříž Zlín Region Czechia 5 Brno-venkov South Moravian Region Czechia 31 Kutná Hora Central Bohemian Region Czechia 6 Bruntál Moravian-Silesian Region Czechia 32 Liberec Liberec Region Czechia 7 Břeclav South Moravian Region Czechia 33 Litoměřice Ústí nad Labem Region Czechia 8 Česká Lípa Liberec Region Czechia 34 Louny Ústí nad Labem Region Czechia 9 České Budějovice South Bohemian Region Czechia 35 Mělník Central Bohemian Region Czechia 10 Český Krumlov South Bohemian Region Czechia 36 Mladá Boleslav Central Bohemian Region Czechia 11 Děčín Ústí nad Labem Region Czechia 37 Most Ústí nad Labem Region Czechia 12 Domažlice Plzeň Region Czechia 38 Náchod Hradec Králové Region Czechia 13 Frýdek-Místek Moravian-Silesian Region Czechia 39 Nový Jičín Moravian-Silesian Region Czechia 14 Havlíčkův Brod Vysočina Region Czechia 40 Nymburk Central Bohemian Region Czechia 15 Hodonín South Moravian Region Czechia 41 Olomouc Olomouc Region Czechia 16 Hradec Králové Hradec Králové Region Czechia 42 Opava Moravian-Silesian Region Czechia 17 Cheb Karlovy Vary Region Czechia 43 Ostrava-město Moravian-Silesian Region Czechia 18 Chomutov Ústí nad Labem Region Czechia 44 Pardubice Pardubice Region Czechia 19 Chrudim Pardubice Region Czechia 45 Pelhřimov Vysočina Region Czechia 20 Jablonec nad Nisou Liberec Region Czechia 46 Písek South Bohemian Region Czechia 21 Jeseník Olomouc Region Czechia 47 Plzeň-jih Plzeň Region Czechia 22 Jičín Hradec Králové Region Czechia 48 Plzeň-město Plzeň Region Czechia 23 Jihlava Vysočina Region Czechia 49 Plzeň-sever Plzeň Region Czechia 24 Jindřichův Hradec South Bohemian Region Czechia 50 Praha* Prague Czechia 25 Karlovy Vary Karlovy Vary Region Czechia 51 Praha-východ Central Bohemian Region Czechia Pro-Europeanness from a subnational perspective 191
Order Order District Region Country District Region Country Number Number 52 Praha-západ Central Bohemian Region Czechia 102 Liptovský Mikuláš Žilina Region Slovakia 53 Prachatice South Bohemian Region Czechia 103 Lučenec Banská Bystrica Region Slovakia 54 Prostějov Olomouc Region Czechia 104 Malacky Bratislava Region Slovakia 55 Přerov Olomouc Region Czechia 105 Martin Žilina Region Slovakia 56 Příbram Central Bohemian Region Czechia 106 Medzilaborce Prešov Region Slovakia 57 Rakovník Central Bohemian Region Czechia 107 Michalovce Košice Region Slovakia 58 Rokycany Plzeň Region Czechia 108 Myjava Trenčín Region Slovakia Rychnov nad 109 Námestovo Žilina Region Slovakia 59 Hradec Králové Region Czechia Kněžnou 110 Nitra Nitra Region Slovakia 60 Semily Liberec Region Czechia Nové Mesto 111 Trenčín Region Slovakia 61 Sokolov Karlovy Vary Region Czechia nad Váhom 62 Strakonice South Bohemian Region Czechia 112 Nové Zámky Nitra Region Slovakia 63 Svitavy Pardubice Region Czechia 113 Partizánske Trenčín Region Slovakia 64 Šumperk Olomouc Region Czechia 114 Pezinok Bratislava Region Slovakia 65 Tábor South Bohemian Region Czechia 115 Piešťany Trnava Region Slovakia 66 Tachov Plzeň Region Czechia 116 Poltár Banská Bystrica Region Slovakia 67 Teplice Ústí nad Labem Region Czechia 117 Poprad Prešov Region Slovakia 68 Trutnov Hradec Králové Region Czechia 118 Považská Bystrica Trenčín Region Slovakia 69 Třebíč Vysočina Region Czechia 119 Prešov Prešov Region Slovakia 70 Uherské Hradiště Zlín Region Czechia 120 Prievidza Trenčín Region Slovakia 71 Ústí nad Labem Ústí nad Labem Region Czechia 121 Púchov Trenčín Region Slovakia 72 Ústí nad Orlicí Pardubice Region Czechia 122 Revúca Banská Bystrica Region Slovakia 73 Vsetín Zlín Region Czechia 123 Rimavská Sobota Banská Bystrica Region Slovakia 74 Vyškov South Moravian Region Czechia 124 Rožňava Košice Region Slovakia 75 Zlín Zlín Region Czechia 125 Ružomberok Žilina Region Slovakia 76 Znojmo South Moravian Region Czechia 126 Sabinov Prešov Region Slovakia 77 Žďár nad Sázavou Vysočina Region Czechia 127 Senec Bratislava Region Slovakia Bánovce nad 128 Senica Trnava Region Slovakia 78 Trenčín Region Slovakia Bebravou 129 Skalica Trnava Region Slovakia 79 Banská Bystrica Banská Bystrica Region Slovakia 130 Snina Prešov Region Slovakia 80 Banská Štiavnica Banská Bystrica Region Slovakia 131 Sobrance Košice Region Slovakia 81 Bardejov Prešov Region Slovakia 132 Spišská Nová Ves Košice Region Slovakia 82 Bratislava* Bratislava Region Slovakia 133 Stará Ľubovňa Prešov Region Slovakia 83 Brezno Banská Bystrica Region Slovakia 134 Stropkov Prešov Region Slovakia 84 Bytča Žilina Region Slovakia 135 Svidník Prešov Region Slovakia 85 Čadca Žilina Region Slovakia 136 Šaľa Nitra Region Slovakia 86 Detva Banská Bystrica Region Slovakia 137 Topoľčany Nitra Region Slovakia 87 Dolný Kubín Žilina Region Slovakia 138 Trebišov Košice Region Slovakia 88 Dunajská Streda Trnava Region Slovakia 139 Trenčín Trenčín Region Slovakia 89 Galanta Trnava Region Slovakia 140 Trnava Trnava Region Slovakia 90 Gelnica Košice Region Slovakia 141 Turčianske Teplice Žilina Region Slovakia 91 Hlohovec Trnava Region Slovakia 142 Tvrdošín Žilina Region Slovakia 92 Humenné Prešov Region Slovakia 143 Veľký Krtíš Banská Bystrica Region Slovakia 93 Ilava Trenčín Region Slovakia 144 Vranov nad Topľou Prešov Region Slovakia 94 Kežmarok Prešov Region Slovakia 145 Zlaté Moravce Nitra Region Slovakia 95 Komárno Nitra Region Slovakia 146 Zvolen Banská Bystrica Region Slovakia 96 Košice – okolie Košice Region Slovakia 147 Žarnovica Banská Bystrica Region Slovakia 97 Košice* Košice Region Slovakia 148 Žiar nad Hronom Banská Bystrica Region Slovakia 98 Krupina Banská Bystrica Region Slovakia 149 Žilina Žilina Region Slovakia 99 Kysucké Nové Mesto Žilina Region Slovakia 100 Levice Nitra Region Slovakia Notes: * whole city as a one district for this purpose Source: Czech Statististical Office, Statistical Office of the Slovak Republic 101 Levoča Prešov Region Slovakia (2019). 192 Martin Plešivčák highest population density, while the Medzilaborce Highest values also applied to regions using EU agri- 2) has the lowest popula- cultural subsidies and promoting a policy of guaran- tion density. teeing the rights of ethnic minorities (especially the district (27 inhabitants/km 4.2 Results and Findings typical of an approach toward the values of Chris- Trnava Region and the Nitra Region in Slovakia), or 4.2.1 Regions - Pro-Europeanness Index esttian values Democracy were recorded(the Zlín Regionin the Czechin Czechia regions and with the - aPrešov peripheral Region position in Slovakia). geographically On the contrary, and socio-eco the low- nomically (Ústí nad Labem Region and Karlovy Vary At the regional level (Fig. 3), the highest pro-Europe Region), with increased support for the far-left or far- anness index values (ranging from 0.700 to 0.906) right parties and the Slovak region with a traditionally votesreflecting for Eurosceptic five variables parties in total in (votes the 2004–2019 for accession Euro in- peanthe 2003 Parliament European elections) Union weremembership recorded referendum, by Bratisla- Region). - egalitarian, etatist and nationalist electorate (Trenčín 4.2.2 Districts almostva followed achieved by other also four by the regions capital from region Slovakia of Czechia, (Trna va, Nitra, Košice and Prešov). A level of 0.700 was pro-Europeanness index at the district level. Bystrica and Žilina) reached the values of index from We came up with other interesting findings for the Prague. Other two regions from Slovakia (Banská Pro-Europeanness Index obtained by two regions from Czechia (Zlín and South - 0.500 to 0.600. The values between 0.400–0.500 were ences (Fig. 4). In the districts with the highest values contrary, the lowest values of the index (0.029–0.294) ofThe the resulting index, the index Slovakian shows considerable districts absolutely spatial differdomi- wereMoravia) reported and the by thelast Ústí Slovak nad region,Labem Region,Trenčín. the On Kar the- lovy Vary Region and the Moravian-Silesian Region valuesnate, with of index the first are being located Dunajská exclusively Streda in (0.994) the south and- Czechia. second Komárno (0.903). Other districts with high in Czechia.Considering Of theregional last twelveperspective, regions, one all of arethe highfrom- est spatial concentrations were recorded also in inwest the of agro-sector Slovakia. These and adistricts, strong Hungarianwith geographical minor- regions of Prague, Bratislava and Košice, i.e. centres ity.proximity In the tocase the of capital, Czechia, enjoy the ahighest significant values position were that have been more successful in the post-socialist transformation, with higher economic performance, and Praha-východ) and the city district of the second - reached by districts of Prague (Praha-západ, Praha ulation with higher education and socio-economic the contrary, districts with the lowest values (below status,localization supporting, of significant in general, foreign liberal investment, political parties. a pop 0.200)largest arecity located of the country,in west Bohemia Brno (Brno-město). in the Ústí nad On
1.0
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0.0
Zlín Nitra Žilina Plzeň Trnava Košice Prešov Prague Trenčín Liberec Bratislava Pardubice Vysočina Olomouc Karlovy Vary Hradec Králové Banská Bystrica South Moravian South Bohemian Ústí nad Labem Central Bohemian Moravian-Silesian
Fig. 3 NUTS 3 regions of Czechia and Slovakia by Index of Pro-Europeaness. Note: Light greyindicates the Czech regions, dark gray indicates the Slovak regions. Source: Czech Statistical Office, Statistical Office of the Slovak Republic (2019), author’s research. Pro-Europeanness from a subnational perspective 193
Fig. 4 Districts of Czechia and Slovakia by Index of Pro-Europeaness. Source: Czech Statistical Office, Statistical Office of the Slovak Republic (2018), author’s research.
Labem Region (Most 0.144, Louny 0.156, Chomutov “districts of significantly pro-European orientation” is comprised of 19 districts, 0.161) and the Karlovy Vary Region (Sokolov 0.190). The second cluster, named 0.164 and Teplice 0.171), the Plzeň Region (Tachov concentration of these districts can be found in the areaof which located all are in thelocated southwest in Slovakia corner again. of the The country largest This group of districts is completed by couple from (Bratislava, Šaľa, groupMoravia of –districts Bruntál is (0.151, characterized Moravian-Silesian by a relatively Region) high the Prešov Region in the northeast (a compact belt unemploymentand Znojmo (0.168, rate andSouth relatively Moravian low Region). wages, withThis Galanta, Senec, Nové Zámky) and in support for left-wing, far-left or far-right parties (Hav- Ľubovňa, Bardejov, Sabinov and Prešov) continuing - toof thethe sixKošice districts Region – Poprad, in the Kežmarok,east of Slovakia Levoča, (Košice Stará kia, the districts with the lowest values are located in - thelík and northwestern Voda 2016; corner Maškarinec of the 2017, country 2019). in the In SlovaŽilina est towns of Slovakia (Bratislava, Košice and Prešov). Čadca Comparedand Trebišov). to the This national cluster average also includes (Fig. 6), thisthree catego larg- of long-term support for the values of egalitarianism, ry declared a strong support for accession in the EU Region (Kysucké Nové Mesto, and Bytča) typical pp) and lower support for Euro- 2012). pp). etatism and nationalism (Plešivčák 2011; Madleňák referendum (+12.71 “dis- 4.2.3 Clusters trictssceptic of parties mildly inpro-European the EU elections orientation” (−7.99 . It consists - The third category of districts is described as gories based on the values of the pro-Europeanness - Byindex using (Fig. the 5). cluster method, we derived five cate of 30 spatial units, the larger part of which (23) is “districts of markedly thelocated southwest in Slovakia. of Slovakia The main with concentration seven districts is recog over- pro-European orientation”, consists of only two spatial nized in the regions of Bratislava, Trnava and Nitra in The first category, named - - ing to this cluster are located in southwest Slovakia all (Pezinok, Malacky, Trnava, Senica, Nitra, Piešťany units, both located in Slovakia. The districts belong ofand the Hlohovec). Žilina The second compact area can be iden position of agriculture and a strong Hungarian minor- tified in the north of the country in the eastern part- (Dunajskaity. Compared Streda to the and national Komárno), mean with(Fig. 6),a traditional this group tration is comprehensively Region (Dolný Kubín, complemented Tvrdošín, byMartin, four recorded strong support for country’s accession to Liptovský Mikuláš and Ružomberok). This concen the EU (+14.06 pp) and very low support for Euros- pp). adjacent districts, from the Banská Bystrica Region (Banská Bystrica, Zvolen, Veľký Krtíš and Lučenec). ceptic parties (−16.92 The third concentration of districts of this type is 194 Martin Plešivčák
Fig. 5 Categories of districts in Czechia and Slovakia clustered by Index of Pro-Europeaness. Source: Czech Statistical Office, Statistical Office of the Slovak Republic (2019), author’s research.
located in the eastern part of Slovakia in the regions - ed in the Ústí category, with 22 (88%) from Czechia. Most are locat of Banská Bystrica (Revúca), Košice (Rožňava, Spišská nad Labem Region (7 of 7), theÚstí Plzeň nad ofNová districts Ves and of mildlyMihcalovce) pro-European and Prešov orientation (Snina, Vranov can be LabemRegion Region (4 of 7)and and the theMoravian-Silesian Moravian-Silesian Region, Region have nad Topľou and Svidník). In Czechia, the only cluster certain(3 of 7). socio-economic Most of them, problems mainly from (relatively the high unemployment, and low wages), favouring left-wing, found in the capital region (districts of Prague). One- far-left and recently protest parties to some extent. district is located in the east of the Vysočina Region (Žďár nad Sázavou), three in Moravia in South Moravi northwest of Czechia, encompassing the regions of an Region (Brno-město) and Zlín Region (Zlín, Vsetín). The most visible concentrationÚstí is locatedÚstí in nad the (Fig.Two largest6) slightly cities above-average of Czechia (Prague support and for Brno)EU acces fell- - sioninto (+2.44 this category. pp) and lower This groupsupport is for characterized the Eurosceptic by tov),Liberec Central (Česká Bohemia Lípa), (Rakovník nad Labem and (Děčín, Kladno), Kar- pp). Labem, Teplice, Litoměřice, Most, Louny and Chomu numerousparties in EU of electionsall clusters, (−4.06 we labelled “transitional” regionslovy Vary of (Sokolov Moravia and and Cheb) Silesia, and there Plzeň are (Domažlice, six such districtsThe fourth group of districts, which is the most districts,Tachov, Plzeň-sever namely Znojmo, and Rokycany). Vyškov (South In the Moravianhistorical
. Cluster analysis marked 73 spatial units- - (almost a half of all districts), of which 48 (66%) are vakia,Region), districts Přerov of (Olomouc this type Region),are located Jeseník, exclusively Bruntál in easternlocated inBohemia. Czechia. In The Slovakia, largest the concentration most compact of areadis theand northwestKarviná (Moravian-Silesian of the territory inRegion). the Žilina Within Region Slo tricts of this type can be identified in the central and- Čadca, - isized situated by (Fig. in 6) all slightly territory below of Trenčín average Region support and for adja EU for(Bytča, its traditional Kysucké support Nové of nationalist Mesto), whereparties, theyand cent districts. Thispp category) and mildly of districts higher is support character for theformed values a compact of etatism concentration. and egalitarianism. This area Compared is known to Eurosceptic parties in EU elections (+2.15 pp). the national average (Fig. 6), this cluster of districts is accession (−3.51 characterized by a markedly below average approv- by the notable degree of Euroscepticism (“districts pp) and an outstanding withThe mild last or group significant consists elements of districts of Euroscepticism” characterized). support for the Eurosceptic parties in EU elections (+8.25al for EU pp ).accession (−6.44
Of the total number of 149 districts, 25 fall into this Pro-Europeanness from a subnational perspective 195
20.0
15.0
10.0
5.0
“YES” in EU Referendum 2003 0.0 Districts of markedly Districts of significantly Districts of mildly “Transitional” Districts with mild pro-European pro-European pro-European Districts or significant elements Support for Eurosceptic parties orientation orientation orientation of Euro-scepticism (in EU Elections 2004–2019) –5.0
–10.0
–15.0
–20.0
Fig. 6 Categories of districts in Czechia and Slovakia clustered by Index of Pro-Europeaness – selected EU electoral characteristics. Note: Difference from average value for Czechia and Slovakia, difference is measured by percentage points. Source: Czech Statistical Office, Statistical Office of the Slovak Republic (2019), author’s research.
5. Discussion the left (Bale 2010). Nevertheless, in Czechia, the
Some regions of Czechia and Slovakia are character- in the case of elections to the European Parliament, ized by increased turnout and low level of Euroscep- thoughfar-left KSČMcurrently has of been declining enjoying trend. significant Economic support prob- ticism in the case of European issues (referendum, lems and migratory pressures are causing uncertain- elections), while others are more passive in voting ty in Europe and create a breeding ground for popu- lists (Charvát - recently addressed the causes of support for far-right peanor support parliamentary Eurosceptic elections parties toinclude a great disappoint extent. The- 2007). Quite a number of authors have mentmost frequentwith politics, reasons non-alignment for not participating with the electoralin Euro program of any of the parties, and a lack of political increaseparties in in Czechia Euroscepticism and Slovakia can (e.g. be linkedKluknavská to ongoing 2012, globalization2013; Gregor (Salo2015; 2014). Mikuš It et has al. created2016). Ina groupgeneral, of interest in voting in European elections as well as sup- - portcommitment for Eurosceptics, as such (Greffet it is also 2007). important In the to debate mention on nerable to the current liberalization, worrying about the voter’s relationship to the idea of Euro-citizenship their“bereaved” economic who future are losing and growing certainty cultural and feeling diversity. vul (Frognier 2000), the project of European integration, Eurosceptics does not necessarily have to be a result Europeanism, the degree of awareness related to the of populism, but also as a legitimate part of the politi- (Blon- - zation (Salo 2014). It is said about the new European role of the EU and the benefits of membership cleavagecal arena ofthat social represents polarization those based “bereaved” on the by existence globali fromdel et EUal. 1997)membership, and the visibilitye.g. through of the structural EU in ordinary funds, of various life and material opportunities (in our havelife (Irwin a higher 1995). tendency Voters to ofparticipate regions benefiting in the European more research, especially socio-economic status, civil and Parliament elections and support Euro-optimist par- minority rights) perceived by different actors with different interests depending on the process of terri- voter (non)participation were investigated by Linek ties (Jesuit 2003). In the case of Czechia, reasons for It turns out that the current wave of Euroscepti- torial integration (Bartolini 2007). cism(2013), or populismin Slovakia is bybetter Gyárfášová understood (2019). by the far-right 6. Conclusions - kia by ĽSNS). - erssubjects who (inpreviously Czechia inalmost particular unreservedly by SPD and supported in Slova They are able to attract manual work Considering the NUTS 3 regions of both countries (14 in Czechia and 8 in Slovakia) in terms of the final index 196 Martin Plešivčák of pro-Europeanness and regional cleavage, the high- est values were achieved by the regions of the largest found mainly in the Ústí nad Labem Region, the Kar- cities (Prague in the Czechia, Bratislava and Košice in lovyare from Vary the Region Czechia. and Thethe highestMoravian-Silesian concentration Region, was Slovakia) typical of a more educated, urban popula- - tion with a higher socio-economic status, more eco- - lowi.e. in wages), regions with with increasing relatively support significant for radical socio-eco (left- cialist transformation, with a higher concentration of wingnomic or difficulties right-wing) (relatively and protest high (anti-system) unemployment parties. and large,nomically especially efficient foreign and moreinvestment, successful and ain predomi post-so- Within Slovakia, districts of this type are located in the northwest of the territory in the Žilina Region second case of a pronounced tendency towards the Čadca, - EUnant project right-wing is represented (or central/liberal) by regions usingelectorate. EU subsi The- ally support nationalist parties and espouse values of dies under its largest agricultural policy, also depend- etatism(Bytča, and egalitarianism. Kysucké Nové Mesto), which tradition ent on the supranational policy securing the rights of In general, ideas of European integration and Euro-optimism as such in Czechia and Slovakia are the Nitra Region in Slovakia) or known for values of more common among the urban electorates (Prague, ethnic minorities (especially the Trnava Region and Brno, Bratislava, Košice), areas with a higher con- - lowestChristian values Democracy of the (thepro-Europeanness Zlín Region in Czechiaindex were and cant agricultural production (southwest of Slovakia), registeredthe Prešov inRegion the socio-economically in Slovakia). On the and contrary, geograph the- whichcentration in this of acase particular is probably ethnic related group, to with the signifi status ically peripheral regions of Czechia, with increased of the EU as a guarantor (higher instance for pro- support for far-left, far-right and protest parties (the tection) of civil and minority rights, and a provider Ústí nad Labem Region, the Karlovy Vary Region and of agro-subsidies. In the case of Euroscepticism, the the Moravian-Silesian Region) and parts of Slovakia Czech districts and regions prevail, especially from with a population traditionally close to the values the peripheral northwest and the Moravian-Silesian border areas. In this context, relatively important Region). socio-economic problems (in comparison with the of egalitarian, etatism and nationalism (the Trenčín national average high unemployment and low wages) theAt top the two district categories level (149(“districts districts of markedly in total, 77 pro-Eu from- the EU (membership) for them (or at least in the ropeanCzechia orientation”and 72 from andSlovakia), (“districts it was of found significantly that in formcan be of mentioned. a penalty for Part their of thefailure electorate to solve can them) “blame” and pro-European orientation”), with the highest values thus, on a practical level, can prefer populist (radical, of the index of pro-Europeanness, the Slovak districts anti-system) and Eurosceptic parties. In the case of (mainly from the southwest) dominated over Czech Slovakia, regions located to the northwest without - category (“districts with mild or significant elements ditional vote for (ultra)nationalists (ĽSNS, formerly ofdistricts. Euroscepticism” On the contrary,) came from 88% Czechia, districts with of the the abso last- votingsignificant for SNS), socio-economic egalitarian problems, and etatist-minded but with the politi tra- lute lowest values being registered for districts from Ústí nad Labem Region (northwestern Bohemia). - cal movements (SMER-SD, formerly voting for HZDS), gories of districts across the countries based on the ofwere previous shown studies as least from pro-European the Czech (Pink oriented. 2012; Voda The By using the cluster method, we derived five cate 2015;obtained Kostelecký results are et al.in 2016;accordance Koubek with 2019) the andfindings Slo- “districts of mark- edly,values significantly of the pro-Europeanness or mildly pro-European index. The orientation” first three) “evidently pro-European” groups ( socio-economicvak literature (e.g. causes Krivý affecting et al. 1996; the spatialPlešivčák distribu 2011;- - tionMadleňák of election 2012; results. Przybyla 2019) on the historical and terizedconsists by of a 51 clear spatial support units, for 86% accession of them to locatedthe Euro in- - peanSlovakia. Union The and districts low support of these for clusters Eurosceptic are characparties in European Parliament elections. Several districts of This study provides new insights into the “geog ofraphy regions of pro-Europeanness” and districts of two over countries a relatively that in long the in Czechia, and Bratislava, Košice and Prešov in Slo- pastperiod formed of time a (2003–2019),single state, applying at the sub-national a methodology level vakia),this type while are othersurban are(districts located of in Prague, southwest Brno-město Slovakia not used before to assess the territorial context of EU - - tion for policy- and decision-makers on the regions (Senec,duction Dunajskaand a strong Streda, concentration Komárno and of the Nové Hungarian Zámky) inintegration which EU support. assistance This should work alsobe targeted provides to informa sustain where there is relatively significant agricultural pro (restore) the meaningfulness of both the idea of Euro- pean integration and EU membership in those parts of minority. On the opposite side, there are districts with the countries that are currently most critical to the EU supportmild or significant for Eurosceptic elements parties of Euroscepticism, when EU elections with project, and thus to stop encouraging Euroscepticism significantly lower support for EU accession relevant across the EU, starting with its partial regions. take place. Of the 25 spatial units in this category, 22 Pro-Europeanness from a subnational perspective 197
Acknowledgements Goodwin, M. (2011): Right Response: Understanding and Countering Populist Extremism in Europe. London, Chatham House. - This paper was prepared with the support provided Greffet,Macmillan, F. (2007): Basingstoke, Abstention. 1–6. In: Encyclopaedia of by research grants of Slovak Research and Develop European Election. Y. Déloye, M. Bruter. Palgrave ment Agency APVV-17-0079 “Population Analysis and Forecast of the Slovak Republic in Time-horizon 2080: Gregor, K. (2015): Kto sú Kotlebovi voliči? Preskupovanie Identification and Modelling the Impacts on Society in G voličov v, Banskobystrickom kraji v rokoch 2009 až 2014. Different Spatial Scales”. Sociológia – Slovak Sociological Review 47(3), 235–252. References yárfášová O. (2007): Jak volili Slováci ve volbách do EP v roce 2004 v kontextu veřejného vnímáníústav Evropské AV ČR, unie. In: Volby do Evropského parlamentu 2004. L. Linek, J. Outlý, G. Tóka, A. Batory. Sociologický Bale,from T. (2010): the Populist If You Radical Can’t Beat Right Them in Western Join Them? Europe. Praha, 279–292. Explaining Social Democratic Responses to the Challenge Gyárfášová,K. Falasca, O. M. (2019): Grusell Public’s and L. Nord.perception Mittuniversitetet, of the EU and Sundsvall,turnout in 56–56.the EP election. In: Euroflection. Eds. N. Bolin, Political Studies 58(3), 410–426, https://doi.org /10.1111/j.1467-9248.2009.00783.x. širšie Bartolini, S. (2007): Cleavages. In: Encyclopaedia of Gyárfášová, O. Velšic. M.členstva (2004): Volebné správanie European Elections. Ed. Y.Šaradín, Déloye, M. Bruter. Palgrave slovenských voličov v prvých eurovoľbách a Macmillan, Basingstoke, 43–51. súvislosti vnímania v EÚ.. Bratislava, Inštitút pre Baun,and M., the Durr, EU Referendum. J., Marek, D., Journal ofP. Common(2006): The Market verejné otázky. StudiesEuropeanization 44(2), 249–280, of Czech https://doi.org/10.1111 Politics: The Political Parties Hastie,Prediction, T., Tibshirani, Second R. Edition Friedman, (Springer J (2016): Series The in ElementsStatistics). /j.1468-5965.2006.00622.x. Berlin,of Statistical Springer. Learning: Data Mining, Inference, and Blondel, J., Sinnott, R., Svensson, P. and Voter participation. European Journal for Political (1997): Representation Havlík, V. (2008): Jak je měřit? Přístupy a metody analýzy Havlík,postojů V. (2010): politických České stran k evropské integraci. Central Research 32(2), 243–272, https://doi.org/10.1023 European Political Studies Review 10(4), 340–369. Sage./A:1006832419724. Brno, Masarykova univerzita. politické strany a evropská Blumer, J. G. Šaradín,(1983): Communicating to Voters. London, integrace: Evropeizace, evropanství, euroskepticismus?České České republice. řádu. Bradová,Ed. P. Šaradín. E., P. (2004): Volební kampaň. In: Volby Havlík, V., Hoskovec, L. (2009): Krajské volby v Budge,do Evropského I. (2001): Mapping parlamentu Policy v Preferences: Estimates for zastupitelstevrepublice v kontextu v letech konceptu 2000–2008. voleb European druhého Electoral Periplum, Olomouc, 179–207. Analýza vybraných aspektů voleb do krajských
Parties, Electors and Governments (1945–1998). Oxford, Studies 4(1), 22–47. Oxford University Press. Havlík,politologický V., Kaniok, ústav P. (2006): Masarykovy Euroskepticismus univerzity. a země Cabada, L. (2010): Volby do Evropského parlamentu jako střední a východní Evropy. Brno, Mezinárodní volby druhého řádu: reflexe voleb v nových členských and Re-Alignment of Party Politics in the Czech Republic. zemích EU ze středovýchodní. (2005): Evropy. Explaining Working vote Papers Havlík,Acta V.,Politologica Voda, P. (2016): 8(2), 119–144. The Rise of New Political Parties Fakulty mezinárodních vztahů 4(11), 5–26. Carrubba,elections: C., evidenceTimpone. from R Europe. Comparative Political switching across first-order and second order Havlík,of the V., Czech Vykoupilová, Republic. H. Communist (2008): Two and dimensions Post-Communist of the Europeanization of election programs: The case Studies 38(3), 260–281, https://doi.org/10.1177 /0010414004272693. Studies 41(2), 163–187, https://doi.org/10.1016 Czech Statistical Office. (2019): Elections Database /j.postcomstud.2008.03.007. C(2003–2019).. (2011): Individual Prague, and Czech contextual Statistical variation Office. in EU Henderson,Euroscepticism, K. (2008): Vol. The 1 Case Slovak Studies Republic: and Country Eurosceptics Surveys. De Vries, C., van der Brug, W., van Egmond, information. M., van Electoral der Eijk, and Phoney Europhiles. In: The Party Politics of
/j.electstud.2010.09.022.issue voting: The role of political Hix,Eds. S., MarshA. Szczerbiak, M. and P. Taggart. Oxford University Studies 30(1), 16–28, https://doi.org/10.1016 UnderstandingPress, Oxford, 277–294. European Parliament Elections. (2007): Punishment or Protest? Fiala, P., Hloušek, V., Pitrová, M., Pšeja, P. (2006): časopis Evropeizace politických stran a zájmových skupin: Hix,The S., MarshJournal, M. of (2011):Politics 69(2),Second-order 495–510, effects https://doi. plus pan- základní problémy a směry analýzy. Politologickýélectorale: Europeanorg//10.1111/j.1468-2508.2007.00546.x. political swings: An analysis of European 13(1), 3–26. élections Frodnier, A.-P., (2000): Identité et participationélections pour une approche européenne des Hloušek,Parliament V., Pšeja, elections P. (2009): across Europeanization time. Electoral ofStudies Political 30(1), européennes.P. Perrineau, C. In: Ysmal, Le vote C. Pressesdes Quinye. de Sciences, Les Paris, Parties4–15, https://doi.org/10.1016/j.electstud.2010.09.017. and the Party System in the Czech Republic. européennes du 13 juin 1999. Ed. G. Grunberg,
75–94. Journal of Communist Studies and Transition Politics 198 Martin Plešivčák
25(4), 513–539, https://doi.org/10.1080 Krivý, V., Feglová, V., Balko, D. (1996): Slovensko a jeho 2009:/13523270903310902. Prostor pro evropeizaci politických stran ve Ladrech,regióny: R. (2002):sociokultúrne Europeanization súvislosti volebného and Political správania. Parties: Hricová, H. (2009): Česká republika. In: Eurovolby Bratislava, Nadácia Médiá.
středovýchdoní Evropě. L. Cabada, V. Hloušek. Aleš 4002.Towards a Framework for Analysis. Party Politics 8(4), Making:Čeněk, Plzeň, Methods 41–66. and Applications. New York, Springer- 389–403, https://doi.org/10.1177/135406880200800 Hwang, C. L., Yoon, K. (1981): Multiple Attribute Decision volební účasti v České republice v letech 1990–2010. Šárovec Linek,Brno, L. Centrum(2013): Kam pro sestudium ztratili demokracie voliči? Vysvětlení a kultury. vývoje Verlag, https://doi.org/10.1007/978-3-642-48318-9. , A. Hynčica, P., , D. (2018): Slovenské politické strany ústav optikou konceptu novosti. Acta Fakulty filozofické Linek,AV ČR.L., Outlý, J., Tóka, G., Batory (2007): Volby do Západočeské univerzity v Plzni 10(2–3), 7–34, Evropského parlamentu 2004. Praha, Sociologický Paris,https://doi.org/10.24132/actaff.2018.10.2-3.1 Publications de la Sorbonne. Alignments: Cross–national Perspectives. London, Charlot, M. (1986): The June 1984 European Elections. Lipset,Collier-Macmillan. S. M., Rokkan, S. (1967): Party Systems and Voter
Charvát, J. (2007): Současný politický extremismus Parliament. Basingstoke, Palgrave Macmillan. campaigna radikalismus. for the Praha, elections Portál. for the European Parliament Lodge, J. (2010): The 2009 elections to the European Irwin G. (1995): Second-order or third-rate?: Issues in the Manokaran,Madleňák, T. E., (2012): Subhashini, Regionálna S., Senthilvel, diferenciácia S., volebného 1994. Electoral Studies 14(2), 183–199, https://doi.org Muruganandham,správania na Slovensku R., Ravichandran, (1998–2010). K. Bratislava,(2011): Veda. Electoral/10.1016/0261-3794(95)96843-7. Politics: Participation in National and European Application of multi criteria decision making tools Jesuit,Contests D. (2003): in the The 1990s. Regional European Dynamics Union of Politics European 4(2),
2001. and validation with optimization techniquecase Kaniok,139–164, P., Havlík, https://doi.org/10.1177/146511650300400 V. (2016): Populism and euroscepticism 112–115.study using TOPSIS, ANN & SAW. International Journal of Management & Business Studies 1(3),
in the Czech Republic: Meeting friends or passing by? Maškarinec, P. (2017): A Spatial Analysis of Czech Romanian Journal of European Affairs 16(2), 20–35. Parliamentary Elections, 2006–2013. Europe-Asia Klíma, M., Outlý, J. (2010):Čeněk. Volby do Evropského parlamentu Studies 69(3), 426–457, https://doi.org/10.1080 2009. V vybrané otázky přípravy, průběhu a výsledků parties/09668136.2017.1313962. in Czech parliamentary elections between voleb. Plzeň, Aleš Maškarinec, P. (2019): The rise of new populist political Kluknavská,časopis pro A. výzkum (2012): radikalismu,Krajne pravicové extremismu strany v a terorismu parlamentných voľbách 2012 naSlovensku. Rexter – 2010 and 2017: the geography of party replacement. Štúra Euroasian Geography and Economics 60(5), 511–547, 10(1), 1–35. https://doi.org/10.1080/15387216.2019.1691928. Kluknavská, A. (2013): Od k parazitom: Tematická Mikuš, R., Gurňák, D. (2016): Rómska otázka ako jeden adaptácia krajnej pravice v parlamentných voľbách Cassoviensisz mobilizačných 10(1), faktorov 29–46. volebnej podpory krajnej na Slovensku. Politologický časopis 20(3), 258–281, pravice na Slovensku a v Maďarsku., A. (2016): Geographia Analýza Euroscepticism:https://doi.org/10.5817/PC2013-3-258. Party Positions on European Integration Kopecký, P., Mudde, C. (2002): The Two Sides of Mikuš, R., Gurňák, D., Máriássyová volebnej podpory Mariána Kotlebu ako reprezentanta in East Central Europe. European Union Politics 3(3), krajnej pravice v krajských voľbách 2013. Sociológia 297–326, https://doi.org/10.1177/146511650200Čermák, 48(1), 48–70. Bernard,3003002. J., Šimon, Opricovic, S., Tzeng, G. H. (2004): Compromise solution Kostelecký, T., Mikešová, R., Poláková, M., D., Researchby MCDM 156(2), methods: 445–455, a comparative https://doi.org/10.1016 analysis of VIKOR Česku M. (2016): Geografie výsledků ústav and TOPSIS. European Journal of Operational AVparlamentních ČR. voleb: prostorové vzorce volebního Pink, M. (2012): Volební mapy České chování v 1992–2013. Praha,českých Sociologický /S0377-2217(03)00020-1. pro studium demokracie a kultury. a Slovenské republiky Koubek, J. 2019. Volební chování velkoměst v po roce 1993: vzorce, trendy, proměny. Brno, Centrum historické perspektivě se zaměřením na výsledky levice. ElectionsPraha, Masarykova in the Czech demokratická Republic: A akademie.Light at the End of the Plešivčák, M. (2011). Regionálny obraz korelačnej závislosti Kovář, J. (2013): The Supply Side of Second-Order 5(1),medzi 2–16. volebnými preferenciami a vybranými sociálnymi charakteristikami populácie Slovenska. Regionální studia Tunnel? Romanian Journal of European Affairs 13(1), 57–82. Plešivčák, M. (2015): First Two Elections to the European Kovář, J. (2014): Europeizace volebních programů pro volby Parliament in Slovakia: Turnout and Results in Regional do EP. Politické vedy 17(3), 32–67. Perspective. Evropská volební studia 10(1), 14–32. Kovář,Slovakia J., Kovář, Compared. K. (2014): Perspectives Representation on European of Women Politics Przybyla, V. (2019): Vliv vybraných demografickýchčasopis in Second-order Elections: The Czech Republic and charakteristik na volební chování na Slovensku v meziválečném období. Geografický 71(2), and Society 15(1), 1–18, https://doi.org/10.1080/1570 181–199, https://doi.org/10.31577/geogrcas.2019 5854.2013.793533. .71.2.10. Pro-Europeanness from a subnational perspective 199
Reif, K., Schmitt, H. (1980): Nine second-order national elections: a conceptual framework for the analysis of European election results. European Journal of Political Computer Modelling 45(7–8), 801–813, https://doi.org /10.1016/j.mcm.2006.03.023. Statistical Office of the Slovak Republic (2019): Elections RobertResearch Schuman 8(1), Foundation 3–44, https://doi.org/10.1111 (2004): Survey on Evaluation and Referenda Database (2003–2019). Bratislava, of/j.1475-6765.1980.tb00737.x. competences at the European and national levels Statistical Office of the Slovak Republic. Šaradín, P. (2008): Teorie voleb druhého řádu a možnosti Salo, S. (2014): Eurosceptics in the 2014 EP Elections. jejich aplikace v České republice. Olomouc, Univerzita Protest(in EU member Parties Mobilizedstates), 28–30 on National April 2004. Cleveages Between EuroscepticismPalackého v Olomouci. in EU Member and Candidate States. Globalization Winners and Losers. Helsinki, Finnish Taggart,Brighton, P., Szczerbiak, Sussex European A. (2002): Institute. The Party Politics of Institute of International Affairs. České Voda, P. (2015): Jaká je role postkomunismu? Volební Shih, H. S., Shyurb, H. J., Lee, E. S. (2007): An extension of ageografie kultury. a Rakouské republiky v letech TOPSIS for group decision making. Mathematical and 1990–2013. Brno, Centrum pro studium demokracie 200 Original Article Which subjects contribute to the teaching of cross-curricular topic Environmental Education at elementary schools in Czechia? Tomáš Matějček1,2,3,*, Jan Bartoš1, Silvie R. Kučerová3
1 Charles University, Faculty of Science, Department of Social Geography and Regional Development, Czechia 2 Charles University, Faculty of Science, Department of Phyisical Geography and Geoecology, Czechia 3 J. E. Purkyně University in Ústí nad Labem, Faculty of Science, Department of Geography, Czechia * Corresponding author: [email protected]
ABSTRACT This article summarizes the results of the research focused on the realization of the cross-curricular subject Environmental Educa- tion (CCSEE) at elementary schools (pupils’ age 6–15 years) in Czechia. The introduction of cross-curricular subjects into the Czech educational system is linked to curricular reform and it has been implemented in Czech schools since 2007. CCSEE is one of the six currently implemented cross-curricular topics. The main objective of the present study is to determine which school subjects are involved in its implementation. The study was conducted through an internet questionnaire and responses were received from 640 schools. Data were processed by basic statistical methods. A school typology depending on the subjects involved in implementing EE was developed with the help of cluster analysis. The research shows that EE is implemented through most subjects, but their representation varies considerably for individual schools.
KEYWORDS curricula; elementary school; geographical education; environmental education; cross-curricular subject; Czechia
Received: 11 July 2019 Accepted: 15 July 2020 Published online: 21 October 2020
Matějček, T., Bartoš, J., Kučerová, S. R. (2020): Which subjects contribute to the teaching of cross-curricular topic Environmental Education at elementary schools in Czechia?. AUC Geographica 55(2), 200–209 https://doi.org/10.14712/23361980.2020.14 © 2020 The Authors. This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0). Teaching of cross-curricular topic Environmental Education at elementary schools 201
1. Introduction (including CCSEE) were then elaborated in detail later (see Pastorová et al. 2011; Činčera 2011). One of the changes brought to the Czech schools In practice, CCSEE can be implemented in schools by Framework Educational Programs (FEPs) was either by integrating it into the educational content of the introduction of cross-curricular subjects. The existing school subjects, by creating a separate school cross-curricular subject Environmental Education subject or through project teaching (cross-curricu- (CCSEE) is one of them. The aim of the research pre- lar projects, project days). The individual forms may be combined with each other. The same situation is involved in the implementation of CCSEE in schools in Slovakia (Kelcová 2009). The opposite example andsented to whathere extent,was to findand outwhether which it school is possible subjects to trace are certain typical “models” of CCSEE implementation in Environemtal Curriculum gives examples of the imple- terms of the involvement of individual school sub- mentationis Great Britain of EE in (specifically individual subjects England), (Green where 2018). The jects. According to our professional orientation, we Thus, all teachers in the school may theoretically were especially interested in the role of geography in participate in the implementation, but the situation in the individual schools may be quite different and the The development of EE in Czechia is clearly sum- actual situation has not yet been closely monitored. fulfillingmarized bythe Máchal objectives (2000) of EE. and Činčera (2013a, 2014). This paper summarizes the results of research con- Putting into the international context is discussed in ducted at Czech elementary schools (primary schools) more detail in Činčera (2013b). The roots of EE can and focused on determining the degree of involve- be found in the activities of volunteer organizations already in the interwar period of Czechoslovakia. of the EE, respectively. implementation of CCSEE. These activities were followed by organizations work- mentIn theof individual past, this issuesubjects was in only fulfilling partially the monitored, objectives ing with children and youth in the 1970s (especial- ly the Czech Union for Nature Conservation and the (2009). A wide range of aspects of EE goals imple- - mentationas part of theat schools analysis was carried also addressed out by Daňková by a detailed et al. ronmental education (Máchal 2000) were later estab- study of Činčera et al. (2016), however, the rate of lishedBrontosaurus from these Movement). and other The organizations. first centers of envi involvement of individual subjects was not real- However, the implementation of EE into formal ized in mentioned study. Finding out which subjects (school) education was delayed in comparation with are involved in the implementation of EE can show some countries of northern and western Europe. whether EE really has a cross-curricular character At present, the basic framework for EE at schools is based on the State Program of EE in the Czech Repub- EE can be realized in various forms and in vari- lic (approved in 2000). When planning EE implemen- ousand schoolthus fulfills subjects. its potential. The potential for implementing tation goals, schools are based on valid curricular documents (especially from FEP or school education- al programs – see below), from the Methodological butCCSEE in different objectives, degrees. introduced The interdisciplinary by FEPs (see conceptJeřábek Guideline of the Ministry of Education, Youth and ofand EE Tupý prevails 2007), (Aikens have practically et al. 2016). all Numberschool subjects, of con- Sports on ensuring EE, from the regional concept of crete examples of linking environmental issues with EE and their action programs, and from the analysis of other subjects, including fewer common ones, was described in the literature (see below). designated at individual schools and is responsible for In general, Godemann (2008) deals with integra- specific school conditions. The CCSEE coordinator is tion and transdisciplinary concepts of environmental authors mention that its evaluation is very important issues in a comprehensive way, which also summariz- forthe fulfillmentthe development of EE objectives of EE (Verma according and toDhull FEP. 2017; More es the main principles of working with information in such a teaching approach. The importance of an 2004; Nam 1995). interdisciplinary approach to EE is also highlighted SsoziCCSEE 2012; got Fergusoninto Czech 2008; curriculum Grodzińska-Jurczak documents in 2007, in connection with the implementation of FEP The traditional is linking of EE with nature science - education,by Jančaříková especially (2009). with biology, physical geogra- ulum previously used. CCSEE is one of six currently phy and chemistry (Mwendwa 2017; Florentina and (Jeřábekimplemented and cross-curricularTupý 2007), which subjects. replaced However, the curric in the context of actual FEP revision, the future of cross-cur- vá 2009). Aikens and McKenzie (2016) also state that ricular topics is uncertain and actualy discussed. mostBarbu of 2015; the topics Ryplová used and in Reháková environmental 2011; educationŘezníčko According to Činčera (2005), the inclusion of CCSEE belong to the natural sciences, but recent studies as a cross-curricular subject represents a major shift are beginning to address the social sciences. Educa- in its understanding in Czechia. Thus, EE started to tion in these school subjects enables pupils to know be understood as a real cross-subject issue that inte- principles of natural processes, introduces them to grates both the natural and human sciences. The the diversity of nature, and to understanding of the expected outcomes for cross-curricular subjects human activity impact on the natural environment. 202 Tomáš Matějček, Jan Bartoš, Silvie R. Kučerová
Yet some dichotomy can be observed in this respect The questionnaire was created according to the prin- as well. While nature science education was mainly ciples for quantitative research (Gavora 2010; Chrás- ka 2007). andmotivated strengthen by the competitiveness need to create arounda sufficiently the mid-20th strong the basic informations about the respondent (length century,scientific EE and which technical emerges base in to the accelerate 1960s as innovationa response of practice,The first partsex, ofapprobation), the questionnaire the second was focused part was on to the environmental crisis is in favor of the aim is focused on the implementation of CCSEE. The ques- to develop the environmental literacy necessary to tionnaire was sent out by a pilot survey (around understand the broader contexts from which these 20 respondents), after this phase some items were problems have arisen and are solved within them - (Wals et al. 2014). ry schools to which the questionnaire was sent out The possibilities of interconnection of EE with (3203),changed we or receivedclarified. responses From the totalfrom number640 respondents, of prima mathematics and physics, including concrete exam- after removing a few incomplete answers (return of about 20%). The headmasters of the schools listed in (2006) and Palivec (2013). The possibilities of inte- the Atlas of Education database (http://atlasskolstvi gratingples, are EE presented and social by sciences Sýkora (2007), in the curriculum Melichar et are al. .cz) were asked to send a request for forwarding to addressed by the example of Nigeria by Adedayo and the EE school coordinator. Olawepo (1997), Ferstl and Parkan (2007) summa- The length of teaching practice addressed by rize the possibilities of linking to history teaching. the EE coordinators varied from 1 year to 54 years. Numerous suggestions for linking environmen- The average length of practice was 19.1 years, most tal issues with language and literature teaching was respondents were women (86%). The most frequent provided by Bowers (2010), which emphasizes the - importance of using appropriate concepts in teach- phy, chemistry, mathematics and physical education. ing environmental topics and introduces misconcep- Oftenqualifications repeated of were respondents different combinationswere biology, ofgeogra these tions that may result from the use of inexact terms. school subjects. Possible reasons for children’s concepts and miscon- The main part of the questionnaire consisted of ceptions discusses Pavlátová (2019). Kubrická and items aimed at determining the rate of participation of individual school subjects in the implementation use of environmental topics for teaching English. The of CCSEE. For each school subject, respondents chose possibilitiesHromádka (2015) of linking provided EE with specific language examples and literature of the one of the following options: not involved or the sub- teaching are mentioned by Howard (2010), Lustyan- ject is not taught at our school (0), very little (1), mod- tie (2015) and Soetaert et al. (1996). The importance of linking different forms of artis- tic activities with EE is dealt with by Dielman (2013). erate (2), significantly (3), a core subject for achieving- Navrátil (2012) presents on concrete examples the quencyCCSEE objectives of individual (4). Aresponses coefficient in was each assigned category to each and theoption average (see above).was calculated. This coefficient Thus it was multiplied found out the how fre of key competences of EE according to FEP, through each subject participates in the realization of CCSEE. artisticpossibilities activities. of fulfilling the goals and development The total dataset was divided into two parts, the Various examples how to utilize works of art in geography and EE are also presented by Parkinson 1 of elementary school (n = 153) and the second part - withfirst with answers answers of teachers of teachers from who complete teach only (nine-year) at Stage elementary schools – with both Stage 1 and Stage 2 (2009),(2010), SánchezVočadlová (2013), (2009) Quigley or Kučera et al. (2012). (2014) Haloshow (n = 487). chathe possibilities(2008), Trojanová of developing (2009), Řezníčkothe skillsvá of a acquiring Boháček geographical or environmental information from incorporation of the EE topics into the curriculum, image sources or photographs. Several options for it wasTo findnecessary out a to typology choose ofthe schools proper accordingmethod. The to integrating environmental topics can also be found in metod would divide the objects (i.e. individual schools or better the answers of respondents from individual
(McNaughtonmusic education 2004). (Campos Integration 2013; with Váňová physical et al. educa 2007;- of particular school subjects that comprise EE topics Jurmu 2005) or in connection with drama education andschools) second into according categories to first intensity according of presence to composition of EE (Dechano and Shelley 2004). topics (i.e. extent of the EE curriculum) within these tion can be realized primarily through field activities subjects. Therefore, the multidimensional statistical method of hierarchical clustering was found as the 2. Research methodology most suitable for application. The method enables to divide the objects into categories according to mutual Data collection was carried out by anonymous on-line both similarity and dissimilarity of their characteris- questionnaire, which was addressed to the EE coor- tics. The analysis was conducted in the statistical soft- dinators at most of all elementary schools in Czechia. ware SPSS. Teaching of cross-curricular topic Environmental Education at elementary schools 203
The clustering of non-standardized variables was 3. Results of Research conducted. The variables were not standardized because all of them represent the same type of respon- The degree of involvement of individual school sub- dents’ answers, originating from the same time peri- jects in the implementation of EE objectives at prima- od, therefore they don’t vary in their values. During ry schools with Stage 1 only, is shown in Figure 1. The the hierarchical clustering the method of Average predominance of science-related subjects is evident, linkage between groups was applied to obtain max- - imal similarity within the groups together with the ject with mainly geographical and historical content) maximal dissimilarity between groups. The linkage of but the role of fine arts and homeland studies (a sub the variables (value of their distance) was measured The situation in complete (nine-year) elementa- with utilization of Pearson correlation intervals. Their ryis also schools significant. is shown in Figure 2. Biology is the most utilization ensures that the structural similarity of the important subject in these schools, but geography, answers is preferred – in this case the proportion of which is the second most important subject in this frequency in appreciation of EE topics between indi- vidual school subjects by the respondents. The educational objectives of CCSEE can be ful- Since the number of input variables (i.e. number of respect, also plays a significant role. school subjects) was too wide for such type of analy- also through a special separate school subject. One of sis, several groups of school subjects were created. At thefilled questions not only inin thealready questionnaire existing school survey subjects, was there but- Stage 1, the appreciation of basic biology (originally fore focused on using this option. Results show that přírodověda), homeland studies (originally vlastivě- it is used by 114 schools (18%). In about half of the da) and elementary teaching (originally prvouka) was cases the title of the subject contains the word ecology observed. At Stage 2 biology, geography and health or ecological - education were distinguished separately. The other tion of the subject of natural history or a practically subjects were grouped into: science subjects (phys- conceived subject. In the focused remaining mainly cases on it the is arealization modifica ics, chemistry, mathematics), languages (mother lan- guage, foreign language), humane science subjects etc. (history, civics), artistic and practical subjects (music, of scientificA simple experiments, typology of research-orientedschools was based teaching, on the
The number of 3 clusters was selected as the most respondents. The aim of this typology is to try to representativefine arts, physical number education). of clusters in the dataset of classifycontribution schools to accordingCCSEE goals to curriculum fulfilling, assessedstrategies byof teachers at Stage 1 and 5 clusters in the dataset implementation EE goals. of teachers at complete elementary schools. The clus- Applying multivariate analytical statistical meth- ters were tested about their independency at 95% ods (cluster analysis) it was possible to distinguish 3 different clusters of elementary schools with the means by method One Way ANOVA. Stage 1 only according to the strategy of integrating confidence interval through comparison of their
Basic Biology
Elementary Teaching
Fine Arts
Homeland Studies
Mother Language
Health Education