ACTA UNIVERSITATIS CAROLINAE
AUC GEOGRAPHICA 54 2/2019
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, 2019 ISSN 0300-5402 (Print) ISSN 2336-1980 (Online) Original Article 117
Territorial Impact Assessment – European context and the case of Czechia Štěpán Nosek*
Charles University, Faculty of Science, Department of Social Geography and Regional Development, Czechia * Corresponding author: [email protected]
ABSTRACT The aim of this paper is to critically review recent EU level debates on territorial impact assessment, which serves as a tool to improve the understanding of uneven territorial impacts of the EU sectoral policies. The paper also seeks to elicit (1) which Europe- an countries employ territorial impact assessment when designing various national policies and (2) how this tool is used in different governance environments. Particular attention is paid to the case of Czechia. The paper elaborates upon the state-of-the-art tools used on a national level and analyses the motivation of actors on regional and local levels to use such tools in their decision-making process. The research shows that EU member states that employ this tool can be categorized into four groups, for instance, depend- ing on whether they enrich other impact assessments tools by territorial aspects, or whether they focus coherence of regional development strategies and sectoral strategies. Territorial impact assessment is not implemented in Czechia, although recent years have witnessed an upsurge of interest concerning territorial aspects of investments.
KEYWORDS Territirorial Impact Assessment; European Union; Czechia
Received: 27 April 2019 Accepted: 4 July 2019 Published online: 20 August 2019
Nosek, Š. (2019): Territorial Impact Assessment – European context and the case of Czechia. AUC Geographica 54(2), 117–128 https://doi.org/10.14712/23361980.2019.1 © 2019 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). 118 Štěpán Nosek
1. Introduction with special attention to TIA. Chapter 3 analyses the current state of the art in selected European coun- For a long time, there has been a quest to understand tries in using TIA or similar tools. Chapter 4 focuses spatial or territorial consequences of public policies on Czechia and explains how TIA or similar tools are (Nijkamp and Van Pelt 1989). Some public policies are employed. Chapter 5 summarizes the main conclu- focused directly on increasing the development poten- sions and presents several suggestions for further tial of certain regions (regional policy). Other public development of TIA in the Czech context. policies are not directly focused on regions, but their From a methodological standpoint, the paper is impact is often greater compared to explicit region- based mostly on content analysis of available docu- al policies. For decades, there have been attempts to ments dealing with the topic of policy evaluation and develop tools which would enable to assess or to mea- (territorial) impact assessment in particular coun- sure a spatial impact of public policies. tries. This source of information has been comple- mented with questionnaires sent to academics and to assess the impact of policies or projects with a focus on Inenvironmental the 1970s, first and tools economic emerged impacts in the (e.g. USA Caldwell in order the study was to elaborate on the current situation in 1988; Bond and Pope 2012; Francis 1975). However, TIApublic use officials in all EU in countries.several countries. Thus, respondents The ambition were of geographers have been traditionally interested in asked an open question, i.e. to describe use of TIA in more complex tools assessing all aspects of territorial particular EU member states. Combination of those development. Discussions at the EU level regarding the two sources (content analysis and questionnaires) territorial dimension and territorial cohesion paved enabled to divide EU member states into four cate- the way for more complex tools enabling a holistic gories (see Chapter 3 and Table 2). Nonetheless, rel- approach in the evaluation of impact on territorial evant information has not been found in case of par- development (e.g. Medeiros 2017). ticular countries. The analysis concerning Czechia is based on this assessment (TIA) is a distinctly EU policy assessment author’s experience gained at the Czech Ministry closelyAccording connected to Fischer to the et al. EU territorial (2015), territorial cohesion impact agen- of Regional Development and at the International da and “does not make much sense in, e.g. the USA Advisory Centre for Municipalities, which imple- where people are used to move to areas of greater mented impact assessment techniques in the city of economic opportunity” (p. 9). The Territorial Agenda of the European Union 2020 (European Commission 2011) states that territorial cohesion “enables equal Litoměřice. opportunities for citizens and enterprises, wherever 2. The Role of Impact Assessment Tools they are located” (p. 3). Following the addition of ter- in Policy Coordination ritorial cohesion among strategic goals of the Europe- an Union in 2009, it seems to be even more important It is only natural that since the addition of territo- to evaluate potentially uneven territorial impacts of rial cohesion among the European Union’s goals, it particular projects, strategies, or policies. has become even more pressing to assess territorial Territorial Impact Assessment is a tool which impacts of sectoral policies that do not explicitly fol- appeared in a practical debate in 1989 with the launch of the European Spatial Development Perspective particular regional policies in different national con- (ESDP; Zonneveld and Waterhout 2009). textslow regional are, the policy Territorial goals – no Impact matter Assessment what the method goals of- - ology has been developed as a tool to measure poten- ically review recent developments regarding impact assessmentThe aims tools of this at thepaper EU arelevel twofold: and, secondly, firstly, toto anal crit- Marot 2011). yse awareness about these tools and the current state tialTIA territorial might effectsbe understood of sectoral as policies a new (Golobič addition and to a vast array of impact assessment tools that have Czechia. Czechia was selected for a case study due to been applied since the 1970s. Impact assessment, theof their fact thatimplementation the author of in this one paper EU member has followed state – in the which is mostly used and embedded in legislation of discussion and subsequent application of territorial many countries all over the world focuses primari- ly on economic and environmental dimensions (see years. Therefore, the paper seeks primarily to answer Table 1). However, proper evaluation of impacts on theseimpact questions: assessment in this country over the past five territorial development should include other dimen- (Q1) How do European countries deal with the sions as well. Stutz and Warf (2012) state that, “in challenge of evaluating territorial impacts of their conventional usage development is a synonym for policies? economic growth” (p. 268). Nevertheless, territorial (Q2) What is the current state of Territorial Impact development should be understood in a more com- Assessment implementation in Czechia? plex way, Medeiros (2017) concludes that territorial The paper unfolds as follows: Chapter 2 presents - an overview of different types of impact assessment ic competitiveness, social cohesion, environmental development should follow five main goals: econom Territorial Impact Assessment – European context and the case of Czechia 119
Tab. 1 Types of Impact Assessment.
Impact assessment focused on Complex Impact Assessments Economic dimension Environmental dimension Social dimension Examples Cost-benefit analysis Environmental Impact Poverty Impact Assessment Territorial Impact Assessment Assessment Regulatory Impact Health Impact Assessment Impact Assessment (European Union) Assessment Strategic Impact Assessment Source: Author sustainability, sound processes of territorial gover- assessments embedded in the law in other advanced countries. However, according to Nijkamp and Van or territorial articulation. TIA is then supposed to Pelt (1989), the OECD collected a number of case encompassnance, and all efficient those dimensions. processes of spatial planning studies of urban impact analyses in Sweden, Canada, or France. 2.1 Impact assessment – definition and origins In the late 1980s, the term strategic environmen- The purpose of impact assessment is to assess poten- tal assessment was coined in an interim report to tial impacts of certain actions before they are imple- the European Commission (Fungisland Tetlow and mented (Hayes 2017). Impact assessment is thus Hanusch 2012). Strategic environmental assessment supposed to evaluate what might happen as a direct of policies and strategies has been clearly developed consequence of a particular project, strategy, or a poli- upon the practice of Environmental Impact Assess- cy. Today, there are more than 40 types of IA discussed ment (EIA) focused on the level of projects. However, whereas EIA is primarily concerned with how a cer- most frequently used ones are Environmental Impact tain proposed development could affect the environ- Assessments,in literature (Morrison-SaundersStrategic Environmental et al. Impact 2014). Assess The- ment, Strategic Environmental Assessment (SEA) was ment, or Regulatory Impact Assessments, which are also supposed to analyse potential effects of alterna- embedded in law in most of countries. tive developments early in the decision-making pro- - cess (Fungisland Tetlow and Hanusch 2012).
EnvironmentalThe first impact Impact assessments Statement focused (EIS), whichon environ came 2.2 Territorial Impact Assessment of projects, intomental force impacts: on January the first 1, 1970 legally in thebinding USA IA(Bond was andthe policies and future trends Pope 2012). This tool was aimed at federal actions Since the addition of territorial cohesion among the goals of the European Union in 2009, the importance of evaluating potential territorial impacts of policies, usedthat mightfor large significantly projects such affect as dams the quality or nuclear of human power strategies or particular projects became yet more plantsenvironment but it was (Francis also used 1975). for assessmentAt first, the of tool smaller was relevant, at least in the European context. Neverthe- projects in urban environment relatively soon (Fran- less, policy-evaluation tools were developed already cis 1975). in 1980s as an integral part of the Structural Funds The system of impact analysis in the United States intervention process (Medeiros 2017). Between 1995 was strengthened in the 1980s, when the Urban and and 1999, the European Commission (EC) produced Community Impact Analysis (UCIA) was developed a detailed evaluation methodology under the MEANS as a consequence of Jimmy Carter’s presidential cam- programme, which inter alia clearly distinguished - - icy. The main aim of the impact assessment program nected to direct measurable consequences, the latter waspaign to promise ensure that to develop there will the be first no actionsexplicit contradicturban pol- focusedbetween on “results” long-term and consequences. “impacts”. The Medeiros first one (2017) con ing the urban policy (Hack and Langendorf 1980). analysed “the MEANS collection” and concluded that Glickman (1980) summarizes the methodology as “the notion of territorial impacts is entirely absent” being focused on assessing policies rather than proj- and that “no concrete references are made to the need ects and being focused on places rather than people for a more holistic impact assessment” (p. 149). UCIA was being focused on assessing macroeconomic Another set of documents dealing with TIA were impacts. Evers (2011) understands the Urban Impact the outcomes stemming from the EVALSED pro- Analysis as a clear predecessor to the TIA, since “it sought gramme concluded in 2008 (updated in 2013). to assess the impact of non-urban policies (region- According to Medeiros (2017), EVALSED took territo- al, local or even international) on urban areas” (p. 8). rial dimension of policy evaluation more into account. The implementation of the Urban Impact Analysis sparked a debate concerning potential unintended Assessment (IA) guidelines only in 2005. Surprisingly, spatial effects of public policies. With the exception theThe guidelines European did Commission not include released any reference its first toImpact “ter- of the United States, there were no impact analyses or ritorial impacts”; instead they have focused on the 120 Štěpán Nosek assessment of economic, social, and environmental introduction of the TIA Quick Check Tool. This interac- impacts in a way that is complementary to SEA, and tive web-based tool was introduced by ESPON in 2015 EIA (Tscherning, König, Birthe, Helming and Sieber and enables to measure territorial impacts at the level 2007). At the same time, during public consultation of the NUTS 3 regions. The TIA Quick Check Tool Is to the IA procedure, the EC published guidelines to based on the ESPON ARTS methodology and was pre- assess regional and local impacts, which could be pared by the Austrian Institute for Spatial Planning. However, according to Medeiros (2017), the TIA quick procedure (Medeiros 2017). consideredAs a result as theof allfirst these EC attempt programs to establishand activities, a TIA crucial components of territorial impact evaluation. the awareness about TIA, at least on the European Medeiroscheck tool (2017) is too simplifiedcalls for a and more does robust not include(and rele all- level, was growing. The European Commission doc- vant) evaluation technique, despite the fact that such techniques demand more resources and time. on Economic Social and Territorial Cohesion” includ- Overall, one can argue that TIA moved from a tool umented a statement “Investing that in “both Europe’s policies Future – Fifth with and without Report to be used for assessment of territorial impact of proj- ects (Austria, Germany) to a tool used for evaluation assessment of territorial impact” (European Commis- of potential territorial impacts of policies. There have sionan explicit 2010, p.spatial 195). dimension could benefit from an been even efforts to analyse territorial impacts of However, at the beginning (in the 1990s and future trends (for instance Böhme and Lüer 2017). the 2000s), TIA was conceived rather as an ex-post Broadly speaking, the methodology of TIA is sim- assessment tool of the EU policies such as the Com- ilar to those of EIA and SEA, or other impact assess- mon Agricultural Policy or Transport and Trans-Euro- ment tools. All of these methodologies “explore the pean Network Policies and was based on quantitative causal links between proposed actions and impacts” - The European Spatial Development Perspective models (Fischer et al. 2015). as(Perdicoulis prescriptive et al. regarding 2016, p. 42).which EIA, procedures SEA, or Cost should Ben development of a TIA methodology as one of key tasks beefit followed Analyses and (CBA) which methodologies aspects should can bebe considered addressed forAction the PlanEuropean adopted Spatial in Tampere Planning in Observation1999 defined Net the- - posed to be more holistic, and should also take into account(Fischer et al.,governance 2015). TIA,arrangements. on the other According hand, is sup to towork support (ESPON) numerous (Fischer research et al. 2015). projects. ESPON Thus, enabled it was Evers (2011), the emergence of TIA methodology has theto involve ESPON experts programme in the field (established of TIA and in subsequently 2002) that also been caused by the need “to gain information on - policy effects within an increasingly fragmented pub- odologies. According to Medeiros (2017), the ESPON lic sector” (p. 76). programmetriggered the funded development more than of first 20 rigorousanalyses TIAdedicated meth Medeiros (2017) sees TIA as a tool with a signif- icant potential to replace the EIA and SEA, since all but only a few of them could be regarded as TIA tools. existing TIA methodologies take into account the to the identification of impact of EU sectoral policies, environmental dimension. However, it might be into several period which demonstrates gradual shift a challenging task, given the path-dependency of pub- ofFischer TIA design. et al. (2015) divide ESPON methodologies lic policies. Consequently, when compared to its pre- decessors focusing on impacts in a narrow sense (EIA, use of quantitative models and ex-post assessment of SEA, CBA), TIA enables to evaluate a potential impact particularFirst period European (2004–2006) policies (e.g.was Commoncharacterized Agricul by- in a more complex way (Evers 2011). tural Policy). During this period for instance STIMA (Spatial Telecommunications IMpact Assessment) assessment tools and concludes that TIA is “the most methodology has been established (see ESPON 2004). Medeiros (2014) summarizes briefly all impact take into consideration all the aspects of territorial 2010) were still based on quantitative models, but developmentdifficult and complete (socio-economic, IA procedure, environmental, since it needs gov to- focusedMethodologies mostly on prepared ex-ante assessment. in second period This was(2008– the ernance and spatial organization)” (p. 198). However, case for instance for TEQUILA tool/model which according to Medeiros (2015), the territorial develop- included all dimensions and components of the con- ment does not equal the sum of socio-economic and cept of territorial cohesion (Medeiros 2017), but the environmental impacts. Medeiros (2018 and 2019) also - suggests to apply TIA to evaluate spatial planning pro- cesses and (Medeiros 2019) and sectoral policies, for resultsqualitative were models difficult and to more read participatory (Medeiros 2017). approach. Spe instance cross-border cooperation (Medeiros 2018). cificThird was periodARTS methodology (since 2012) whichis connected admitted especial use of- According to Medeiros (2015), the economic dimension of policy evaluation has always been omni- present because of the “strong position of economists methods.ly with EATIA EATIA methodology is still based which on ex-ante brings assessment. significant simplification of TIA procedure and use of qualitative Medeiros (2015, p. 9) also explains that complex IA in this particular scientific field of policy evaluation”. Another step in a quest for simplification was the Territorial Impact Assessment – European context and the case of Czechia 121 tools started to be tested after some of the EU stra- Regarding TIA implementation EU countries could tegic documents explicitly recognized the “territorial be distinguished into several categories which are not dimension as a key element of policy intervention”. mutually exclusive. First category consists of German Medeiros (2015, p. 9) suggests that there is a main- speaking countries (chapter 3.1). Second category stream vision that policies should contribute primar- represents countries where focus on compliance ily to the economic development. Nevertheless, he with regional or sustainable development strategy is calls for a more holistic approach that would take into required when sectoral policies are designed (chap- account several other dimensions of development ter 3.2). Countries in third category (chapter 3.3.) such as social cohesion, environmental sustainability, enrich traditional impact assessment tools by terri- territorial governance and spatial planning. Similar- torial aspects. Countries in fourth category (chapter ly, Zonneveld and Waterhout (2009) argue that TIA 3.4) elaborate Territorial Impact Assessment ad-hoc “has to focus on the impact of policy and contextual in order to assess impact assessment of large projects, developments on the spatial organization and spatial or European directives/policies. position of spaces and places” (p. 3). 3.1 TIA as a traditional tool in German speaking countries 3. TIA implementation in EU member states Tools similar to TIA are used traditionally mostly in German-speaking countries. The roots of TIA can be TIA at the EU level certainly represents a valuable traced back to tools used traditionally in Austria, Ger- tool that helps to understand unintended territori- many, and Switzerland. Even the fact that the term ter- al impacts of EU policies or new legislature. Never- ritorial impact assessment is directly translated from its German equivalent demonstrates that the roots of the tool can be traced in the German-speaking world Therefore,theless, such TIA a toolat the at theEU EU-levellevel should cannot be reflectfollowed all (Othengrafen and Cornett 2013; Healy 2001). In Aus- byspecificities national TIAsat the (Zonneveld national, regional, and Waterhout or local 2009).levels. tria and Germany, TIA is used extensively for ex-ante However, according to Medeiros (2017), no EU mem- evaluation of major projects. ber state or other country has encoded a mandatory Particularly in Austria, the use of TIA can be traced use of TIA. Consequently, the TIA procedures could be back to 1959, thus even predating the introduction labelled as an “EU experiment” (Medeiros 2017). of EIA (Healy 2001). In German-speaking countries, In European countries, different techniques are competences of different levels of public administra- employed to ensure that sectoral (non-regional) tion (national, regional, and local) are set out very interventions are in line with the principles of region- - ous factors, including the governance system, admin- theclearly federal (Tosics systems et al. 2010). of Germanic This holds countries true for result spatial in al policy (no matter how defined). Obviously, numer aplanning strong regional as well. level According of planning. to Tosics et al. (2010), way (territorial) impact assessment is implemented Particular authors even declare that TIA is embed- (Meulemanistrative cultures, 2015). and Meuleman political settings,(2015) suggests influence that the ded in Austrian, German, or Swiss Law (e.g. Zonne- the governance approach may vary: it can be more veld and Waterhout 2009). This was, nevertheless, market-driven in Anglo-Saxon countries, more net- dismissed by Dallhammer (2016) who argues that work-driven in northwest Europe, and more hier- “Raumverträglichkeitsprüfung” (which is a tool shar- archical in Central, Eastern, and Southern Europe. ing similarities with TIA) is proceeded on volun- Each governance system implies certain weaknesses, tary basis and only in particular federal states (e.g. - Carinthia). en model, the risk of losing focus in a network-driv- Paradoxically, Austria does not have a very strong ensuch approach, as prioritization or rigidity of efficiency in the hierarchical in a market-driv model explicit regional policy (Polvevari and Michie 2011). (Meuleman 2015). The Austrian Conference on Spatial Planning has Within the EATIA project, there were some, albeit been established only in 1971 in order to enhance limited, efforts to disseminate the European TIA meth- cross-sectoral coordination. The body comprises odology to particular EU member states (Slovenia, Por- representatives of the state (Bund), federal states (Länder), and municipalities. Each ten years the Aus- of the conclusions stemming from the dissemination trian Spatial Development Concept is adopted to coor- tugal, the United Kingdom) (Fischer et al. 2015). One- dinate federal policies with spatial impact and poli- ible enough to be adaptable to different policy-making cies implemented by Länder (Polvevari and Michie of results was a suggestion that the TIA concept is flex 2011). Similarly, in Germany, sectoral policies are (e.g. Golobic and Marot 2011) and Portuguese (e.g. coordinated through the Regional Joint Task (Polve- traditions (Fischer et al. 2015). Although Slovenian vari and Michie 2011). Germany evaluates ex-ante impacts of trans- PortugalMedeiros TIA 2015) procedures academics are contributedstill rather unknown significantly con- port infrastructure or large retail projects “to verify ceptto the (Marot research 2015; in the Medeiros field of 2015).TIA, both in Slovenia and whether these are in line with the aims and objectives 122 Štěpán Nosek
2013, p. 13). The process (called Raumordnungsver- evaluates impacts on services, infrastructure, busi- offahren) official usually planning takes policies” around (Othengrafen three months. and AccordCornett- ness,impacts employment, on rural areas environment (Hague 2010). etc. Rural (Department proofing ing to Healy (2001), the German procedure is focused for Environment, Food and Rural Affairs 2017). In on testing the conformity of new projects and new 2011, IA Toolkit was published by the Department for regional plans with existing plans. Bundesländers Business, Innovation and Skills (2011). In the British (sixteen units) are in charge of this process. The use of impact assessments in the United Kingdom 3.2 TIA as tool to align territorial priorities with hascontext, been IA discussed is based inmostly connection on a cost-benefit with Brexit, analysis. but no sectoral and regional strategies impact assessment of Brexit has been conducted so Such strategy is relevant for instance for Switzerland, far. This was later criticized by Fischer and Sykes (2017), who suggested that impact assessment tak- states that the implementation of the document is ing into account not just economic impact but a wider coordinatedwhere Sustainable by the Development interdepartmental Strategy 2016–2019 Sustaina- scale of impact should be conducted in case of Brex- ble Development Committee, which is supposed to it, because it is clear that it is going to have uneven ensure the incorporation of sustainable development regional impacts. principles into sectoral policies (Swiss Federal Coun- cil 2016). 3.4 Ad-hoc use of TIA Similar strategies are employed in Ireland, where Such TIAs were applied in Netherlands, Poland, or the compliance of strategies at both local and regional Lithuania. In the Netherlands, there is no standard levels with the Irish National Spatial Strategy is a key TIA procedure. Large projects are usually assessed requirement (Hague 2010), and in Estonia (Ministry of Regional Development 2018) where compliance with the most important elements of the Estonian analysisvia a cost-benefit (Evers 2016). analysis. Nonetheless, Currently, thethere Netherlands are efforts experiencedto include a social an ad-hoc dimension implementation via a social of cost-benefit TIA in con- assessed while preparing new regulations, strategic nection with possible regional implications of the policyRegional documents, Development and grant Strategy schemes for 2014–2020 supported by is EU’s policies. According to Evers (2011), the Dutch the EU Structural Funds. Ministry of Housing, Spatial Planning and the Envi- ronment requested the Environmental Assessment dimension of strategic documents have been pre- Agency to conduct TIA as a reaction to the European paredIn Poland in order specific to bridge document the gap entitledbetween Territorialterritorial territorial cohesion policy and to the energy-climate and sectoral strategic documents (Szlachta, Zaucha package. Evers (2011) stresses that “given the two and Komornicki 2017). Similar tool with focus on very different points of the policy process, the studies operational programmes have been invoked in Cze- are very different in character in terms of methodol- chia (see chapter 4). ogy and conclusions” (p. 98). The TIA focusing on the European territorial cohesion policy was affected by 3.3 TIA as an enrichment of traditionally used tools This strategy is relevant for Finland, Belgium, or contrast, the energy-climate package had clearly stip- Denmark. Those countries use standard EIA, SEA, not having a clear definition of territorial cohesion. In or Regulatory Impact Assessment (RIA) procedures policy actions in order to meet targets. Ad hoc impact and “upgrade” them, at least in some cases in order assessmentsulated goals areand undertaken it was thus asambition well for oflarge TIA projects to find to assess wider development processes. Finnish EIA (Poland), or after implementation of the EU program- aims to evaluate not only the impact on human health ming period (Lithuania). and environmental aspects (soil, water, air, climate etc.), but also the impact on “the urban structure, buildings, landscape, townscape and cultural herit- 4. Summary age”. (Ministry of Environment Finland 2009). Simi- larly, in Wallonia the EIA procedure was adjusted in Table 2 shows that there is a rather limited use of order to encompass certain aspects of TIA (Hague Territorial Impact Assessment or similar tools across 2001; Böhme and Lüer 2017). Denmark held discus- European countries. Nevertheless, according to sever- sions regarding a holistic approach within the EIA al studies (e.g. Dabrowski 2012; Raagmaa and Stead procedure currently, as well (Kørnøv, Christensen and 2015), a certain trend towards Europeanization of Nielsen 2005). In another cases the RIA is “upgraded”. public policies can be observed especially in Central This is the case for Belgium and its Flemish Region, where RIA includes the Sustainable Development planning and regional policy. This might be seen, for Impact Assessment (OECD 2010; Huge and Waas instance,and Eastern in Partnership Europe, especially Agreements in the signedfields of between spatial 2011). the European Commission and particular member - states, which include an analysis of tools to be used by the member states to enhance the territorial approach The United Kingdom has specific impact assess ment called rural proofing, which evaluates potential Territorial Impact Assessment – European context and the case of Czechia 123
Tab. 2 Typology of TIA use in EU member states and Switzerland.
Type of TIA use Description Countries
TIA used extensively for ex-ante impact TIA encoded in the law Germany, Austria assessment of large infrastructure projects TIA as tool to align of territorial priorities Compliance of policies, plans, or strategies with sectoral and regional strategies compared with main regional development Estonia, Ireland, Switzerland, Poland strategies TIA as an enrichment of traditionally used Territorial aspects more or less incorporated Belgium, Finland, Greece, Denmark, Great Britain tools into existing evaluation techniques TIA used “when needed” such as before Ad-hoc TIA implementation of European directives and Netherlands policies into national contexts Mostly Central and European countries Slovakia, Slovenia, Hungary, Lithuania, Latvia, No TIA procedure declare no use of standard TIA procedures Romania, Bulgaria, Croatia, Czechia, Portugal During desk research no relevant information Sweden, Malta, Cyprus, Spain, Italy, France, No relevant information founded regarding TIA in seven EU countries has been Luxembourg found Source: Author in implementation of the EU cohesion policy (Nosek 2017). Ministry of Industry and Trade or Ministry of Envi- The position of (spatial) planning in Central and ronment.often considered As a consequence, less influential certain in ministriescomparison see with the territorial dimension concept as a threat and inter- mostly because planning of any kind was associated vention into its competencies. withEastern the European communist countries regime in by the both 1990s policy-makers was difficult, Czechia is one of the countries which do not use and population in general (Raagmaa and Stead 2015). TIA methodology at the national level. However, an There is no doubt that an EU membership push- interest in the tool has been growing in recent years, es member states to adopt changes in the strategic similarly to other Central and Eastern European planning process (Dabrowski 2012), including an countries and certain measures have been taken to application of more sophisticated methods in impact assessment. Therefore, it can be expected that some sort of mapping of territorial impact will become adjust sectoral policies in order to reflect territorial more embedded in spatial planning processes in the specificities. near future. 5.1 TIA at the policy level in 2014–2020 Namely, some Central and Eastern European coun- programming period tries plan to apply impact assessments tools over the According to the Partnership Agreement between next few years. This is the case, for instance, in Roma- the European Commission and Czechia concluded in nia, which expects that the methodology for mea- 2014, the territorial dimension is in the Czech context suring territorial impact of sectoral policies will be understood as a possibility to concentrate resources approved by the Government following the approval from the European Structural and Investments Funds of the Territorial Development Strategy (Ministry of Regional Development 2018). competitiveness and moderating territorial dispari- ties(ESIF) (social, in specific economic, types environmental). of territories supporting both In Czechia, a relatively comprehensive approach 5. The case of Czechia to incorporate the territorial dimension into the pro- - - oped. In particular, a special document called the ticular Czech ministries, the Ministry of Regional Nationalgrammes Document within the for period Territorial 2014–2020 Dimension was (NDTD; devel DevelopmentAccording to theis obliged law defining to coordinate responsibilities other ministries of par Ministry of Regional Development 2015) was pre- in order to implement regional policy. Such setting might seem favourable for pushing sectoral ministries paredThe with incentive the aim for itsto existenceidentify, whichcame implicitly specific goals from Nonetheless, the actual promotion of regional theshould European be preferred Commission, in specific which types pushed of territories. for more toaspects reflect into regional sectoral specifics. strategies and policy-making territorially sensitive implementation of European process cannot be taken for granted. Such promotion Structural and Investment Funds. This document also demands numerous discussions and eventually also served as a “substitute” for Czech NUTS II regions - which implemented its own regional operational pro- ly lacking. The Ministry of Regional Development is political support – which is in Czech context frequent grammes in 2007–2013 programming period. 124 Štěpán Nosek
stipulating the territorial dimension within nation- al grant schemes or national sectoral policies. For ruralOverall, areas, for eight socially types excluded of specific localities, territories for solving have example, even the Regional Development Strategy labourbeen identified market problems,in the NDTD for (e.g. enhancement for urban areas,of entre for- preneurship etc.). NDTD could be understood as a kind of ex-ante Czech2014–2020 Regional does Development not stipulate Strategyhow sectoral only policies identi- TIA. This document indicates which measures will are supposed to reflect territorial specificities. The- most likely have uneven territorial impacts, similar- ed powers as “economically weak” and without any ly to concepts described in the previous chapter and fies micro-regions of 57 municipalities with extend- applied by ESPON or by particular countries (e.g. gions should be taken into account within sectoral Austria or Germany). However, in contrast to ESPON policies.further specification states that these micro-re models, in Czechia, no data analysis has been per- Each year, progress evaluation in implementa- tion of the NDTD is to be prepared. This evaluation have the territorial dimension. should be primarily based on an analysis of particular formed in order to decide which specific goals should calls for proposal and on an analysis of allocation of the state will distribute via the territorial dimension. funds. So far, evaluation has been prepared for years ResourcesThe NDTD used also for specifies addressing the thevolume territorial of funds dimen that- 2015, 2016, and 2017. According to the evaluations sion include funds distributed through tools designed to support urban development (i.e. Integrated Terri- torial Investments and Integrated Territorial Devel- Analysesfor the period of spatial 2015–2017, distribution the of calls allocation were of in funds most opment Plans) and tools aimed at rural development havecases been targeted limited at territories by poor quality specified of available by the NDTD. data- (Community-led Local Development). The total allo- bases, in which a large proportion of projects did not cation is expected to be approximately 8 billion EUR, include information concerning their localization at which is almost one third of the Czech allocation from the municipal level. Nonetheless, available data sug- the EU European Structural and Investment Funds. gest that the share of the Cohesion policy expected by The very existence of the NDTD might be seen as NDTD to be spend on the territorial dimension has progress in promoting the territorial dimension into been achieved. sectoral policy making. However, the NDTD was the The analysis of individual calls for proposals raised several important questions. First, the pivotal ques- and suffers from several weaknesses. The NDTD was tion was which calls should be declared as taking into first document of its kind implemented in Czechia account the territorial dimension. The NDTD did not - include any rigorous methodology in order to assess arguablyritorial dimension created too was late, thus since not it implementedwas finalized wellinto calls with regard to its territorial dimension. Eventu- afterthe operational the finalization programmes, of operational but programmes.it was described Ter ally, four types of possible scenarios have been identi- ex-post, based on discussion with responsible minis- tries/managing authorities. The NDTD does not cover national grant schemes fied by the Ministry of Regional Development: - and currently, there is no alternative document ry (for instance, on a selected NUTS 2 region). – The call is focused exclusively on a specific territo
Tab. 3 Territorial dimension in operational programmes in Czechia in 2014–2020.
Level of Territorial Operational Programme Description of Territorial Dimension Dimension Specific calls for certain type of regions (e.g. regions with socially excluded Integrated Regional OP High localities); Funds in certain calls distributed in accordance with regional needs. OP Entrepreneurship and Innovation for Regions with high unemployment rate or economically weak regions High Competitiveness benefited in selection process. Projects contributing to regional development priorities supported OP Transport Medium primarily. Most of the activities within Community-led Local Development financed via Rural Development Programme Medium Rural Development Programme. Local and regional action plans in education define local/regional priorities OP Research, Innovation and Education Medium in elementary schools and high schools. In certain cases specific territorial dimension applied – for instance localities OP Environment Low specified in flood plans. OP Employment Low Except of labour market policy. OP Prague – Growth pole Not Relevant Focused only on the Capital City of Prague. Source: Author Territorial Impact Assessment – European context and the case of Czechia 125
It might be argued that “territorial” stakeholders aim to ensure allocation of funds in future policies – The call explicitly declares the distribution of funds implementing the strategy. In their understanding, regionsamong a (e.g. specific those type with of regions/municipalities.an above-average rate of more territoriality equals higher chances to get fund- – unemployment).The call sets up the criteria favouring particular ed. Ministries, on other hand, aim to avoid, from their point of view, unnecessary intervention into their pol- that are concentrated in certain regions/territories icies. In such context, the role of Ministry of Regional – (i.e.The implicitcall is focused territorial on addressing dimension). specific problems Development might be seen as of a mediator.
When applying such criteria, it appeared that most 5.3 TIA at the Project level calls that were expected to have a territorial dimen- sion indeed incorporated it in one way or another into in Czechia at the level of projects. Finally, in 2016, the the calls’ designs. MinistryUntil 2016, of Regionalthere was Development no certified publishedTIA methodology project The Ministry of Regional Development itself man- impact assessment methodology prepared in cooper- ation with regional policy experts. Following the pub- Regional Development Programme, which is seen as lication of this methodology, the Ministry of Regional ages in 2014–2020 programming period Integrated Development plans to start a project (to be imple- important aspects of regional development were howevera flagship implemented programme forby otherregions (sectoral) and cities. ministries, Certain as well. bemented implementing in 2019–2020) the methodology aimed at raising on all awareness hierarchical at Table 3 summarizes operational programmes the methodology and to train public officials who will also a plan to develop a software tool to guide poten- territorial dimension expected by National Document tiallevels – national, users through regional,the methodology. and municipal. There is forimplemented Territorial in Dimension. 2014–2020 programming period and This TIA methodology has been inspired by simi- lar documents based on the initiative of the Europe- 5.2 Territorial dimension expectations in post-2020 an Commission (especially EC 2009), which includes period a list of potential impacts that should be taken into Discussions regarding implementation of post-2020 account when evaluating possible impacts of EU pol- programming period are already underway. Following icies. Therefore, the Czech methodology designed at the project level has been developed on the basis of large decrease of European funds allocation in Czechia the EU methodology, which, however, targeted the lev- issignificant expected growth in post-2020 of Czech period. economy Thus, in recentthe Regional years, el of policies. The TIA procedure is divided into two steps. In the territorial dimension within particular national sec- toralDevelopment policies inStrategy addition post-2020 to the territorial aims to dimensiondefine the the project proposal and decide whether the project embedded in operational programmes supported via first step, the evaluator must assess overall quality of the European Structural and Investment Funds. Discussion with bodies responsible for sectoral employment,will have a significant environment, territorial transport, impact market, or not. quality The policies again shows scepticism of sectoral minis- ofpotential life, and impact public is services. evaluated The in evaluatorthe following must fields: also tries regarding the role of Ministry of Regional Devel- assess what the level of impact is (average, important, opment. Promotion of regionally sensitive policies and fundamental). is often understood as a unnecessary intervention Project proposals that are considered as those with a potential socioeconomic impact are subject to fur- process. ther analysis. The subsequent procedure of the TIA whichRegional put artificial Development element Strategy into the post-2020 policy-making aims differs based on the scale of the potential impact in
Prior to 2016, TIA or similar methodologies were comparisonto elaborate to only its predecessor. topics with identifiedThus, key territorialquestions onlythe above-listed exceptionally fields. applied in Czech regions and/or appearingdimension andduring is supposed most of theto be discussions more specific are: with (1) municipalities. One of these exceptions was the city which topics have territorial dimension (i.e. needs - territorially sensitive solutions), and for those which sentatives decided to establish a procedure analog- icalof Litoměřice to TIA in order in Northern to measure Bohemia, the possible whose impact repre of project proposals supported by the city from the motivationshave, (2) how of itparticular can be defined stakeholders. in particular While fields?repre- economic, social, and environmental point of view. sentativesApparently, of the rural debate areas, is cities often and affected regions by haveconflicting a ten- The assessment is undertaken by “Sustainable Devel- dency to promote territoriality as much as possible, opment Strategic Team” consisting of the mayor or his/her deputy, members of the department of stra- (and respective policies) are spatially blind and no tegic development and the department of sustain- ministerial officers tend to declare that certain fields able development, and members of other relevant significant regional adjustments are needed. 126 Štěpán Nosek departments (depending on type of the project). The The second question focused solely on the recent resulting assessment is subsequently presented to development of TIA in Czechia. In this country, in the political representation and forms the basis for order to encompass the territorial dimension within the decision-making process. This process is sim- ilar to the one employed in Austrian Carinthia, as programming period, the National Document for Ter- described in Chapter 3. However, city representatives ritorialthe operational Dimension programmes has been approved for the EU by 2014–2020the govern- are not obliged to follow recommendations of this assessment. some shortcomings, it forced the sectoral ministries toment. take It the is the territorial first document dimension of itsinto kind account and, despiteby one way or another during implementation of operation- 6. Conclusions al programmes they are responsible for. Nevertheless, a large variety of interpretations of the term territorial Studies in the sphere of regional policy and region- dimension have been observed during the implemen- al development have emphasized uneven territorial tation phase of particular operational programmes. impacts of non-regional public policies for a long time. - Over the past two decades, the search for a sophis- ment has been recently accomplished in Czechia, as in ticated tool that would enable to analyse territorial 2016Finally, a TIA at methodology the project atlevel, the aproject significant level hasachieve been impacts of proposed interventions has already found published by the government. This came just in time several potential solutions. to assist some more active cities and regions in their Nevertheless, the ambitions behind the search for attempt to design an impact assessment approach to their own projects. Thus, the recently published TIA the one hand, the tool is supposed to be sophisticated methodology can stimulate stakeholders at the local enougha suitable to TIA take methodology into account seem all potential to be conflicting. territorial On and regional level to deal seriously with uneven ter- impacts, while, on the other hand, it should be sim- ritorial impacts of various public interventions. The extent to which this potential to deliver more effective in the sphere of regional development at all levels of publicple enough administration. to serve the This purposes is hardly of realistic. public officials Never- regional actors remains to be seen theless, the paper suggests that the territorial impact andConsequently, efficient public a policiesthorough will analysis be used ofby advanceslocal and assessment is a tool that is developing quite dynam- - ically both at the EU level and at the level of member ject to further research. states. made in the field in various countries should be sub The process of developing TIA methodologies is still under-researched. The aim of this paper was to References contribute to research of TIA via an analysis of the extent and form to which the EU member states imple- Bond A., Pope J. (2012): The state of the art of impact ment the TIA methodology within their policy-making assessment in 2012. Impact Assessment and Project https://doi.org/10.1080/14615517 Czechia, which recently witnessed an upsurge in the .2012.669140. Appraisal 30(1), 1–4, processes. Specific attention was paid to the case of (2011): How to strengthen the territorial dimension of the national level. Böhme,Europe K., 2020Doucet, and P., the Komornicki, EU Cohesion T., Zaucha, Policy. Report J., Świątek, based D. interest of public officials in the TIA methodology at on Territorial Agenda 2020. Retrieved March 1, 2018, on the current state of affairs in European countries http://ec.europa.eu/regional_policy/sources/docgener regardingThe first the question use of posedTerritorial in this Impact paper wasAssessment focused /studies/pdf/challenges2020/2011_territorial_ or a similar tool. It can be concluded that while the dimension_eu2020.pdf. European Commission has been a clear leader in Böhme, K., Lüer, Ch. (2017): Assessing territorial impacts of developing new (and in the past decade largely sim- future trends. In: Medeiros, E.: Uncovering the Territorial Dimension of European Union Cohesion Policy. Cohesion, EU member states have not paid similar attention Development, Impact Assessment and Cooperation, toplified) measuring methods territorial of territorial impacts impact despite assessment, the recent the https://doi.org/10.4324/9781315169743-11. incorporation of territorial cohesion among the Euro- Caldwell, L. K. (1988): Environmental Impact Analysis (EIA): Origins, Evolution, and Future Directions. Impact pean Union goals. From author’s overview follows https://doi.org/10.1080 Germany and Austria lead the way towards sound TIA /07349165.1988.9725648. implementation. Several other countries aim at intro- Dabrowski,Assessment M. (2012):6(3–4), Shallow75–83, or deep Europeanisation? ducing at least simple methods to measure uneven The uneven impact of EU cohesion policy on the regional territorial impacts of public policies (Estonia, Roma- and local authorities in Poland. Environment and nia). This process might be seen as another proof of Europeanization of public policies in Central and https://doi.org/10.1068/c1164r. Eastern European countries (Dabrowski 2012; Raag- Dallhammer,Planning C: E. Government (2016): Personal and Policy communication, 30(4), 730–745, maa and Stead 2015). 9 February, 2016. Territorial Impact Assessment – European context and the case of Czechia 127
Department for Business, Innovation and Skills (2011): Francis, M. (1975): Urban Impact Assessment and Impact Assessment Toolkit. Retrieved May 6, 2018, Community Involvement: The Case of the John Fitzgerald www.legislationline.org/documents/id/16803. Kennedy Library. Environment and Behavior 7(3), Department for Environment, Food and Rural Affairs https://doi.org/10.1177/001391657500700307. Fungisland Tetlow, M., Hanusch, M. (2012): Strategic impacts of policies on rural areas. Retrieved May 5, 2018, environmental373–404, assessment: the state of the art. Impact (2017):https://www.gov.uk/government/uploads/system Rural Proofing. Practical guidance to assess https:// doi.org/10.1080/14615517.2012.666400. -guidance.pdf. Glickman,Assessment N. J. (1980): and Project Urban Appraisal Impact Analysis:30(1), 15–24, Premises, ESPON/uploads/attachment_data/file/600450/rural-proofing (2004): ESPON 2.1.1: Territorial Impact of EU Promises, Procedures, and Problems. Built Environment Transport and TEN policies. ESPON (2013): ESPON ARTS: Assessment of Regional and Golobic, M., Marot, N. (2011): Territorial impact Territorial Sensitivity. Final Report. assessment:6(2), 84–91. Integrating territorial aspects in sectoral European Commission (2009): Impact Assessment Guidelines. Retrieved April 9, 2018, http://ec.europa.eu https://doi.org/10.1016/j.evalprogplan.2011.02.009. /smart-regulation/impact/commission_guidelines Hack,policies. G., Langendorf, Evaluation R. and (1980): Program Lessons Planning from 34,Urban 163–173, /docs/iag_2009_en.pdf. Impact Assessment. Environmental Impact Assessment European Commission (2010): Fifth Report on Economic, https://doi.org/10.1016/S0195 -9255(80)80103-X. future. Retrieved April 8, 2018, http://ec.europa.eu Healy,Review A. (2001): 1(3), 251–266, Keywork Phrase: Territorial Impact /regional_policy/en/information/publicationsSocial and Territorial Cohesion – Investing in Europe’s Assessment. Retrieved June 18, 2017, http://www .esprid.org/keyphrases/30.pdf. -territorial-cohesion-investing-in-europe-s-future. Huge, J., Waas, T. (2011): Converging impact assessment European/reports/2010/fifth-report-on-economic-social-and Commission (2011): Territorial Agenda discourses for sustainable development: the case of of the European Union 2020. Retrieved April 7, 2018, Flanders, Belgium. Environment, Development, and http://ec.europa.eu/regional_policy/en/information https://doi.org/10.1007 /publications/communications/2011/territorial /s10668-010-9279-z. -agenda-of-the-european-union-2020. Kørnøv,Sustainability L., Christensen, 13(3), 607–626,P., Nielsen, E. H. (2005): Mission European Council of Town Planners and Committee on impossible: does environmental impact assessment in Spatial Development (2001): Proceedings of the one-day Denmark secure a holistic approach to the environment? conference on Territorial impact assessment 26 October 2001 Louvain-la-Neuve, Belgium. Retrived March 14, https://doi.org/10.3152/147154605781765427. 2018, http://www.ectp-ceu.eu/images/stories/download Marot,Impact N. (2015): Assessment On the and maximal Project Appraisalshortest path 23(4), in 303–314, /cp20030515.pdf. a connected component in V2V. Performance Evaluation Evers, D. (2011): Territorial Impact Assessment: a critical https://doi.org/10.1016/j.peva examination of current practice. PBL Netherlands .2015.09.003. Environmental Assesment Agency. Medeiros,94, 25–42, E. (2014): Assessing territorial impacts of the EU Evers. D. (2016): Personal communication, 18, January, cohesion policy at the regional level: the case of Algarve. 2016. Impact Assessment and Project Appraisal 32(3), Hayes, A. C. (2017): What is impact assessment? Some https://doi.org/10.1080/14615517.2014 .915134. posthumously by Adrian C. Hayes. Impact Assessment Medeiros,198–212, E. (2015): personal communication, October 8, personal reflections C. P. Wolf (1933–2015),https://doi.org edited 2015. /10.1080/14615517.2017.1322812. Medeiros, E. (2015): Territorial Impact Assessment and Hague and (2010): Project EuropeanAppraisal Perspectives35(3), 186–199, on Territorial Impact Public Policies: The Case of Portugal and the EU. Assessment. A Background paper for the ESPON EATIA Medeiros, E. (2017): From simple to relevant TIA tools Project. Retrieved March 11, 2018, http://www.rtpi.org for European policies. In: Medeiros, E.: Uncovering The .uk/media/5994/European-Perspectives-on-Territorial Territorial Dimension of European Union Cohesion -Impact-Assessment-April-2010.pdf. Policy. Cohesion, Development, Impact Assessment and Fischer, T. B. (2014): Impact Assessment: there can be Cooperation, https://doi.org/10.4324/9781315169743. strength in diversity! Impact Assessment and Project Medeiros E. (2018) Focusing on Cross-Border Territorial https://doi.org/10.1080 Impacts. In: Medeiros E. (eds) European Territorial /14615517.2013.872844. Cooperation. The Urban Book Series. Springer, Cham, Appraisal 32(1), 9–10, https://doi.org/10.1007/978-3-319-74887-0. P., Waterhout, B., Perdicoulis, A. (2015): Territorial Medeiros, E. (2019) Spatial Planning, Territorial Fischer, T. B., Sykes, O., Gore, T., Marot, N., Golobič, M., Pinho, Development and Territorial Impact Assessment. Journal The Emergence of a New Policy Assessment Instrument. https://doi.org Impact Assessment of European Draft Directives –https://doi /10.1177/0885412219831375. .org/10.1080/09654313.2013.868292. Meuleman,of Planning L. (2015): Literature Owl 34(2), meets 171–182, beehive: how impact Fischer,European T. B., PlanningSykes, O. (2017):Studies 23(3),Impact 433–451, Assessments: assessment and governance relate. Impact Assessment What Impact Assessments? And what is an Impact https://doi.org Assessment? Retrieved March 10, 2018, https://news /10.1080/14615517.2014.956436. .liverpool.ac.uk/2017/12/08/impact-assessments-what Ministry and Project of Environment Appraisal Finland33(1), 4–15, (2009): Act on -impact-assessments-and-what-is-an-impact-assessment.
Environmental Impact Assessment Procedure, Unofficial 128 Štěpán Nosek
translation. Retrieved August 24, 2017, https://www . The (in)coherence of Cohesion policy. European Policy Ministry of Regional Development (2015): National Polverari,Research L., Center,Michie, No.R. (2011): 78, University Complementarity of Strathclyde or conflict? 40 document.finlex.fi/en/laki/kaannokset/1994/en19940468.pdf for Territorial Dimension. George Street Glasgow G1 1QE, http://www.eprc-strath. Ministry of Regional Development (2018): Questionnaire eu/public/dam/jcr:23283c26.../EPRP%2078.pdf. among National Territorial Cohesion Contact Points. Raagmaa, G., Stead, D. (2015): Spatial Planning in the Baltic Morrison-Saunders, A., Pope, J., Gunn, J. A. E., Bond, A., States: Impacts of European Policies IN Raagma, G., Stead, Retief, F. (2014): Strengthening impact assessment: a call D. (eds.). Impacts of European Territorial Policies in for integration and focus. Impact Assessment and Project https://doi.org/10.1080/14615517 .2013.872841. Stutz,the F.P.,Baltic Warf, States. B. (2012): London: The Routledge – Taylor World Economy: & Geography, Francis Nijkamp,Appraisal P., van 32(1), Pelt, 2–8, M. (1983): Spatial Impact analysis for Business,Group, 1–9. Development. Hoboken, NJ: Pearson developing countries. A framework and a case study, Education. Serie Research Memoranda 0014, VU University Szlachta J., Zaucha, J., Kormonicki, T. (2017): Polish Amsterdam, Faculty of Economics, Business development policy and its territorial dimension. In.: J. Administration and Econometrics. Bradley, J.Zaucha (eds.): Territorial Cohesion: A missing Nosek, Š. (2017): Territorial cohesion storylines in link between economic growth and welfare. Lessons
https://doi.org/10.1080/09654313 Swiss Federal Council (2016): Sustainable Development .2017.13490792014–2020 Cohesion. Policy. European PlanningSstudies from the Baltic tiger, Gdansk, 49–72. OECD25(12), (2010): 2157–2174, Better regulation in Europe. Retrieved May www.eda.admin.ch/agenda2030/en/home.html. 17, 2018, http://www.oecd.org/gov/regulatory-policy Tosics,Strategy I., Szemzö, 2016–2019. H., Illés, Retrieved D., Gerheis, April A., 23, Lalenis, 2018, K., https:// /betterregulationineuropebelgium.htm. Kalergis, D. (2010): National spatial planning policies OECD (2011): Territorial Reviews Switzerland. Retrieved and governance typology, www.plurel.org/images May 15, 2018, https://read.oecd-ilibrary.org/urban-rural /D221.pdf. -and-regional-development/oecd-territorial-reviews Tscherning, K., König, H., Birthe, S., Helming, K., Sieber, S. -switzerland-2011_9789264092723-en#page1. (2007): Ex ante impact assessments (IA) in the European Othengrafen, F., Cornett, A. P. (2013): A Critical Assessment Commission: An overview. In: Helming, K., Marta- of the Added Value of Territorial Cohesion. Refereed Pérez-Soba, M., Tabbush (eds.): Sustainability Impact article No. 53, October, 2013, European Journal of Assessment of Land Use Changes. New York: Springer, Spatial Development, https://archive.nordregio.se/en https://doi.org/10.1007/978-3-540-78648 /News/A-Critical-Assessment-of-the-Added-Value-of -1_3. -Territorial-Cohesion/index.html. Zonneveld,17–33, W., Waterhout, B. (2009): EU Territorial Impact Perdicoúlis, A., Batista, L., Pinho, P. (2016): Logical chains Assessment: Under what Conditions? Final Report, in territorial impact assessment. Environmental Impact Delft University of Technology, Delft, http://citeseerx.ist https://doi.org /10.1016/j.eiar.2015.11.006. . Assessment Review 57, 46–52, .psu.edu/viewdoc/download?doi=10.1.1.615.647&rep =rep1&type=pdf Original Article 129
Influence of neotectonics on land surface evolution in the upper part of the Blue Nile Basin (Ethiopia): findings from a DEM Michal Kusák1,*, Vít Vilímek2, Jozef Minár3
1 Institute of Rock Structure and Mechanics of the Czech Academy of Sciences, Czech Republic 2 Charles University, Faculty of Science, Department of Physical Geography and Geoecology, Czech Republic 3 Comenius University in Bratislava, Faculty of Natural Sciences, Department of Physical Geography and Geoecology, Slovak Republic * Corresponding author: [email protected]
ABSTRACT The morphometric analysis of lineaments, valleys and signs of erosion taken from a digital elevation model (DEM) made it possible to not only confirm most of the conclusions of the morphotectonic development of the Blue Nile Basin from the previously pub- lished results of structural, petrological, tectonic and geochronological analyses, but also to expand our knowledge by applying several new hypotheses. The relative age of the morpholineaments of particular directions was estimated from the character of topographic profiles. Faults, lineaments and valleys are predominantly oriented in a direction compatible to the published concepts of the tectonic development of the area. Overall, the most abundant NE-SW and NNE-SSW lines reflect a change of extension from a NW-SE to WNW-ESE direction during the Pliocene, in relation to the creation and development of the Main Ethiopian Rift (MER). This is confirmed by a more developed character of the valleys and less pronounced erosion activity of the NE-SW oriented valleys contrary to the deeper narrower NNE-SSW valleys characterised by downward and headward erosion in the second direction. The most pronounced morphological manifestations of the E-W extension of the MER and western Afar during the Quaternary are confined to the borders of the MER. The directions of the Pre-Neogene rift structures to the NW-SE and WNW-ESE are compatible with the oldest elements of the current landscape and with the relict fragments of the valley network on the SE boundary of the upper Blue Nile Basin, which could have been drained across current shoulders of the MER to the S and E before the Late Miocene.
KEYWORDS landscape evolution; neotectonics; river piracy; lineaments; DEM; Main Ethiopian Rift; Ethiopian Highlands
Received: 17 April 2019 Accepted: 28 June 2019 Published online: 26 August 2019
Kusák, M., Vilímek, V., Minár, J. (2019): Influence of neotectonics on land surface evolution in the upper part of the Blue Nile Basin (Ethiopia): findings from a DEM. AUC Geographica 54(2), 129–151 https://doi.org/10.14712/23361980.2019.13
© 2019 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). 130 Michal Kusák, Vít Vilímek, Jozef Minár
1. Introduction the upper part of the Blue Nile Basin. Lineaments are linearly arranged elements of the GaniThe et al. main 2007, objectives 2009; ofWolela this work 2010). are We to: focused 1) analyse on the morpholineaments, faults and valley networks of the Blue Nile Basin and their interrelationships; landscape – lineara potential zone of brittle sections fractures of a valley, of bedrock ridges and/ or orstraight the concentration sections of slopes – should of erosion processes be considered with anas and the geomorphological evolution, in order to determine2) study the which relation elements between are thepotential morpholineaments zones of brit- - tle fractures of bedrock and serve as concentrations mentsinfluence can on give the an geomorphological insight into landscape evolution evolution, of an for erosion processes; 3) identify areas with the most thusarea (Abdullahproviding informationet al. 2010). onAn tectonicanalysis activityof the linea over dynamic changes of valley networks. large areas, which is particularly useful for areas with Several works deal with the tectonic and volca- nic history of the upper part of the Blue Nile Basin Highlands. Recently, the lineaments are increasingly usedlimited not field only access as a surface(Ehlen 2004),expression like theof individualEthiopian faults (fracture zones) but their statistics is used also (Abebe et al. 1998; Chorowicz et al. 1998; Wolfenden paidet al. to2004; the Pikrelated et al. geomorphic 1998, 2003; development. Kieffer et al. 2004;Local morphotectonicGani et al. 2009); and however, morphodynamic less attention aspects has beenhave to define/confirm the character (directions) of the mainly been investigated (Ayalew and Yamagishi palaeomorphotectonicIn this paper we aim stress to fields understand (Minár and relations Sládek 2008;between Koronovskya landforms andet al. the 2014; tectonic Šilhavý structures et al. 2016). of the - 2018).2003; Ismail Works and by Abdelsalam Gani and Abdelsalam 2012; Kusák (2006) et al. 2016; and - Mäerker et al 2016; Kycl et al. 2017; Gani and Neupne eneEthiopian and Abdelsalam Highlands, 2005; as the Gani region and Abdelsalam is strongly 2006;influ the selected aspects of geomorphic development. enced by tectonics (Kazmin 1975; Pik et al. 2003; Bey Gani et al. (2007) deal most comprehensively with
Fig. 1 Map of the Blue Nile basin and adjacent area, Ethiopian Highlands. Subregions particularly analyzed in the following text: 1 – the Jemma River; 2 – Beshlo Wenz River; 3 – Volcanoes Mt. Choke and Mt. Gish; 4 – Lake Tana; 5 – Didessa River and Baro River; A – Guder River and Muger River. Influence of neotectonics on land surface evolution 131
They refuse the older idea that present physiography landscape area at a lower elevation and with gentle of the elevated plateau exists since the Oligocene (e.g. valley slopes is situated in the western part of the of incision rate of the Blue Nile initiates ca. 10 Ma ago. with a low vertical division (Figure 1). ThisPik et al. paper 2003) presents and they evidence argue from that thesatellite rapid data increase that studyThe area modern and thetopography Blue Nile of flows the area through formed the Tertia basin- - ry and Quaternary volcanism and tectonics. The uplift ing morphotectonic activity of the territory, specifying supportdirections young and space (Late impactMiocene – Quaternary) of palaeomorphotectonic increas ledof the to Ethiopianthe formation Highlands of faults (during and cracksthe last (Kazmin 29 Ma) and opening of the MER (during the last 18 Ma) has fields on the land surface development. 2. The study area: geomorphological and1975; Abdelsalam Pik et al. 2003; 2012). Beyene and Abdelsalam 2006; and geological settings GaniThe et al. plateau 2007; in Gani the easternet al. 2009; part Wolela of the study2010; areaIsmail is built by various types of Cenozoic volcanic rocks and The study area includes the Blue Nile Basin between several Cenozoic volcanoes tower above it. Erosion Lake Tana and the Roseires water reservoir by the Blue Nile and its tributaries exposed Jurassic 2) and the adjacent area (Figure 1). The and Crateceous sedimentary complexes in the valleys Blue Nile Basin consists of a large plateau that is incised into the plateau. The western part of the study inclined(175,400 km from the border of Afar and the MER (more area is built by a Precambrian crystalline complex and volcanics (Figure 2; Mangesha et al. 1996). The Ethiopian Highlands are marked by a num- than 3,000 m a.s.l.) to the W and NW (approximately- ber of faults increasing from the west to the east, bolian2,000 m Quaternary a.s.l.). Huge volcanoes Neogene areshield concentrated volcanoes towerto the SSWup to of 2,000 m Lake Tana above reach the severalplateau, hundreds and smaller of metersstrom Rift (Chorowitz 2005). The oldest Permo-Trias and Crateceous-Palaeogenewhere they border the MER – part rift systems of in the the East study African area tectono-thermal uplift of the rift shoulders and a swell - centredat most (Kiefferon the Lake et al. Tana 2004). basin Pliocene are responsible and Quaternary for the gor 2015). The main phase of extension in the MER have a NW-SE orientation (Gani et al. 2009; MacGre Ismail and Abdelsalam 2012; MacGregor 2015; Gani and the majority of fault activity occurred in the last andrecent Neupne elevation 2018). of the Uplift plateau caused (Chorowicz the rivers et al. to cut 1998; into started approximately 11 Ma (Ukstins et al. 2002)- the bedrock and divide the area by canyons and deep gor 2015). The NW-SE Late Miocene extension formed valleys with steep slopes linked to faults in the east- 5–9 Ma, and particularly in the last 1–2 Ma (MacGre
NE(NNE)-trending faults (Wolfenden et al. 2004; Gani ern part of the study area (Gani et al. 2009). An older et al. 2009). E-W and NNE-SSW Quaternary extensions
Fig. 2 Geological map of the Blue Nile basin and adjacent area, Ethiopian Highlands (scale 1 : 250,000; after Mangesha et al. 1996). 132 Michal Kusák, Vít Vilímek, Jozef Minár opening the MER are related to the development of mega-lineamentsand lineaments) for an analysis of the E-trending transverse faults. N-ESE- and NW-trending main directions of the landscape forming processes. extensional structures are connected with the NE-SW The faults in the Ethiopian Highlands were taken from Quaternary extension in southern Afar and the E-W the geological map of Ethiopia (scale 1 : 250,000; west of Addis Ababa, an elongated Yerer-Tullu Wellel A valley network model for the Ethiopian High- volcano-tectonicextension in western lineament Afar (Gani is probably et al. 2009). an inherent To the landsMangesha was et al.extracted 1996). from SRTM DEM by following structure, which has decoupled the Northern Ethiopi- method of Jenson and Domingue (1988) and using an Plateau and the MER as a transtensional dominant- ArcHydro tools in ArcGIS 10.2.: Fill Flow Direction ly oblique E-W to ENE-WSW structure since the Late Flow accumulation (ESRI 2015), i.e. the valley net- -
Miocene (Abebe et al. 1998). sink-likework was artefacts identified caused by analysing by noise the and flow inaccuracies accumula 3. Materials and methods intion the after original a filling DEM function dataset. that Thewas thresholdapplied to removeof con- tributing area for the valley network generation was 3.1 Digital Elevation Model (DEM) and satellite 1,000 pixels (0.9 km2). For the analysis of the azimuth images of stream channels, we converted the valley network In this paper we use the SRTM DEM 1 Arc-Second Glob- raster into a vector format. The azimuths of megalineaments, lineaments, and September 2014. Shaded relief images with varying faults were determined in terms of the orientation directionsal version of(resolution illumination ≈ 30 were m), createdwhich was in ArcMap released 10.2 in of lines to the coordinate system and the azimuth of (ESRI 2015) using the hillshade tool (Spatial Analyst stream channels was also determined as the orienta- Toolbox Surface toolset). Various elevation angles tion of the stream channels (i.e. parts of valleys from of the illumination source (30°, 35°, 40° and 45°) and - four different values of illumination azimuth 1) from ence; or parts of valleys between two valley connec- tionthe valley points) heads to the to thecoordinate first point system. of the Thevalley azimuths conflu (solar azimuth of 45°), 3) from the east (solar azimuth were illustrated by rose diagrams, which are divided ofthe 90°), north and (solar 4) from azimuth the south-eastof 0°), 2) from (solar the azimuthnorth-east of into 72 intervals of 5° (360° in total). The numbers of 135°), were used in order to assure the independency lines in the rose diagrams were multiplied (weighted) of the results from the direction of illumination. This procedure maximizes the visualization of a surface for better observation and graphical display, to allow for by their length according to Belisario et al. (1999), the extraction of lineaments by using shaded-relief Ciotoli et al. (2003) and Kusák et al. (2016). - process by changing the illumination azimuth (Smith ed Moreover,by different several colours types (e.g. ofgreen, profiles red) were to determine created: the1) valley erosion cross-sections intensity and which(if possible) are in the figures relative depict age; - 3.2and Analysis Clark 2005; of linear Ahmed features et al. 2017). In the selected areas of the Ethiopian Highlands we MER2) long margin profiles to inreveal the landscape its block structures. to reveal the possibil mapped the linearly arranged elements of the land- ities of river piracy; 3) profiles perpendicularly to the scape, for example linear sections of valleys, ridges or straight sections of slopes. These linear landforms and 4. Results structures were manually mapped from: a) digital ele- vation model (hillshade images of various elevation 4.1 Aggregated characteristics angles of the illumination source and solar azimuths), In the eastern part of the Blue Nile Basin the old and b) satellite images of Google Earth (pixel size of 15 m, scale 1 : 100,000), and subsequent overlapping they recede in favour of expanding deep valleys. The of the layers (according to Wladis 1999; Ehlen 2004). fragmentsplateaus are of thestrongly plateaus eroded consist by offluvial Paleocene-Oligo erosion and- Many of these long linear structures are parallel cene-Miocene basalts and tuffs and deeply cut valleys to each other and the processes forming the current are also located there. In the western part of the Blue landscape are concentrated there. In order to deter- Nile Basin, Precambrian rocks are exposed (Figure 2). mine whether certain parallel linear structures are of A total of 161 mega-lineaments(with an aver- the same age and were created by the same process age length of 60 km), 1,264 faults (with an average (tectonic activity and erosion), the parallel mega-lin- length of 21 km), 3,429 lineaments (with an average eamentswere marked with the same colour (red, length of 10 km) and 81,137 valleys segments (with an green, blue, etc.). The shortest length of the mega-lin- average length of 2.3 km) were mapped and analysed in eamentswas 20 km. Shorter linear structures (length the upper part of the Blue Nile Basin. The results were < 20 km) are called lineaments and thousands of them then compared using azimuths and rose diagrams are seen on the Ethiopian Highlands. We mapped the (Figure 3). The analysis revealed that the mega-lin- linearly arranged elements of the landscape (the eaments(Figure 3A,B) and faults (Figure 3C, D) Influence of neotectonics on land surface evolution 133 are strongly correlated and concentrated into two orange lines), NNW-SSE (azimuth 170°; violet lines) directions: azimuth 20° (NNE-SSW orientation; red and N-S (azimuth 0°, black lines). Not all these direc- lines in Figure 3A) and azimuth 40° (NE-SW orienta- tion; green colour in Figure 3A). The third-most domi- nant direction being E-W (azimuth 90°) is clearly visi- tions are clearly reflected in aggregated faults but can have4.2 Local a fault scale reflection analysis at a local level. ble in the lineaments. It has not a significant reflection order to represent various types of landscapes in the ofin mappedthis direction. faults, Smallerhowever concentrations Abebe et al. (1998) of linear and The sub-regions 1–5 from Figure 1 were chosen in Gani et al. (2009) also mentioned important faults- (covering the area of 18,800 km2) and 2 (11,500 km2) tions: ENE-WSW (azimuth 80°; blue lines), WNW-ESE upperare on partthe boundary of the Blue of Nile the RiverEthiopian Basin. Highlands Sub-regions 1 close (azimuthfeatures were 115°; also brown identified lines), NW-SEin the following (azimuth direc145°; to the MER. Sub-region 3 (20,500 km2) is a typical area
Fig. 3 A) Map of megalineaments; B) rose diagram of megalineaments; C) map of faults; D) rose diagram of fault azimuths; E) map of lineaments; F) rose diagram of lineament azimuths. 134 Michal Kusák, Vít Vilímek, Jozef Minár adjacent to the central volcano, while sub-region 4 type of relief, and could represent an old structural (21,000 km2) represents the Lake Tana area. Sub-re- gion 5 (35,000 km2) in the southwest represents evolution in the Ethiopian Highlands. a highly-eroded part (i.e. the lowest part) of the Ethi- predisposition, which influenced the valley network opian Highlands. 4.2.1 The Jemma River It is evident that the NE-SW direction (green lines) The local scale analysis of sub-region 1 in the catch- is important in the sense of the valley evolution. This ment of the Jemma River revealed that apart from the is true across all of the sub-regions, regardless of the importance of the NE-SW orientation (azimuth 40°;
Fig. 4 The sub-region 1 (Jemma River Basin sub-region): A) map of the valley network and valley border; B) map of megalineaments, faults and lineaments; C) rose diagram of valley network azimuths; D) rose diagram of mega-lineamentsazimuths; E) rose diagram of fault azimuths; F) rose diagram of lineament azimuths. Influence of neotectonics on land surface evolution 135
Fig. 5 Profiles in sub-region 1 (Jemma River Basin). Profiles 1–21 are valley cross-sections across green lines and profiles 22–33 are valley cross-sections across red lines; A – longitudinal profile No. 34 along a single red line (a site with a tendency for river piracy). 136 Michal Kusák, Vít Vilímek, Jozef Minár green line), the NNE-SSW orientation (azimuth 20°; diagrams (Figure 4). The shapes of the red line val- ley cross-sections suggest more intensive downward network evolution and lineament depiction. The rose erosion compared to the green line valley cross-sec- diagramred line) forplays the a faults more suggestssignificant tectonic role in predisposi the valley- tions, which suggests that the red line direction is tion. Both directions (green and red) have a similar younger (Figure 5). One site with a tendency to river azimuth (40° and 20°, respectively); however, they are sharply bounded and distinguished in the rose red lines (Figure 5A). The headward erosion is strong piracy was also positively identified along one of the
Fig. 6 Sub-region 2 (Beshlo Wenz River Basin sub-region): A) map of the valley network and valley border; B) map of megalineaments, faults and lineaments; C) rose diagram of valley network azimuths; D) rose diagram of mega-lineamentsazimuths; E) rose diagram of fault azimuths; F) rose diagram of lineament azimuths. Influence of neotectonics on land surface evolution 137
Fig. 7 Profiles in the Beshlo Wenz River Basin sub-region. Note: Profiles 1–10 are cross-sectional profiles directed across the valley marked by green lines; and 11–23 are cross-sectional profiles directed across the valley marked by red lines. Profiles 24–35 are cross-sectional profiles directed across the valley marked by an orange line. 138 Michal Kusák, Vít Vilímek, Jozef Minár from both sides of the two different valleys and river is well pronounced (Figure 6). The Afar area adjacent piracy can be expected to happen soon because the to the Blue Nile Basin, is also included in this sub-re- remnant of the original relief created by the basalt gion. The Rift has a NNW-SSE direction herein. Bor- layer is only 1.3 km wide. It is clear that the shape of the cross-sectional pro- and megalineaments. The rose diagram for the valley - networkder faults is (violet rather colour) complex are because reflected the in western lineaments part sitions), including intensity and type of erosion, dis- is similar to sub-region 1 (Jemma River Basin), while file is influenced by several processes (or predispo will reveal the general trend in valley evolution (if the facing the MER. lithologycharge and is lithology;homogeneous), nevertheless, which canthe setbe comparedof profiles theIn eastern accordance part bears with significantsub-region features 1, the valleys of the inslope the to other sub-regions. Valley cross-sections along the Beshlo Wenz catchment lying along green lines, are again more open and wide, while those running along the red lines are markedly incised (Figure 7). The only green lines (profiles 1–21 in Figure 5) are wider and narrowmore open and comparedtighter. This to is the regardless profiles ofalong the positionthe red to the rift, Figure 7), which are also very narrow like lines (especially profiles 26–33), which are more thoseexceptions connected are profiles to the red 3 through lines. The 6, andarea 10next (closest to the of one direction are similar in the upper, middle and border of the MER is under very intensive young evo- lowerof the partsvalley of in the the valley river system.network system i.e. profiles lution (strong downward erosion) regardless of the type of valley (red or green lines). 4.2.2 Beshlo Wenz River Several tectonic blocks evolved on the east-fac- The valley network predisposition of the landscape ing slope of the rift, which are asymmetrical in their evolution in the Beshlo Wenz Basin (sub-region 2) is similar to the Jemma River Basin. However, the faults do not have an azimuth of 20°, but the NNW-SSE ori- strongcross-sectional headward profile erosion and retainalong thesmall steep, lakes east-facing (Figure 8). entation (azimuth 170°; violet colour on the graphs) The drainage in a NNW-SSE orientation is modified by
fault slope. In our fieldwork, we identified the case of
Fig. 8 Cross-sectional profiles oriented towards the MER. Tectonic blocks, which are asymmetrical in the cross-sectional profiles evolved here and retain Lake Hayk (A), Lake Hardibo (B) and Lake Mybar (C). Influence of neotectonics on land surface evolution 139
by tectonics with an azimuth of 170° (NNW-SSE east nowadays because a small river reached the lake orientation) depicted by a violet colour in Figure Lakeshore Mybar due to very(Figure intensive 8C) – the headward lake is erosiondrained and to the 6 and 8, while in the surrounding of the adjacent former drainage to SSE was abandoned. Lake Hayk (Figure 8A) both NNE-SSW and NNW-SSE In this locality the tectonic predisposition is rep- resented by faults with an azimuth of 20°. This is also respectively). the area with highest frequency of faults. Neverthe- trending faults were identified (red and violet lines, less, 30 km to the north the tectonic predisposition 4.2.3 Volcanoes Mt. Choke and Mt. Gish is changing in connection with the main structures The situation around the volcanoes of Mt. Choke of the rift. Lake Hardibo (Figure 8B) is predisposed (4,100 m a.s.l.) and Mt. Gish (2,890 m a.s.l.) is rather
Fig. 9 Sub-region 3 (Mt. Choke and Mt. Gish volcano sub-region): A) map of the valley network and valley border; B) map of mega- lineaments, faults and lineaments; C) rose diagram of valley network azimuths; D) rose diagram of mega-lineaments azimuths; E) rose diagram of fault azimuths; F) rose diagram of lineament azimuths. 140 Michal Kusák, Vít Vilímek, Jozef Minár
the Ethiopian Highlands due to the NE-SW trending volcanoes in this sub-region (radial network with preferencessimple. The towardsvalley system certain reflects directions); the position nevertheless, of the (Figure 10). Position on the boundary between the all other rose diagrams show NE-SW trending faults, mega-lineaments identified in all of the sub-regions - ence of the younger NNE-SSW tectonics (Figure 10). parallelLT and HT to provinces this NE-SW of flood oriented basalts valley, (Pik et al.a younger, 1998; lineamentsThe long andvalley mega-lineamentswithout system between the volcanoes any influ of deeperKieffer et al.and narrower 2004) could valley also has play been a role.created However, on the Mt. Choke and Mt. Gish is probably a part of an old- lower slopes of Mt. Choke, which were dated by Kief- er river network, which is pronounced in this part of
fer et al. (2004) to an age of 22.4 Ma. This direction
Fig. 10 Profiles in the Mt. Choke and Mt. Gish volcano sub-region. Note: Profile 1 along the single green line (an area with a tendency for river piracy); Profiles 2–15 are cross-sectional profiles directed across the valley and marked by green lines; Profiles 16–22 are cross-sectional profiles directed across the valley and marked by brown lines. Influence of neotectonics on land surface evolution 141
Fig. 11 Sub-region 4 (Lake Tana sub-region): A) map of the valley network and valley border; B) map of megalineaments, faults and lineaments; C) rose diagram of valley network azimuths; D) rose diagram of mega-lineamentsazimuths; E) rose diagram of fault azimuths (E); rose diagram of lineament azimuths (F). 142 Michal Kusák, Vít Vilímek, Jozef Minár seems to be dominant in the area of Mt. Choke and Another example of a tendency for river piracy was Mt. Gish, even though the river network should have a radial character around the volcanos. If we consid- original landscape was left before the rivers from the er the fact that the NE-SW valleys (green lines) also identified here (Figure 10A). Less than 1.5 km of the cut 10.7 Ma old basalts, then they have to be younger than the basalts itself. SW or NE came across the old lava flows. Here we can presume expansion of a SW valley due to the profile (Fig. 10, profile 1).
Fig. 12 Longitudinal profiles in the Lake Tana sub-region: Longitudinal profiles 1 and 2 show the situation for rivers stretching WNW-ESE; Longitudinal profile 3 shows the situation for rivers stretching SW-NE. Influence of neotectonics on land surface evolution 143
Fig. 13 Sub-region 5 (Didessa River and Baro River): A) map of the valley network and valley border; B) map of megalineaments, faults and lineaments; C) rose diagram of valley network azimuths; D) rose diagram of mega-lineamentsazimuths; E) rose diagram of fault azimuths (E); rose diagram of lineament azimuths (F). 144 Michal Kusák, Vít Vilímek, Jozef Minár
Fig. 14 Profiles in the Didessa River and Baro River sub-region. Note: Profiles 1–4 are cross-sectional profiles directed across the valley marked by a brown line; and profiles 5–8 are cross-sectional profiles directed across the valley marked by an orange line. Profiles 9–16 are cross-sectional profiles directed across the valley marked by a yellow line. Profiles 17–24 are cross-sectional profiles directed across the valley and ridge marked by a red line. Influence of neotectonics on land surface evolution 145
4.2.4 Lake Tana lake; however, the headward erosion has to remove - by tectonics, namely in a NNE-SSW direction (azimuth - ofThe 20°). valley The network valleys seems trending to be WNW-ESE only partly (azimuth influenced of uationa 200 m in highthe SW crest of the (escarpment) lake shows, riverforming piracy the taking lake 115°) do not show any connection to the mapped placeshore inin athe SW west to NE (Sembroni direction et al. (along 2016). the greenAnother lines) sit the Lake Tana Basin. Nevertheless, both directions are long. Headward erosion along this green line is more markedfaults and because are significant they represent for the a serious area to potential the west for of probableand the remnant, because which it is predisposedmust be crossed, by the is only fault. 3.8 km The - recent erosional activity can be shown in a NNE-SSW directionriver piracy (red (see lines), the paragraph like in other below). sub-regions No significant closer (2,890landscape m a.s.l.). on the lakeshore is also flatter and is creat ed Currently,by a 150-m the high area remnant of the of Blue lava flowNile frombasin Mt. covers Gish is not well pronounced in this part of the Ethiopian 175,400 km2 (between Lake Tana and the Roseires Highlands.to the MER. Nevertheless, It seems that NNE-SSW the influence faults ofexist the in MER this water reservoir). However, when Lake Tana will change the drainage to the SW, the cut off part of the megalineaments. Blue Nile Basin (98,600 km2) will not be supplied by area,The and whole they of influence the Lake both Tana the Basin, lineaments which is and tecton the- water from Lake Tana, but only by precipitation, which is not evenly distributed during the year. The annual drains into the Blue Nile to the SE. Rivers on the west- ernically margin predisposed of the basin (e.g. areChorowicz approaching et al. the1998), lake now due 80% of this precipitation falls from July to September to headward erosion from two different directions (Klimadiagrammerainfall in the Ethiopian weltweit Highlands 2016). is 1,800 mm and show the situation for rivers stretching WNW-ESE 4.2.5 Didessa River and Baro River (brown(Figure 12).lines Longitudinal on the rose diagram profiles in1 andFigure 2 in 11). Fig. The 12 Many of the valleys (Figure 13) follow either a NW-SE headwaters are approximately 10 km away from the or NE-SW direction where the valley bottoms are
Fig. 15 Earthquake epicentres in the Jemma River Basin (A) and in the Beshlo Wenz Basin (B) (Earthquake Track 2017). 146 Michal Kusák, Vít Vilímek, Jozef Minár wider, which suggests an older river network. If we (NNW-SSE orientation) depicted by a violet colour. across brown and orange lines (azimuth 115° and catchment. The average depth of the hypocentres is 145°,check respectively),the cross-sectional the shapes profiles resemble (1–8) that wide stretch and 20This km, is duewhich to indicatesa logical influencerather shallow of the earthquakesMER in this open valleys (Figure 14). In addition, they are rath- and the average magnitude is 5.5. er deep, which suggests a long period of evolution. A slightly different situation can be observed on the 4.4 Features of Guder River and Muger River - network reorganisation sion is already visible on the valley bottoms. An oppo- Features of drainage network reorganization along profiles crossing the yellow lines, whereby young ero lineaments of various types can be detected in var- ious parts of the Blue Nile Basin. The drainage net- Fig.site 14),situation which to are these again valley very profilesdeep, young can andbe found narrow in works of the Guder and Muger catchments (Figure (similarthe valleys to Figuredepicted 7). by The the area red of lines the (profilesDidessa and21–24 Baro in 1A, Figure 16) are organised dominantly along diago- rivers seems to be older in the sense of the evolution nal green and orange lines; however, the most active of the valley network with young erosion in the direc- recent river incision is along axial lines (yellow, violet, tion of the red lines (azimuth 20°). black and red, see Figure 16). This activity is evident mainly in lower parts of both watersheds as well as in 4.3 Analysis of earthquake epicentres the adjacent part of the Blue Nile valley. On the oth- A total of 12 earthquakes have occurred in and around er hand, valleys along diagonal lines are wider, more the Jemma River Basin since 1961 (Figure 15A; Earth- denuded and are partly abandoned. A set of highly ele- quake Track 2017). Most of the epicentres relate to “red lines” (eight events), only three to “green lines”. In the Beshlo Wenz Basin (Figure 15B) the epicen- onvated the depressions southern watershed (2,600–3,000 to the m NWa.s.l.) of elongated Addis Aba in- tres are located on faults with an azimuth of 170° a NNW-SSE direction (orange lines) can be identified
ba. Their very flat and wide bottoms (too large for the
Fig. 16 The Sub-region A (Guder River and Muger River basin sub-region): A) Guder River and Muger River networks organised dominantly along diagonal green and orange lines; B) Profiles of Guga river (1), terrace (2) and plateau (3); C) lineaments of recent drainage network, lineaments indicating paleo-drainage network and other lineaments. Influence of neotectonics on land surface evolution 147 recent small streams) with differing descent (partly to 5. Discussion the S and partly to the N) suggest that they belong to an ancient drainage network that might have drained Landscape analysis from SRTM DEM data is of great the area towards the SE. More detailed features of river network reorgan- et al. 2015). SRTM DEM data are often used for three typesinterest of to morphometric geoscientists (e.g. analysis: Martino discrete et al. 2009; elevation Haider Guder River (Figure 16B). The river recently rose function f onisation the bordercan be seenscarp in of the the longitudinal E-W oriented profile Yerer-Tul of the- - lu Wellel lineament, below one of the above-men- ful for performing (x, y, z), hill a morphostructural shading as well as analysis flow direction of the tioned paleo-valleys at an elevation of approximately upperand flow Blue accumulation. Nile Basin and All its three surroundings. were also very use 3,000 m a.s.l., where only several hundred horizontal Hill shading offers many possibilities for interpret- meters are missing in order to capture the stream of ing detailed surface structures that are often com- the paleo-valley by headward erosion. Leaving the pletely imperceptible using simple DEM visualization (e.g. Kennelly 2008). Hill shading allowed us to per- lineament without any notable incision. At an ele- form a much more accurate delineation of lineaments vationscarp, theof ~1,800 river flows m a.s.l. along it turns the toN theboundary NW, where of the it and borders of erosion than mapping from certain stepped canyon with a small decrease in gradient. The the overlapping of layers (the digital elevation mod- flows across a denuded border scarp in a 0.5 km deep el;satellite hillshade images images only of (Kusákvarious et al.elevation 2016). angles However, of the lineament to the east of Lake Chomen and then before illumination source and the solar azimuths; satellite profile only becomes steeper when it crosses the W-E- images with a pixel size of 15 m placed on the virtu- mately ~950 m. Along with the typical character of al globe of Google Earth) proved to be the most use- it flows into the Blue Nile at an elevation of approxi- ful when checking the manually mapped landscape shapes (following Wladis 1999; Ehlen 2004; Gani and the past.upper The part general of the trendprofile, of severalthe plateau other inclination morpho Abdelsalam 2006). Morphometric analysis based on (fromlogical N features to S) has point been to reversed reverse byflow a strong of the upliftriver inof a calculation of the potential drainage of each pixel the border of the Yerer-Tullu Wellel lineament (Figure of an SRTM DEM shows the shape, density and hier- 16); therefore, natural drainage to the S could have existed in the past. The strongly eroded, wide and of structure and erosion processes on the formation branched middle part of the watershed contrasts with ofarchy the oflandscape valley networks, (sensu which reflect the influence the narrow lower part closed between well preserved 2016), similarly to an evaluation of various topo- remnants of the volcanic plateau in the north. On the Stoddart 1997; Kusák et al. contrary, the uplifted plateau is markedly denuded in effectively used in the central part of the area by Gani the south. andgraphic Abdelsalam profiles. (2006) A set ofbut all our of workthese brings analyses a differ was- Probably the most interesting feature is the char- ent type of information. acter of a wide terrace (2 in Figure 16B), which is The character of faults, their direction and space directly incised by the river. The elevation of the differentiation as well as time of formation and reac- terrace rises from the mouth (~1,350 m a.s.l.) only tivation, are important for the reconstruction of the approximately 20 km (~1,450 m a.s.l.) to the nar- morphotectonic development. However, direct inves- rowest part of the valley bounded by remnants of tigation of faults is limited to outcrops, boreholes the plateau. Subsequently (from 20 km to 70 km), the terrace is reversed to recent river as for W-E ele- Therefore, they are not mapped equally in detail on vated belt at Lake Chomen lineament. Nevertheless, geologicaland geophysical maps profiles,and the distinctionwhich is time of demanding.their origin a small paleo-valley on the southern border sloping and activation is limited. On the other hand, linea-
Further to the S, the river is not incised and prob- uniformly mapped with a certain probability on the to the S is also an indication of reverse flow here. basisments of reflecting a DEM. Mega-lineaments the geological fractures are most can probably only be the point where a general rise in the elevation of the connected with important faults, but because of their plateauably flows begins down (Figure to the 16).terrace The (~1,300former waterm a.s.l.) shed to limited number they are not as suitable for statistical divide, approximately 20 km from the recent mouth, analysis as regular lineaments. The stages of develop- ment of structurally determined landforms can also can be inferred. Considering all other indices, drain- - ageand towardsat least athe 50 SSE, km mainlylong reversed along diagonal flow of lines,the river can ence of morphodynamic phenomena (erosion, land- be assumed before the uplift of the border of the Yer- slides,be estimated earthquakes) from topographic (Gani and Neupne profiles 2018). and the Finally, pres er-Tullu Wellel lineament. Therefore, reorganization the combined information from all of these sources of the drainage network in relation to the current sit- can be compared to the results of various geological uation took place later by headward erosion, mainly and geochronological analyses to obtain a compre- along the axial lines. hensive picture of the morphotectonic development. 148 Michal Kusák, Vít Vilímek, Jozef Minár
Two basic groups of structural features can be dis- Pliocene tectono-thermal uplift of the MER shoulders tinguished based on our analysis of the upper part of and relative subsidence of the Tana basin (Chorowicz the Blue Nile basin: a) Axial (including black, yellow, violet, blue and red that in this initial stage the upper part of the Blue lines) around N-S and E-W directions. Nileet al. could 1998; have Sembroni also drained et al. 2016). along the We ancient can speculate NW-SE b) Diagonal (including green, orange and brown lines) (orange) or WNW-ESE (brown) structures to the area around NE-SW and NW-SE directions. of the recent Afar depression. An indication of this may The difference between the red and green lines be the mature erosional pattern of the E and SE part (representatives of the basic groups mostly con- of the Blue Nile Basin in comparison with the more tained in the faults and lineaments) is associated homogeneous (younger) erosional pattern of the W with a transition from a 130°E-oriented extension and SW parts (Gani and Abdelsalam, 2006) and also to a 105°E-oriented extension of the MER sometime indications of diagonally-elongated palaeovalleys in between 6.6 to 3 Ma ago, which caused a change in the recent watershed around the city of Dessie. the normal fault orientations from predominantly The spatial distribution of the particular types - of the studied linear features can be linked with the time-space character of the background tectonic N35°EBlue Nile to N10°EBasin, (Wolfendeni.e. its formation et al. 2004).as a NW-trending The diago events. Sub-region 1 lies closely westward of the MER, nal features also reflect yet older development of the and sub-region 2 directly on the border of the Afar mainly respond to Quaternary development, which depression. Sub-regions 3 and 4 are far from the rift rift (Gani et al. 2009; MacGregor 2015). Axial features- boundary and the eastern part of sub-region 5 is close sions related to an oblique opening of the MER and to the Rift (Figure 1). The investigated (coloured) developmentis influenced ofby E-trending E-W and NNE-SSW-orientedtransverse faults as wellexten as directions represent only 44% of the possible drain- age courses. Nevertheless, 62% (sub-region 4) to 87% (sub-region 2) of the drainage agrees with the anal- E-W-oriented extensions in western Afar (Gani et al. 2009).The Landslidesyoung Quaternary linked with character axial directions of the majority confirm Quaternary activity too (Kycl et al. 2017). ysed directions. This not only confirms the significant- They are characterized by prevailing deep and narrow aryinfluence of the ofMER tectonics and in particularon the drainage the Afar network, depression. but valleys,of the axial fault lines scarps, was earthquake also confirmed activity by andour aanalyses. tenden- Axialalso an directions increase dominatein this influence in sub-regions around thenear bound or on cy for recent river piracy. On the contrary, the diag- the rift boundary (sub-regions 1, 2 and 5), whereas onal lines are generally linked with wider and more open valleys with a higher tendency for meandering the last two sub-regions. or braided channels, a dominant abundance of higher diagonalThe brown directions direction are comparably is mainly about significant the SW only part in of sub-regions 3 and 5 with the most mature relief features. Reorganisation of the drainage networks of and exposed pre-Tertiary bedrock, which suggests its theorder Guder streams and Mugerand significant rivers from “well-developed” the diagonal to relief axial ancient character. The orange direction has a similar directions supports this dichotomy. character, but its important representation in sub-re- The development of lineaments and river networks gion 2 and the eastern part of sub-region 5 suggest is generally accelerated by the uplift of the area. The initial Paleogene surface of the volcanic plateau of the green direction is generally the most represented in Ethiopian Highlands probably formed at elevation sig- allneotectonic of the linear reactivation features, mainly (Abebe in et al. sub-regions 1998). The1, 3 - and 4. Transformation of tectonic failures into the ing until the Late Miocene (MacGregor 2015). Indeed, drainage network indicates the oldest neotectonic nificantly lower than now, and then had been decreas The red lines are second in the order of areal Pike et al. (2003) suggest that erosion initiated in the importanceevent of regional (after significance, the green lines) i.e. opening and are of theconnect MER.- Blueof the Nile Blue canyon Nile up as to early the Late as 25–29 Miocene, Ma whereas ago, but rapid Gani - et al.erosion (2009) would evaluated have started only aonly very 6 gentle Ma ago. incision Generally, rate nounced in the inner part of the Ethiopian Highlands. the neotectonic uplift of the territory can be connected Theed with black the and stress violet field lines of the have MER, only but a limitedare less areal pro - (their higher importance was revealed only in sub-re- andwith Abdelsalamthe opening 2012)of the MERbut uplift (11 – 6 of Mathe ago Lake accord Tana gioninfluence 2), while on the they adjacent are the parts youngest of the andMER have and Affarbeen sub-regioning to various probably authors – c.f. began earlier, Gani et al. i.e. before 2009; the Ismail end active for the short period of time. The yellow lines - are located in the south of the upper part of the Blue nospheric mantle intruded the lithosphere as a result Nile Basin and they are primarily a result of activation of plumethe mid-Tertiary activity, inducing flood volcanism thermal uplift when (Chorowicz the asthe of the Yerer-Tullu Wellel lineament in the Late Mio- cene. Nevertheless, recessing the lower parts of their - tionet al. 4.4), 1998). which Further was lateran initial blocked radial by drainageLate Miocene network and activity increases in an E-W direction. could have formed including outflow to the S (see sec cross profiles means that the Quaternary tectonic Influence of neotectonics on land surface evolution 149
The Ethiopian Highlands are one of the most tec- The beginning of the neotectonic stage (during which the modern topography and drainage network formed) can be connected to the opening of the MER tonicallyerosion and influenced river piracy areas are in thevery world active (Beyene there. Theand during the Late Miocene. The most wide-spread NE-SW Abdelsalamvalley networks 2005; are Gani young et al. and 2007, water 2009), divides so change river (green) linear features parallel to MER represent the position easily. Thanks to the steeper slopes and a more humid climate inside the MER, the water extension in the whole territory. Reactivation of some courses (tributaries of the Awash River) have more NW-SElong term (orange) morphotectonic linear features influence mainly near of an the NW-SE MER erosive power and the Jemma and Beshlo Wenz riv- boundary match with activity of the offset faults. Devel- er basins lose parts of their catchments as the water oped (wider, more open) valleys and more ‘mature relief’ are typical for the leading green direction. Pliocene counter clockwise rotation of the MER ofdivide river is piracy moving would slightly probably towards be seenthe west in the (Kusák Lake Tanaet al., Basin. 2016). However, the most significant effect the generation of NNE-SSW (red) linear features. The Currently, the Lake Tana Basin is in the focus of stress field led (mainly in the MER surrounding) to uplift of the whole territory, especially the rift shoul- - stress field transition was linked also with intensive many scientific disciplines in terms of: 1) geology toward S or E. Continuing relative subsidence of the 2014);and geochronology 3) chemical (Praveanalysis et al. of water2016); and 2) reproducsediment Lakeders Tana that definitelybasin, headward blocked erosion any existing and river outflow piracy ing historic sediment concentrations (Kabaa et al. could contribute to spreading of the Blue Nile basin - too. Morphological features such as incised narrow (Ahrens et al. 2016); 4) assessing the implications of valleys, indices of river piracy, valley asymmetry and ofwater the aboveharvesting stated intensification research activities on upstream–down and many other the generally ‘young age relief’ characteristic for the authorsstream ecosystem write about services the Lake (Dile Tana et al. Basin 2016); as aetc.. single All red direction support this interpretation. unit, but do not include the processes taking place in N-S (black and violet) linear features represent the the surroundings. River piracy approaches the Lake youngest Quaternary development of the territory Tana Basin from SW to NE (along the green lines) and connected with E-W extension in MER and western Afar. They are morphologically least developed and long (Figure 12). The potential change of drainage of represented by huge fault slopes with tectonic blocks Lakethe remnant, Tana to which the SW has (instead to be surpassed, of to NE) is willonly greatly 3.8 km and active youngest upper parts of valleys. Orthogo- nal features with E-W direction are dominant along large area. older (Late Miocene) Yerer-Tullu Wellel volcanotec- influence the hydrological conditions over a rather tonic lineament, however they show evident features
6. Conclusion to river piracy. Accelerating Quaternary uplift creat- edof recentrecent reactivationvery high relief in cross mainly profiles along andthe boundarytendency The results of simple morphostructural analysis from with MER and Afar depression and so conditions for a DEM can be interpreted in high concordance with intensive headward erosion from the outside the Blue the present results of stratigraphic, structural, petro- Nile basin that could lead to its reduction by river logic, tectonic and geochronologic investigations of piracy in the near future. the Ethiopian Highlands. The azimuths of lineaments, faults and river networks combined with erosional features, cross-sectional and longitudinal valley pro- Acknowledgements stages of topography development. The study was supported by the GACR project “Mass filesThe reveal NW-SE a change (orange in thelines) dominant and WNW-ESE stress fields (brown and wasting and erosion as an indicator of morphotecton- lines) linear features correspond with older stress ic activity in the Ethiopian Highlands based on remote - sensing approaches” (P209/12/J068) and by the gene rifts and transposed onto the volcanic plateau Grant Agency of Charles University in Prague, Czech fields organizing the formation of ancient pre-Neo Republic “Morphostructural analysis of Ethiopian determined a palaeo-drainage system that could par- Highlands based on remote sensing” (1436314). Jozef during the first phases of its rise. They probably shoulders and are well pronounced in sub-regions Development Agency under the contract No. APVV-15 withtially uncoveredflow to the southpre-volcanic and east bedrock. across the These recent direc rift- Minár-0054. The was SRTM supported DEM was by the provided Slovak by Research the U.S. Geo and- tions are characterized by extremely braided valleys logical Survey. This work was carried out thanks to with locally brained channels, small remnants of vol- support of the long-term conceptual development canic plateau and features of well-developed relief’. research organization RVO: 67985891. We wish to They are typical for the oldest part of Blue Nile and thank Craig M. Hampson and Frederick Rooks for lan- Didesa valley. guage editing. 150 Michal Kusák, Vít Vilímek, Jozef Minár
References ESRI Inc. ArcMap™. Ver. 10.1. Copyright 2015. [online]. [downloaded 01/01/2015]. www.esri.com. Abdullah, A., Akhir, J.M., Abdullah, I. (2010): Automatic Gani, N., D., Abdelsalam, M. G. (2006): Remote sensing Mapping of Lineaments Using Shaded Relief Images analysis of the Gorge of the Nile, Ethiopia with emphasis Derived from Digital Elevation Model (DEMs) in the on Dejen-Gohatsion region. Journal of African Earth https://doi.org/10.1016 /j.jafrearsci.2005.10.007. Abebe,Maran – Sungi T., Mazzarini, Lembing F., Innocenti, Area. Malaysia. F., Manetti, EJGE P. (1998):15, Gani,Science N. D., 44,Gani, 145–150, M. R., Abdelsalam, M. G. (2007): Blue Nile The949–957. Yerer-Tullu Wellel volcanotectonic lineament; incision on the Ethiopian Plateau: Pulsed plateau growth, a transtensional structure in central Ethiopia and the Pliocene uplift, and hominin evolution. GSA Today 17, associated magmatic aktivity. Journal of African Earth https://doi.org/10.1130/GSAT01709A.1. https://doi.org/10.1016 Gani, N. D., Abdelsalam, M. G., Gera, S., Gani, M. R. (2009): /S0899-5362(97)001413. Stratigraphic4–11, and structural evolution of the Blue Nile Ahrens,Sciences L., Gashaw, 26(1), 135–150, H., Sjoholm, M., Gebrehiwot S. G., Basin, Northweastern Ethiopian Plateau. Geological Getahun, A., Derbe, E., Bishop, K., Åkerblom, S. (2016): https://doi.org/10.1002 /gj.1127. Gani,Journal N. D., 44(1),Neupane, 30–56, P. C. (2018): Understanding transient Tana,Poly- andEthiopia perfluoroalkylated and implications substances for human (PFASs) exposure. in landscape of the Ethiopian Plateau in relation to mantle water, sediment and fish musclehttps://doi.org/10.1016 tissue from Lake https:// /j.chemosphere.2016.09.007. doi.org/10.1002/gj.2903. Ayalew,Chemosphere L., Yamagishi, 165, H.352–357, (2003): Slope failures in the Googledynamics. Earth, Geological Google Inc. Journal Available 53(1), from 371–385, URL. [online]. Blue Nile basin, as seen from landscape evolution [downloaded 01/01/2015]. http://www.google.cz/intl https:// /cs/earth doi.org/10.1016/S0169-555X(03)00085-0. Ayman,perspective. A. A., Abdelkareem, Geomorphology M., Asran, 57(1–2), A. M., 95–116, Mahran, T. M. (2017): Geomorphic and lithologic characteristics of Haider,multi-parameter V. L., Kropáček, assessment J., Dunkl, ofI., Wagner,digital elevation B., and Von models. Wadi Feiran basin, southern Sinai, Egypt, using remote EarthEynatten, Surface H. (2015): Processes Identification and Landforms of peneplains 40(11), by https://doi.org/10.1002/esp.3729. https://doi.org/10.1007 Ismail, E. H., Abdelsalam, M. G. (2012): Morpho-tectonic /s12040-017-0861-8.sensing and field investigations. Journal of Earth System analysis1477–1492, of the Tekeze River and the Blue Nile drainage Belisario,Science F., 126(85), Del Monte, 1–25, M., Fredi, P., Funiciello, R., Palmieri, systems on the Northwestern Plateau, Ethiopia. Journal E. L., Salvini, F. (1999): Azimuthal analysis of stream https://doi.org /10.1016/j.jafrearsci.2012.04.005. Jensonof African S. K., Domingue, Earth Sciences J. O. (1988):69, 34–47, Extracting Topographic Beyene,orientations A., Abdelsalam, to define M. regional G. (2005): tectonic Tectonics lines. of Zeitschrift the Afar Structure from Digital Elevation Data for Geographic Depression:für Geomorphologie A review Supplementary and synthesis. JournalIssues 118, of African 41–63. Information System Analysis. Photogrammetric https://doi.org /10.1016/j.jafrearsci.2005.03.003. http://citeseerx.ist.psu.edu/viewdoc/download Chorowitz,Earth Sciences J. (2005): 41(1–2), The East 41–59, African rift system. Journal of ?doi=10.1.1.138.6487andrep=rep1andtype=pdf.Engineering and Remote Sensing 54(11), 1593–1600, https://doi Kaba, E., Philpot, W., Steenhuis, T. (2014): Evaluating .org/10.1016/j.jafrearsci.2005.07.019. suitability of MODIS-Terra images for reproducing Chorowicz,African Earth J., Collet, Sciences B., Bonavia, 43(1–3), F. F.,379–410, Mohr, P., Parrot, J. F., historic sediment concentrations in water bodies: Lake Korme, T. (1998): The Tana basin, Ethiopia: Intra-plateau Tana, Ethiopia. International Journal of Applied Earth
https://doi.org/10.1016/S0040-1951(98) doi.org/10.1016/j.jag.2013.08.001. 00128-0.uplift, rifting and subsidence. Tectonophysics 295(3–4), Kazmin,Observation V. (1975): and Geological Geoinformation Map of 26, Ethiopia, 286–297, Geological https:// CiotolI,351–367, G., Seta, M. D., Del Monte, M., Fredi., P., Lombardi, Survey of Ethiopia, Addis Ababa, Ethiopia. S., Palmieri, E. L., Pugliese, F. (2003): Morphological and Kennelly, P. (2008): Terrain maps displaying hill-shading geochemical evidence of neotectonics in the volcanic area of Monti Vulsini (Latium, Italy). Quaternary https://doi.org/10.1016/j.geomorph.2008.05.046. https://doi.org Kieffer,with B.,curvature. Arndt, N., Geomorphology Lapierre, H., Bastien, 102(3–4), F., Bosch, 567–577, D., /10.1016/S1040-6182(02)00093-9. Pecher, A., Yirgu, G., Ayalew, D., Weis, D., Jerram, D. A., Dile,International Y., T., Karlberg, 101–102, L., Daggupati, 103–113, P., Srinivasan, R., Wiberg, Keller, F., Meugniot, C. (2004): Flood and shield basalts D., Rockströmd, R. (2016): Assessing the implications from Ethiopia: Magmas from the African Superswell. https://doi.org downstream ecosystem services: A case study in the /10.1093/petrology/egg112. Lakeof water Tana harvesting basin. Science intensification of the Total on Environment upstream- KlimadiagrammeJournal of Petrology weltweit. 45(4), [online]. 793–834, [downloaded https://doi.org/10.1016 10/10/2016]. www.klimadiagramme.de. /j.scitotenv.2015.10.065. Koronovskya, N. V., Bryantsevaa G. V., Goncharova M. A., 542(part A), 22–35, Naimarka A. A, Kopaev A. V. (2014): Lineaments, Encyclopaedia of geomorphology. Routledge, London, Planetary Jointing, and the Regmatic System: Main Points Ehlen, J. (2004). Lineation. In Goudie, A. S. (ed.) et al.: of the Phenomena and Terminology. Geotectonics 48(2), Earthquake Track [online]. [downloaded 10/01/2017]. https://doi.org/10.1134 www.earthquaketrack.com.623–624. /S0016852114020058. 151–162, Influence of neotectonics on land surface evolution 151
(U/Th)/He thermochronometry: implications for uplift Aspects of the evolution of the Jemma River basin in the and erosion of the Afar plume dome. Earth and Planetary Kusák,Ethiopian M., Kropáček, Highlands J., SchillacI, determined C., Vilímek, by means V. of(2016): a relief https://doi.org analysis based on Digital Elevation Model (DEM), /10.1016/S0012-821X(03)00457-6. Pravea,Science A. R., Letters Batesa, 215(1–2), C. R., Donaldsona, 73–88, C. H., Tolandb, H., Condonc, D. J., Markd, D., Rauba, T. D. (2016): Geology Ethiopia. Geografia Fisica e Dinamica Quaternaria 39, and geochronology of the Tana Basin, Ethiopia: LIP of37–50, automatically https://doi.org/10.4461/GFDQ.2016.39.4. and manually extracted lineaments from volcanism, super eruptions and Eocene-Oligocene Kusák,DEM M., and Krbcová, geologically K. (2017): mapped Analysis tectonic of thefaults relationship around environmental change. Earth and Planetary Science the Main Ethiopian Rift and the Ethiopian Highlands, https://doi.org/10.1016 https://doi.org /j.epsl.2016.03.009. /10.14712/23361980.2017.1. Sembroni,Letters 443,A., Faccenna, 1–8, C., Becker, T. W., Molin, P., Abebe, B. Ethiopia. AUC Geographica 52(1), 5–17, (2016): Long-term, deep-mantle support of the Ethiopia- Mišurec, J., Eshetu, H., Haile, E. T., Alemayehu, L., Goslar, T. https://doi Kycl,(2017): P., Rapprich, Tectonic V., control Verner, ofK., complex Novotný, slope J., Hroch, failures T., in the .org/10.1002/2015TC004000. Ameka River Valley (Lower Gibe Area, central Ethiopia): Smith,Yemen M. J.,Plateau. Clark, C.Tectonics D. (2005): 35(2), Methods 469–488, for the Implications for landslide formation. Geomorphology visualisation of digital elevation models for landform https://doi.org/10.1016/j.geomorph mapping. Earth Surface Processes and Landforms 30(7), .2017.03.020. https://doi.org/10.1002/esp.1210. Macgregor,288, 175–187, D. (2015): History of the development of the Stoddart, D. R. (ed.) (1997): Process and form in East African Rift System: A series of interpreted maps geomorphology.885–900, Routledge, London. through time. Journal of African Earth Sciences 101, https://doi.org/10.1016 artefacts resistant method for automatic lineament /j.jafrearsci.2014.09.016. Šilhavý,extraction J., Minár, using J., Mentlík,Multi-Hillshade P., Sládek, Hierarchic J. (2016): Clustering A new 232–252, https://doi.org/10.1016/j.cageo.2016.03.015. Mäerker M., Hochschild V., Maca V., Vilímek V. (2016): Ukstins,(MHHC). I. A., Computers Renne, P. R. and Wolfenden, Geosciences E., BakeR, 92, 9–20, J., Ayalew, D., Stochastic assessment of landslides and debrishttps:// flows in Menzies, M. (2002): Matching conjugate volcanic rifted doi.org/10.4461/GFDQ.2016.39.5the Jemma basin, Blue Nile, Central. Ethiopia. Geografia margins: 40Ar39Ar chrono-stratigraphy Mangesha,Fisica e DinamicaT., Chernet, Quaternaria T., Haro, W. 39(1), (1996): 51–58, Geological Map of Ethiopia (1 : 250,000). Geological Survey of Ethiopia: and Yemen. Earth and Planetary Science Letters Addis Ababa, Ethiopia. of pre and syn-rift bimodalhttps://doi.org/10.1016 flood volcanism in Ethiopia Martino, C., Nico, G., Schiattarella, M. (2009): Quantitative /S0012-821X(02)005253. analysis of InSAR digital elevation models for U.S.198(3–4), Geological 289–306, Survey [online]. [downloaded 10/01/2017]. www.usgs.gov. southern Italy. Earth Surface Processes and Landforms Wladis, D. (1999): Automatic Lineament Detection Using identificationhttps://doi.org/10.1002/esp.1681. of areas with different tectonic activities in Digital Elevation Models with Second Derivative Filters. Photogrammetric Engineering and Remote Sensing 34(1), 3–15, https://pdfs.semanticscholar.org Minár,Western J., Sládek, Carpathians J. (2008): (Slovakia). Morphological Zeitschrift network für as /3c9e/4f7863002102ead6d390954ff4d1225e1f6d an indicator of a morphotectonichttps://doi.org field in the central .pdf65(4),. 453–458, /10.1127/0372-8854/2009/0053S3-0023. Wolela, A. (2010): Diagenetic evolution of the Ansian- Pik,Geomorphologie R., Deniel, C., Coulon, 53(2), C., 23–29, Yirgu, G., Hoffmann, C., Ayalew, Pliensbachian Adigrat Sandstone, Blue Nile Basin,
https://doi.org/10.1016/j.jafrearsci.2009.05.005. types.D. (1998): Journal The of northwestern Volcanology Ethiopianand Geothermal Plateau Research flood Wolfenden,Ethiopia, E.,Journal Ebiger, of C., African Yirgu, Earth G., Deino, Science, A., Ayalew, 6 (1), 29–42, D. basalts: classificationhttps://doi.org/10.1016 and spatial distribution of magma (2004): Evolution of the northern Main Ethiopian Rift: /S0377-0273(97)00073-5. birth of a triple junction. Earth and Planetary Science Pik,81(1–2), R., Marty, 91–111, B., Carignan, J., Lavé, J. (2003): Stability of the https://doi.org/10.1016 Upper Nile drainage network (Ethiopia) deduces from /j.epsl.2004.04.022. Letters 224(1–2), 213–228, 152 Original Article
Is the Second Demographic Transition a useful framework for understanding the spatial patterns of fertility change in Serbia at the beginning of the 21st century? Vladimir Nikitović1, Daniela Arsenović2, Aleksandar Sekulić3, Branislav Bajat3,*
1 Institute of Social Sciences, Belgrade, Serbia 2 University of Novi Sad, Faculty of Sciences, Serbia 3 University of Belgrade, Faculty of Civil Engineering, Serbia * Corresponding author: [email protected]
ABSTRACT Gaps in comprehension of demographic change in the region of ex-Yugoslavia after 1990, caused by a lack of reliable data series, frequent change of borders, and distinctive historical and cultural tradition in comparison to other post-communist societies, moti- vated us to contribute to the understanding of the spatial diffusion of recent profound fertility changes in South-Eastern Europe. We analysed changes in the spatial pattern and distribution of typical fertility indicators of the second demographic transition at the sub-national level in Serbia in order to find out whether these demographic shifts could be interpreted to be similar to those in Central and Eastern Europe. We found that differences in economic, historical, and cultural development between sub-regions of the country strongly affect spatial patterns of fertility change. Also, this paper suggests that the sub-regions forerunners of the first demographic transition could be considered as the cores of diffusion for the second demographic transition.
KEYWORDS South-East Europe; second demographic transition; sub-national fertility patterns; spatial autocorrelation; Serbia
Received: 7 March 2019 Accepted: 1 August 2019 Published online: 20 September 2019
Nikitović, V., Arsenović, D., Sekulić, A., Bajat, B. (2019): Is the Second Demographic Transition a useful framework for understanding the spatial patterns of fertility change in Serbia at the beginning of the 21st century?. AUC Geographica 54(2), 152–167 https://doi.org/10.14712/23361980.2019.14
© 2019 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). Second Demographic Transition and the spatial patterns of fertility change in Serbia 153
1. Introduction cultural tradition in comparison to the other former socialist countries (Sobotka 2003). In the European context, the second demograph- Moreover, divergent trends in the onset and tem- ic transition (SDT) was studied through various po of the fertility transition had been noted across socio-demographic processes including: fertility sub-regions of the former Yugoslavia during the decline, mortality impact on fertility change, marital socialist period, which was distinctive from the rela- postponement, non-marital fertility, partner relation- ships, changing household forms etc. Besides report- - ing general remarks on the very concept of the second gionaltive homogeneity differences of inEastern demographic Europe (Nikitović change, which2011, demographic transition as regards to its ideational were2016; well Josipovič documented 2016; Lerch in Yugoslav 2018). The demographic strong sub-re lit- and cultural drivers and universality issues (Coleman 2004; Esping-Andersen, Billari 2015), certain doubts of fertility transition in the region (Slovenia, Croatia, are often expressed with respect to the applicability anderature Serbia (Josipovič but excluding 2016), Kosovo). highlighted Declines early in adopters fertility of this theoretical concept outside the region where below replacement level started in Croatia and Serbia it was originally recognized as a process (Lesthaeghe, (excluding Kosovo) approximately at the same time as Neidert 2006). As one of the SDT proponents recently stressed, there is growing evidence that phenome- 2016). This peculiarity was acknowledged to a limited na such as sub-replacement fertility, rising shares of extentit started in newerin the vanguard studies (Frejka, countries Gietel-Basten of SDT (Nikitović 2016; cohabitations, and non-marital fertility are spreading worldwide (Lesthaeghe 2014). Still, some authors partly can be attributed to data issues generated by associate these factors with the rise of poverty formu- theNikitović, frequent Bajat, change Blagojević of national 2016; Lerchboundaries 2018), inwhich the lated through the ‘pattern of disadvantage’ (cohabi- region. tation as a characteristic of lower social stratum); When analysing the post-1990 period, Sardon therefore, expressing doubts that SDT mechanisms of (2001) noticed that strong sub-regional differences - of the former Yugoslavia with respect to other former ern Europe (Perelli-Harris, Gerber 2011). Further- socialist countries could be hidden if the levels of indi- more,cultural divergence shifts influence in trends demographic of SDT fertility changes indicators in East cators are examined as a whole, especially when ana- between and within European countries is evident lysing nuptiality levels and non-marital fertility. Addi- - ghe (2014) has recently suggested that SDT should to the Eastern bloc variant were often marked as a dis- be(Billari, rather Liefbroer conceived 2010). as a narrative For that for reason, understanding Lesthae tionally, specificities of Yugoslav socialism in relation- the profound cultural change that will sooner or later induce a variety of demographic changes thus imply- tinctive feature that can induce specific demograph ing that cross- and within-country differences in lev- concludeic outcomes that in the the patterns socialist of context change (Rašević typical for 2004). SDT els of SDT indicators are inevitable. couldHowever, not thesebe relevant specificities for explaining have led some fertility authors chang to- In this regard, exploring whether a fertility change es in the region of former Yugoslavia (Kuhar 2009); typical for SDT has been spreading from its source whereas, others argue that SDT has been blocked in towards the East seems to be a relevant research top- ic. The interpretation of abrupt and fast demographic Despite opinions that the concept of the second changes that followed the fall of the Iron curtain in demographicthis region since transition 1990 (Bobić, is limited Vukelić in 2011). Balkan coun- European post-socialist countries (the CEE term often tries, recent macro-level evidence suggests that Ser- comprises all former socialist countries including the bia1 may have begun to follow the post-socialist pace states successors of Yugoslavia) using the framework of this transition after a period of stagnation in the of SDT has been a frequent research topic in this cen- - tury (Philipov, Kohler 2001; Sobotka 2003, 2008, fore aim to examine if there is enough evidence for determining1990s (Nikitović, spatial Bajat, diffusion Blagojević of fertility 2016). change We there typi- - cal for SDT in Serbia at the beginning of this century. 2011;er 2010; Sobotka, Kurek 2011a; Zeman, Perelli-Harris, Kantorova 2003; Gerber Rašević 2011; In other words, this study aims to contribute to the 2004;Botev Muresan2012; Walford, 2007; Hoem Kurek et al. 2016). 2009; Yet, Kuhar, literature Reit understanding of the spatial diffusion of recent pro- explaining demographic changes in CEE countries found fertility changes occurring throughout South- after 1990 usually do not fully cover the region of for- East Europe. mer Yugoslavia, which is mainly due to the lack of reli- Although the proponents of SDT have no doubts able data series and frequent changes of borders ever when it comes to spatial spreading of this transition since Yugoslavia’s dissolution. In addition, the region from North-West to South-East Europe since the of former Yugoslavia also had a distinctive history or 1990s, a recent comprehensive study on non-marital
1 In this paper, it does not include the disputed territory of Kosovo, which unilaterally proclaimed independence of Serbia in 2008. Kosovo is currently recognized as an independent state by 108 of 193 UN member states and 23 out of 28 European Union member states. 154 Vladimir Nikitović, Daniela Arsenović, Aleksandar Sekulić, Branislav Bajat fertility (a typical SDT indicator) across sub-national regions of Europe showed that its spatial diffusion is researchers highlighted the role of macroeconomic not so straightforward and suggested that state bor- factorsThe explanations while others of thisthought relationship of economic differed – some and politi- ders are still more relevant than regional ones with - respect to non-marital fertility (Klusener, Perelli-Har- graphic transition in this region (Philipov, Kohler 2001:cal transitions 38; Kurek as a2011b). specific Anyway,part of the after second 1990, demo the from the aforementioned study showed that a more detailedris, Sanchez spatial Gassen analysis 2013). would However, shed somemore findingslight on demographic development inside the group of former understanding recent spatial patterns of fertility socialistdiversification countries in terms is obvious of social, due economic, to different and speed even change in South-Eastern Europe. For that reason, we of societal transformation (Botev 2012: 69). There- focus on the spatial pattern and distribution of select- fore, some authors suggest that the process of second ed SDT indicators at the municipality level in Serbia demographic transition ‘is not a unitary movement based on spatial autocorrelation analysis in order that reached all the countries in Central and Eastern Europe roughly at the same time and had the same be interpreted within the framework of the second features throughout, no more than it was in Western demographicto find out whether transition. recent In addition,demographic we aim shifts to discan- - marriage formation risk in selected countries in Cen- sity in comparison to countries of former communist tralEurope’ and (HoemEastern et al. Europe 2009: and 250). found They that investigateda particular regimes.cuss these findings in the context of regional diver drop in this indicator started in Hungary and Bulgaria Geographic differentials in diffusion of innovations after the early 1980s and in Russia and Romania half- a-decade and full decade later, respectively. According to Sobotka (2011: 262), a strong proved to be significant for sub-national differences decrease in period fertility rates have been accompa- wein the expect onset that of an the analysis first demographic using lower spatialtransition levels in nied and partly caused by a shift towards a later tim- willSerbia help (Nikitović, us to understand Bajat, Blagojević the spatial 2016). pattern Therefore, of the current fertility transition in Serbia and whether it with family transformations that resulted in the rise is similar to the patterns observed in CEE countries. ofing non-marital of first births. births. These However, changes in relationwere also to followedthe peri- The analysis is focused on common SDT indicators in od before 1990, when the share of non-marital births terms of fertility aspects of the process: total fertility did not exceed 10 per cent, the CEE region is cur- rate (TFR), mean age of women at childbearing (MAC), rently characterised by diversity of this typical SDT and percentage of births outside marriage (BOM).
(justindicator – from above 20 per Nordic cent). values in Bulgaria (more than 2. Overview of the Fertility Change in Serbia 50 per cent) to half the West Europe values in Poland in the Context of the Post-Socialist Societies 2.1 Change in Fertility Rates in Serbia, 1950–2017 The general pattern of changes in the period total fer- Fertility changes that could be interpreted in the tility rate in Serbia since 1950 was similar to those framework of SDT lagged by at least 20 years in observed in most European countries: post-war baby post-socialist countries when compared with the boom followed by a decrease to the replacement level, forerunners of the process in North-Western Europe. and then to the sub-replacement level (Figure 1). The During the socialist period, CEE countries were char- main difference in Serbia is that the baby boom came acterized by different demographic development, par- to its end already in the late 1950s. The decline in fer- ticularly concerning fertility, in relation to the rest of tility was faster and ran deeper than in most Euro- Europe (Frejka, Gietel-Basten 2016). ‘A demographic pean countries in spite of the same general determi- distinctiveness of European former Communist coun- nants of fertility decline, such as the adoption of new tries as a relatively homogeneous region had gradu- norms and values and growing levels of female labour ally developed between the mid-1960s and the mid- force participation (Kupiszewski, Kupiszewska, Niki- 1980s’ (Sobotka 2003: 453). In the 1970s and 1980s, institutional and cultural factors that jointly sustained below the replacement level by 1971. Several factors the regime of universal and early reproduction under statedtović 2012). in literature As a result, could TFR help was one already to comprehend 15 per cent state socialism protected Eastern Europe from the fer- such a distinctive pattern. Those often include early tility decline experienced in all other parts of Europe. liberalization (from 1952 to 1969) of women’s right Then, following the collapse of communist regimes, within a decade the CEE region experienced the low- was typical for former socialist countries (Sardon est fertility rates in Europe (Sobotka 2003, 2011). 2001;to abortion Sobotka, (Rašević Zeman, and Kantorova Sedlecky 2003; 2009: Frejka 385), 2008),which A sharp decrease in TFR in Eastern Europe dur- and structural factors such as fast secularization and ing the 1990s was usually associated with political, industrialization during the period of Yugoslavia social and economic distortions during the period. implying abrupt migration from rural to urban areas Second Demographic Transition and the spatial patterns of fertility change in Serbia 155
Fig. 1 Total fertility rate (double-line) and 5-year age-specific rates (solid, dashed and dotted lines) in Serbia excluding Kosovo, 1950–2017. Notes: Time series after the census year of 2011 does not include the effect of international migration. TFR – Total fertility rate (the mean number of children that would be born alive to a woman during her lifetime if she were to pass through her childbearing years conforming to the fertility rates by age of a given year, and surviving); ASFR – age-specific fertility rates (the number of live births occurring during a given year per 1,000 women of reproductive age classified in five-year age groups). Source: Statistical Office of the Republic of Serbia (Demographic Yearbook 2017; documentation tables). associated with housing issues in urban areas. Yugo- belonged to the lowest fertility group have recently slav society, unlike socialist states in CEE at the time, experienced a rise in TFR (Latvia, Russia, and Slove- was quite open to the impact of Western values in the nia have even exceeded 1.5), which should be expect- sense of satisfying individual needs to a much great- ed considering the pace of ‘postponement transition’ er extent, which arguably contributed to an increase (Goldstein, Sobotka, Jasilioniene 2009). in the economic and psychological cost of children. Figure 1 demonstrates a transition in the age pat- tern of fertility in Serbia, which occurred over the last population policy during the period were introduced 60 years. Very high baby-boom rates at ages older than Furthermore, no efficient measures in the domain of 20 rapidly declined during the late 1950s and early The subsequent sharp fall in the TFR in Serbia 1960s and kept principally stable during the 1970s began(Rašević in 2004:16).the late 1980s, similar to Slovenia and Croa- and 1980s. Only the fertility rates of adolescents expe- tia. These three republics of former Yugoslavia were rienced an increase and relatively high levels during recognised as the ‘early starters’ sub-region in terms all the period of socialist Yugoslavia. This is usual- ly explained by the moral primacy of patriarchy in a steeper fall could have been expected due to wars spite of the new social institutions introduced by the andof fertility institutional transition crisis (Nikitovićduring the 2017). 1990s, Althoughthe total socialist regime (Lerch 2018). However, since the late fertility rate in Serbia has not experienced the low- 1980s, the fertility rates of females younger than 25 est levels recently reported in Europe. Moreover, the began to decline sharply, while those of females older war ambience could be a reason that contributed to than 30 started to increase, actually at a slower pace. prolonged consistency of traditional patriarchal val- This was particularly prominent for the age groups ues associated with family and childbearing (Sardon total fertility rate in Serbia has oscillated between 1.4 20–24 and 30–34, as they swapped their positions in and2001; 1.48 Rašević since 2004;2005. PetrovićMost former 2011; socialist Lerch 2018). states thatThe theterms last of 60 contribution years. Moreover, to the totalthe difference fertility – the between latter group currently has a higher rate for the first time in 156 Vladimir Nikitović, Daniela Arsenović, Aleksandar Sekulić, Branislav Bajat these two groups regarding their contribution to total cent in 1950 to 26.3 per cent in 2017. However, the fertility rate has been increasing since 2009. The high- position of Serbia regarding the share of non-marital births has shifted downwards in the European con- text during the last 60 years. Since the late 1950s and est fertility rates in the twenty-first century refer to during the 1960s the share was stable around 12%, thewomen current aged age 25–29 pattern (despite that is the quite recent similar drop) to the and one is placing Serbia just behind Iceland and Austria in 1960 observedclosely followed in Poland by women(Walford, aged Kurek 30–34, 2015). which Although forms (12.2 per cent), and behind Iceland, Sweden, Estonia, the indication of ‘the postponement of childbearing for opportunistic reasons and on account of structur- 2010). Given the universality of marriage and child- al impediments’ could have been noticed during the bearing,and Austria the in rise 1970 of the(11.8 share per cent)at the (Penev, time was Stanković due to clear onset in postponing childbirths in Serbia seems caused by the termination of the short-lived baby- toperiod be associated of former withYugoslavia the downfall (Rašević of 2004: the socialist 9), the booma significant period (Figuredrop in 2), the much number shorter of withmarital respect births to system, thus resembling the trends observed in oth- most European countries (Kupiszewski, Kupiszewska, er post-socialist societies. Consequently, the average age of childbearing increased from 25.9 in 1991 to The share of non-marital births was almost the sameNikitović in 1968 2012). and 1988, while in between it dropped. half a year than in Slovenia, Estonia and Croatia, and by This was induced primarily by the emergence of the almost29.8 years a year in 2017. from Thethe EU latter average, value whileis lower it is by similar about baby-boom echo generations, which accordingly to that in Latvia, Hungary and Poland (Eurostat 2019). affected the number of marital births and lowered
2.2 The Share of Non-Marital Births in Serbia, of the indicator started in the late 1980s as the end 1950–2017 ofthe the share socialist of non-marital system was births. approaching, The significant a trend thatrise was similar to what was observed in most of former The share of non-marital births is one of the typical socialist countries. However, the share of non-mari- - tal births in Serbia in 1990 (13.1 per cent) was just iour associated with the second demographic tran- above those in the Netherlands and Belgium and indicatorssition. This reflecting indicator changes raised in demographicSerbia from 8.0 behav per ahead of most former socialist societies. In 2017,
Fig. 2 Number of livebirths by legal marital status (primary axis) and the percentage of extra-marital births (secondary axis) in Serbia excluding Kosovo, 1950–2017. Note: Time series after the census year of 2011 does not include the effect of international migration. Share of extra-marital births – the proportion of live births outside marriage. Source: Statistical Office of the Republic of Serbia (Demographic Yearbook 2017; documentation tables). Second Demographic Transition and the spatial patterns of fertility change in Serbia 157
Tab. 1 Distribution of municipalities (LAU-2) according to the categories of relevant input data in the census years of 2002 and 2011.
Frequency of LAU-2 Number of Frequency of LAU-2 Frequency of LAU-2 Total population Number of women aged size livebirths 2002 2011 15–49 2002 2011 2002 2011 <10,000 6 11 <2,000 5 13 <100 11 36
10,000–19,999 50 56 2,000–4,999 60 64 100–199 49 44
20,000–49,999 60 51 5,000–9,999 40 37 200–499 50 45
50,000–100,000 27 23 10,000–20,000 30 26 500–1000 30 16
>100,000 18 20 >20,000 26 21 >1000 21 20
Average / LAU 46,571 44,639 Average / LAU 11,238 10,033 Average / LAU 488 418 Note: Data on the number of livebirths refer to the three-year average over the census years. Source: Statistical Office of the Republic of Serbia (Census of Population, Households and Dwellings in the Republic of Serbia, 2002 and 2011; documentation tables). it was positioned much lower, leaving behind, apart - from Switzerland, Azerbaijan, Cyprus, Greece, and Turkey, only countries from Eastern Europe (Poland, liveStatistical births Officeoutside of marriage.the Republic Each of Serbia.indicator Legal was mar cal- Belarus, Ukraine, Moldova) and the successor states ital status was used as a criterion for identification of of the former Yugoslavia except Slovenia (Eurostat 2019). The share of non-marital births in Serbia culated for the last two census years – 2002 and 2011, as a three-year average over the periods 2001–2003 pre-marital cohabitations that are more typical for calculationsand 2010–2012, were respectively, made following in order the administrative to smooth out peoplein the in twenty-first disadvantaged century economic resulted positions, mainly such from as divisionannual fluctuations in 2002 as noin thedata number for 2002 of werelivebirths. available The lone mothers including those in extended multiple for the new six municipalities that had been estab- - lished by the time of the 2011 census. Consequently, the results of the analyses in the paper were distrib- generation families (Penev, Stanković 2010; Stank uted across the total of 161 municipalities in Serbia ović 2014). It suggests that the recent rise of the excluding Kosovo. Out of that number, 42 municipal- indicator does not reflect spreading of cohabitation ities have city status (16 of them are actually includ- post-socialistas an alternative societies, to marriage where (Petrovićthe diffusion 2011; of cohabBobić,- ed within the City of Belgrade) based on the admin- itationVukelić has 2011). been This rather is in slow line (Muresanwith findings 2007, from Sobotka most istrative criterion, which comprise 63.2 per cent of 2008). the total population in the country as to the 2011 The following spatial analysis of the SDT indicators census. Table 1 shows the basic descriptive statistics (TFR, MAC and BOM) will show that a variety of oth- of the relevant input data across municipalities in er factors need to be considered in Serbia in order to the census years of 2002 and 2011. The interactive gain an understanding of their substantially diverse map depicting the values of calculated SDT indica- distribution at sub-regional levels. tors on the municipality level (LAU-2) is available in HTML format at the Web page http://osgl.grf .bg.ac.rs/static/materials/demography/tfr_mac_bom 3. Data and Methods .html
3.1 Data sets and territorial disposition 3.2 Global and Local Spatial Autocorrelation Indices
The spatial analyses were conducted at the local The global Moran’s I statistics (O’Sullivan, Unwin administrative unit (LAU) 2 level, which refers to 2003) and local indicators of spatial association municipalities in Serbia (Figure 3). The values of the (LISA) (Anselin 1995) were calculated in order to outline spatial autocorrelation globally and to identify the municipalities’ borders. These analyses are exclu- spatial clusters with respect to percentage difference analysed indicators were matched to GIS shapefiles of- in the SDT fertility indicators for 2002 and 2011. tion estimates for the years out of intercensal periods As global Moran’s I index recapitulates autocorre- dosively not referred include international to census years migration. as the official In this populaway, we lation values over the entire study area, inherent spa- wanted to get the values of the indicators as realistic tial autocorrelation is necessary to calculate indices as possible, which is particularly important in case of typical emigration municipalities. clustering around an individual location. like LISA in order to assess significant local spatial- tion (ISA) (ESRI 2012) was conducted to determine Data on live births and population classified In the first step, Incremental Spatial Autocorrela across five-age groups and sex were provided by the 158 Vladimir Nikitović, Daniela Arsenović, Aleksandar Sekulić, Branislav Bajat
Fig. 3 The map of Serbia with local administrative units (LAU-2). Second Demographic Transition and the spatial patterns of fertility change in Serbia 159 the distance associated with peak clustering for all could be more typical for highly urbanized and more SDT indicators. The value obtained from ISA was then developed sub-regions of the country. In both cases used as a distance threshold or radius for determin- this correlation suggests that cohabitation is still not ing proximity weights (O’Sullivan, Unwin 2003) when common alternative to marriage as it was noted in calculating global and local spatial autocorrelation recent studies on fertility and union formation in Ser- indices. determined among countries of ex-Yugoslavia exclud- bia. The similar findings on this SDT indicator may be 4. Spatial Analysis of Fertility Change in 2014). The following spatial analysis across munici- Serbia, 2002–2011 palitiesing Slovenia in Serbia as well provides as in Poland a deeper and insight Slovakia in sub-na (Bobić- tional differences in fertility. The analysis of the three fertility indicators connect- The two methods used aimed at checking wheth- ed with SDT at the level of local administrative units er there is a spatial spread of SDT fertility indicators refers to the main objective of the study. Table 2 sum- throughout Serbia. One implies a descriptive analy- marizes the changes of these indicators across munic- sis of the spatial differences of each indicator across ipalities in Serbia from 2002 to 2011. the country and their dynamics between the two last census years, 2002 and 2011, while the other relies than in urban areas in this period. As expected, the on spatial autocorrelation indices calculated upon gapThe in MACdrop hasin TFR even was extended significantly between higher urban in rural and relative changes (%) over the intercensal period rural areas. Although the higher rise in BOM could be expected in rural than in urban areas in accord- ance with the recent changes in the European context 4.12002–2011. Total Fertility Rate (Walford, Kurek 2015), a BOM in rural municipali- ties is 30 per cent higher than in the cities of Serbia According to the 2002 census, TFR was below the according to the 2011 Census. This peculiarity, which replacement level in almost all of the 161 municipal- is analysed in more details across municipalities, ities in Serbia, as only 8 of them reached 2.1 or above - (Figure 4(a)). Five of these eight municipalities were ical for countries where this process was originally predominantly populated by Muslims, the religious described.might point to the specific forms of SDT that are atyp group who were the last in the country to enter the The next introductory step of the spatial analysis - involves inspecting whether there is a correlation ing three municipalities belong to the group with the between each of the two pairs of indicators distrib- highestfinal stage share of ofdemographic persons living transition. abroad for The more remain than uted across municipalities during the two census a year. Unlike 2011, the TFR for 2002 included the births by this group of emigrants indicating that the was found indicating that the mean age of mothers TFR was actually lower than 2.1 in these municipali- atyears, childbirth 2002 and is negatively2011. The onlyassociated significant with correlation the share Due to a further decline of TFR, only three munici- palitiesties (Rašević, with predominantly Penev 2010). Muslim populations were of births outside marriage, with Pearson coefficient- above 2.1 in 2011, while ten municipalities dropped tionof correlation pattern includes raising betweentwo types 2002 of association, (r = −0.35, butp < .001) and 2011 (r = −0.44, p < .001). This correla- cent of municipalities (21 of 161) did not decline in ties where MAC is lower while BOM is higher and the TFRto just between around the1 (0.99–1.09) census 2002 (Figure and 2011. 4(b)). TheOnly highest 13 per otherwith opposite to municipalities directions – one where refersMAC is to higher municipali while increase (mainly lower than 10 per cent) took place in the northern, most urbanized, part of the country developed and mainly rural areas, while the latter one including municipalities in the area of the two largest BOM is lower. The first one could be expected in less
Tab. 2 Quantitative analysis of the second demographic transition indicators in urban and rural local administrative units in Serbia, 2002–2011.
Total fertility rate Mean age of women at childbirth Births outside marriage (%)
Urban Rural Urban Rural Urban Rural
year 2002 1.54 1.65 27.37 26.18 20.46 23.69
year 2011 1.44 1.35 28.90 27.16 23.27 30.29
% difference −6.51* −17.79** 5.59** 3.73** 13.75* 27.87** ** Significant atp = 0.001; * significant atp = 0.05 using Mann–Whitney test. Note: A municipality is classified into urban or rural depending on its administrative status according to the current administrative division of local spatial units (cities and municipalities). Source: Statistical Office of the Republic of Serbia (documentation tables). 160 Vladimir Nikitović, Daniela Arsenović, Aleksandar Sekulić, Branislav Bajat
Fig. 4 Total Fertility Rate in Serbia at LAU-2 level: (a) 2002, (b) 2011, (c) percentage difference for the period 2002–2011, (d) LISA cluster map of percentage difference for the period 2002–2011. Second Demographic Transition and the spatial patterns of fertility change in Serbia 161
- 4.2 Mean Age of Women at Childbirth tal of the province of Vojvodina). This result indicates thecities – Belgrade effect of postponing (the capital) births andto later Novi ages Sad (more (the capi typ- The highest values of MAC in 2011 are recorded in ical in urban areas), rather than the rise in the cohort the largest urban centres of Serbia, almost exclu- sively in those located in north and west parts of the A quarter of all municipalities had experienced fertilitya decrease (Rašević in TFR 2015). between 20 and 40 per cent from Serbia, while the lowest values refer to the east (Fig- 2002 to 2011 (Figure 4(c)). Most of these munic- urecountry, 5(b)). except The numberfor Niš – the of municipalities largest city in with southern MAC ipalities are in eastern and southern Serbia, the least developed regions in the country according to 60) between 2002 and 2011. The minimum values of higher than 28 years increased five times (from 12 to website of the National Agency for Regional Devel- area of the country) remained almost the same after opment.the 2014 Theofficial observed scale of decline development during listed the periodon the aMAC decade (observed (24.4 yearsin the insoutheast – the 2002 compared least with developed 24.7 in was even larger than 40 per cent in several munic- 2011), while maximum MAC (recorded in the histor- ipalities in the northeast of this area, which is dis- ic centre of the Belgrade City) increased from 30.4 to tinguished by the highest share of population living - - ities did not report a rise of MAC between 2002 and ference in methodology when the rates for 2002 and 2011,32.2 years clearly (Figures indicating 5(a–c)). the transition Only 9 of 161in the municipal peak fer- 2011abroad are (Penev, compared Predojević-Despić (the latter one 2012). did notGiven include a dif tility towards older ages throughout the country. The births by emigrants), this strong reduction could not largest increase in MAC (between 6 and 10 percent) during the 2000s occurred mainly in municipalities Penev 2010). with smaller urban centres (most of them located in onlyGlobal be attributed Moran’s toI = the 0.342 actual (Z =fertility 11.824, fall p (Rašević,< 0.001) Central Serbia) and those where relatively lower val- ues of MAC had previously prevailed. indicates moderately clustered spatial patterns with Global Moran’s I = 0.144 (Z = 5.338, p < 0.001) for for the TFR (threshold radius h = 70 km) in Serbia- - pronounced clustered spatial patterns with high sta- torshigh assortedstatistical by significance. different types The ofLISA spatial cluster correlation: map des the MAC (threshold radius h = 140 km) suggests less- 1)ignates High-High locations association, with significant a municipality local Moran with indica high tistics for the mean age at childbearing, there is almost TFR value has neighbouring municipalities with high tisticala continuous significance. High-High As regards cluster to connecting local Moran’s the Iarea sta TFR values (positive spatial correlation); 2) Low- of the Belgrade City and Central Serbia (Figure 5(d)). Low association, a municipality with low TFR value The distribution of all four types of clusters could be has neighbouring municipalities with low TFR values (positive spatial correlation); 3) Low-High associ- areas where MAC was previously very low indicating ation, a municipality with low TFR value has neigh- generalexplained spatial by two diffusion processes. of demographic The first one innovation. refers to bouring municipalities with high TFR values (neg- The other process is associated with recent trends in ative spatial correlation); 4) High-Low association, urbanisation (mostly pronounced along important a municipality with high TFR value has neighbouring road corridors in Serbia) indicating a higher mean age municipalities with low TFR values (negative spatial correlation). 2015). The local Moran’s I statistics singled out a big com- at childbearing in urban than in rural areas (Rašević pact High-High cluster in TFR denoting a metropolitan 4.3 Live Births Outside Marriage area that includes the two largest cities (Belgrade and Novi Sad) located in northern Serbia (Figure 4(d)). Eastern Serbia differs quite distinctively from the rest of the country in terms of the percentage of in TFR between 2002 and 2011 in the area. It might be births outside marriage. The BOM values higher than explainedThe cluster, by primarily, both the development reflects a quite of housing small increase associ- 30 per cent were recorded in almost all municipali- ated with a recent concentration of industry particu- ties in eastern Serbia, where BOM is even higher than larly in suburbs of the area induced by the proximity 40 per cent in 11 municipalities for the year 2002. to important highways, and the ongoing process of On the other hand, BOM is lower than 20 per cent in - 71 of 161 municipalities in the country (Figure 6(a)). In 2011, the observed spatial polarization became insouth-to-north East Serbia indicatesinternal migration a decrease (Nikitović, in TFR, whichPredoje is even more pronounced (Figure 6(b)). The increase partlyvić-Despić, caused Marinković by the changes 2015). in A methodology big Low-Low between cluster of BOM values in eastern Serbia is noticeable when - compared to stagnating rate in southwestern Serbia rect effects of continuous international emigration (9 municipalities with BOM less than 10 per cent). the 2002 and 2011 census, but also reflects the indi- Generally, the increase of BOM values and BOM range mann, Musil, Gruber 2014). throughout the country are obvious in the period of from this area (Penev, Predojević-Despić 2012; Fass 162 Vladimir Nikitović, Daniela Arsenović, Aleksandar Sekulić, Branislav Bajat
Fig. 5 Mean age of women at childbirth in Serbia at LAU-2 level: (a) 2002, (b) 2011, (c) percentage difference for the period 2002–2011, (d) LISA cluster map of percentage difference for the period 2002–2011. Second Demographic Transition and the spatial patterns of fertility change in Serbia 163
Fig. 6 Share of Live Births Outside Marriage in Serbia at LAU-2 level: (a) 2002, (b) 2011, (c) percentage difference for the period 2002–2011 (%), (d) LISA cluster map of percentage difference for the period 2002–2011. 164 Vladimir Nikitović, Daniela Arsenović, Aleksandar Sekulić, Branislav Bajat
national level, however special focus was maintained below 20 per cent dropped to 50, while the num- at the level of municipalities. ber2002–2011. of those Thewith number a share of of municipalities births outside with marriage BOM Concerning the national level, the total fertility rate higher than 40 per cent reached 29, of which 13 had and mean age at childbearing in Serbia follow paths - similar to those in CEE countries. Even though the last ities experienced a decline of BOM between 2002 sharp decrease of TFR began a little earlier than in anda share 2011, higher while than the 50 strongest per cent. drops Only were20 municipal localized CEE countries, it did not reach the lowest values as in the southwestern area characterized by the lowest reported in CEE. This can be explained through a not values of the indicator. A total of 23 municipalities so sharp transition from state socialism to new soci- experienced a rise of BOM higher than 50 per cent etal conditions and stubborn subsistence of tradition- (Figure 6(c)). al values associated with family formation during the Global Moran’s I = 0.078 (Z = 2.818, p < 0.006) for - 2011; Lerch 2018). The latter reason could also have tain degree of clustered spatial patterns with statis- ancivil impact wars inon the the 1990s much (Petrovićlower rate 2011; of non-marital Bobić, Vukelić fer- the BOM (threshold radius h = 52 km) indicates a cer tility in Serbia than in most of CEE countries, which of BOM of around or below 20 percent was recorded is almost the same as in Poland, and higher only than intical the significance. largest cities It of should the country be noted regardless that the of share their those in former Soviet republics and Croatia (Eurostat 2019). of them in eastern Serbia) experienced BOM higher thangeographic 30 per location, cent. The while local only Moran’s five urban I statistics centres (allfor the analysed fertility indicators with respect to the BOM denoted two clusters that represent two dis- The three general findings concerning changes of tinctively different areas (east and west in relation to decade of this century: a) the decline in fertility has the main road corridor) with regards to the percent- diffusedtype of municipality from urban into Serbiarural areas, were b)pointed the gap in inthe mean first age difference in the share of births outside marriage age at childbearing between urban and rural areas between 2002 and 2011 (Figure 6(d)). One distinctive has widened, and c) the share of non-marital births area refers to the largest urban centres in the more - developed west sub-region, including central urban as. The diffusion of values and ideas associated with zones of the Belgrade City, characterized by the low- lowbecame fertility significantly norms from higher urban in rural centres than to in peripheralurban are est share of BOM and its change over the period. The areas of the country was already well noted during other distinctive area refers to municipalities in the less developed east sub-region with the highest share Also, the pronounced postponing of births in cities isthe in firstaccordance demographic with the transition recent changes (Watkins in post-so 1990).- observation could be considered as a paradox. How- cialist societies (Walford, Kurek 2015). However, ever,and highest it was notedincrease that of ‘cohabitations BOM. On the first and glance,extra mar this- such a discrepancy between rural and urban areas in ital births in Serbia are more common among actors terms of non-marital fertility might be unexpected if we neglect that non-marital births in Serbia are most- - ly reported among women from lower social strata, comingat the lower structural end ofrisks, the although, stratification without ladder ideational … [as a result of] … rationally developed strategies in over 78). In addition, East Serbia could be considered as analysis,which are as often the pivotal located analytic in rural tool, areas denoted (Bobić that 2014). the changes typical to the theory of SDT’ (Petrović 2011: metropolitanMore specific area, results which of includes the autocorrelation the two largest spatial cit- demographic change since the mid-nineteenth centu- ies located in the northern Serbia, is the only one with a specific sub-region given the distinctiveness of its at the beginning of this century. At the same time, ry (Knežević 2013). spatiala mild butdiffusion significant of the increase rise in mean of TFR age in at the childbear country- 5. Discussion and Conclusions ing was registered throughout the country but was mostly pronounced along important road corridors The objective of this study was to determine wheth- connecting bigger towns. This indicates that the pro- er recent spatial patterns of fertility change in Serbia, cess is particularly associated with recent trends in which we assumed to be associated with the second urbanisation. Our results generally coincide with the demographic transition, agree with those observed in post-socialist societies of Europe. Whereas most of transformation of reproductive behaviour began ear- the literature exploring whether and how SDT taking lierfindings in the from largest Czechia city districtsindicating and that districts ‘this stage with in eco the- place in former Yugoslavia refers to the national lev- nomic centres.’ As such, the spatial patterns of fertility el and to aspects of family change, we wanted to see change demonstrate different speed of adaptation to if sub-national geographic differentials in diffusion the changes in reproductive conditions (Šprocha, Šid- of demographic innovations measured by key fertil- lo 2016: 228). In this respect, East Serbia is clustered ity indicators associated with SDT could shed more in terms of a decrease in TFR, which could be to some lights on the issue. We deployed our analyses from the extent associated with indirect effects of continuous, Second Demographic Transition and the spatial patterns of fertility change in Serbia 165 long-lasting international emigration from that area. are not so dividing with respect to SDT indicators. For The diffusion of the increase in the proportion of example, the rise in non-marital fertility rates in Bul- non-marital births from its core in East Serbia was garia resembles the one in East Serbia. evident with respect to spatial statistics on the per- On a national scale, Serbia and Poland are very centage difference of this indicator between 2002 and similar in terms of the level and family settings of 2011. The two distinctive opposite clusters (East and non-marital fertility. However, only a closer look at Southwest Serbia) suggest that not only the tempo of the sub-national patterns of the diffusion of this indi- this diffusion was rather slow, but the persistence of traditional values associated with family formation is very strong in the Southwest, which is recognized as andcator religious reveals differencesthat ‘diversified of the levelPoland’s of extramarital population’ the sub-region that was the last to enter fertility tran- births in the spatial arrangement reflects cultural the diffusion of this SDT indicator in Serbia is clearly - guided(Kurek 2011a: 393).by differences Likewise, in historical our results and cultural suggest herthat- sitionlysed indicators(Nikitović, Bajat,of SDT, Blagojević particularly 2016). as to the share itage between sub-regions of the country. Therefore, of non-maritalSpecificities fertility,of East Serbia belong that to the relate outcomes to the anathat the spatial pattern of the share of non-marital births, might be unexpected to some extent. Furthermore, as an indicator of ideational change, shows that typi- this sub-region stands as an exception with respect cal SDT drivers may not be the decisive determinant for recent fertility change in Serbia. The fact that BOM division between West and East Europe, which seems is not only the result of modern cohabitations but is to bethe still Hajnal relevant line thatfor differentials reflects the in nuptiality non-marital regime fer- more attributed to a rise in unstable partnerships that tility between sub-regions in Europe (Walford, Kurek end up with single mothers of poor socio-economic 2015). East Serbia was the innovating sub-region as status supports this conclusion. Furthermore, the slow rise of BOM since the beginning of this century nineteenth century). Also, it could be an innovating along with the highest and increasing level of post- areafor the regarding first fertility the transition second transition, in the country similarly (the late to ponement of births in urban and most developed what Lesthaeghe and Neels (2002) concluded explor- sub-regions of the country questioning ideational ing the spatial diffusion of SDT in Belgium and France. roots of the change assumed by SDT. The societal anomie, slow transition to market economy and sta- East Serbia were atypical in a way similar to what ble political conditions typical for the post-socialist The drivers of the first demographic transition in- defer childbearing in case of the higher educated or thecould Vlach be the population, current specificity. the fertility In this transition sub-region, already rec tocountries, adapt their may strategies also influence associated individual with decisions family for to- startedognized inby the the secondgreat influence half of the of the nineteenth cultural model century of mation in case of persons in disadvantaged socio-eco- with the acceptance of birth control through the ‘one- nomic position (Perelli-Harris, Gerber 2011). However, as Van de Kaa (2004: 8) pointed in his of the transition was agrarian overpopulation due to paper revisiting the concept of SDT, ‘it is not of any thechild extensive system’. farming The most of predominantlycommonly identified rural popula cause- - tions, which increased after the end of an immigration ent cultural endowments do not rapidly converge to that was induced by the enlargement of the territory agreat standard significance pattern’. if countries Also, it is and possible regions that with societies differ of Serbia. On the other hand, there were also indica- tions that the trigger was endemic diseases such as progression of SDT, as suggested by Sobotka (2008). Inlike other the Serbian words, society the pattern just follow of advantage a specific and path the of pattern of disadvantage could be just two opposite bettersyphilis perceived (Knežević if 2013).it is considered from a longer his- sides of the same coin. This seems to be reasonable toricalTherefore, and socio-cultural the specificity perspective, of East Serbia which might high be- - ing decisions of individuals concerning family forma- understanding current fertility change in the sub-re- because both modern and traditional values influenc lights the importance of specific ethnic differences for- ders on the diffusion of fertility changes associated slowtion were transition noticed to atwhat the some same authors time in referSerbia to (Bobić, as the withgion. SDT Furthermore, could be questioned the influence in this of part national of Europe. bor genderVukelić egalitarian2011). It could family be normsalso interpreted (Esping-Andersen, as a very A recent study (Klusener, Perelli-Harris, Sanchez Gas- Billari 2015). The similar trend was noted in Romania sen 2013) pointed out that most prominent nation- indicating that the SDT process might be in progress al dividing lines with respect to non-marital fertility as most of the ‘threshold levels’ of the onset of SDT after 1990 shifted from West and Central to Southeast were surpassed with exception of cohabitation, ulti- - mate celibacy, and modern contraception (Muresan itations regarding different levels of regional decom- 2007: 65). Moreover, our study shows that the tran- positionEurope. The particularly finding was between probably formerdriven by Yugoslavia study lim sition between traditional and modern concerning and its neighbours. The spatial analysis in this paper the demographic regime might have its geographical showed that national boundaries in Eastern Europe dimension in Serbia, too. 166 Vladimir Nikitović, Daniela Arsenović, Aleksandar Sekulić, Branislav Bajat
Finally, as Thornton and Philipov (2009) noticed, Hoem, J. M., Kostova, D., Jasilioniene, A., Muresan, C. (2009): the ideational change behind SDT is important in Traces of the second demographic transition in four explaining changes in marriage, cohabitation and selected countries in Central and Eastern Europe: Union childbearing in CEE; particularly for the belief that formation as a demographic manifestation. European adopting modern family systems will help to produce Journal of Population 25(3), 23 55, https://doi .org/10.1007/s10680-009-9177-y. modern political and economic accomplishments as 9–2 they are interdependent. If one assumes that Serbia a demographic crossway: 25 years after the collapse of is on the course of SDT, it seems that further diffusion Josipovič,Yugoslavia. D. (2016): Stanovništvo The post-Yugoslav 54(1), 1 space0, https://doi on of the ideational changes throughout the region could .org/10.2298/STNV160415006J. be closely related to the process of EU enlargement Klusener, S., Perelli-Harris, B., Sanchez5–4 Gassen, N. (2013): Spatial aspects of the rise of nonmarital fertility across include the whole region of the former Yugoslavia in Europe since 1960: The role of states and regions in furtherto the East. research For thaton sub-national reason, it will patterns be beneficial of fertility to shaping patterns of change. European Journal change. Currently, the issues on availability and quali- of Population 29(2), 13 65, http://dx.doi ty of relevant demographic data are a serious obstacle .org/10.1007/s10680-012-9278-x. to accomplish this. 7–1 Knežević,Srbije. [PhDA. (2013): diss.]. Istorijskodemografske Beograd: Univerzitet u Beogradu. Kuhar,i etnodemografske M., Reiter, H. (2010): osnove Transformation razvitka stanovništva Istočne References and demographic change in the ex-Yugoslav
Anselin, L. (1995): Local indicators of spatial association- perspectives on reproductive trends. Annales for Istrian http://dx.doi andcountries – materialist, Mediterranean Studies, idealist Series and Historia institutionalist et Sociologia .org/10.1111/j.1538-4632.1995.tb00338.x. 20(1), 1 6. Frejka,LISA. T., Geographical Gietel-Basten, Analysis S. (2016): 272, Fertility 93–115, and family policies in Central and Eastern Europe after 1990. Impact of3–2 Demographic and Migration Flows on Comparative Population Studies 41(1), 6, http:// Kupiszewski,Serbia. Belgrade, M., Kupiszewska, International D., OrganizationNikitović, V. (2012): for dx.doi.org/10.12765/CPoS-2016-03en. Billari, C. F. (2005): Europe and its fertility:3–5 from low to Kurek, S. (2011a): Population changes in Poland: A second lowest low. National Institute Economic Review 194(1), demographicMigration – Mission transition. to Serbia. Procedia Social and Behavioral 5 3, https://doi Sciences 19, 38 96, http://dx.doi.org/10.1016/j .org/10.1177/0027950105061496. .sbspro.2011.05.146. Billari,6–7 C. F., Kohler, H. P. (2004): Patterns of low and Kurek, S. (2011b):9–3 Double transitions? Regional patterns lowest-low fertility in Europe. Population Studies 58(2), 16 76, http://dx.doi.org/10.1080/00324720420002 Series B, Human Geography 93(2), 16 84, https://doi 13695. .org/10.1111/j.1468-0467.2011.00367.x.of population ageing in Poland. Geografiska Annaler: 1–1 Lerch, M. (2018): Fertility and union formation3–1 during transformation in Serbia. Sociologija i prostor 52(2), crisis and societal consolidation in the Western Balkans. Bobić,20 M. 25,(2014): https://doi.org/10.5673/sip.52.2.5 Specific (blocked) post-transition. union Population Studies 72(2), 21 34, https://doi.org /10.1080/00324728.2017.1412492. de-blocked?7–2 Sociologija 53(2), 14 76, https://doi Lesthaeghe, R. (2014): The second7–2 demographic transition: Bobić,.org/10.2298/SOC1102149B M., Vukelić, J. (2011): Second. demographic transition A concise overview of its development. Proceedings of Bongaarts, J. (2002): The end of the fertility9–1 transition the National Academy of Sciences 111(51), 1811 in the developed world. Population and Development 18115, http://dx.doi.org/10.1073/pnas.1420441111. Review 28(3), 419-433, http://dx.doi 2– .org/10.1111/j.1728-4457.2002.00419.x. second demographic transition: An interpretation Botev, N. (2012): Population ageing in Central and Eastern Lesthaeghe,of the spatial R., Neels, continuity K. (2002): of demographic From the first innovation to the Europe and its demographic and social context. in France, Belgium and Switzerland. European Journal European Journal of Ageing 9(1), 6 9, http://dx.doi of Population 18(4), 32 60, https://doi .org/10.1007/s10433-012-0217-9. .org/10.1023/A:1021125800070. ESRI (2012): Spatial Statistics Tools, ArcMap9–7 10.1. Lesthaeghe, R., Neidert, L. 5–3(2006): The second Redlands, CA, Environmental Systems Research Institute. demographic transition in the United States: Exception Fassmann, H., Musil, E., Gruber, K. (2014): Dynamic or textbook example? Population and Development historical analysis of longer term migratory, labour Review 32(4), 66 98, http://dx.doi market and human capital processes in the SEEMIG .org/10.1111/j.1728-4457.2006.00146.x. region. SEEMIG Working papers No. 3. Budapest, Muresan, C. (2007): 9–6Family dynamics in pre and post- Hungarian Demographic Research Institute. transition Romania: a life table description. MPIDR http://www.seemig.eu/downloads/outputs Working Paper WP 2007-18. Rostock, Max Planck /SEEMIGWorkingPapers3.pdf. Institute for Demographic Research. https://www Goldstein, R. J., Sobotka, T., Jasilioniene, A. (2009): .demogr.mpg.de/papers/working/wp-2007-018.pdf. The end of lowest-low fertility? Population and Development Review 35(4), 66 99, http://dx.doi Serbian fertility. Stanovništvo 49(2), 7 9, https://doi .org/10.1111/j.1728-4457.2009.00304.x. Nikitović,.org/10.2298/STNV1102073N V. (2011): Functional data. analysis in forecasting 3–6 3–8 Second Demographic Transition and the spatial patterns of fertility change in Serbia 167
of the communist regimes. Population 58(4/5), region of former Yugoslavia. Stanovništvo 54(2), 2 8, 45 86, http://dx.doi.org/10.1007 Nikitović,https://doi.org/10.2298/STNV161115009N V. (2016): Long-term effects of low fertility. in the /s10433-012-0217-9. 7–5 Sobotka,1–4 T. (2008): The diverse faces of the second demographic processes in the region of former demographic transition in Europe. Demographic Nikitović, V. (2017): The impact of migration on Research 19, 17 24, http://dx.doi.org/10.4054 Understanding of Contemporary Migration: Causes, /DemRes.2008.19.8. Yugoslavia. In: Bobić, M., Janković, S. Towards Sobotka, T. (2011):1–2 Fertility in Central and Eastern Europe of Philosophy, University of Belgrade and Serbian after 1989: collapse and gradual recovery. Historical SociologicalConsequences, Society, Policies, 61-81. Reflections. Belgrade: Faculty Social Research 36(2), 24 96, http://dx.doi .org/10.12759/hsr.36.2011.2.246-296. of recent demographic trends in Serbia 196 010. Sobotka, T., Zeman, K., Kantorova,6–2 V. (2003): Nikitović, V., Bajat, B., Blagojević,43. D. (2016): Spatial patterns Demographic shifts in the Czech Republic after 1989: 1–2 A second demographic transition view. European Geografie 121(4), 521–5 Journal of Population 19(3), 24 77, https://doi Nikitović, V., Predojević-Despić, J., Marinković, I. (2015): .org/10.1023/A:1024913321935. Migracijeza statistiku, stanovništva. 9 28. In: Nikitović, V. Populacija Srbije 9–2 O’sullivan,početkom D., 21.Unwin, veka. D. Beograd, J. (2003): Republički Geographical zavod Information support: Sociodemographic research of nonmarital Analysis. New 7–1Jersey, John Wiley and Sons. Stanković,parenthood B. (2014): in Serbia. Lone Stanovništvo mothers and 52(1), their 5network6, https://doi.org/10.2298/STNV1401055S. emigration out of Serbia: Three ‘hot’ emigration zones. Šrocha, B., Šidlo, L. (2016): Spatial differentiation5–7 and Penev,Stanovništvo G., Predojević-Despić, 50(2), 3 4, J. (2012):https://doi.org/10.2298 Spatial aspects of fertility postponement transition in Czechia. AUC /STNV1202035P Geographica 51(2), 21 33, https://doi 5–6 .org/10.14712/23361980.2016.18. extramarital births in Serbia in the second half of Van de Kaa, D. J. (1987): Europe’s7–2 second demographic Penev,the 20thG., Stanković, and at the B. beginning (2010): Characteristics of the 21st century. of transition. Population Bulletin 42(1), 9. Stanovništvo 48(1), 3, https://doi.org/10.2298 Van de Kaa, D. J. (2004): Is the second demographic /STNV1002001P. transition a useful research concept: Questions1–5 Perelli-Harris, B., Gerber,1–2 T. P. (2011): Nonmarital and answers. Vienna Yearbook of Population childbearing in Russia: Second demographic transition Research 2, 0, http://dx.doi.org/10.1553 or pattern of disadvantage? Demography 48(1), /populationyearbook2004s4. 31 42, https://doi.org/10.1007/s13524-010 Walford, N., Kurek,4–1 S. (2016): Outworking of the second -0001-4. demographic transition: National trends and regional 7–3 patterns of fertility change in Poland, and England family patterns in post-socialist countries: Delayed, and Wales, 2002-2012. Population, Space and Place Petrović, M. (2011): Changes of marital behaviour and 22(6), 50 25, http://dx.doi.org/10.1002/psp.1936. Stanovništvo 49(1), 5 8, https://doi.org/10.2298 Watkins, S. C. (1990): From local to national /STNV1101053Pincomplete or specific. second demographic transition? communities:8–5 The transformation of demographic Philipov, D., Kohler, H.-P. 3–7(2001): Tempo effects in the regimes in Western Europe, 187 960. Population fertility decline in Eastern Europe: Evidence from and Development Review 16(2), 24 72, https://doi Bulgaria, the Czech Republic, Hungary, Poland, and .org/10.2307/1971590. 0–1 Russia. European Journal of Population 17(1), 3 0, 1–2 https://doi.org/10.1023/A:1010740213453. 7–6 Data sources 1990s. Stanovništvo 42( ), 7, https://doi Rašević,.org/10.2298/STNV0404007R M. (2004): Fertility trends. in Serbia during the Fertility indicators. Luxembourg, Eurostat Data Base. 1–4 7–2 (cited on May 13, 2019). http://appsso.eurostat. ec.europa.eu Rašević, M. (2015): Fertilitet ženskog stanovništva. In: Nikitović, V. Populacija Srbije početkom 21. veka. OsnovniBeograd, demografski, Republički zavod ekonomski za statistiku, i socijalni 74-95. pokazatelji of/nui/show.do?dataset=demo_find&lang=en. the Republic of Serbia. http://publikacije.stat.gov.rs Rašević,relevantni M., Penev, za populacionu G. (2010): politiku.Opštine republikeBeograd, UNFPASrbije. and Demographic/G2018/PdfE/G201814015.pdf. Yearbook 2017. Belgrade, Statistical Office Društvo demografa Srbije. Population by Age and Sex, Book2, 2002 Census of Population, Households and Dwellings in the Republic in Serbia. European Journal of Contraception and Rašević,Reproductive M, Sedlecky, Health K. (2009):Care 14(6), The 38abortion90, issuehttp://dx.doi of Serbia (In Serbian). http://publikacije.stat.gov.rs .org/10.3109/13625180903215422. /G2002/Pdf/G20024002.pdf.of Serbia. Belgrade, Statistical Office of the Republic Sardon, J-P. (2001): Demographic change5–3 in the Balkans Population by Age and Sex, Book2, 2011 Census since the end of the 1980s. Population: An English of Population, Households and Dwellings in the Republic Selection 13(2), 4 0. Sobotka, T. (2003): Re-emerging diversity: rapid fertility of Serbia. http://publikacije.stat.gov.rs/G2012/Pdf changes in Central9–7 and Eastern Europe after the collapse /G20124002.pdf.of Serbia. Belgrade, Statistical Office of the Republic 168 Original Article
Effect of rural transport infrastructure on the intensification of purchased input use for major food crop production: the case of smallholder farmers in Horro Guduru Wollega Zone, Western Ethiopia Sileshi Tamene*, Tebarek Lika Megento
Department of Geography and Environmental Studies, Addis Ababa University, Ethiopia * Corresponding author: [email protected]
ABSTRACT This paper examines the effect of rural transport on smallholder farmers’ purchased input use. A random sample of 500 respondents was selected and relevant data was collected. Descriptive, correlation, and regression statistics were used to analyze the data. The multiple linear regression analysis revealed that farmers’ purchased input use was found to be significantly and negatively related to distance to major market, distance to all weather road, distance to farm plot, transport cost, and size of land holding. In contrast, farmers’ purchased input use was found to be significantly and positively related to family size, off farm income, membership in a cooperative, being in Horro district, having animal cart, and access to good road. Further, the results of hierarchical multiple regression showed that approximately 82% of the total variation in purchased input use can be explained by the linear combina- tion of all independent variables. Furthermore, the result showed that rural transport infrastructure-related variables, as a set, contributed 13.3% to the prediction of farmers’ purchased input use over and above the remaining predictors. The results suggest that improving the rural road infrastructure and access to rural transportation services is vital in encouraging farmers’ purchased input use.
KEYWORDS hierarchical regression; Horro Guduru; intensification; purchased input use; rural transport
Received: 13 January 2019 Accepted: 10 September 2019 Published online: 19 December 2019
Tamene, S., Megento, T. L. (2019): Effect of rural transport infrastructure on the intensification of purchased input use for major food crop production: the case of smallholder farmers in Horro Guduru Wollega Zone, Western Ethiopia. AUC Geographica 54(2), 168–181 https://doi.org/10.14712/23361980.2019.15
© 2019 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). Effect of rural transport on smallholder farmers’ purchased input use 169
1. Introduction 2011; Weir and Knight 2004). Past agricultural research in the country also focused on the impact of improved agricultural technologies on smallholder remains the greatest challenge facing smallholder farm income and its implication for poverty reduc- farmersThe achievement in developing of major countries. food Thecrop poor intensification state of the - - tionAlthough strategies these (Hailu past 2014; studies Katungi provided et al. useful 2011; infor Sala- rural transportation network and inefficient logis- mationmi et al. on2010). the trends, patterns, and determinants of ticsport continue infrastructure to hinder challenges agricultural are inadequateintensifications and agricultural input adoption, rigorous assessments of inpoor Africa conditions (Delaney of the et al. rural 2017). road The network key rural and limited trans rural transport constraints on the use of purchased availability of vehicles, which has led to an increase (variable) inputs for major food crop production, in the cost of transportation further affecting agricul- as a measure of household commercialization from the input side have rarely been studied in Ethiopia. condition of transport infrastructure, the effect of A better understanding of rural transport constraints geographictural input pricesisolation, (Salami high et al. transport 2010). costs, The very and poortime that hinders smallholder farmers’ participation in lost to roadblocks, can completely undermine the agricultural input markets as a buyer is therefore - important for designing promising pro-poor agri- cultural and transport policies that could stimulate returnsLike into manyinvestments other sub-Saharan in crop intensification African countries, practic the use of modern agricultural inputs and increase small-scalees (Delaney agricultureet al. 2017). is the most important sector small-scale agricultural production. This study aims for achieving sustainable household food security in Ethiopia. It accounts for 80% of the working popula- information to empirically address the relationship tion, 86% of the total foreign exchange earnings, and to fill this knowledge gap and provide quantitative- 48% of gross domestic product (GDP) (Worku 2011). An increasing use of modern agricultural inputs such ofbetween this research smallholder are supposed farmers’ to food be usedcrop byintensifica different as chemical fertilizers and improved seeds remains stakeholderstion and rural involved transport in infrastructure.rural transport The policy, findings agri- one of the best hopes for greater agricultural produc- cultural land use planning, and sustainable food crop tion and productivity of rural Ethiopia, where more production strategies. The overall objective of this study was to investi- However, lack of enhanced supply and promotion of gate rural transport constraints of major food crop improvedthan 80% seeds,of the populationorganic and lives inorganic (Minten fertilizers, et al. 2013). and low level of an irrigation system are major obstacles to sustain the agricultural production in the country to:intensification (1) examine strategies relationship by small-scalebetween proximity farmers toin all-weatherEthiopia. The roads specific and total objectives values of purchasedthis study inputwere As a landlocked country with largely non-naviga- use among smallholder farmers; (2) determine the (Elias et al. 2006). extent to which a combination of rural transport the performance of the Ethiopian economy. In Ethi- infrastructure, institutional factors, resource endow- bleopia rivers, immature road ruraltransport transport plays and a significant other key rolephysi in- ment, and physical factors predict smallholder farm- cal infrastructure have led to high transport costs for ers’ purchased agricultural input use; (3) identify agricultural products to the market as well as of farm the extent to which rural transport infrastructure inputs, reducing farmers’ competitiveness (Fufa and (distance to major market, distance to all-weather Hassan 2006; Lulit 2012). This could be a disincen- road, distance to farm plot, transport cost, mode of tive to the use of productivity-enhancing agricultural transport and road conditions) predicts smallholder inputs and therefore could discourage smallholder farmers’ purchased agricultural input use, controlling for the effects of demographic, institutional, resource Most of the literature on small-scale agriculture endowment, and physical factors. farmer’smarket participation major food crop concerns intensification. only the output side of marketing production (Arethun and Bhatta 2012; 2. Theoretical underpinnings of the study However, the sustainable marketing of smallholders alsoBekele requires et al. 2010; integration Gebremedhin into the and input Hoekstra markets 2008). (Pin- This study aims to look at the nature and extent of gali and Rosegrant 1995). To bridge the gap in the lit- rural transport infrastructure and its effect on mod- erature on the marketing participation of rural house- ern agricultural input use among smallholder farm- hold on the input side, we analyze the rural transport ers of Horro Guduru Wollega Zone, Western Ethiopia. determinants of purchased modern agricultural There are several theories that attempted to explain inputs use. In Ethiopia, there is relatively large litera- how rural transport infrastructure investment can ture dealing with the role of technological innovation bring about economic growth and development and diffusion in increasing agricultural productivity Kachallah 2017; Margarian 2011; Roland-holst 2009). (Banerjee et al. 2012; Didenko et al. 2017; Jelilov and and intensification (Elias et al. 2013; Katungi et al. 170 Sileshi Tamene, Tebarek Lika Megento
Of these multiple theories, “the theory of induced 3. Data and Methods technical and institutional change” was used as rel- evant theoretical perspective to design the research 3.1 Selection and description of study site questions that this study is based on. In addition, this The study was conducted in Horro Guduru Wollega theoretical perspective was used for organizing and Zone, western Ethiopia. This Zone lies between Lat- induced technical and institutional change provides interpreting the findings of this study. The theory of itude 9°10′2 (CSA N and 2011; 9°50′ Tamene N and andLongitude Megento 36°00′ 2017). E will philosophically, epistemologically, methodologi- andShambu 36°50′ is theE (Figure capital 1). town It has of thea total zone land and area found of callythe structure and analytically to define approached. how this particular research 8,097 km Economic historians are increasingly drawing on of CSA (2011), this zone had a total population of the theory of induced technical change in attempt- 641,575314 km West of which of Addis 50.09% Ababa. are According male and to 49.91% the report are ing to interpret differential patterns of productivi- female. According to the same source, about 89% of ty growth among countries and over time (Ruttan the population lives in the rural areas driving their 2008; Ruttan and Hayami 1984). Agricultural econ- livelihoods from agriculture. omists like Hayami and Ruttan (1993, p. 6) argue by saying that in agriculture, changes in the ‘relative resource endowments’, especially land and labor, induce a derived demand for technological innova- tions to facilitate the replacement of relatively less scarce and cheap factors for more scarce and expen- sive ones. For example, when labor is in short sup- ply, there is a tendency for capital in the form of labor-saving machinery to be substituted for human labor. Whereas, in a land-scarce economy, yield-in- creasing and land enhancing inputs such as fertilizers and improved seeds are substituted for land which in turn depends on the agricultural input market condi- tions (Hayami and Ruttan 1985, 1993; Ruttan 2008). Moreover, roughly 22 years ago, Ruttan (1996, p. 54) in one of his seminal papers addressed that induced technical change acts to make the ‘scare factor more abundant’. Induced technical and institutional change theo- Fig. 1 Location map of Horro Guduru Wollega Zone. ry offers a theoretical understanding appropriate to Source: Adapted Finance and Economic Development Bureau examine the complex and dynamic relations between of Oromia, 2016. rural transport access and smallholder farmers’ pur- chased input use. When we examine the appropri- The average annual temperature in the study area ateness of this theory for understanding the nature is 22.1 °C, with an average minimum of 13 °C and an of rural transport and purchased input use in the average altitude of Horro Guduru Wollega Zone rang- effort of promoting the rural agricultural econo- esaverage from 860maximum to 2657 of meters30 °C (Beyene above sea et al. level 2015). (Beyene The study area, two related realities emerge: first, the farmers needs technical change and innovations to mainstay in the study area with notable food crops transformmy and bringing the most maximum common benefits and tedious to smallholder traditional includinget al. 2015). wheat Mixed (Triticum crop-livestock aestivum agriculture), barley (Horde is the- rural transport mechanisms-human porterage (head, um vulgare), teff (Eragrostis tef), maize (Zea mays), shoulder, and back-loading) to improved rural trans- pulses (Vicia faba, Pisum sativum) and cash crops like port means-pack animals and animal drawn carts. sesame (Sesamum orientale), niger (Guizotia abyssi- Such rural transport improvements brought about nica), and linseed (Linum usitatissimum) (CSA 2014). by rural transport innovations can enhance the use Even though few independent variables used in of purchased input use among smallholder farmers. our previously published article (Tamene and Megen- Second, this theory also acknowledges that improved to 2017) are also used in current manuscript, a great- production technologies (fertilizer and improved er number of independent variables are newly intro- seed) as well as improved farm management prac- duced to the current manuscript to emphasize on the tices (credit and extension institutions) can play in originality of this paper. Furthermore, due the follow- replacing relatively scarce resources like land and ing reasons each manuscript has a distinct focus and labor. purpose. First, each manuscript addressed different Effect of rural transport on smallholder farmers’ purchased input use 171 research questions. Second, each manuscript studies software program version 20. To determine wheth- the data from completely different angles. Third, each manuscript used different relevant literature. Fourth, to the prediction of smallholder farmers’ purchased they both differ in their analytic methods, interpreta- inputer rural use, transport over and infrastructure above the variance added predictedsignificantly by tions, and conclusions. Lastly, the dependent variables demographic, household resource endowment, insti- in both manuscripts are conceptually different and - empirically not related. Therefore, the manuscripts bles, the independent variables were entered in to should be considered independently. twotutional separate infrastructure, blocks of the and regression location specificanalysis varia(hier-
3.2 Study Design, Sampling, Data Collection, variables except rural transport infrastructure relat- and Analysis edarchical variables method). (variables In the of interest)first block, were all includedindependent into The study utilized a smallholder household based the regression analysis. At the second step, the new cross-sectional quantitative survey design using independent variables (rural transport infrastructure a structured questioner with face to face interview. related variables) were added, and all of the independ- Descriptive and analytical cross-sectional micro-lev- el data have been used to estimate the effect of rural the independent set to see the unique contribution of transport infrastructure on smallholder farmers’ pur- ruralent variables transport (those infrastructure entered at related first block) variables remain in the in chased input use. Horro Guduru Wollega Zone was - ing corridors of Ethiopia. On the contrary, the existing intensification3.3 Description of purchasedof Variables agricultural input use. ruralidentified road astransport one of the infrastructure potential cereal in the crop zone produc is not Our variables of interest (both dependent and inde- satisfactory to support smallholder agricultural pro- pendent) used in this study and their levels of meas- duction system including the purchases of necessary urement are shown in table 1. These variables are inputs (Tamene and Megento 2017). Therefore, Horro Guduru Wollega Zone was purposively selected. independent variables on the purchased input use as A multistage simple random sampling technique asupposed dependent to variable.capture the influence of the potential was sequentially employed to select four districts from Horro Guduru Wollega Zone, four rural kebeles 3.3.1 Rural transport infrastructure (RKs) from each district, and rural farm households It is assumed that a well-functioning rural trans- from each rural kebele-Kebele is the smallest admin- a random selection of four districts from the nine smallholderportation infrastructure farmers. Smallholder is significant level determinants purchased districtsistrative of unit Horro in Ethiopia.Guduru Wollega The first Zone. stage As involvesa result, inputof the use form for and major pace food of cropfood production crop intensification (Y) is there of- four districts (Hababo Guduru, Horro, Amuru and Abe fore modeled as a function of smallholder farmers’ Dondoro) were selected (Figure 1). The second stage access to rural transport infrastructures and services. involves the random selection of four RKs from each Access to rural transport infrastructures and servic- of the four districts making a total of 16 RKs. es, captured as an average distance to major market, - average distance to farm plot, average distance to all tion of farm households with farmland size of 0.25 ha weather-road, are expected to be negatively correlat- andThe above third from and each final RK. stage The waslist ofthe farm random households selec ed to the total values of purchased input use (Bekele in each RK was compiled with the assistance of the extension agents and RK managers. According to Gray road is considered crucial to improving access to agriculturalet al. 2010). inputAvailability markets, of resultinggood quality in greater rural access use of productivity-enhancing modern agricultural inputs et al. (2007) suggestion, the researcher used a 95%- confidenceably in the numberlevel to determine of smallholder the sample farmers size they for con this- of Transportation for local-level transport services specifictain and study. hence It to is avoid usual bias, that probability RKs may vary proportional consider in(Jayne rural et al. areas 2003). is also Ownership expected of to Intermediate promote increased Means to size (PPS) was employed. Thus, 500 smallholder use of purchased modern inputs. Hence, household farmers from the four districts were sampled for the level purchased input use is modeled as a function of study. distance to major market, distance to farm plot, dis- The household survey was conducted from Feb- tance to all weather-road, transport cost, ownership ruary to June 2016, which followed shortly after the of Intermediate Means of Transportation and road main season (Meher) harvest. Interviews were con- quality. ducted in places convenient to farmers either at home 3.3.2 Demographic factors Statistical analysis such as descriptive statistics, To capture the effects of demographic factors on correlation,or in the field. and multiple regression with a hierar- smallholder farmers’ purchased input use, we used the age of household head and family size. It was Statistical Package for Social Sciences (IBM SPSS) assumed that older households tend to be more chical model specification was performed using the 172 Sileshi Tamene, Tebarek Lika Megento
Tab. 1 Summary of variables used and their measurements.
Variable name and type Description Measurement Expected sign A dependent variable indicating the total amount of birr spent for Purchased input use Ethiopian Birr* purchasing agricultural inputs for major food crops Age of household head Age of household head Number of years − (continuous) Family size (continuous) Total family size of the household head Number + Size of land holding (continuous) Area of farm land owned by the household Hectare + Livestock ownership Total livestock unit Total livestock ownership of the household −/+ (continuous) (TLU) Off farm income (continuous) Income earned from non-agricultural activities Ethiopian Birr* + Membership in a cooperative Being a member of an agricultural cooperatives 1 = Yes; 0 = No + (Dummy) Extension visit monthly 1 = Monthly visit; Frequency of extension visit + (Dummy) 0 = Twice in a year Level of annual rainfall Amount of annual precipitation Millimeter (mm) −/+ (continuous) Distance to major market Distance travelled by the household to reach the nearest major Kilometer − (continuous) market Distance to farm plot Average farm plot distance from the homestead Kilometer − (continuous) Distance to all weather road Distance travelled by the household to reach the nearest all Kilometer − (continuous) weather road Transport cost incurred to move 100 kg of agricultural input Birr per 100 kg Transport cost (continuous) − over 1 km per km 1 = Animal cart; Animal cart (Dummy) Transport mode used + 0 = Headloading 1 = Good road; Good road (Dummy) Road quality + 0 = Bad road 1 = Highland; High land (Dummy) Agro-ecology type −/+ 0 = Otherwise 1 = Horro district; Horro (Dummy) District type − 0 = Otherwise Source: Own construction. * During data collection 1 USD equals 23.73 Ethiopian birr. market-oriented and have higher agricultural input 3.3.4 Institutional infrastructure variables market participation. The number of members in Apart from household resource endowment and the household is also assumed to be very important demographic factors several other factors such as the in the use of purchased input use since the use of institutional frameworks affect the demand for and such productivity-enhancing inputs are said to be consumption of purchased farm inputs such chemical more labor-intensive than conventional subsistence fertilizers (DAP and urea), improved seeds (maize and farming. wheat) and pesticides. Therefore, the third impor- tant set factors affecting purchased input use among 3.3.3 Household resource endowment smallholder farmers in rural areas include the insti- Off-farm income, livestock ownership and size of tutional infrastructure variables such as membership landholding increase the welfare of farmers because they help farmers in getting the required input need- 2009) and contact with agricultural extension agents (Belayin a cooperative 2015). (Hellin et al. 2009; Markelova et al. throughed for agricultural cash from production off-farm activities (Abdullah and et al. livestock 2017). 3.3.5 Physical and location specific factors salesPurchased (Christiaensen agricultural 2017). inputs Therefore, are mainly household financed Several environmental variables were hypothesized resource endowment (size of land holding, livestock to encourage/discourage farmers to invest in pur- ownership, and off farm income) are hypothesized to chased agricultural input use for major food crop increasingly recognized as a resource that can signif- production. These include the amount of annual rainfall received and agro-ecology. The total values among smallholder farmers. of purchased input use and physical factors (amount icantly influence the quantity of purchased input use Effect of rural transport on smallholder farmers’ purchased input use 173
of annual rainfall received and agro-ecology where Y 0 1X1 2X2 16X16 a household belongs) are related. It is expected that, on average, smallholder farmers in midland Here = β in + theβ case + β of both + … +models: β + ε (2) Y more purchased inputs. To capture the differences in inputs (Ethiopian birr), purchasedagro-ecology input and usewith among sufficient study rainfall districts, tend district to use 0is is the the financialconstant costterm or or monetarythe intercept, value of variable dummy was considered. β respective independent variables 1–16 3.3.6 Rural transport infrastructure βX1 is ageare ofthe household regression head coefficients (years), associated with It is assumed that a well-functioning rural trans- X2 is family size (number), X3 is size of land holding (hectare), X4 is livestock ownership (TLU), smallholderportation infrastructure farmers. Smallholder is significant level determinants purchased X5 is off farm income (Ethiopian birr), inputof the useform for and major pace food of food crop crop production intensification is there of- X6 is membership in a cooperative (dummy), fore modeled as a function of smallholder farmers’ X7 is extension visit monthly (dummy), access to rural transport infrastructures and services. X8 is level of annual rainfall, Access to rural transport infrastructures and servic- X9 is high land (dummy), es, captured as an average distance to major market, X10 is Horro (district dummy), average distance to farm plot, average distance to all X11 is distance to major market (km), weather-road, are expected to be negatively correlat- X12 is distance to all weather road (km), ed to the total values of purchased input use (Bekele X13 is distance to farm plot (km), X14 is transport cost (Ethiopian birr), road is also considered crucial to improving access to X15 is animal cart (Dummy), agriculturalet al. 2010). inputAvailability markets, of resultinggood quality in greater rural access use of X16 is good road (Dummy), productivity-enhancing modern agricultural inputs - of Transportation for local-level transport services in ε isrelationship. the random error component reflecting the dif rural(Jayne areas et al. is 2003). also expected Ownership to promoteof Intermediate increased Means use ference between the observed and fitted linear of purchased modern inputs. The above described independent variables were and computing the parameter estimates some pre- entered in to into a hierarchical linear regression dictorAfter variables fitting a arehierarchical omitted since linear they regression are much model less analysis through sequential block-enter approach so important or most likely affect the explanatory power as to predict smallholder farmers’ purchased input of the model if included. use. The choice of how to include independent varia- bles was determined by researchers based the overall - 4. Results cal linear regression analysis, demographic variables, purposehousehold of theresource analysis. endowment In the first variables,block of hierarchi institu- 4.1 Preliminary analyses tional infrastructure variables, and physical and loca- There are a number of assumptions that must be met for multiple linear regression model to be reliable of independent variables. (Osborne and Waters 2002). Preliminary analysis to tion specific variables were entered as the first block ensure the non-violation of the assumptions of nor- Y 0 1X1 2X2 10X10 mality (Kim 2015; Miot 2017; Yap and Sim 2011), linearity (Osborne and Waters 2002), homoscedas- = βSince + β we are + β also +interested … + β in + examining ε the effect(1) ticity and multicollinearity (Daoud 2017; Friday and of rural transport infrastructure variables on small- - holder farmers’ purchased input use after demo- pleted prior to the analysis. The preliminary analysis graphic variables, household resource endowment revealedEmenonye that 2012; these Imdadullah assumptions et al. were 2016) not were seriously com variables, institutional infrastructure variables, and violated. - trolled for, we entered rural transport infrastructure 4.2 Descriptive statistics variablesphysical and (variables location specificof interest) variables in the have subsequent been con Descriptive statistics were used to summarize data. blocks of independent variables in the hierarchical Table 2 shows a descriptive analysis of continuous linear regression analysis (block 2). Therefore, to see variables by presenting numerical facts about the if the rural transport infrastructure variables predict quantitative dataset. The total number of observa- smallholder farmers’ purchased input use above and tions (n) was 500 smallholder farmers in four dis- beyond the effect of the controls and to test if succes- tricts of Horro Guduru Wollega Zone, Western Ethio- pia. The mean score for household size was 6.48, with model was developed. a standard deviation of approximately 3.31 points, sive model fit better than previous one the following 174 Sileshi Tamene, Tebarek Lika Megento
Tab. 2 Descriptive statistics for continuous variables (n = 500).
Minimum Maximum Mean Std. Deviation Skewness Kurtosis Variable Statistic Statistic Statistic Statistic Statistic Statistic Age of household head 17 64 34.29 13.02 .56 −.77 Family size 2 12 6.48 3.31 .19 −1.29 Size of land holding .4 6.0 2.41 1.18 .45 −.48 Livestock ownership (TLU) .55 11.51 3.26 1.64 .95 1.59 Off farm income 0 9100 2416.30 2730 .89 −.68 Distance to major market 5 35 20.18 7.87 −.15 −1.16 Distance to farm plot .4 13.0 6.41 3.38 −.08 −1.07 Distance to all weather road 0 27 12.97 6.76 −.16 −.79 Transport cost .50 3.20 2.18 .78 −.53 −.74 Purchased input use 717.50 37644.43 10096.92 9357 1.44 1.09 Source: Compiled from field data, 2016.
4.3 Correlation analysis persons (CSA 2017). - figuresThe averagewhich are distance above to the the national nearest average major market of 4.6 - uesThe of strongest purchased negative input use significant was with the Pearson proximity prod to smallholders in the study area, while the average dis- all-weatheruct-moment roads: correlation r coefficientp for< .001. the totalWhereas val is about 20.18 km, indicating poor market access of the strongest positive Pearson product-moment cor- skewness and the kurtosis indices are very small and (498) = −.772, falltance within to all-weather the acceptable roads rangeis 12.97 km. which Theindicates values that for input use was with off-farm income: r prelation coefficient for the total values of purchased- cases or outliers (see Table 2). mination (r2 (498) = .654, theThe variables average most purchased likely do input not (chemical include influential fertilizer, weather < .001 road(Table explains 3). The 59.6% simple of coefficientthe variation of in deter total improved seed, and herbicides) value used for major values of purchased = 0.596) input indicated use by smallholder that distance farmers. to all food crop production (maize, wheat, and teff) is ETB Approximately, the other 40% of the total variance 10096.92. The result also indicates that on average between distance to all weather road and purchased smallholder farmers get about ETB 2416.30 income input use remains unexplained. This value of simple from off-farm employment. A household on average operates about 2.41 ha (SD relationship between distance to all weather road andcoefficient purchased of determination input use in the shows absence the ofbinary other linear inde- (CSA 2015; Tamene and Megento = 1.18), 2017). a result Finally, which on is pendent variables. Thus distance to all weather road average,two times households greater than incur the national a 2.18 averageETB to transportof 1.14 ha purchased input use. center to home. seems to explain a significant amount of variation in one kg of farm input for 100 km from input market
Tab. 3 Pearson correlations of purchased input use (dependent variable) and independent variables.
Independent Variables Correlation coefficient (r) P value Age of household head (years) −.419** .000 Family size (no.) .571** .000 Size of land holding (ha) −.416** .000 Livestock ownership (TLU) .263** .000 Off farm income (ETB) .654** .000 Distance to major market (km) −.745** .000 Distance to farm plot (km) −.467** .000 Distance to all weather road (km) −.772** .000 Transport cost (km) −.727** .000 Level of annual rainfall (mm) .114* .011 n = 500, ** P < .001 (2-tailed), * P < .05 (2-tailed) Source: Compiled from field data, 2016. Effect of rural transport on smallholder farmers’ purchased input use 175
4.4 Hierarchical multiple regression model and disturbance term. The adjusted R2 value is A multiple regression analysis was conducted to 0.816, which is very close to the multiple R2, indicat- examine whether rural transport infrastructure, ing that we shouldn’t worry too much about whether institutional factors, resource endowment, and phys- we are using too many variables in the model. farmer’s purchased input use. The results of multiple 4.4.1 Rural transport infrastructure regressionical factors analysis could significantly using all sixteen predict independent smallholder var- iables are summarized in table 4 of model 2. A strong displayed in table 4 showed the relative importance relationship between purchased input use and the ofThe each estimated predictor unstandardized in the model. regression Among rural coefficients trans- independent variables was observed and the model port infrastructure variables, transport cost was the strongest negative predictor as indicated by its (R2 F p < .001). was a significant predictor of purchased inputR2) valuse- p < .001) followed by good road dum- ue of = .822, 0.822, indicates(16, 483) = 138.95, that approximately 82.2% of the estimated unstandardizedp < .001) and regression animal cart coefficient, dummy totalThe variation multiple in coefficient purchased of input determination use can be account( - (β = −3707.88,p < .05). The unstandardized regres- ed for by the linear combination of explanatory var- my (β = 1832.57, iables. But the remaining 17.8% of the variance has to(β = 1090.49, major market, distance to all weather road and been attributed to other variables not included in the distancesion coefficients to farm for plot the on association one hand between and smallholder distance
Tab. 4 Hierarchical egression for variables predicting purchased input use (n = 500).
Unstandardized Coefficients Standardized Coefficients Model t Sig. B Std. Error Beta (Constant) 22750.229 11824.431 1.924 .055 Age of household head −98.328 23.189 −.137 −4.240 .000 Family size 651.186 102.507 .230 6.353 .000 Size of land holding −2572.812 288.302 −.325 −8.924 .000 Livestock ownership (TLU) 577.724 205.493 .101 2.811 .005 1 Off farm income 1.240 .113 .362 10.989 .000 Membership in a cooperative 2801.652 602.726 .145 4.648 .000 Extension visit monthly_Dummy 571.666 542.847 .027 1.053 .293 Level of annual rainfall (mm) −8.864 7.432 −.038 −1.193 .234 High land_Dummy 2078.012 725.521 .100 2.864 .004 Horro_Dummy 1046.011 611.352 .044 1.711 .088 (Constant) 22594.157 10233.971 2.208 .028 Age of household head −7.988 18.552 −.011 −.431 .667 Family size 189.180 82.571 .067 2.291 .022 Size of land holding −1070.479 246.050 −.135 −4.351 .000 Livestock ownership (TLU) 38.776 168.814 .007 .230 .818 Off farm income .812 .092 .237 8.848 .000 Membership in a cooperative 1457.595 467.319 .076 3.119 .002 Extension visit monthly_Dummy 722.301 424.473 .035 1.702 .089 2 Level of annual rainfall (mm) 1.000 6.519 .004 .153 .878 High land_Dummy −269.070 574.761 −.013 −.468 .640 Horro_Dummy 3379.733 492.646 .142 6.860 .000 Distance to major market −160.575 46.674 −.135 −3.440 .001 Distance to all weather road −316.427 50.336 −.229 −6.286 .000 Distance to farm plot −215.478 61.132 −.078 −3.525 .000 Transport cost −3707.879 380.257 −.310 −9.751 .000 Animal cart_Dummy 1090.493 548.139 .044 1.989 .047 Good road_Dummy 1832.569 457.855 .096 4.003 .000 a. Dependent Variable: Purchased input use Note: R2 = .689 for model 1, R2 = .822 for model 2 and 2 = .133. Source: Compiled from field data, 2016. ΔR 176 Sileshi Tamene, Tebarek Lika Megento farmers’ purchased input use on the other hand are 5. Discussion - - 5.1 Family size and purchased input use −160.57,cients are −316.43 46.67, 50.34 and −215.48 and 61.13 respectively; respectively. the asso The results of hierarchical multiple regression ciated standard errors for these regression coeffi revealed a positive relationship between family size 4.4.2 Household demographics and purchased input use indicating that large families The results of the regression model showed that spend more on purchased input use as compared to - - ings that explain how the number of family mem- p family size is an important factor identified to influ small families. There are many prior research find ence purchased input use (β = 189.18, t(485) = 2.29, smallholder farmers. For example, Nambiro (2008) 4.4.3= .02). Household resource endowment bers’ influences purchased input use decisions of the proportion of farm allocated to the cultivation t p < .001. Off- offound improved a significant hybrid positive maize seed.effect Another of family study size onby farmSize of income land holding is positively significantly associated predicted with total purchased values Kamara (2004) and Perz (2003) found availability ofinput purchased use, β = −1070.48, input use such (485) = −4.35, that, holding everything of family labor as a precondition for greater use of else constant, for each additional ETB off-farm chemical fertilizers, pesticides, and improved seeds. income, the total values of purchased input use is pre- Therefore, labor-augmenting technologies are impor- dicted to increase by .81 ETB, and this association is tant to encourage small families so that they can take part in input market participation. As a result, this statistically4.4.4 Institutional significant infrastructure (p < .001). variables that agricultural equipment designed for use in small As can be seen in table 4, membership in a coopera- farmfinding plots supports make inducedit feasible technical for farmers innovation to shift theory from labor-intensive practices to higher-yielding mecha- indicating smallholder farmers who are members of nized practices (Hayami and Ruttan 1993). tivecooperative had a significant associations positive were expectedregression to coefficient, invest ETB 1457.59 more than the nonmembers, after controlling 5.2 Household-level asset variables and purchased for the other variables in the model, and this result is input use p 5.2.1 Landholding 4.4.5statistically Location significant specific ( factors/district = .002). dummy - The value associated with being a farmer living in cant relationship between farm size and smallholder Horro district is ETB 3379.73 (adjusting for the oth- farmer’sThis study purchased found a inputnegative use. and Past statistically studies have signifi found a mixed result. For instance, the largest share of house- dummy variable is both positive and statistically sig- holds renting mechanization is more likely related to er variables in the model), and the coefficient on this . In
4.4.6nificant The (p unique < .001). contribution of rural transport infra- rolelarge of farming increased size farmland (Diao et al. values 2016; on Ma the et al. level 2018) of sup- structure related variables to purchased input use plementaryaddition, Ma feedet al. used(2018) for investigated dairy production the determinant. Kiplimo The percent of the variability in the purchased input use that can be accounted for by all the predictors a negative relationship between a continuous reduc- tionand inNgeno farm (2016)size (farm and fragmentation) Hung et al.(2007) and farm also house found- and hence is an excellent model. The change in var- together is 82.2%. ThisR is2 a significant contribution (2015) found the use of seed and fertilizer technolo- gieshold tolevel be inputscale-neutral, use. In contrary which tois thesetheir usefindings, intensity FAO ianceF(6, 483) accounted forp < (Δ .001.) wasIn this equal case, to the .689 − .822 percent- does not depend on farm size. FAO further underlined age= .133, of variabilitywhich was in significantly purchased inputdifferent use accountedfrom zero that, since fertilizer is ‘a land augmenting input’, small- for went up= 59.76, from 68.9% to 82.2%. The F change asso- holders use it intensively, probably to substitute for ciated with R2 - showing that adding rural transport infrastructure ing supporting the theory of induced technical inno- variables to the change model of .133 increases is statistically the model’s significant pre- vationland (p. which 10). Our advocates finding goesthe needin line to with substitute FAO’s find the dictive capacity. This is to mean that rural transport relatively less scarce and cheap factors of production infrastructure explained additional13.3% of the vari- (fertilizer) for more scarce and expensive ones (land) ance in purchased input use, after controlling for the (Hayami and Ruttan 1985;Hayami and Ruttan 1993). possible effects of potential confounding variables 2 p < 0.001]. The 5.2.2 Off-farm income β), intercept, The study revealed that income from off-farm sourc- [∆R = 0.133, ∆F (6, 483) = 59.764, Beta), - forunstandardized the full model regression are reported coefficients in table 5.( mentary inputs such as fertilizer, pesticides, and high and the standardized regression coefficients ( es positively influenced the application of comple Effect of rural transport on smallholder farmers’ purchased input use 177 yielding seed varieties. This is because off-farm income was an important source for smallholder farmers to increase their market access to agricultural inputs to greaterThis finding access also to corroboratesurban markets the increases ideas of theNin-Pratt inten- intensify production. Findings of this study are in line sity(2016) of input and Katungi use and et al. productivity (2011), who in the suggested rural sector. that off-farm income earning opportunities drive small- 5.4.2 Distance to farm plot with those of Dahal et al. (2007), who concluded that As expected, the regression result showed that plot (2004) also found that farmers’ access to adequate and distance from the homestead has a negative and sig- holders towards agricultural intensification. Kamara their use of modern agricultural inputs. Moreover, in input use. It may be the case therefore that the more sustainable off-farm income has a significant effect on- remotenificant therelationship farm plot with from total farmer’s values residence, of purchased the lesser would be the probability of purchased agricul- atheir source seminal of income work, for Lim-Applegate Australian farm et al. families. (2002) point tural input utilization. This result is in line with the ed out the significance of off-farm employment as 5.3 Cooperative membership and purchased input - use findings of a great deal of the previous work in this Membership in agricultural cooperatives is among the field.distance For from instance the homestead Hailu et al. and (2014) probability found ofa statischem- variables that determine the propensity of smallhold- icaltically fertilizer significant adoption negative decision. relationship By using between a dataset plot er farmer’s participation in the agricultural input mar- ket. In this research, it was expected to have a positive located adjacent to the homestead are more likely to - from Ghana, Kotu et al. (2017) also found that plots- es than the more distant ones. ofinfluence. participation Accordingly, in agricultural membership input marketing. in farmer Therecoop adopt sustainable agricultural intensification practic areeratives several was points found that to significantly help us to maintain influence this the view.level 5.4.3 Distance to all weather road Firstly, agricultural cooperatives play a pivotal role in The current study found a negative and statistically subsidizing fertilizer and seed distribution. Secondly, participating in farmer organizations has the poten- road and total values of purchased input use. It is evi- tial to secure better prices for produce. And thirdly, dentsignificant that the association cost of transport between is distance determined to all by weather taking they also play a key role in improving farmer’s access account of road roughness and seasonality. All weather to technical advice. The results thus obtained are road was reported from earlier studies to be an impor- compatible with previous studies. For instance, in tant variable which explains variations in purchased - input use. For example, in Madagascar, Ninnin (1997) ed that membership to cooperative society was the reported that dry season fares were less costly than their research findings, Birachi et al. (2011) indicat wet season fares. By using data drawn from longitudi- among food crop farmers in Burundi. Furthermore, nal Ethiopian Rural Household Surveys, Wondemu and significant driver of agricultural commercialization Weiss (2012) also reported that improving the class relationship between participation in pools (cooper- of rural roads to a degree that allows all-weather road atives)Carrer and et al. the (2018) adoption found of forward significant contracts and positive among access sharply increases average household income. citrus growers in the State of Sao Paulo, Brazil. They further established that with the equal level of farmland ownership, having paved road access allows 5.4 Rural transport infrastructure and purchased a smallholder farmer to generate 82% higher income input use than would be the case with poor access road. Using cross-sectional data, Beshir (2014) also examined the 5.4.1 Distance to major market factors that affect the probability of improved forage Access to agricultural input markets is expected seeds adoption in two districts of South Wollo zone, to be negatively correlated with the total values of - purchased input use. It is, therefore, hypothesized cant relationship between distance to all weather road that reduced distance to major market will positive- andEthiopia. the probability He found aof negative adoption and and statistically intensity of signifi use of ly affect smallholder farmers’ purchased input use. improved forage seeds. association between distance to the nearest major 5.4.4 Transport cost marketAs expected, and total there values was of a purchasednegative and input significant use. The - ers with reduced physical distance to input markets totalThe coefficientvalues of purchased of transport input cost use (seeincurred Table had 4). This the havecurrent a higher results probability confirmed theof using hypothesis modern that agricul farm- isexpected because, negative on the whole, sign and it has significant been established effect on that the tural inputs than those who are remote (Hailu and there is a strong correlation between transport cost Fana 2017). This is partly because the costs of obtain- incurred and the ability of smallholder farmers to ing agricultural inputs such increase more quick- purchase and use modern agricultural inputs. This ly with increased distance to input supply centers.
finding is consistent with previous studies of Kotu 178 Sileshi Tamene, Tebarek Lika Megento
variables are controlled, a good road condition, as input markets characterized by high transaction and opposed to bad road condition, will result in an addi- transportet al. (2017), costs, who in reported Ghana, farmersthat because mostly of inefficientpay high- tional ETB 1832.569 investment in purchased input use by smallholder farmers. This current study fur- potassium (NPK) fertilizer. In another study in Aus- ther showed that a good road system from input mar- tralia,er than Freebairn official prices (2003), for reportednitrogen, thatphosphorus Australians and ket to the farms would allow easy and timely access who live in remote rural areas will incur additional to purchased inputs. This might also partly explain transport costs to get access to some services offered why good-quality roads (paved roads) provide a good in large urban areas. Similarly, according to Wonde- mu and Weiss (2012), high input prices due to lack rural economy. of infrastructure, such as underdeveloped rural road stimulusA study to by farm Quan profitability (2009) showed and productivity that good physiof the- networks have led to high transport costs for farm cal connections to input markets are a fundamental inputs, thus holding back farmers’ demand for pur- enabler for smallholder farmer’s purchased input chased input use. use. He further showed that good road access being paramount for farm input commercialization through 5.4.5 Mode of transport new technologies that increased the yield of basic The survey showed a marked variation in the pur- food crops. Another recent study in Kenya and Tan- chased input use among smallholder farmers by type of transport mode owned and used. This research poor accessibility challenges that smallholder farm- found that an improvement in the mode of rural erszania experience by Bradbury in getting et al. (2017) agricultural explicitly inputs revealed from the transport use from head-loading to animal cart will market to the farms. In their examination of the trans- result in an additional 1090.493 birr investment in port costs and access constraints for well connected purchased input use by smallholder farmers. The and remote rural farmers of Kenya and Tanzania, they supply of agricultural inputs was expected to increase found that smallholders who are linked by a network substantially with the increased probability of mod- maintained and mostly impassable during the rainy used can seriously affect access to farm inputs. As seasonsof unclassified, are less likelyearth toaccess use of tracks purchased that inputsare poorly than aal result, shift. The with overall respect efficiency to transport of the mode transport owned mode and those who are connected to good road networks link- used, modern input use was expected to be higher ing farming areas to major markets. for smallholder farmers who owned animal cart as compared to those who use the various methods of human porterage (head, shoulder, and back-loading). 6. Conclusion and policy implications the most deriving factor for input market participa- It is incontestable that rural road connectivity and tionOur findingwhere the reveals difference that ownership in purchased of pack input animals use was is rural transport services are among the key compo- seen among those who own and not. nents for rural development, as it promotes access to economic and social services, stimulating the demand access to an improved mode of transport that helps to and consumption of purchased agricultural inputs movePast farm empirical inputs from findings input suggesteddelivery center that to limited home- that in turn enhance production and productivity of stead and from homestead to farm remains a major the farmers. To this end, this study revealed that size challenge for smallholder farmers (Hine 2004). of land holding, distance to major market, distance According to Zewdie (2015), those households who to all-weather road, distance to farm plot, transport face binding transport constraints may be unsuccess- ful to afford the maximum desired levels of input use. purchased input use. Whereas family size, off-farm The traditional mode of transport like human porter- income,cost have and a membershipsignificant negative in a cooperative relationship are found with age negatively affected the level of input use for agri- cultural production in Nigeria (Akramov 2009; Orak- farmers’ purchased input use. to beTwo significant important and policy positively implications related toemerged smallholder from as fertilizer and improved seeds is very low among the study: First, to free smallholder farmers from thosewue et al. smallholders 2015). Yet, who the do use not of ownimproved transport inputs, animals such a vicious cycle of subsistence production, policy as compared to those who own the same. In Ethiopia reforms in the area of rural infrastructure, access to too, transport mode choice is said to be an increas- input markets and to credit facilities must be the cen- ingly important area in getting access to agricultural tral government and local government’s rural devel- inputs for smallholder farmers (Kassa 2014). need a close policy follow-up so as to enhance produc- 5.4.6 Road Condition tionopment and top productivity priority. Second, of farmers. input Hence, use intensification a policy-mix Linking smallholder farmers with a good road net- that can increase smallholder farmers’ off-farm income is desirable as income is a critical predictor of to and utilization of purchased farm inputs. If all other improved seed and fertilizer use. work was found to be positive and significant in access Effect of rural transport on smallholder farmers’ purchased input use 179
7. Limitations and future research directions Beshir, H. (2014). Factors affecting the adoption and intensity of use of improved forages in North East First and foremost, the most obvious limitation of this Highlands of Ethiopia. American Journal of Experimental study is its cross-sectional design. This limitation calls upon future research to adopt a longitudinal (time /AJEA/2014/5481. Beyene,Agriculture B., Hundie, 4(1), D., 12–27, Gobena, https://doi.org/10.9734 G. (2015): Assessment series) study approach to robustly capture the impact on dairy production system and its constraints in of public infrastructure investment towards determin- Horoguduru Wollega Zone, Western Ethiopia. Science, ing the exact role of rural transport infrastructure in purchased input use. Second, the environmental costs https://doi.org/10.4314/star.v4i2.28. of modern agricultural input use were overlooked by Birachi,Technology E. A., Ochieng, and Arts J., Research Wozemba, Journal D., Ruraduma, 4(2), 215–221, C., this study; therefore, future research concerns should consider a more in-depth analysis of the impact of farmers’ bean production and supply to market in modern agricultural input use on environmental qual- Niyuhire, M. C. (2011): Factors influencing smallholder ity. To measure smallholder farmers’ purchased input Bradbury, A., Hine, J., Njenga, P., Otto, A., Muhia, G., Willilo, S. (2017):Burundi. Evaluation African Crop of the Science effect Journal of road 19(4), condition 335–342. on the of variable inputs was used. Purchased input use in quality of agricultural produce (Report No. RAF2109A). thisuse intensity,study is therefore the total financial refers to cost the (monetaryquantity of value) mon- United Kingdom: Cardno Emerging Markets (UK) Ltd. Carrer, M. J., Silveira, Rodrigo L. F., Filho, H. M. (2019): ey that smallholder farmers spent on major food crop variable inputs (chemical fertilizer, improved seed, from Brazilian citrus growers. Australian Journal of and pesticides) in 2015/16 crop production season. Factors influencing hedging decision: evidence Indeed, in order to estimate the extent of purchased https://doi.org/10.1111/1467-8489.12282. input use we used information on quantity of varia- Agricultural and Resource Economics 63(1), 1–19, ble inputs used and the prices at which they are pur- chased. Since the effect of agricultural input use on Christiaensen,https://doi.org/10.1016/j.foodpol.2017.02.002. L. (2017): Agriculture in Africa – Telling farmers’ productivity is not the concern of this paper, CSA.myths (2011): from Federal facts: ADemocratic synthesis. RepublicFood Policy, of Ethiopia 67, 1–11, it can be another potential area for future study. Central statistical Agency Statistical Abstract, Addis Ababa, Ethiopia. CSA. (2014): Agricultural sample survey: Report on farm management practices for private peasant holdings, 3, References Addis Ababa, Ethiopia. CSA. (2015): Key Findings of the 2014/2015 Agricultural Abdullah, Ahamad, R., Ali, S., Chandio, A. A., Ahmad, Sample Surveys. The Federal Democratic Republic of W., Ilyas, A., Din, I. U. (2017): Determinants of Ethiopia Central Statistical Agency, Addis Ababa, Ethiopia. Commercialization and its impact on the Welfare of CSA. (2017): Ethiopia: Demographic and health survey Smallholder rice Farmers by using Heckman’s two-stage 2016. Central Statistical Agency, Addis Ababa, Ethiopia approach. Journal of the Saudi Society of Agricultural and ICF International, Maryland, USA. Dahal, B. M., Sitaula, B. K., Bajracharya, R. M. (2007): /j.jssas.2017.06.001. Akramov,Sciences K. 18(2), (2009): 224–233, Decentralization, https://doi.org/10.1016 agricultural services and food security in Nepal. Asian Journal of Water, and determinants of input use in Nigeria (IFPRI Sustainable agricultural intensification for livelihood Discussion Paper No. 00941). Washington, D.C., USA: Daoud, J. I. (2017): Multicollinearity and regression International Food Policy Research Institute (IFPRI). Environment and Pollution 5(2), 1–12. Retrieved from http://www.ifpri.org/category https://doi.org/10.1088/1742-6596/949/1/012009. /publication-type/discussion-papers. Delaney,analysis. S., Livingston,Journal of Physics: G., Schonberger, Conference S. (2017). Series 949,Right 1–6, Arethun, T., Bhatta, B. P. (2012): Contribution of rural place, right time: Increasing the effectiveness of roads to access to- and participation in markets: agricultural development support in sub-Saharan Africa Theory and results from Northern Ethiopia. Journal of development support in sub-Saharan Africa. South
//doi.org/10.4236/jtts.2012.22018. https://doi.org/10.1080/10220461.2011.622950. Transportation Technologies 2, 165–174, https: Diao,African X., Silver, Journal J., Takeshima, of International H. (2016): Affairs Agricultural 18(3), 341–365, to transportation infrastructure and economic. Paper mechanization and agricultural transformation Banerjee,presented A., Duflo, at IGC E., Conference, Qian, N. (2012): London, On https: the Road: Access (Background Paper for African Transformation Report //doi.org/10.3386/w17897. 2016: transforming Africa’s agriculture). International Bekele, A., Kassa, B., Legesse, B., Lemma, T. (2010): Effects Food Policy Research Institute. of crop commercial orientation on productivity of Didenko, N., Skripnuk, D., Mirolyubova, O., Radion, M. (2017): smallholder farmers in drought-prone areas of the Analysis of rural areas development of the region using Central Rift Valley of Ethiopia. Ethiopian Journal of the adl-model. Research for Rural Development 2, https://doi.org/10.22616/rrd.23.2017.061. Belay, K. (2015). Agricultural extension in Ethiopia: The case Elias, A., Nohmi, M., Yasunobu, K., Ishida, A. (2013): Effect ofAgricultural participatory Sciences demonstration 20, 16–34. and training extension of agricultural extension program on smallholders’ system. Journal of Social Development in Africa 18(1), farm productivity: Evidence from three peasant associations in the highlands of Ethiopia. African Journal
49–84, https://doi.org/10.4314/jsda.v18i1.23819. 180 Sileshi Tamene, Tebarek Lika Megento
Jayne, T. S. Ã., Govereh, J., Wanzala, M., Demeke, M. (2003): //doi.org/10.5539/jas.v5n8p163. Fertilizer market development: a comparative analysis of FAO.of (2015):Agricultural The economicResearch 6(2),lives of476–487, smallholder https: farmers: An analysis based on household data from nine countries, https://doi.org/10.1016/j.foodpol.2003.08.004. Rome, Italy. Jelilov,Ethiopia, G., Kachallah, Kenya , and M. B.Zambia. (2017): Food The Policy nexus 28, among 293–316, road Freebairn, J. (2003): Economic policy for rural and regional transport and the economic growth in Nigeria. The Australia. Australian Journal of Agricultural and Journal of Middle East and North Africa Sciences 3(9),
//doi.org/10.1111/1467-8489.00220. Kamara, A. B. (2004): The impact of market access on Friday,Resource O. R., EconomicsEmenonye, 47(3),C. (2012): 389–414, The detention https: and input22–29, use https://doi.org/10.12816/0040813. and agricultural productivity: evidence from correction of multicollinearity effects in a multiple regression diagnostics. Elixir International Journal 49, https://doi.org/10.1080/03031853.2004.9523645. Kassa,Machakos B. (2014): district, Assessment Kenya. Agrekon of factors 43(2), affecting 202–218, Fufa, B., Hassan, R. (2006): Determinants of fertilizer use on10108–10112. maize in Eastern Ethiopia: A weighted endogenous farmers of Southern Tigray, Northern Ethiopia (master’s sampling analysis of the extent and intensity of adoption. thesis).agricultural Mekelle production : University, Evidence Mekelle, from Ethiopia. smallholder Katungi, E., Horna, D., Gebeyehu, S., Sperling, L. (2011): /03031853.2006.9523732. Gebremedhin,Agrekon 45(1), B., Hoekstra, 38–49, https://doi.org/10.1080 D. (2008): Market orientation of common bean in Ethiopia: A microeconomic analysis. smallholders in selected grains in Ethiopia: Implications Market access, intensification and productivity of for enhancing commercial transformation of subsistence Kim, N. (2015): Tests based on skewness and kurtosis for African Journal of Agricultural Research 6(2), 476–487. agriculture. IPMS (Improving Productivity and Market multivariate normality. Communications for Statistical Success) of Ethiopian Farmers. Project Working Paper No. 11. ILRI. //doi.org/10.5351/CSAM.2015.22.4.361. Applications and Methods 22(4), 361–375, https: Gray, P. S., Williamson, J. B., Karp, D. A., Dalphin, J. R. (2007): Kiplimo, L. B., Ngeno, V. (2016): Understanding the effect The Research Imagination: An Introduction to Qualitative application of quantile regression-based thick frontier and Quantitative Methods. Cambridge: Cambridge approachof land fragmentation to maize production on farm levelin Kenya. efficiency: Paper Anpresented University Press, https://doi.org/10.1017 at the 5th International Conference of the African /CBO9780511819391. Association of Agricultural Economists, September Hailu, B. K., Abrha, B. K., Weldegiorgis, K. A. (2014): Adoption and impact of agricultural technologies on Kotu, B. H., Alene, A., Manyong, V., Hoeschle-Zeledon, I., farm income: evidence from Southern Tigray, Northern Larbi,23–26, A. 2016, (2017): Addis Adoption Ababa, andEthiopia. impacts of sustainable Ethiopia. International Journal of Food and Agricultural
Hailu, C., Fana, C. (2017): Determinants of Input doi.org/10.1080/14735903.2017.1369619.intensification practices in Ghana. International Journal Economics 2(4), 91–106. Commercialization as Buyers of Agro-chemicals and Lim-Applegate,of Agricultural H., Sustainability Rodriguez, G., 15(5), Olfert, 539–554, R. (2002): https:// improved seed: Evidence from farm households’ of Determinants of non-farm labour participation rates Ambo and Toke Kutaye districts, West Shewa Zone, among farmers in Australia. Australian Journal of Ethiopia. American Research Journal of Agriculture 3(1), https://doi.org/10.1111/1467-8489.00168. .17004. Lulit,Agricultural A. (2012): and Impact Resource of Road Economics on Rural 46(1),Poverty 85–98, Evidence 1–14, https://doi.org/10.21694/2378-9018 Hayami, Y., Ruttan, V. W. (1985): Agricultural development. Form Fifteen Rural Villages in Ethiopia (master’s thesis). Johns Hopkins University Press, Baltimore. Erasmus University Rotterdam, Institute of Social Hayami, Y., Ruttan, V. W. (1993): Induced technical and Studies (ISS), The Hague, the Netherlands. institutional change evaluation and reassessment: two Ma, W., Bicknell, K., Renwick, A. (2019): Feed use chapters. Bulletin Number 93-1, University of Minnesota Economic Development Center, Minnesota, USA. in New Zealand. Australian Journal of Agricultural and Hellin, J., Lundy, M., Meijer, M. (2009): Farmer organization, intensification and technical efficiency of dairy farms collective action and market access in Meso-America. doi.org/10.1111/1467-8489.12283. Margarian,Resource A. Economics (2011): Endogenous 63(1), 20–38, Rural https:// Development: /j.foodpol.2008.10.003. Empowerment or Abandonment? Paper presented at Hine,Food J. (2004): Policy, 34(1),Good Policies 16–22, https://doi.org/10.1016and Practices on Rural the 4th International Summer Conference in Regional Transport in Africa Planning Infrastructure and Services Science, Dresden. (SSATP Working Paper No. 100). Washington, D.C., USA: Markelova, H., Meinzen-dick, R., Hellin, J., Dohrn, S. (2009): The World Bank. Collective action for smallholder market access. Food Hung, P. Van, MacAulay, T. G., Marsh, S. P. (2007): The economics of land fragmentation in the North of /j.foodpol.2008.10.001. Vietnam. Australian Journal of Agricultural and Resource Minten,Policy B., 34(1), Koru, 1–7, B., Stifel, https://doi.org/10.1016 D. (2013): The last mile(s) in
/j.1467-8489.2007.00378.x. Imdadullah,Economics M., 51, Aslam, 195–211, M., Altaf, https://doi.org/10.1111 S. (2016): mctest: doi.org/10.1111/agec.12078.modern input distribution: Pricing, profitability, and An R Package for Detection of Collinearity among Miot,adoption. H. A. (2017): Agr Econ Assessing 44(6), normality629–646, https://of data in clinical
doi.org/10.32614/RJ-2016-062. doi.org/10.1590/1677-5449.041117. Regressors. The R Journal 8(2), 495–505, https:// and experimental trials. Vasc Bras 16(2), 88–91, https:// Effect of rural transport on smallholder farmers’ purchased input use 181
Nambiro, E. (2008): Trends in land use and agricultural Tamene, S., Megento, T. L. (2017): The effect of rural road transport infrastructure on smallholder farmers’ dissertation, Rheinische Friedrich-Wilhelms University, agricultural productivity in Horro Guduru Wollega Nairobi,intensification Kenia). in Retrieved Kakamega, from Western http://hss.ulb.uni-bonn Kenya (Doctoral .de/diss. https://doi.org/10.14712/23361980.2017.7. Negari, K. I. (2017): Compiled Body of Works in Field Vandercasteelen,zone, Western J.,Ethiopia. Tamru, AUCS., Minten, Geographica B., Swinnen, 52(1), J. 79–89, Epidemiology. Ethiopia Field Epidemiology Training (2016): Cities and Agricultural Transformation in Africa: Program(EFETP) Addis Ababa University. Evidence from Ethiopia (Working Paper No. 374/2016). Belgium: LICOS Centre for Institutions and Economic Fertilizer Use. In: Samuel Benin (Ed.), Agricultural Performance, https://doi.org/10.2139/ssrn.2744504. Nin-Pratt, A. (2016): Agricultural Intensification and Weir, S., Knight, J. (2004): Externality effects of education: Dynamics of the adoption and diffusion of an innovation Internationalproductivity in Food Africa : Policy trends, Research patterns, Institute. and in rural Ethiopia. University of Chicago, https:// Ninnin,determinants B. (1997): (pp. Transport 199–246). et Developpement Washington, DC: doi.org/10.1086/423254. A Madagascar. French Co-operation Ministry and Wondemu, K. A., Weiss, J. (2012): Rural Roads and Malagasy Public Works Ministry, INRETS. Development: Evidence from Ethiopia. EJTIR 12(4), Olana, B. T. (2006): People and Dam: Environmental and Socio-economic changes induced by reservoir in Fincha Worku, I. (2011): Road sector development and economic Water shades, Western Ethiopia. PhD Dissertation, growth417–439. in Ethiopia. Ethiopia Support Strategy Program Wageningen University, Netherland. II, International Food Policy Research Institute, Addis Orakwue, C., Umeghalu, I., Ngini, J. (2015): Effects of Road Transport on Agricultural Productivity: A Case Study of Yap, B. W., Sim, C. H. (2011): Comparisons of various types Ayamelum Local Government Area of Anambra State, ofAbaba, normality Ethiopia, tests. 101–146. Journal of Statistical Computation and
Osborne, J. W., Waters, E. (2002): Four assumptions of /00949655.2010.520163. multipleNigeria. Inter regression J Appl thatSci Engr researchers 3(1), 1–4. should always test. Zewdie,Simulation T. D. 81(12),(2015): 2141–2155, Access to https://doi.org/10.1080Credit and the Impact of Credit constraints on Agricultural Productivity Perz, S. G. (2003): Social determinants and land use in Ethiopia: Evidence from Selected Zones of correlatesPractical Assessment, of agricultural Research, technology and Evaluation adoption in 8(2), a forest 1–5. Rural Amhara (master’s thesis). Addis Ababa frontier: A case study in the Brazilian Amazon. Human University, Ethiopia.Salami, A., Kamara, A. B., Brixiova, Z. (2010): Smallholder Agriculture in East Africa: Trends, /A:1022838325166. Constraints and Opportunities. Working Papers Series Pingali,Ecology P. L., 31(1), Rosegrant, 133–165, M. W. https://doi.org/10.1023 (1995): Agricultural No. 105, African Development Bank, Tunis, Tunisia. Tamene, S., Megento, T. L. (2017): The effect of rural road transport infrastructure on smallholder farmers’ //doi.org/10.1016/0306-9192(95)00012-4.commercialization and diversification: processes and agricultural productivity in Horro Guduru Wollega Quan,policies. T. T. (2009): Food Policy Transition 20(3), from 171–185, subsistence https: farming to commercial agriculture in Quang Binh province, Vietnam https://doi.org/10.14712/23361980.2017.7. (doctoral dissertation). Lincoln University, Oxford, USA. Vandercasteelen,zone, Western J.,Ethiopia. Tamru, AUCS., Minten, Geographica B., Swinnen, 52(1), J. 79–89, Roland-holst, D. (2009): Infrastructure as a Catalyst (2016): Cities and Agricultural Transformation in Africa: for Regional Integration, Growth, and Economic Evidence from Ethiopia (Working Paper No. 374/2016). Convergence: Scenario Analysis for Asia. In: F. Zhai (Ed.), Belgium: LICOS Centre for Institutions and Economic From Growth to Convergence: Asia’s Next Two Decades Performance, https://doi.org/10.2139/ssrn.2744504. Weir, S., Knight, J. (2004): Externality effects of education: doi.org/10.1057/9780230250604_4. Dynamics of the adoption and diffusion of an innovation Ruttan,(108−149). V. W. (2008): Palgrave Induced Macmillan technical UK, https://change induced in rural Ethiopia. University of Chicago, https://doi institutional change and mechanism design. Paper .org/10.1086/423254. presented at the 10th International Workshop on Wondemu, K. A., Weiss, J. (2012): Rural Roads and Institutional Economics, Institutions, Technology and Development: Evidence from Ethiopia. EJTIR 12(4), their Roles in Economic Growth, University of Hart. Ruttan, V. W. (1996): Induced Innovation and Path Worku, I. (2011): Road sector development and economic Dependence: A Reassessment with Respect to growth417–439. in Ethiopia. Ethiopia Support Strategy Program Agricultural Development and the Environment. II, International Food Policy Research Institute, Addis
https://doi.org/10.1016/0040-1625(96)00055-8. Yap, B. W., Sim, C. H. (2011): Comparisons of various types Ruttan,Technological V. W., Hayami, Forecasting Y. (1984): and Toward Social Changea theory 53, of 41–59,induced ofAbaba, normality Ethiopia, tests. 101–146. Journal of Statistical Computation and institutional innovation. The Journal of Development /00949655.2010.520163. doi.org/10.1080/00220388408421914. Zewdie,Simulation T. D. (2015): 81(12), Access 2141–2155, to Credit https://doi.org/10.1080 and the Impact Salami,Studies A., Kamara,20(4), 203–223, A. B., Brixiova, https:// Z. (2010): Smallholder of Credit constraints on Agricultural Productivity in Agriculture in East Africa: Trends, Constraints and Ethiopia: Evidence from Selected Zones of Rural Amhara Opportunities. Working Papers Series No. 105, African (master’s thesis). Addis Ababa University, Development Bank, Tunis, Tunisia. Ethiopia. 182 Original Article
The fertility revolution in Zimbabwe with special regards to proximate determinants of fertility Collet Muza*
Charles University, Faculty of Science, Department of Demography and Geodemography, Prague, Czechia * Corresponding author: [email protected]
ABSTRACT The role of proximate determinants in moderating fertility decline is well documented in developing countries. In Zimbabwe, howev- er, there is a limited understanding of the role of proximate determinants on fertility levels and trends. This study aimed to examine the role of proximate determinants of fertility (namely marriages, postpartum infecundity and contraception) using the Bongaarts proximate determinants model. The impact of these determinants is studied on a sample of married women aged 15–49 years’ and corresponding cross-sectional data obtained through the six consecutive Zimbabwe Demographic Health Surveys (ZDHSs) hold in 1988, 1994, 1999, 2005 and 2015. The results reveal that the overall fertility declined from 5.4, 3.8 and 4.0 children per woman observed among 1988, 1999 and 2015 ZDHSs, respectively. This change was caused by the contraceptive inhibitive effect, which correspondingly increased from 3.00 to 4.65 and 6.45 children per woman. The fertility stalling observed in 1999 and after that is caused by postpartum infecundity and marital fertility inhibition which decreased with time. Moreover, contraceptive inhibition effect increased with education, wealth quintiles, and urban residence. In contrast, marital and postpartum infecundity fertility inhibition effects inversely correlate with education, wealth quintiles, and the place of residence. Therefore, to foster further fertil- ity decline to replacement level, policies should promote contraceptive adoption, more extended breastfeeding periods and delay entry into early marriages. Furthermore, women empowerment, especially the promotion of female education to higher education and female employment, could be useful tools to further fertility decline.
KEYWORDS proximate determinants; fertility stalling; decomposition; Bongaarts model
Received: 18 April 2019 Accepted: 25 September 2019 Published online: 9 December 2019
Muza, C. (2019): The fertility revolution in Zimbabwe with special regards to proximate determinants of fertility. AUC Geographica 54(2), 182–193 https://doi.org/10.14712/23361980.2019.16
© 2019 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). The fertility revolution in Zimbabwe 183
1. Introduction - sis-driven fertility decline is when a combination of - high(Mhloyi costs 1988; of living, Müller declining et al. incomes, 2013). In civil contrast, unrest and cri enced either sustained high fertility or fertility stalling persistence droughts in Zimbabwe has forced couples, atA significant midway transition, number of which African has countries put future have prospects experi regardless of their levels of modernization, to adopt of fertility decline into doubt and made demographers their fertility downwards (Mhloyi 1988). For instance, question the demographic transition theory. Accord- the costs related to the education of the children rath- ing to Bongaarts (2008), fertility stalling is when er than the education of mothers themselves have led women to decrease their demand for children. In other decline of 10% or more between two successive sur- aveys country demographic does not and experience health surveys. significant Consequently, fertility from children to parents (Caldwell and Caldwell 2002). because of sustained high fertility and low mortality, words,However, the high Gould costs (2015) of living has reverse argued thatthe wealth crisis-driv flow- the population growth rate in Africa has remained en fertility decline cannot be sustained in the long high at about 2.5% per annum since 1960 (United run. Precisely fertility in Zimbabwe has stalled at Nations 2017). The population doubling time asso- midway transition at about 4 children per woman since 2000 (ZIMSTAT and ICF International 2016). - The question is whether fertility stalled as a result of lionciated from with 1960 such to an 2015 annual (ibid.). growth The challengerate is 28 years. is that changes in background variables or changes in proxi- highThus, population the population growth of Africa has been grew shown from 0.3to negatively to 1.2 bil mate determinants? affect socio-economic and human development (Bird- Although the effects of the proximate determi- sall and Sinding 2001). This is alarming since Africa nants on fertility have been documented in Zimbabwe is already the world’s least developed region (United (Guilkey and Jayne 1997; Letamo and Letamo 2001; Nations 2017). Mhloyi 1986; Muhwava and Muvandi 1994; Siban- In 1956, it was propounded that socioeconomic da 1999), little has been done recently to show the background factors operate through 11 proximate current scenario with regard to proximate determi- nants. Undoubtedly, the socioeconomic circumstances society (Davis and Blake 1956). Further, Bongaarts (1978)determinants revised to andinfluence reduced the Davies level of and fertility Blake’s in any 11 which might also change the proximate determinants. into 7 proximate determinants of fertility. Further- Thus,appear this to studyhave changed seeks to sinceassess significantly the proximate changed, deter- more, Bongaarts (1982) demonstrated, that 96% of minants of fertility in Zimbabwe, incorporating recent fertility levels and variation in developing countries data. The data were pooled from 6 consecutive (1988, is a result of four main proximate determinants; 1994, 1999, 2005 2010, and 2015) Zimbabwe Demo- marriage, postpartum infecundity, abortion and con- graphic and Health Surveys (ZDHSs). This study will traception. The remaining seven variables explain contribute to our knowledge and understanding of a small part of the total fertility rate variation (ibid.). proximate determinants of fertility in Zimbabwe and The challenge is that Africa lacks historical and cur- how they have associated fertility levels and trends. rent complete vital registration systems for conduct- A substantial body of literature has generally ing such analysis (Dyson 2013). However, the seminal agreed that Africa’s early and universal nuptiality con- work of Notestein (1945) has established that the ditions are more conducive to high fertility rates than fertility revolution from natural to controlled fertil- witnessed in Europe’s historical conditions (Chesnais ity regimes is a result of the deliberate adoption of 1992; Coale 1973; Dyson 2013). Undeniably, entry birth control methods (contraception). Unlike Europe into marriage serves as a risk factor to childbearing where women of older reproductive ages used stop- as most of the childbearing was happening in mar- ping reproductive behaviour and natural methods of riages (Bongaarts 1982; Coale 1973). In Zimbabwe, - a recent study by Sayi and Sibanda (2018) based on ods of contraception were used for both child spacing 2014 Zimbabwe Multiple Indicators Cluster Survey andfertility stopping limitation, in younger in Africa, and olderefficient reproductive modern meth ages respectively (Caldwell and Caldwell 2002). Thus, in Africa, the fertility revolution, although still incom- about(MICS) 32% indicates were thatmarried 1 in before4 women 18 agedyears 15–19 of age. years Fur- plete, once started, is more rapid than Europe’s his- ther,were thecurrently gap between married, marriage while among and childbearing 20–49 years is torical fertility transitions. very small as couples seek to strengthen their mar- In Zimbabwe, fertility had declined from TFR 5.4 in riage with children (Chitereka and Nduna 2010). Mar- 1984 to 3.8 in 1999 has been both development and crisis-driven fertility decline (Mhloyi 1988; Müller to fertility levels as remarriages occur early (Mhloyi 1988).ital dissolutions However, ithave is possible been found that nuptialityto be insignificant patterns typical of the Becker’s theory of fertility decline when changed with time and consequently fertility levels. foret al. instance, 2013). Development-driven highly urbanized, educated fertility declineand high is Polygyny is generally high and variable in Africa. incomes women were further ahead in the fertility McDonald (1985) reported that polygyny ranged from transition than disadvantaged socioeconomic groups 10% to 67% in Lesotho and Senegal, respectively. In 184 Collet Muza
Zimbabwe in a recent survey, polygamous unions a later study using 2005/06 DHS data on Zimbabwe have been found to be 8% and 16% in urban and rural Mturi and Kembo (2011) reveal that high contracep- areas, respectively (ZIMSTAT and ICF International tion prevalence rate was the most important factor of 2016). Studies on the effects of polygamy on fertility fertility decline even during periods of socioeconomic offer contradictory evidence. On the one hand, since development and socioeconomic crisis. The question women in the polygamous union share time of hus- is why fertility in Zimbabwe has stalled at midway band, they have less exposure to sex than women in transition given such high contraceptive prevalence monogamous relationships and correspondingly have rate? South Africa (TFR 2.3) and Botswana (TFR 2.6) less risky to pregnancy. Moreover, women in polyg- with similar contraceptive prevalence levels and they amous unions have been noted to breastfeed longer have lower fertility approaching replacement lev- than women in monogamous unions (Sayi and Siban- el fertility (United Nations 2017). This is important da 2018; Mhloyi 1988). On the other hand, Mhloyi given that studies have revealed contraception prev- (ibid.) argued that women in polygamous unions’ alence of 75% is associated with replacement level demand for children is higher than women in monog- fertility (Mahjabeen and Khan 2011). amous unions, as more children offer them security, It is also possible that contraceptives such as con- respect and access to husband’s wealth. The question doms might be used for HIV prevention but not par- is, has the polygyny patterns changed with inevitable modernization? Furthermore, although variable and decreas- contributoryity-specific fertility effects limitationof the proximate behaviour determinants (Terceira ing with modernization, Africa’s long postpartum varyet al. positively2003). Nonetheless, with women studies empowerment, have revealed i.e. wom the- abstinences were strategies ensuring the survival en’s education, wealth quintiles, urban and rural res- of already born children through prolonged birth idence. In Zambia, secondary education and urban spacing rather than deliberate birth control methods residence and wealth were established to be positive- (ibid.). Caldwell and Caldwell (2002) have argued ly related to higher relative fertility inhibition effect that child spacing is embedded in the African cul- in marriage and contraception whilst postpartum ture, and high levels of contraceptive use might be for infecundability was inversely correlated (Chola and spacing and not stopping childbearing. Nonetheless, it can be hypothesised that postpartum infecundity has decreased with modernization. In the absence of KemboMichelo 2011; 2016). Sibanda Similar 1999). findings were made in Uganda alternative contraceptive adoption or increments, fer- (Rutaremwa et al. 2015) and in Zimbabwe (Mturi and tility can increase. Several studies have employed the Bongaarts 2. Data and methods framework for the analysis of the fertility differen- tials, levels and transitions in Africa (Chola and Mich- 2.1 Sources of data This paper utilizes pooled data from 6 consecutive 2015; Mturi and Kembo 2011; Sibanda 1999; Bon- ZDHSs conducted in 1988, 1994, 1999, 2005, 2010 gaartselo 2016; 1978). Majumder A recent and analysis Ram 2015; of countries Rutaremwa in the et al. ear- and 2015. In each survey year, a nationally represent- ly stages of fertility transition namely Namibia (Pal- amuleni 2017); Zambia (Chola and Michelo 2016), which were conducted under the collaboration of Zimbabweative survey National of ever-married Statistics women Agency of (ZIMSTAT)age 15–49 1996), and found the dual importance of marriage and ICF International. The ZDHSs provides nation- andUganda postpartum (Rutaremwa infecundability et al. 2015), Malawias most (Palamuleni important ally representative data on basic health social and predictors of fertility outcomes. The contribution of demographic indicators. The study utilized this data contraception was least although it was increasing rapidly. assessing the contribution of marriage, contraception - andin order postpartum to fit the infecundability aggregate fertility in Zimbabwe.model, thereby The cant fertility decline, it has been established that the fertility estimates were also disaggregated by a num- contributionIn African countriesof contraception that have is theexperienced most important signifi ber of selected socioeconomic variables, namely resi- and amenable to fertility reduction (Finlay, Mejia- dence, education and wealth quintiles. Guevara and Akachi 2018; Majunder and Ram 2015; Mturi and Kembo 2011; Sibanda 1999). Sibanda 2.2 Sampling (1999) looking at the relative contribution of prox- The sample sizes of the interviewed women age imate determinants in Zimbabwe and Kenya using two consecutive surveys for each country show that in Zimbabwe contraception is the most important Sampling(15–49) were was selecteddone using based a two-stage on a master cluster sampling sam- factor in fertility decline in younger and middle ages. plingplan which process. was Initially, provided clusters by ZIMSTAT are selected(1988–2015). from In Kenya, with slightly higher fertility than Zimba- a list of clusters obtained from the master sampling bwe, it was found that postpartum infecundity was plan provided by ZIMSTAT, followed by a section of the most important factor followed by marriage. In households from each cluster. The data obtained were The fertility revolution in Zimbabwe 185
2.4 Index of marriage (Cm) considered adequate to enable analysis and com- The index of marriage (Cm) measures the inhibit- stratified by rural and urban areas. The samples are ing effects of marriage on fertility in the population. of socioeconomic and demographic locus that could The lower the proportions of married, the higher the guideparisons fertility that andwould population be useful policy in the interventions identification in inhibiting effects of marriage and the inverse is true. Zimbabwe. since inhibiting effects of marriage are marital age 2.3 Analytical framework distributionHowever, age-specific sensitive. Childbearingmarital proportions is greatest are in used the Bongaarts’ proximate determinants model (1978, 1982) was applied for analyzing proximate deter- - minants of fertility from the six successive ZDHSs tioncentral is the age assumption distribution that years. only Marriage women inis definedmarriages as named above. Bongaarts and Potter (1983) devel- areformal exposed or consensus to the risk marriages. of childbearing. Implicit in The this ZDHSsdefini oped a technique to quantify the impact of four proxi- mate determinants on fertility, namely marriage, con- cohabitation as a marital union in Zimbabwe. This traception, abortion, and postpartum infecundability. wasdone not in the1999, case 2005, in 1988 2010 and and 1994 2015 since has no classified informa- They assume that the total fecundity rate (TF) of all tion was collected on this variable. In order to make women is the same, but their real reproductive per- proximate determinants. The mechanism of the mod- estimatedthe ZDHSs asmarital follows: definition comparable, we take the elformance is summarized is modified by relating by the the above fertility four measures mentioned to earlier definition used in 1988 and 1994. The index is ( ) ( ) the proximate determinants. The basic model equa- = (Bongaarts 1982), ( ) tions are: ∑ ∑ (a) TFR = Cm × Ci × Cc × Ca × TF (Bongaarts 1982), where m(a (b) TM = Ci × Cc × Ca × TF (Chola and Michelo 2016), married women. This characteristic is computed by (c) TN = Ci × TF (ibid.), dividing the) number is age-specific of married proportions women by currently the total number of women in the same age group. Also, g(a) where TFR is the number of births, a woman would have at the end of her reproductive years if she were by dividing the births of a particular age group (from marriedis age-specific women) marital by the fertility number rates. of married It is computedwomen in rates throughout the reproductive period. Total mari- the same age group. talto bear fertility children rate ( TMat the) is prevailingthe number age-specific of births a fertilitywoman would have at the end of her reproductive year if she 2.5 Index of contraception C( c) The index of contraception measure the inhibiting marital fertility rates and remain married during the effects of modern methods of contraception on fertil- entirewere to reproductive bear children period at the (Bongaarts prevailing 1987).age-specific Total ity in a population which also varies with the preva- natural fertility rate (TN) is observed under condi- lence and use-effectiveness of contraception used by tions in which contraception and abortions are elimi- couples in the reproductive age groups. The higher nated (Bongaarts 1982). Bongaarts noted that whilst the level of contraception in a population, the higher TFR, TM, and TN vary in many populations, total the inhibiting effect of contraception and vice versa. fecundity is constant in all populations. Total natural The index of contraception is estimated using the fecundity rate (TF) index is estimated as follows: following:
TF = TFR / (Cm × Ci × Cc × Ca) (Bongaarts 1982), Cc u × e (Bongaarts and Potter 1983), where, Cm, Ci, Cc, and Ca are indices of marriage, where = 1 – u 1.08 is the × average proportion of married wom- postpartum infecundability, and contraception and en currently using contraception; is average contra- induced abortion, respectively. The indices can only - take values between 0 and 1. Where there is no fer- rility correction factor (represents the adjustment tility inhibiting effect of a given intermediate fertili- forceptive the facteffectiveness. that women The do coefficient not use contraception 1.08 is the ste if ty variable, the corresponding index equals one, and they know they are sterile). The parameter e values when the fertility inhibition is complete the index proposed for particular contraceptives are as follows: equal 0 (Bongaarts 1982). Abortion is illegal in Zim- pill 0.90, IUD 0.95, injection 0.99, sterilization 1.00, babwe (ZIMSTAT and ICF International 2016; Siban- others 0.70 (Bongaarts and Potter 1983). da 1999), and there is limited and unreliable infor- mation. Therefore, the index of abortion in this study 2.6 Index of postpartum infecundability (Ci) will be 1. The contribution of abortion is regarded as The index of postpartum infecundability is a meas- ure of the inhibiting effect of breastfeeding or sexu- the measure of proximate variables. al abstinence on fertility in a population (Bongaarts insignificant. The indices can be approximated from 186 Collet Muza
1982). The index of postpartum infecundability in the 2016). Moreover, abortion data in the ZDHSs include model is estimated using the effect of breastfeeding (lactation amenorrhea) or postpartum abstinence. to isolate abortion data. The index of abortion is esti- The index was calculated as follows: matedstillbirths using and the miscarriages. following formula: Therefore, it is difficult
C1 = 20/18.5 + i (Bongaarts 1982), Ca = TFR /TFR + b × TA = TFR/TFR + 0.4 × (1 + U) × TA (Bongaarts and Potter 1983) where the symbol i represents the average duration where U is contraceptive prevalence use, b is the postpartum ovulation (menses), in this research, the average number of births averted per induced abor- indexin months of postpartum of infecundability infecundability from birth was to estimated the first tion and b = 0.4(1 + U); b is 0.4 when U = 0 and b = 0.8 using the mean duration of breastfeeding. when U = 1.0. TA is total abortion (an average num- ber of induced abortion per woman at the end of the 2.7 Index of abortion C( a) reproductive period if induced abortions remain at The index of abortion measures the inhibiting effect of prevailing levels throughout the reproductive period). abortion on fertility in a population. In this research, Then Ca = 1.0 if TA = 0. Therefore, the index of abor- the index of abortion was set 1.0 due to lack of data. tion in this study is equal to . Abortion is illegal in Zimbabwe excerpt for health Based on the studies of historical populations with and legal reasons (ZIMSTAT and ICF International the highest recorded fertility, Bongaarts recommends
Tab. 1 Percentage distribution of respondent women, 15–49 years by selected characteristics (weighted), 1988–2015, Zimbabwe.
1988 1994 1999 2005 2010 2015 Children born 1–3 34.9% 38.1% 44.7% 47.6% 52.3% 50.2% 4–5 16.7% 15.0% 12.9% 14.0% 14.6% 17.0% 6+ 19.9% 17.5% 12.6% 9.0% 6.7% 6.1% Never Marital status 27.0% 26.9% 27.7% 27.0% 24.0% 25.2% married Married 62.9% 61.8% 56.3% 56.3% 59.4% 58.7% Divorced 7.6% 7.8% 3.5% 4.5% 3.7% 5.0% Widowed 2.5% 3.5% 4.2% 7.5% 6.1% 4.4% Age group 15–19 24.3% 24.0% 24.5% 24.2% 21.2% 22.1% 20–24 20.0% 20.7% 21.9% 21.9% 20.1% 17.0% 25–29 16.2% 14.9% 17.5% 16.5% 18.4% 16.6% 30–34 14.0% 14.2% 11.3% 13.7% 14.1% 16.3% 35–39 11.0% 10.8% 10.8% 9.4% 11.5% 12.4% 40–44 7.6% 8.7% 7.9% 7.8% 8.0% 9.7% 45–49 6.9% 6.6% 6.1% 6.6% 6.8% 5.8% Residence Urban 33.5% 32.2% 38.6% 39.3% 38.7% 38.5% Rural 66.5% 67.8% 61.4% 60.7% 61.3% 61.5% Education No education 13.5% 11.1% 6.7% 4.3% 2.3% 1.35 Primary 55.9% 47.3% 40.2% 32.6% 28.0% 25.8% Secondary 29.7% 40.0% 50.2% 60.15 65.1% 65.6% Higher 0.9% 1.6% 2.8% 3.0% 4.6% 7.3% Wealth quintile Poorest – 18.2% 17.1% 17.4% 16.9% 17.1% Poorer – 16.1% 17.3% 16.8% 17.4% 17.0% Middle – 18.9% 18.4% 17.4% 18.3% 17.6% Richer – 21.8% 22.3% 22.5% 22.6% 23.2% Richest 25.0% 24.9% 25.9% 24.8% 25.1% Modern (any) 36.1% 42.2% 50.4% 58.4% 57.3% 65.8% Family Planning (married) Traditional 7.0% 4.3% 2.8% 1.4% 1.1% 1.0% Mean duration of breastfeeding Months 18.9 18.8 19.0 18.7 17.5 18.1 Total (n) 4,201 6,128 5,907 8,907 9,171 9,955 – missing values Source: ZDHSs The fertility revolution in Zimbabwe 187 using 15.3 as the maximum number of births per 1% to 7% from 1988 to 2015 respectively. The dis- woman. This is referred to as the Total fecundity rate tribution of the respondents by age groups clearly (TF) (Bongaarts 1978; 1982). The value is the theo- demonstrate that percentage distribution diminished retical number that a woman would have if she were with an increase in age for the respective periods. The 1988 ZDHS did not collect data on wealth quintiles. use contraceptives and did not abort any pregnan- Nonetheless, for the later ZDHSs periods with avail- cies.to continuously Multiplying all married the indices from together age 15–44, by the did total not able data, the wealth quintiles stayed relatively the fecundity rate of 15.3 produces the predicted TFR for same with respondents in the lowest wealth quintile the population. The predicted TFR will normally dif- slightly decreasing from 18% to 17% and a simulta- fer from the observed TFR because of underreport- neous increase in the fourth wealth quintile from 22 ing of births; misreporting of behaviours measured to 23% in 1988 and 2015 respectively. by the indices or omission of proximate factors that Table 1 also shows the distribution of the study by help determine fertility levels in the population under selected key proximate determinants of fertility. The study. The Bongaart’s model was also applied for cal- percentage currently married and single declined culation of the proximate determinants for selected from 63%, 56%, 59% and 27%, 24%, 25% in 1988, background characteristics as the place of residence, 1999 and 2015 respectively. The currently married education and wealth quintiles. percentages were considered as proportion currently married in computing the indices for marriage. Mean duration of breastfeeding from ZDHSs reports was, 3. Results on average 18.5 months from 1988 to 2015, respec- tively, and this is what was used in the calculations The individual characteristics of the respondents for to represent the average duration of the postpartum six successive ZDHSs conducted between 1988-2015 infecundability. With regard to contraception, there are shown in table 1. The results show that the per- has been an increase in modern contraception use by centages living in rural areas decreased from 67% in married women from 36% in 1988 to 66% by 2015 1988 to 61% by 2015. The results reveal that wom- and a simultaneous decrease in traditional methods en with no education and primary education were from 7% to 1% from 1988 to 2015 respectively. The reduced from 14% to 1% and 56% to 26% while above-mentioned information was used in the cal- women with secondary education simultaneously culation of fecundity. As highlighted in the equation increased from 30% to 66% between 1988 and 2015 above, the observed fertility levels in a population are respectively. The percentage of women with higher a product of the relationship of proximate determi- education is still very low, although it increased from nants and the maximum biological level of fertility.
20 T R marriage postpartum infecunda ilit 18
1
n 2,9
a ,15 m 14 ,91 o , , 7 ,02
12