Introduction

Geoheritage research is an increasingly important aspect of Earth Science education (Barasoain and Azanza 2018; Justice, 2018; Reynard et al. 2018; Zecha and Regelous 2018; AlRayyan et al. 2019; Lee et al. 2019; Migon and Pijet-Migon 2019; Rapprich et al. 2019), nature conservation (Brilha, 2002; Henriques et al. 2011; Németh and Moufti 2017; Escorihuela 2018; Gordon, 2018; Melendez 2018; Prosser 2018; Brocx and Semeniuk 2019) and geotourism (Hose 1997; Hose 2006; Newsome and Dowling 2006; Hose et al. 2011; Hose 2012a; Hose 2012b; Lim 2014; Németh et al. 2017; Dowling and Newsome 2018; Erfurt 2018; Newsome and Dowling 2018; Rozycka and Migon 2018; Suzuki and Takagi 2018; Vasiljevic et al. 2018; Beraaouz et al. 2019; Guo and Chung 2019; Maghsoudi et al. 2019). Identifying the geoheritage values of the natural environment is an important aspect of this research (Reynard and Brilha 2018). Methods used to define the value of geological sites include systematic quantification of their main (e.g. scientific, educational, aesthetic, protection etc) and additional (e.g. functional touristic etc) values (Brilha 2016; Sa dos Santos et al. 2016; Brilha 2018; Brilha et al. 2018; Forte et al. 2018). While there is no general agreement on the best practice for doing this, such valorisation defines the fundamental roles in land-use planning, development of various geotourism activities or geoeducation opportunities. The ultimate goal of the work is to develop thematic geoconservation sites with recognised local, regional or international (global) geoheritage values that will promote Earth Science and the fundamental role of the geosystem in understanding Planet Earth (Reynard et al. 2011b; Brocx and Semeniuk, 2015; Chakraborty et al. 2015; de Vries et al. 2018; Gray, 2018a). The formal geoconservation aspect of such projects is  to establish local geoparks and promote geotourism (Joyce 2010b; Zeng 2014; Errami et al. 2015; Warowna et al. 2016; Németh and Moufti 2017; Dowling 2018; Gray 2018b; Beraaouz et al. 2019; Mikhailenko and Ruban 2019). Depending on the significance of the geosites identified during the research some features may be sufficiently important as to warrant the establishment of an UNESCO Global Geopark. UNESCO Global Geoparks are protected and managed regions with geological and geomorphological features of global significance (Eder 1999; Keever and Zouros 2005; Dowling 2011; Petrovic et al. 2013; Turner 2013; Henriques and Brilha 2017; Justice 2018; UNESCO,2018). Achieving UNESCO Global Geopark status requires scientific expertise and the involvement of local communities, including indigenous peoples and cultures. Visitor experiences are the key to provide real, hands on information about Earth history and the processes that shape the non-living natural environment to promote geoheritage (Clary 2018; Dowling 2018; Farsani et al. 2018; Lagally 2018). The geosites need to be destinations keenly sought by visitors either as individual tourist activities or more structured ways, such as guided visits (Migoń, 2018; Newsome and Dowling 2018). Their scientific significance must be clearly articulated (Brilha 2016; Brilha 2018). To ensure this the geosites need to be presented in an appealing and informative manner which can take many years of work (Hose 2006; Migoń 2018). Defining the geoheritage values is the key and is best achieved through ongoing research and dissemination of the findings to the wider community by publication in local media and discussion or workshop opportunities to share knowledge with the local population. Geoheritage research will often define the scientific values of geosites and these findings can be distributed through academic journals (Brilha 2016; Brilha, 2018). The geodiversity of sites is also an important factor that warrants inclusion in the research (Gray 2018a; b; Reynard et al., 2018).

 

Geoheritage values tend to be biased toward the tourism values of geosites in particular the visual landscape (Panizza 2001; Reynard et al. 2007; De Waele and Melis 2009; Ielenicz 2009; Ilies and Josan 2009; Reynard et al. 2009; Erhartic 2010; Joyce 2010b; Pereira and Pereira 2010; Coratza et al. 2011; Doniz-Paez et al. 2011; Miccadei et al. 2011; Pellitero et al. 2011; Reynard et al. 2011a; Serrano et al. 2011; Harmon and Viles, 2013; Reynard and Coratza, 2013; Zglobicki and Baran-Zglobicka, 2013; Reynard et al. 2018). However, research is concluding that geoheritage itself has absolute values that are not based solely on morphology of the landscape and its visual aspects  so often (Brilha 2002; Brilha 2016; Brilha 2018). Tourism driven geoheritage is too often based on the preconception of what visitors wish to experience at a geological site. However, for many sites the geological value is not limited to the visual scene (Dowling and Newsome 2018). When the focus is primarily on geotourism, especially where there is government support for tourism (Hayes and Lovelock 2017), tourist hot spots develop and tourism routes tend to build around a network of popular sites while less well known sites, some of which are off the beaten track, are rarely promoted and not considered as valuable  (Pearce 2000; Morrow and Mowatt 2015; Migon and Pijet-Migon 2016). Recently there appears to have been shift away from focusing only on popular geosites  to studying less well known regions where new  sustainable tourism developments could be established (O’Connor, 2008; Gerner et al. 2009; Joyce 2010a; Bitschene and Schueller, 2011; Boley et al. 2011; Burlando et al. 2011; Cutler 2011; Dowling 2011; Farsani et al. 2011; Gordon, 2011; Grant 2011; Letos and Muntele 2011). The recent COVID-19 pandemic and the resultant absence of international tourists to New Zealand has demonstrated that “hot spot” driven tourism with its narrow view of the geoheritage values is unsustainable long term. Other crisis situations such as various natural disasters, for example the eruption of Whakaari White Island on 9 December 2019, demonstrated the impact such an event can have and how  recovery, if any, can take a very long time (Orchiston and Higham 2016; Tucker et al. 2017). In addition, global changes, such as COVID-19 and the accompanying cessation of international visitors, can alter the tourism portfolio dramatically. A profit centered approach to attracting tourist visitors to a few select destinations is unsustainable long term. To avoid fragility there is a need for a more comprehensive portfolio of sites with a more even distribution pattern that appeal to the local tourist market as well as international visitors (Stewart et al. 2016; Vikneswaran et al. 2016; Fang et al. 2018). COVID-19 has had a shock effect on tourism (Assaf and Scuderi 2020; Ateljevic 2020; Baum and Hai Nguyen Thi 2020; Benjamin et al. 2020). It has long been known that the numbers of international visitors fluctuate widely and the countries they come from changes quickly (Bozzato et al. 2020; Brouder et al. 2020; Chang et al. 2020; Edwards 2020). COVID-19 has reinforced that a static tourism model is not sustainable long term (Balli and Tsui 2016; Bremner and Wikitera 2016; Divisekera 2016; Nair et al. 2016). An alternative paradigm, “slow tourism” (Caffyn, 2012; Guiver and McGrath 2015; Weaver, 2016; Agyeiwaah et al. 2017; Pike 2017), is a better fit to the fundamental goals of geotourism away from the main-stream tourist networks (Gravis et al. 2020).

 

Locations with strong geological heritage in regions that are considered to be outside of major tourism hot spots are commonly disadvantaged. They tend to be remote and not well-promoted despite having high geoheritage values. Post COVID-19 will provide an opportunity to develop local programs focussing on sites with valuable scientific values that are also attractive to tourists (both domestic and international), recreation, conservation or education (Farzanegan et al. 2020; Fong et al. 2020; Galvani et al. 2020). There is also opportunity to include topical subjects such as climate change, sea level rise and geological processes. This is behind a new trend in geoeducation to share research findings to the wider community both within and outside of the formal education systems (Fang et al. 2018).

 Understanding global changes and the geological processes in a New Zealand context will raise awareness about our ever-evolving landscape.  This includes New Zealand’s rising landmass, creation of new coastal regions such as around Kaikoura and the neotectonic processes that occur along a convergent plate margin responsible for earthquakes and tsunamis (Mortimer et al. 2014; Mortimer et al. 2017).

 

In this paper we explore the geological record of rapidly rising terrestrial environments, sea level changes and an accretionary prism that has been active since early Miocene. It is represented by an interbedded succession of mudstone, siltstone and sandstone that are generally considered to be “unattractive” (Palmer and Németh 2018). We will demonstrate the visual attractiveness of the succession, that stretches a few hundreds of kilometres along the eastern part of the North Island of New Zealand, lies with their genesis and significance to the neotectonic development of Zealandia along an active subduction zone.

Geological setting

Wairarapa region is located within the southern part of the East Coast Basin, in southeastern North Island, New Zealand (Fig. 1). It is flanked to the west by the North Island Axial Ranges which separate Wairarapa from Whanganui Basin to the west (Lee and Begg 2002; Mortimer et al. 2014). The region is the onshore part of the deformation front associated with the Hikurangi subduction margin between the Pacific and Australian lithospheric plates (Fig. 1) that has been active since the early Miocene (c. 25Ma)(Cole and Lewis 1981; Nicol et al. 2007; Sutherland et al. 2009). The active subduction margin is located 65-125km east and southeast of the Wairarapa coast (Lee and Begg 2002; Mortimer et al. 2014).

 

Figure 1. Geotectonic situation of the Te Riu-a-Māui/Zealandia dominantly submerged continent. The Te Waipounamu/South Island and Te Ika-a-Māui/North Island are the two largest landmass that together form the mainland of Aotearoa/New Zealand. Dashed rectangle show the Lower North Island of New Zealand presented on figure 2. Abbreviations: Mac – Macquarie Island; Cam – Motu Ihupuku/Campbell Island; Auc – Auckland Island; Ant – Antipodes Islands; Ste – Rakiura/Steward Island; Cha – Rekohu/Wharekauri/Chatham Islands; AF – Alpine Fault; How – Lord Howe Island; Nor – Norfolk Island; Rao – Raoul Island. Bathymetry data is from GNS and NIWA.

 

The geology of Wairarapa is complicated by deformation resulting from uplift of west-dipping thrust faults and congruent folding and the presence of later NE-SW orientated strike-slip faulting (Lee and Begg 2002). This deformation, associated with shortening along the Hikurangi margin, most likely commenced around the inception of the subduction margin c. 25Ma, that intensified in the early Pliocene and continues today (Ilies and Josan 2009).

 

The North Island Axial Ranges are composed of four lithologically distinct sedimentary units collectively known as the Torlesse composite terrane (Mortimer et al. 2014). They are well exposed in the Manawatu Gorge and crop out in upthrown fault blocks farther east (Lee and Begg 2002). These Torlesse rocks consist of intensely deformed, indurated marine sandstones and mudstones of Mesozoic age and are the oldest rocks found in Wairarapa (Fig. 2). The North Island Axial Ranges are a Quaternary feature of the region. Uplift of the ranges began well north of the Manawatu Gorge about 2.4Ma (Rees et al. 2018). At that time a shallow sea extended across the entire region now known as Whanganui and East Coast basins. In the early Quaternary northern and southern landmasses formed dividing the western and eastern regions of the southern North Island (Rees et al. 2018). The region likely became an archipelago with small islands fringed with limestone producing a paleoenvironment that evident in today’s geological architecture of the region.

 

Figure 2. The geological setting of lower Te Ika-a-Māui/North Island. Major towns are marked with red dots. Rock units are ordered by their ages. Rocks older than Cretaceous form the basement that were part of Gondwana. The so-called Axial Ranges, geographically named as Rimutaka, Tararua and Ruahine Ranges dominantly composed of Mezosoic greywackes part of the Torlesse Composite Terrain. East from the Axial Ranges along the axis of the Wairarapa region forearc basin sediments from Pliocene and younger rock assemblages. Along the coastline faulted and folded successions of early Cenozoic to Miocene rocks are part of the accretionary prism. In this succession siliciclastic sediments dominate the exposed rocks among the Lower Miocene thick flysch successions provide the surface rocks as mud and sandstones. Cretaceous and early Cenozoic rocks form the exposed axis of folded rocks. Key locations along the coastline marked with yellow rectangles. Topographic data is from the LINZ NZ 8m Digital Elevation Model (2012).

 

Figure 3. Geological map of the Castlepoint region. Geological data is from the QMap 1:250,000 scale Wairarapa sheet. Topographic data is from GoogleEarth Pro satellite imagery. Geosites mentioned in the text labeled by black rectangles; a – Waipouri’s Mark shore platform with turbidite succession of Whakataki Formation; b – inland exposure of turbidite succession of Whakataki Formation; c – coastal cliffs of the turbidite succession of the Whakataki Formation;  d – coastal cliffs of the turbidite succession of the Whakataki Formation with thick Quaternary cover beds; e – Mataikona coastline with tectonic melange and exotic rocks as olistoliths in the Whakataki Formation; f – reef of the Castlepoint section; g – the castle of the Castlepoint section, h – exposed steeply inclined succession of the Whangai Shale. Little needle-like signs show the bedding attitudes with a number referring to the bed angle. Topographic data is from the LINZ NZ 8m Digital Elevation Model (2012).

 

Wairarapa is located on the overriding Australian Plate of the Hikurangi margin (Fig. 1). Pre-Miocene the region lay in a passive margin tectonic setting where the Cretaceous to Oligocene stratigraphic succession consisted initially of coarse clastic sediments (Fig. 2) that were succeeded by progressively finer grained mudstones and marls (Neef 1995). The tectonic regime changed in the early Miocene with the inception of the plate boundary resulting in regional uplift involving reverse/thrust faulting and congruent folding (Neef 1997; 2000) (Figs. 2 & 3). While much of the area continued to occupy a continental shelf to slope position along the eastern sea border of Zealandia some areas experienced shallowing conditions with terrestrial deposits recorded in places. Where the tectonic regime resulted in the subsidence of basement the increased accommodation space together with a ready sediment supply into the basin resulted in olistostrome, breccia and turbidite sandstone and mudstone deposits (Neef 1992) dominating the Miocene basin fill (Lee and Begg 2002). There was a gradual emergence of Wairarapa from early Pliocene and the development of local marine basins flanked by uplifted blocks led to laterally restricted depositional facies that created a mosaic of marine mudstones and bioclastic limestone deposits over much of the northern Wairarapa (Fig. 3) (Palmer and Németh 2018). A major new sediment type appeared about 1.8Ma with the arrival of calc-alkaline volcanic deposits from the central North Island. Primary air-fall and reworked volcanic ash from some of the largest central North Island volcanic events are found throughout the Wairarapa providing valuable marker beds and chrono-horizons within the succession (Pillans et al. 1993; Shane et al. 1995; Wilson 2001). This volcanism marks the southward propagation of Central Volcanic Zone with the mega-eruption events that occurred from that time. Lee & Begg (2002) suggest this volcanism was a response to reactivation of west to southwest directed subduction along the plate boundary. By the Pleistocene glacio-eustatic sea level fluctuations were a global phenomenon and are recorded in the Quaternary landscapes and sediments of Wairarapa. Wairarapa was never glaciated except for one known small valley glacier in the heart of the axial ranges (Lee and Begg 2002). A distinctive and widespread landscape feature is river terraces which begin their development during cold climates when the tree line is lowered, and areas of high elevation are subjected to intensive erosion. The resultant gravels choke up the river channels and during floods thick successions of fluvial gravels are spread out as sheets across the floodplain.  Known as aggradational gravels these surfaces become incised when the river cuts its channel into the gravels leaving much of the floodplain elevated forming a landscape feature known as a river terrace.  Preservation of these terraces is enhanced by regional uplift that progressively elevates the terraces beyond the reach of the river. During cold climate, glacial, periods the floodplain, barren of vegetation, is subject to strong winds that entrain silt sized material depositing it on the terraces adjacent to the floodplain. This wind-blown silt is known as loess. Where loess preservation is good the age of the terraces can be determined using the number of loess units in the cover beds. Loess preservation is not good everywhere so alternative dating methods such as TL, OSL (Litchfield and Rieser 2005) and tephras (Rees et al. 2018) are used. Marine terraces form near the coast when the sea level is high during interglacial periods. The uplift occurring in Southern Wairarapa means the marine terraces are preserved as they become elevated beyond the reach of future sea-level high-stands.

Geoheritage Values of “mudstone” Country

Mudstones or sandstones are common rock types associated with an accretionary prism or wedge (Meschede 2014). While they are not the most attractive rock types their geological context and features introduce fundamental principles of geology. The rocks are clearly bedded and when the bedding extends into the distance it captures people’s attention and they start to ask why? This is the perfect opportunity to introduce visitors (tourists and locals) to the principles of stratigraphy through visitor experiences (Hose 2006; Lim 2014). The first observation will be the bedding and that each bed (or layer) looks different from that above and below it.  The first principle is that of Original Horizontality (This is based on the premise that sediment particles deposited in water under the influence of gravity form essentially horizontal sediment layers). The fact you can see the layers means these beds are no longer horizontal. First challenge question: how can that be?  This is an opportunity to explain how when layers of rocks are no longer horizontal rocks it means they have been tectonically disturbed after deposition and lithification. But before we talk about tectonics and how it might move rocks let’s first consider which rocks in the layer might be oldest and which are youngest. This is the Law of Superposition and it applies to an undeformed succession of layers of sedimentary rocks. The youngest beds are always at the top and the oldest at the bottom of the succession. Reading features preserved within the layers it is possible to distinguish the bottom and top of the layers. The layers of rock can be seen extending away into the distance and that introduces the Principle of Lateral Continuity. In simple terms, and under normal conditions, deposits originally extend in all directions until they thin, pinch out, or terminate along the edge of the depositional basin. In places, you may observe a fault or intrusion (e.g., igneous) cut through the rock layers. When this happens, the feature cutting through is younger than the rock it cuts; this is called the Law of Cross-Cutting Relationships. Sometimes a layer may include an inclusion of another rock type. These inclusions are older than the rock that contains them (Principle of Inclusion). In any succession of sedimentary rocks, there are always “breaks in the rock record”. These are called unconformities and they result from interruptions to sedimentation by erosion or non-deposition of sediments. Where the layers of rocks have been tilted or deformed between depositional episodes an angular unconformity can form. Here the tilted or folded sedimentary rocks are overlain by more flat-lying strata. The thick successions of bedded rock such as those present in Wairarapa are ideal for showing superb examples of stratigraphic features and unconformities. Even better is when the rocks contain fossils. These preserved remains of ancient organisms are a record of past life on Earth. Finding fossils is the ideal opportunity to open the conversation about how similar or different the organism is from those living today and introduce the concept of Faunal Succession.

 

 

Figure 4. A) Steeply bedded spectacular shore platform of the turbidite succession of the Whakataki Formation at Waipouri’s Mark. Youngening of the section is marked by an arrow. B) Typical Bouma sequence at Waipouri’s Mark with Tb, Tc and Td-e units. C) The Castlepoint view toward the south showing steeply bedded limestone as capping unit. D) The Castlepoint from the south.

 

 

The repetitive nature of the sandstone – mudstone beds reinforces the concept that geological events take place at regular intervals and with similar intensity over long periods of time (Fig. 4A). Sedimentary features characteristic of processes operating during sediment transport (such as low versus high sediment concentration, cohesive versus non-cohesive materials, and low versus high intensity/sediment flux) make such successions very valuable for introducing the concept of subaqueous sediment density flows in geology (Talling et al. 2012) to visitors. Furthermore, the accessibility of the locations makes them useful destinations for educational fieldtrips.

 

Subaqueous sediment density flows are one of the main processes responsible for producing a thick pile of sediments. The sediment clasts are eroded from a landmass and transported across the continental shelf to deep-water bathyal regions. The moving sediment can either be confined within an under-water canyon or unconfined as broad sediment fans that accumulate where there is a change slope such as happens at the foot of the continental slope where it meets the abyssal plain. Extensive, multilobed deposits of sediment build out into the basin forming fans. Such fans and their sedimentary deposits are commonly called flysch and they composed of alternating layers of mudstone, sandstone and other coarse clastic rocks.

 

Flysch successions are often the result of deposition from turbidity currents, gravity driven flow of sediment laden water, and may be referred to as turbidites (Fig. 4B). Turbidity currents are an integral part of transporting sediment from the terrestrial to deep marine environments in what has been described as the sediment conveyor belt. Turbidites are the product of gravity driven sediment flow that form a Bouma sequence (either complete or partial) that have recognized bedding features relating to the flow conditions at the time of deposition (Bouma 1964; Bouma 2000b, a; Bouma and Rozman 2000; Stelting et al. 2000; Bouma 2004; Bouma and Ravenne 2004; Talling et al. 2012). Classical Bouma sequences are preserved in Wairarapa flysch deposits.

 “Mudstone countries” are common features along convergent plate margins where long lasting subduction processes transport copious amounts of sediment from eroding landmasses into the deep marine environment. In such settings, the complete to incomplete Bouma sequences explicitly demonstrate the ruling sedimentary regime that governed the formation of the rocks.

 Well known flysch successions are common along the convergent orogenic belts of the Africa – Eurasia where they provide a geographically continuous, well exposed successions ideal for studying the processes associated with accretionary prism formation along numerous subduction fronts since the Palaeozoic (Robinson et al. 1995; Boiano 1997; Kurz et al. 1998; Nemcok et al. 1998; Baas et al. 2000; Bielik et al. 2004; Sylvester and Lowe 2004; Braid et al. 2010; Forte et al. 2010; D’Errico et al. 2014; Kovac et al. 2017; Wang et al. 2018; Lapcik 2019; Tet’ak et al. 2019). Over time, active accretionary wedges have been uplifted by compressional tectonic events to become part of a newly emerging landmass. Here are the ideal places to study thrust fold belts where flysch deposits are exposed in complex and tectonically active zones.

 Such regions are huge in geoheritage value. The common sandstone and mudstone rocks they host are ideal for educating the public in fundamental geological processes. In addition, the layered nature of the commonly steeply inclined strata produces a visually attractive array of bedded rocks. The layers are fascinating and when presented with appropriate explanations the stratigraphy of the rocks becomes an interesting and challenging task.

 Sharing the geological heritage of geologically significant regions is the premise behind the establishment and promotion of dedicated UNESCO Global Geoparks that have been established based on their geological entity. The Basque Coast UNESCO Global Geopark in Spain is one of the world’s finest examples of geopark established for geoheritage value of the exposed, in places steeply dipping, bedded mudstones and sandstones (Poch and Llordes 2011; Henriques et al. 2019). The geological value of the geopark is centred around the various flysch successions deposited in the upper Triassic (215 million years old) to Middle Eocene (45 million years old) during the opening and closure of the Bay of Biscay (Tenkate and Sprenger 1993; Alegret and Ortiz 2010; Giannetti 2010). The present-day coastline exposes various stages of this major tectonic event recorded in the accumulation of nearly 5000 m thick flysch deposit. The exposed rock records one of the most recent mass extinctions in the Earth’s history commonly referred to as the K/Pg or P/E extinction event (Zili et al. 2009). This extinction event is inferred to be the result of an impact by a large asteroid striking the Earth some 65.5 million years ago in Chicxulub (Mexico). The Basque geopark includes two Global Stratotypes (GSSP) defined by the IUGS – UNESCO at the base of the Selandian and Thanetian stages (Schmitz et al. 1998; Arenillas et al. 2008; Molina et al. 2009). In 2001 the “Ruta del Flysh” was awarded a regional tourism prize from the Basque Government in appreciation of the importance of the geological heritage aspects of the region and for its role in promoting regional geotourism.

 Such a rich geological heritage for a single well-defined and relatively small territorial geopark is rare. Flysch successions play an important role in the geological heritage aspects of other UNESCO Global Geoparks such as the Muroto UNESCO Global Geopark in Japan. The geopark slogan expresses well the nature of the geology of that geopark; “where the ocean and the land meet for the birth of a new habitable land”. While the flysch succession itself has not been defined and promoted as the backbone of the geopark, its significance is evident in the geology of the region.

The Muroto UNESCO Global Geopark is considered a typical accretionary complex. It is Late Paleogene to Early Neogene and located at the plate boundary between the Eurasian Plate and the Philippine Sea Plate. The sediments of the Nankai Trough are attached to the continental plate to form a mass of built-up sedimentary rock, called an accretionary complex (Moore et al. 2015). Exposed to Cape Muroto, visitors can see and touch the succession which provides a rich geoeducational resource for local schools (Yuhora et al. 2016). The region has also experienced magmatic intrusions that have been metamorphosed and uplifted where they have an important role demonstrating how geological processes can present a present-day hazard in the region. This region has close geological similarities with the East Coast of the North Island of New Zealand from a geoheritage and geological processes perspective.

 Rocks that are part of flysch or accretionary prism succession are common in many other UNESCO Global Geoparks and regional geoparks, however, well-exposed and dedicated localities are rare. The region along the East Coast of the North Island of New Zealand has spectacular sites with superb exposures in a modern subduction system. The same area is a focus for geohazard and neotectonic research as it is subjected to regular earthquakes, active uplift and is at risk of tsunami. The geoheritage value of this region needs promotion through geoeducation and geotourism to preserve its unique qualities. It should be also be included in formal conservation policies.

 The East Coast of the North Island is a location that is considered the haven of locals rather a destination for visitors. Consequently, it is underdeveloped and less-visited when compared to other geologically interesting regions of New Zealand. Access to much of the region is by secondary, often unsealed, roads. Many outcrops are accessed along the spectacular coastline.

 

Figure 5. A) Overview from the top of Castlepoint toward Castle Rock. Note the grey mudstone (Rangiwhakaoma Formation – RF) between the cream coloured limestone (Castlepoint Formation) – CF). Also, note the steep bedding toward west (line). Box mark the location of the major load structures shown of E. B) Overview of Castlepoint from the north. Note the steep contacts between the Rangiwhakaoma (RF) and Castlepoint Formation (CF) marking a partially in situ canyon morphology formed by the Pliocene mudstone of RF filled and partially incorporated into the CF. C) Complex fossil assemblage with broken and fractured shells indicate mass flow transportation from differently sourced material into a submarine canyon. D) Fluid-release features indicate a dynamic sedimentary environment commonly affected by sudden mass emplacement over wet, soft sediments commonly accompanied with syn-sedimentary faults. E) Large load casts and blocks (b) and convolute bedding features with the fluid release (fr) structures suggest sudden and energetic mass movement on the seafloor.

 

One popular destination on the east coast is Castlepoint, about 30km east of Masterton.  Castlepoint is a small settlement with about 100 permanent residents with basic service facilities. It is an end of the road location from where key sites pertinent to understanding New Zealand’s most recent geological history can be accessed. Castlepoint itself is a spectacular landscape feature standing 160 m above the present-day sea level.  Its steep eastern side faces directly out to the Pacific Ocean (Fig. 4C). These young strata are  fault bounded on the landward side from Miocene turbidite rocks (Whakataki Formation) (Lee and Begg 2002). This structural relationship has created a castle-shaped cliff, consisting of the Castle and the Reef, and a partially open laguna (Fig. 4D). The main rock types are accessible cliffs in the area that dip gently westward and consist of limey sandstone and fossiliferous limestone. The succession records marine sedimentation about 2 million years ago (Johnston 1973). The fossil assemblages were thought to record environmental changes from a cold marine environment to warmer conditions within the Pleistocene (Johnston 1973), but a recent study provided a different interpretation of the complex geology of Castlepoint (Buckeridge et al. 2018). The main lithostratigraphic unit cropping out in the region is early Pleistocene (early Nukumaruan) Castlepoint Formation (Fig. 5A). It contains a diverse molluscan, bryozoan and barnacle assemblage, with a few brachiopods, corals, echinoids, crabs and polychaetes, comprising mixed warm- and cool-water taxa and outer shelf and shallow-water taxa (Buckeridge et al. 2018). Early studies of this mixture of warm and cool-water faunal assemblage was interpreted as suggesting some sort of time-progressive change from one climate to another and linked to specific stratigraphic horizons. However, a recent study has demonstrated that most of the fossils are mixed together, they are commonly fragmented (Fig. 5B) and accompanying sedimentary structures within the beds include fluid release structures (Fig. 5C), steep original dips, large load casts (Fig. 5D) and displaced blocks (Buckeridge et al. 2018). All of these are better explained by an energetic and very dynamic sedimentary environment. In the latest interpretation it is suggested Castlepoint consists of a rock succession that formed by gravity sliding into a narrow canyon, which collected variously sourced remains of shells and sedimentary grains or semi-consolidated blobs of sediments. The barnacle- and bryozoan-rich molluscan coquina, named as barnamol and bryomol respectively, that makes up much of the Castlepoint Formation comprises a thanatocoenoses resulting from transportation of shallow-water elements, which accumulated deeper-water taxa in transit, to settle in a shelf-edge canyon environment (Buckeridge et al. 2018). This new model fits very well to recent geosystems on the seafloor around the convergent plate margins in the New Zealand realm (Micallef et al. 2014), hence the Castlepoint geosite itself tells a very important story that can be utilized for public geoeducation. The site could also be used to demonstrate the dynamic sea-floor processes commonly linked to major thrust earthquakes along a convergent margin (Reyners 1998), such as the active Hikurangi subduction zone about 100 km away from the present day location of Castlepoint.

 

To the north of Castlepoint, via a sealed road, is the Whakataki Coast where a 20 km long, broad coastal platform exposes the most spectacular turbidite successions in New Zealand. Most of the sites are exceptionally well exposed during low tide enabling visitors to access an impressive succession of tilted mudstone – sandstone interbeds preserving complete to near-complete Bouma sequences (Fig. 6A). The beds record various dynamic processes that operating on the sea floor when sediment is transported from a landmass to deep marine environment over a period of tens of thousands of years. The succession is impressive and includes various turbidity current channels containing central to marginal facies showing how sediments are dispersed from the shelf region into abyssal depths. Individual beds can be traced along the shore platform and facies changes can be observed directly.  Thickly bedded, channelized Bouma sequences, some complete and some partial,  represent turbidite channels or main arteries, whereas in areas where the units are more regularly bedded and missing Ta (sand) sections of Bouma sequences more likely represent overbank of turbidite facies (Field 2005).

 

Figure 6. A) Steeply bedded turbidite succession just north of the Waipouri’s Mark. Bedding dip direction is perpendicular to the location at Waipouri’s Mark. B) Soft sediment deformation features showing dish structures with uniform creeping direction of the dune-bedded Tc unit of the Bouma sequence at Waipouri’s Mark. C) Soft sediment deformation with convolute bedding at Riversdale Beach. D) Syn-depositional fault and soft-sediment deformation features of the flysch succession at Waipouri’s Mark

 

Thin-bedded turbidites of the early Miocene Whakataki Formation are superbly exposed along the shore platform at Mataikona (Lee and Begg 2002). A study on the Waipouri’s Mark shore platform made a detailed record of the number and thickness of beds and their position utilizing aerial photography, field measurements and some geophysical techniques (Field 2005). The detailed study identified about 360 sandstone beds across a 32 m thick section and three (3) lithofacies associations recognised (Field 2005). The data were integrated into a detailed sedimentology study of the location to establish the time and process needed to form the exposed rocks. The study found the Whakataki Formation accumulated over a 3 million year long time period somewhere between 23.5 and 20.5 Ma (Morgans et al. 1996).

 

Bouma sequences show the variability of sediment input, depositional regime and post-depositional processes that operated.  Dewatering features, such as dish structures or creeping sediments on the gently inclined paleo-seafloor, are examples of post-depositional features found in turbidites (Figs. 6B & C). The classic complete Bouma sequence consists of five units. At the base is an erosional contact overlain by Ta, a sandy, coarse-grained unit which is often normally graded. Tb, a stratified, laminated fine sand unit; Tc, with ripples and wavy laminae; Td, a parallel laminated silt and Te, hemipelagic mud overlay the basal unit. This general sequence varies depending upon whether the turbidity current has a low or high particle concentration. Over the years the original turbidite model  has been refined and many variations on the combination of units comprising the turbidite recognised (Talling et al. 2012). Turbidites are complex and their units cannot easily be accurately generalised nevertheless they are valuable geoeducation resources because of their strongly stratified nature. The geosite at Mataikona is superb with its very regular bedding and clearly recognisable features. Here the thickly bedded, basal sand layer (Ta) is mostly absent. Instead most of the turbidites start with a laminated sand (Tb) and are capped by the thinly laminated fine mud (Te). Studies suggest some 1500 m total thickness of the Whakataki Formation accumulated over about 3 million years (Morgans et al. 1996; Field 2005). Using this rate of accumulation, it’s estimated the 32 m thick succession at Mataikona consists of about 360 sandstone beds (proxy for the number of Bouma sequences). This accumulation rate is consistent with a Bouma sequence-forming turbidite event occurring on average each 177 years (Field 2005). This is an exceptional story to tell the general public as it puts into context the concept of deep time and the resulting rock assemblage.

 

Direct observations of active turbidite fans are uncommon. In regions with high sedimentation rates, such as where the Amazon River enters the Atlantic, it’s estimated that turbidity currents could occur frequently. In a study of a 30 ky record it appears there could have been as little as 2 years between turbidite events (Pirmez et al. 2000; Pirmez and Imran 2003). However, most studies are on older deposits and the inferred time periods for turbidite recurrence is far longer. For example, for the Zumaya flysch on the Basque Coast UNESCO Global Geopark recurrence is more likely 5000 – 7000 years (Tenkate and Sprenger 1993). Interestingly, recurrence of turbidite events in the order of thousands of years is a common conclusion in many older flysch successions. Estimates of recurrence intervals for modern turbidite deposits along the Hikurangi Subduction Margin over the last 16,000 years average at 400 years. These turbidites are believed to have been initiated by major mass movement events such as landslides or earthquakes with a frequency of a few hundred years so it is not surprising that the recurrence of each is strikingly similar r (Pouderoux et al. 2014). Twenty of these turbidites exposed along the shore platform between Mataikona and Castlepoint are correlated with subduction interface earthquakes with a calculated average recurrence of c. 800 yr. This aesthetic succession comprises numerous geosites conveying a fascinating and important geological story that will surely fascinate visitors. Its accessibility and the scientific archive of this outstanding coastal section make it a highly valuable geoheritage region.

 

In addition to the turbidite Bouma sequences syn-depositional faults (Fig. 6D) and various soft sediment deformation features can also be seen at these geosites. Such structures form when recently deposited, unconsolidated marine sediments are subjected to intermittent faulting and associated creep at the time of deposition. The succession exposed along the shore platform shows very clearly the dynamic processes operating on the seafloor of an active accretionary prism along an active convergent plate margin.

Figure 7. A) Tectonized section of the flysch succession at Mataikona Beach. B) Bright white olistoliths from the older Eocene marls transported by turbidites that formed the Whakataki Formation. C) Folded and faulted block of an olistolith at Mataikona Beach D) Folded limy mud with boudinage textures from the Whakataki Formation at Mataikona Beach

 

Another very interesting feature along this coast is a structural zone (Fig. 7A) where sediments are fractured, folded (Figs. 7 B & C), displaced and display a minor metamorphic overprint (Fig. 7D). In places rocks from Early to Mid Cenozoic successions have been caught up in the flysch successions forming small to medium size olistoliths of exotic rocks. These olistostromes provide another geological perspective of the active processes that are operating out into the deep marine environment (Lee and Begg 2002).

 

Figure 8. A) Spectacular coastal section south of Castlepoint exposing flysch successions that are commonly folded and faulted. B) Cretaceous Glenburn Formation forming spectacular coastal erosional features such as Honeycomb Rock near Glenburn Station. C) Axis of an anticline exposing Cretaceous Glenburn Formation in the shore platform of Tora Beach. D) Steeply inclined Paleocene well-bedded muddy limestone (Mungaroa Limestone) as part of the Tora Group forming the White Rock. Along the right side of the view the Ewe Fault inferred to runs as a major thrust fault.

 

To the north of Castlepoint there is a spectacular coastal cliff section where the flysch sediments have a strikingly different bedding attitude to those on the coastal platform (Fig. 8A). Cross cutting faults provide a picture of syn-sedimentary tectonic processes operating as the accretionary prism developed in the region.

 

Farther inland the Miocene flysch successions form steeply dipping bedded mudstone and sandstone hills. These stratigraphically younger units tectonically overlie older late Cretaceous marine successions that were deposited along the Gondwana eastern margin. These are the older Glenburn Sandstone, which forms the indurated, more erosion resistant westward dipping set of hills that crop out near and along the coast. While it is not straight forward distinguishing between the Miocene and Cretaceous sandstone and mudstone successions in inland sections, their co-existence in a narrow belt provides an outstanding example of the longevity of siliciclastic sedimentary basin processes. The persistence of similar units through long periods of geological time together with the challenge of identifying individual units is an excellent geoeducational opportunity highlighting the challenges of research in complex geological settings.

 

The Glenburn Sandstone is best exposure along the coast around Glenburn Station. Here Honeycomb Rock (Fig. 2) forms a spectacular coastal cliff erosion (Fig. 8 B) by the pounding Pacific Ocean has created a characteristic honeycomb-pattern in the hard sandstone rock that is texturally different from the Miocene sandstones (Fig. 8B). The Honeycomb Rock itself emerged as a coastal landscape feature about 6500 years ago. It is a well exposed section of Late Cretaceous Glenburn Formation  comprising quartz-rich sandstone (Lee and Begg 2002).  Farther south at Tora Beach (Fig. 2) it is exposed in a spectacular anticline on the shoreline platform (Fig. 8C). Here younger post-Cretaceous rocks can be traced into a stream valley canyon where the K/T boundary is visible in a spectacular gorge (Laird et al. 2003). This section is unique as it provides an insight to the thickness of sediment that has accumulated since Zealandia detached from eastern Gondwana.

 

Another extensive sedimentary unit  present in the east coast of the North Island is the Whangai Formation, informally referred to Whangai Shale (Neef 1995; Lee and Begg 2002). This sedimentary succession accumulated in a deep in places anoxic marine environment. Not as resistant as the Glenburn Formation, the Whangai Shale tends to erode easily commonly forming narrow gullies and steep scree slopes (Parkner et al. 2006). These features impact on farming in the district and the impact the local economy. Locally small quarries for race metal on farms are common. The soils that develop on these shales are fertile but prone to land sliding.

Figure 9. A) GoogleEarth satellite image of the Tora Beach showing a series of raised beaches.  Along the Pukemuri Stream (middle of image) exposes a complete Late Cretaceous to early Cenozoic successions in a spectacular canyon. Along the Pukemuri Stream in the Lower Member of the Manurewa Formation, the K-T boundary is also exposed adding an extra dimension of the geoheritage value of the location. B) Tora Beach raised beach platforms commonly forming wetlands. C) Ngapotiki fan on a GoogleEarth satellite image. D) Ngapotoki fan demonstrates well the mass movement on steep hillslopes of the Mesozoic basement rocks of greywackes along the Aorangi Range. The fan estimated to start to form about 1800 years ago.

 

Southwards along the coastline features with geoheritage potential include White Rock (Fig. 8D) and the eastern margin of the Aorangi Range where there is outstanding exposure of greywacke basement rocks and huge landslide fans. Along the shore platform regular uplift has preserved an outstanding sequence of raised beaches that have been correlated with major earthquakes (Figs. 9A & B). These raised beach sites, such as those at Tora (Fig. 2), provides a unique neotectonic record of the region. Immediately inland there is archaeological evidence of human occupation of the fertile coastal zones and the tsunami events that almost certainly impacted them. These features are pristine due to the remoteness of area and this gives them added geoheritage value such as young landslide fans along the Aorangi coast (Figs. 9C & D).

 

Geoeducational Potential

The East Coast area is accessible by secondary road and has been the destination of day field trips by researchers and students from Massey University Earth Science programs (Figs 10A, B, C & D). The enormous geoeducation potential is undoubted. Student field trips have explored the complex accretionary prism geology from basic geological principles (Supplementary Material 1). The journey from Palmerston North traverses a dynamic landscape from the Axial Ranges comprising Mesozoic (mostly greywacke-dominated) lithologies (Fig. 10A), through the “flysch” succession dominated by Miocene aged rocks (Fig. 10B) to the older harder Cretaceous successions. Pliocene, fossiliferous limestone deposits can be visited at many locations (Fig. 10C) and the impact of some of the largest volume rhyolitic explosive events from the Central North Island have their record as tephra layers embedded in loess sections. We present a field guide of selected key locations that are readily accessible by school or Earth Science groups (Supplementary Material 1). The fieldtrip also provides excellent examples of the active neotectonic setting of the region, in particular strike slip fault systems bounding the axial ranges and accretionary prism deformation. These features complement the outstanding examples in iconic regions such as along the San Andreas Fault in California.

Massey University Earth science undergraduate field courses have undertaken mapping projects in the East Coast, from Makuri east of Dannevirke to Castlepoint Station (Fig. 10D) and south to Honeycomb Rock and White Rock. Here students tackle complex geological problems involving faulting and folding related to accretionary prism development. The experience and knowledge gained during these field courses have led the authors to realise the enormous geoheritage value of the region and the need for it to be promoted to the general public through geoeducation programmes.

The area around Castlepoint provides superb exposures of layered sedimentary rocks ideal sections, both coastal and inland, illustrating the basic principles of geology (Fig. 10D). Geosites access along coastal shore platform is easy at low tide and here exposures of classic Bouma sequences associated with turbidity current deposits are clearly visible (Fig. 10D). There are locations where textbook style exposures can be seen and these could be developed into dedicated geoeducational sites where the basic concepts of geology could be passed to the general public.

 

Figure 10. A) First year BSc Earth Science fieldtrip visit to the Manawatu Gorge’s greywacke sequences. B) First year BSc Earth Science fieldtrip visit to the Whakataki Formation turbidite sequences near Pongaroa. C) First year BSc Earth Science fieldtrip visit to the Pliocene mudstone sequences accumulated in a forearc basin west of Pongaroa. D) Third year geological mapping exercise for BSc Earth Science graduates along the Mataikona Coast, north of Castlepoint.

 

 

Geotourism Potential

The East Coast area of the lower North Island is not one of New Zealand’s iconic tourism destinations despite outstanding aesthetic vistas accompanied by a wealth of geological history on the development of Zealandia. Sealed roads take the traveller from Masterton to Castlepoint in about 1.5 hours (somewhat longer from the SH1 where most of the tourist travel occurs). The region is clearly detached from the main tourism routes so at present has limited facilities available for visitors. To many locals the region is a holiday destination where fishing and surfing are popular activities. There is an annual horse racing event held on the beach and few tourists venture here.  Local councils have developed some information boards in Castlepoint explaining the geology of the reef and explaining the conservation rules. There is very little information available for independent visitors about what to see and where to see or explore. There is an opportunity to promote the fascinating and important geological processes one can see in the rock record at this convergent plate margin locality.

The rich and growing geological information about the region can provide strong basis to consider investment in geotourism. Geotourism has been originally defined as a “type” of tourism, nowadays it is viewed more as a new “approach” to tourism (Newsome and Dowling 2006). As this view gained more traction conceptual models were proposed to define geotourism that can be undertaken either as “independent visits” or “guided tours” (Newsome and Dowling 2010a, b). As its original definition, geotourism was considered as a geology-based tourism and it has been viewed as a the provision of interpretative and service facilities to let visitors gain knowledge to understand the geology and geomorphology of selected sites more than just their pure aesthetic beauties (Hose 1997). In the most recent definition, geotourism is viewed as a form of natural area tourism that specifically focuses on geology and landscape. It promotes tourism to geosites and the conservation of geodiversity and an understanding of Earth sciences through appreciation and learning. This is achieved through independent visits to geological features, use of geo-trails and viewpoints, guided tours, geoactivities and patronage of geosite visitor centres (Newsome and Dowling 2010a). The area presented here has the potential to provide the fundamental elements for geotourism such as the unique geological history, dramatic viewpoints, and specific geologically and geomorphologically significant locations that are specific to geological processes associated with convergent plate margin environment. In this respect, the Wairarapa Mudstone country needs significant attention and work to document and present geosites through educational panels (both on site and through virtual mediums), geotrails (both on site and virtually), establish visitor centres and develop programs where local communities can be engaged in activities revolving around the rich geology of the region. So far neither the local to regional focused initiative nor the development funds to establish the geotourism infrastructure are in place.  

The region has an adequate secondary (unsealed) roading network. However, many of the roads run through private farmland so  future planning needs to include collaboration with local landowners to provide information about the geoheritage value of their land and to establish a relationship whereby increased visitor numbers will be viewed positively and not as an inconvenience.

The geological aspects of geotourism has recently became more important in the modern definition as a niche type of tourism. Geotourism hence needs to include opportunities to visit specifically and geologically centred geosites either independently or through organised (commonly guided) tours. It should include conservation elements and accurately transfer geological knowledge to visitors (Dowling 2014; Dowling 2015). All the necessary fundamentals for geotourism are present in Wairarapa Mudstone Country and priority needs to be given to its promotion in the next few years.

Based on this research work we plan an open discussion with local communities to establish how this work could be expanded, and its results utilized to grow geotourism in the region.

 

Core of a Future Geopark?

The geological setting, outcrop exposure and accessibility of East Coast region of the North Island of New Zealand has all the elements necessary for a future geopark. Currently, New Zealand has only one application for UNESCO Global Geopark status under consideration. The Waitaki Whitestone Aspiring UNESCO Global Geopark is located in the east coast of the South Island of New Zealand. Its geological heritage is strongly linked to the geological evolution of Zealandia since its separation from Gondwana [https://www.whitestonegeopark.nz/], hence its geological context is different to that of the East Coast region of the North Island of New Zealand.

Convergent plate margins, such as the eastern region of the North Island of New Zealand, with a well-exposed sedimentary succession of a typical flysch environment are not uncommon, but a region where the processes are still on going and the succession records sedimentary environment changes since the early Miocene, and before, with superb, easy to access outcrops in a spectacular coastal scene are rare. In addition, very important cultural heritage locations for New Zealand are present in the region. These include archaeological locations where early Māori settlements flourished and contributed to the geocultural landscape evolution for the region (McFadgen 1980; Leach 1981; Chrisp 1993). The eastern coastal regions of New Zealand, including the Wairarapa,  records many catastrophic geological events including earthquakes and tsunamis (Goff et al., 2000; Berryman et al. 2011; Clark et al. 2011; King et al. 2017), many of them likely responsible for the early Māori settlements to decline and move away (Goff and McFadgen 2001; 2002; 2003; Goff et al. 2004; McFadgen and Goff 2007). In the early colonial period, the region functioned as a major sheep and cattle station area, providing strong link to the global economic assets of the world that has transformed dramatically over the last 150 years (Hedley et al. 2010; Day 2011). Now the region seems an abandoned part of New Zealand in spite its globally significant social geology that can be defined as a complex geoheritage of the region. On the basis of this, it is clear, that the region has all that is needed, from geoheritage perspective, to be promoted and structured into a geopark, that might even reach global status if the local communities and the respected governing councils decide to move toward that plan.

 

Conclusion

The East Coast of North Island, New Zealand, lies about 150km east of SH1, the main North Island tourist route of travellers. The region is home to outstanding geological features consistent with its convergent plate tectonic setting. Set in the accretionary prism the exposed succession includes an exceptional sequence of turbidite deposits. Elsewhere in the world, similarly impressive deposits are the foundation of UNESCO global geoparks, such as on the Basque Coast in Spain. To date New Zealand has fostered few geoparks despite the spectacular setting and fascinating geological history of many areas. The East Coast mudstone country, for example, is comparable to the iconic flysch deposits found in the North American Cordillera, the European Alps, the Pyrenees and the Carpathians, and warrants protection and geopark status.

East Coast of the North Island has spectacular sites with superb exposures in a modern subduction system. Regular earthquakes and the associated tsunami risk have left a record of geohazard events including faults, folds and tsunami deposits.  The area has been the focus for geohazard and neotectonic research but its geoheritage value needs promotion through geotourism and formal conservation policies to preserve its unique qualities.

The geoheritage and cultural value of the East Coast needs to be acknowledged through geotourism and geoeducation so the results of the important geohazard and neotectonic research can be appropriately disseminated to both the local communities and visitors.  The region is accessible, and the uninterrupted views of the landscape and coastline means it offers a unique visitor experience to an area of global geological significance.

 

Acknowledgments

We would like to express our gratitude to Mr Anders Crofoot, the current owner the Castlepoint Station for allowing the Massey University’s Geological Mapping Course to run in its land. We also express our thanks for many excellent students been part of the mapping class in the past few years (many of them returned as tutors after their studies completed), particularly Ilmars Gravis, Pete McGowan, Gareth Burton Jones, Phil Pollock, Michael Morgan, James Ardo, Tessa Mitchell-Anyon, Benjamin Pigott, Pete Royle. John Irons’s help in many first-year field classes at Massey University in Wairarapa was always appreciated. We also thanks for the Journal reviewers’ and Editor comments that made this manuscript better.

 

Supplementary Material 1

Field trip study guide for exploring the Wairarapa accretionary wedge (“mudstone country”) is presented in separate file in the website.

References

Agyeiwaah E, McKercher B & Suntikul W (2017). Identifying core indicators of sustainable tourism: A path forward? Tourism Management Perspectives. 24: 26-33.
 - Alegret L & Ortiz S (2010). The Zumaya section (Spain): record of lower Paleogene benthic foraminifers. Revista Mexicana De Ciencias Geologicas. 27 (3): 477-489.
 - AlRayyan K, Hamarneh C, Sukkar H, Ghaith A & Abu-Jaber N (2019). From Abandoned Mines to a Labyrinth of Knowledge: a Conceptual Design for a Geoheritage Park Museum in Jordan. Geoheritage. 11 (2): 257-270.
 - Arenillas I, Molina E, Ortiz S & Schmitz B (2008). Foraminiferal and delta C-13 isotopic event-stratigraphy across the Danian-Selandian transition at Zumaya (northern Spain): chronostratigraphic implications. Terra Nova 20(1): 38-44.
 - Assaf A & Scuderi R (2020) COVID-19 and the recovery of the tourism industry. Tourism Economics 26 (5): 731-733.
 - Ateljevic I (2020) Transforming the (tourism) world for good and (re)generating the potential ‘new normal’. Tourism Geographies. 22 (3): 467-475
 - Baas J.H, van Dam R.L & Storms J.E.A (2000). Duration of deposition from decelerating high-density turbidity currents. Sedimentary Geology. 136 (1-2): 71-88.
 - Balli F & Tsui W.H.K (2016). Tourism Demand Spillovers between Australia and New Zealand: Evidence from the Partner Countries. Journal of Travel Research .55 (6): 804-812.
 - Barasoain D & Azanza B (2018). Geoheritage and Education: a Practical Example from the Rhinoceros of Toril 3 (Calatayud-Daroca Basin, Spain). Geoheritage. 10 (3): 363-374.
 - Baum T & Hai Nguyen Thi T (2020). Hospitality, tourism, human rights and the impact of COVID-19. International Journal of Contemporary Hospitality Management.32(7): 2397-2407.
 - Benjamin S, Dillette A & Alderman D.H (2020). “We can’t return to normal”: committing to tourism equity in the post-pandemic age. Tourism Geographies. 22 (3): 476-483.
 - Beraaouz M, Macadam J, Bouchaou L, Ikenne M, Ernst R, Tagma T & Masrour M (2019). An Inventory of Geoheritage Sites in the Draa Valley (Morocco): a Contribution to Promotion of Geotourism and Sustainable Development. Geoheritage. 11 (2): 241-255.
 - Berryman K, Ota Y, Miyauchi T, Hull A, Clark K, Ishibashi K, Iso N & Litchfield N (2011). Holocene Paleoseismic History of Upper-Plate Faults in the Southern Hikurangi Subduction Margin, New Zealand, Deduced from Marine Terrace Records. Bulletin of the Seismological Society of America. 101 (5): 2064-2087.
 - Bielik M, Sefara J, Kovac M, Bezak V & Plasienka D (2004) .The Western Carpathians - interaction of Hercynian and Alpine processes. Tectonophysics 393 (1-4): 63-86.
 - Bitschene P & Schueller A (2011). Geo-education and geopark implementation in the Vulkaneifel European Geopark. GSA Field Guide 22: 29-34.
 - Boiano U (1997) Anatomy of a siliciclastic turbidite basin: The Gorgoglione Flysch, Upper Miocene, southern Italy: Physical stratigraphy, sedimentology and sequence-stratigraphic framework. Sedimentary Geology 107 (3-4): 231-262.
 - Boley B.B, Nickerson NP & Bosak K (2011). Measuring Geotourism: Developing and Testing the Geotraveler Tendency Scale (GTS). Journal of Travel Research 50 (5): 567-578
 - Bouma A.H (1964). Turbidites. In Bouma AH & Brouwer A (eds) Developments in Sedimentology. Elsevier, pp 247-256.
 - Bouma A.H (2000a). Coarse-grained and fine-grained turbidite systems as end member models: applicability and dangers. Marine and Petroleum Geology. 17 (2): 137-143.
 - Bouma A.H (2000b). Fine-grained, mud-rich turbidite systems: Model and comparison with coarse-grained, sand-rich systems. Fine-Grained Turbidite Systems. 72 (68): 9-19.
 - Bouma A.H (2004). Key controls on the characteristics of turbidite systems. Geological Society, London, Special Publications. 222(1): 9–22.
 - Bouma A.H & Ravenne C (2004). The Bouma Sequence (1962) and the resurgence of geological interest in the French Maritime Alps (1980s): the influence of the Gres d’Annot in developing ideas of turbidite systems. In Lomas PJ & Lomas SA (eds) Deep-Water Sedimentation in the Alpine Basin of Se France: New Perspectives on the Gres D’annot and Related... (pp 27-38).
 - Bouma A.H & Rozman D.J (2000). Characteristics of fine-grained outer fan fringe turbidite systems. Fine-Grained Turbidite Systems. 72 (68): 291-298
 - Bozzato S, Guadagnoli I & Prosperi M (2020). Redefine tourism model. Reflection: critical issues emerged during Covid-19. Documenti Geografici 1: 529-547.
 - Braid J.A, Murphy J.B & Quesada C (2010). Structural analysis of an accretionary prism in a continental collisional setting, the Late Paleozoic Pulo do Lobo Zone, Southern Iberia. Gondwana Research. 17 (2-3): 422-439
 - Bremner H & Wikitera K-A (2016). Using history for tourism or using tourism for history? Examples from Aotearoa/New Zealand. Journal of Tourism History 8 (3): 260-274.
 
- Brilha J (2002). Geoconservation and protected areas. Environmental Conservation. 29 (3): 273-276.
 - Brilha J (2016). Inventory and Quantitative Assessment of Geosites and Geodiversity Sites: a Review. Geoheritage. 8 (2): 119-134
 - Brilha J (2018). Geoheritage: Inventories and Evaluation. In Reynard E & Brilha J (eds) Geoheritage (pp 69-85). Elsevier,
 - Brilha J, Gray M, Pereira D.I & Pereira P (2018). Geodiversity: An integrative review as a contribution to the sustainable management of the whole of nature. Environmental Science & Policy. 86: 19-28.
 - Brocx M & Semeniuk V (2015). Using the Geoheritage Tool-Kit to Identify Inter-related Geological Features at Various Scales for Designating Geoparks: Case Studies from Western Australia. In Errami E, Brocx M & Semeniuk V (eds) From Geoheritage to Geoparks: Case Studies from Africa and Beyond( pp 245-259).Springer.
 - Brocx M & Semeniuk V (2019). Building Stones Can Be of Geoheritage Significance. Geoheritage. 11 (1): 133-149.
 - Brouder P, Teoh S, Salazar N.B, Mostafanezhad M, Pung J.M, Lapointe D, Higgins Desbiolles F, Haywood M, Hall C.M & Clausen H.B (2020). Reflections and discussions: tourism matters in the new normal post COVID-19. Tourism Geographies. 22 (3): 735-746
 - Buckeridge J.S, Beu A.G & Gordon D.P (2018) Depositional environment of the early Pleistocene Castlepoint Formation, New Zealand: a canyon fill in situ. New Zealand Journal of Geology and Geophysics 61 (4): 524-542.
 
- Burlando M, Firpo M, Queirolo C, Rovere A & Vacchi M (2011). From Geoheritage to Sustainable Development: Strategies and Perspectives in the Beigua Geopark (Italy). Geoheritage. 3 (2): 63-72.
 - Caffyn A (2012). Advocating and implementing slow tourism. Tourism Recreation Research. 37(1): 77-80.
 - Chakraborty A, Cooper M & Chakraborty S (2015). Geosystems as a Framework for Geoconservation: the Case of Japan’s Izu Peninsula Geopark. Geoheritage. 7 (4): 351-363.
 - Chang C-L, McAleer M & Ramos V (2020). A Charter for Sustainable Tourism after COVID-19. Sustainability. 12 (9).
 - Chrisp S (1993). The Maori occupation of Wairarapa - Orthodox and nonorthodox versions. Journal of the Polynesian Society. 102 (1): 39-70
 - Clark K.J, Hayward B.W, Cochran U.A, Grenfell H.R, Hemphill-Haley E, Mildenhall DC, Hemphill-Haley M.A & Wallace LM (2011) Investigating subduction earthquake geology along the southern Hikurangi margin using palaeoenvironmental histories of intertidal inlets. New Zealand Journal of Geology and Geophysics. 54 (3): 255-271.
 - Clary R.M (2018).  Geoheritage and public geoliteracy: opportunities for effective geoscience education within US parks. In Dowling R & Newsome D (eds) Handbook of Geotourism (pp 244-253). Edward Elgar Publishing.
 
- Cole J.W & Lewis K.B (1981). Evolution of the Taupo-Hikurangi subduction system. Tectonophysics. 72 (1-2): 1-21.
 - Coratza P, Bruschi V, Piacentini D, Saliba D & Soldati M (2011). Recognition and Assessment of Geomorphosites in Malta at the Il-Majjistral Nature and History Park. Geoheritage. 3 (3): 175-185.
 - Cutler S.Q (2011). Geotourism: The tourism of geology and landscape. Tourism Management. 32 (1): 205-206
 - D’Errico M, Di Staso A, Fornelli A, Guida D, Micheletti F, Perrone V & Raffaelli G (2014). The Numidian Flysch: a guide formation for the reconstruction of the paleogeography and tectono-sedimentary evolution of southern Apennines. Bulletin De La Societe Geologique De France. 185 (5): 343-356.
 
- Day L (2011). Castlepoint Station. HarperCollins Publishers (New Zealand), p 208.
 - de Vries BvW, Byrne P, Delcamp A, Einarson P, Gogus O, Guilbaud M-N, Hagos M, Harangi S, Jerram D, Matenco L, Mossoux S, Nemeth K, Maghsoudi M, Petronis MS, Rapprich V, Rose WI & Vye E (2018). A global framework for the Earth: putting geological sciences in context. Global and Planetary Change. 171: 293-321.
 - De Waele J & Melis M.T (2009). Geomorphology and geomorphological heritage of the Ifrane-Azrou region (Middle Atlas, Morocco). Environmental Geology. 58 (3): 587-599.
 - Divisekera S (2016). Interdependencies of demand for international air transportation and international tourism. Tourism Economics. 22(6): 1191-1206.
 - Doniz-Paez J, Becerra-Ramirez R, Gonzalez-Cardenas E, Guillen-Martin C & Escobar-Lahoz E (2011). Geomorphosites and geotourism in volcanic landscape: the example of La Corona del Lajial cinder cone (El Hierro, Canary Islands, Spain). GeoJournal of Tourism and Geosites. 8(2): 185-197.
 - Dowling R (2011). Geotourism’s Global Growth. Geoheritage. 3(1): 1-13.
 - Dowling R (2014) Geotourism. In Jafari J & Xiao H (eds) Encyclopedia of Tourism (pp 1-3).Cham: Springer International Publishing.
 - Dowling R (2018)  Geotourism and geoparks. In Dowling R & Newsome D (eds) Handbook of Geotourism (pp 276-291). Edward Elgar Publishing.
 - Dowling R & Newsome N (2018). Geotourism: definition, characteristics and international perspectives. In Dowling R & Newsome D (eds) Handbook of Geotourism (pp 1-22).  Edward Elgar Publishing.
 - Dowling R.K (2015). Geotourism. In Carter C, Garrod B & Low T (eds) The encyclopedia of sustainable tourism (pp 231-232). Oxfordshire: CAB International.
 - Eder W (1999). “UNESCO GEOPARKS” - A new initiative for protection and sustainable development of the Earth’s heritage. Neues Jahrbuch für Geologie und Paläontologie-Abhandlungen. 214 (1-2): 353-358.
 - Edwards R.B (2020). Bubble in, bubble out: Lessons for the COVID-19 recovery and future crises from the Pacific. World Development. 135.
 - Erfurt P (2018). Geotourism development and management in volcanic regions. In Dowling R & Newsome D (eds) Handbook of Geotourism (pp 152-167). Edward Elgar Publishing.
 - Erhartic B (2010) Geomorphosite assessment. Acta Geographica Slovenica-Geografski Zbornik. 50 (2): 296-309.
 - Errami E, Brocx M, Semeniuk V & Ennih N (2015). Geosites, Sites of Special Scientific Interest, and Potential Geoparks in the Anti-Atlas (Morocco). In Errami E, Brocx M & Semeniuk V (eds) From Geoheritage to Geoparks: Case Studies from Africa and Beyond ( pp 57-79).
 
- Escorihuela J (2018). The Role of the Geotouristic Guide in Earth Science Education: Towards a More Critical Society of Land Management. Geoheritage. 10 (2): 301-310
 - Fang Y, Yin J & Wu B (2018). Climate change and tourism: a scientometric analysis using CiteSpace. Journal of Sustainable Tourism. 26 (1): 108-126.
 - Farsani N.T, Coelho C & Costa C (2011). Geotourism and Geoparks as Novel Strategies for Socio-economic Development in Rural Areas. International Journal of Tourism Research. 13 (1): 68-81.
 
- Farsani N.T, de Carvalho CN & Xu K (2018). Education as a key tenet of geotourism. In Dowling R & Newsome D (eds) Handbook of Geotourism (pp 234-243). Edward Elgar Publishing.
 - Farzanegan M.R, Gholipour H.F, Feizi M, Nunkoo R & Andargoli AE (2020). International Tourism and Outbreak of Coronavirus (COVID-19): A Cross-Country Analysis. Journal of Travel Research.
 - Field B.D (2005). Cyclicity in turbidites of the Miocene Whakataki Formation, Castlepoint, North Island, and implications for hydrocarbon reservoir modelling. New Zealand Journal of Geology and Geophysics. 48 (1): 135-146.
 
- Fong L.H.N, Law R & Ye B.H (2020). Outlook of tourism recovery amid an epidemic: Importance of outbreak control by the government. Annals of tourism research: 102951-102951.
 - Forte A.M, Cowgill E, Bernardin T, Kreylos O & Hamann B (2010) Late Cenozoic deformation of the Kura fold-thrust belt, southern Greater Caucasus. Geological Society of America Bulletin 122 (3-4): 465-486.
 - Forte J.P, Brilha J, Pereira DI & Nolasco M (2018). Kernel Density Applied to the Quantitative Assessment of Geodiversity. Geoheritage 10 (2): 205-217.
 - Galvani A, Lew A.A & Perez M.S (2020). COVID-19 is expanding global consciousness and the sustainability of travel and tourism. Tourism Geographies 22 (3): 567-576
 - Gerner D, Rybar P, Engel J & Domaracka L (2009). Geotourizm marketing in Lake Constance’ region. Acta Montanistica Slovaca 14 (2): 197-204.
 - Giannetti A (2010) Influence of climate, sea-level changes and tectonics on ichnoassemblages distribution in a carbonate-dominated, deep-marine environment (Upper Paleocene, Zumaya section). Palaeogeography Palaeoclimatology Palaeoecology 285 (1-2): 104-118
 - Goff J, McFadgen .BG & Chague-Goff C (2004) Sedimentary differences between the 2002 Easter storm and the 15th-century Okoropunga tsunami, southeastern North Island, New Zealand. Marine Geology 204 (1-2): 235-250.
 - Goff J.R & McFadgen B.G (2001). Catastrophic seismic-related events and their impact on prehistoric human occupation, coastal New Zealand. Antiquity. 75 (287): 155-162.
 - Goff J.R & McFadgen B.G (2002). Seismic driving of nationwide changes in geomorphology and prehistoric settlement - a 15th Century New Zealand example. Quaternary Science Reviews. 21 (20-22): 2229-2236.
 - Goff J.R & McFadgen B.G (2003). Large earthquakes and the abandonment of prehistoric coastal settlements in 15th century New Zealand. Geoarchaeology-an International Journal 18 (6): 609-623.
 - Goff J.R, Rouse H.L, Jones S L, Hayward B.W, Cochran U, McLea W, Dickinson WW & Morley MS (2000) Evidence for an earthquake and tsunami about 3100-3400 yr ago, and other catastrophic saltwater inundations recorded in a coastal lagoon, New Zealand. Marine Geology .170 (1-2): 231-249.
 - Gordon J.E (2011). Rediscovering a Sense of Wonder: Geoheritage, Geotourism and Cultural Landscape Experiences. Geoheritage: 1-13
 - Gordon J.E (2018). Geoheritage, Geotourism and the Cultural Landscape: Enhancing the Visitor Experience and Promoting Geoconservation. Geosciences. 8(4):136
 - Grant C (2011). Global geotourism perspectives. Tourism Management. 32 (6): 1486-1487.
 - Gravis I, Nemeth K, Twemlow C & Nemeth B (2020) . The Case for Community-Led Geoheritage and Geoconservation Ventures in Mangere, South Auckland, and Central Otago, New Zealand. Geoheritage. 12 (1).
 - Gray M (2018a). Geodiversity: The Backbone of Geoheritage and Geoconservation. In Reynard E & Brilha J (eds) Geoheritage (pp 13-25). Elsevier.
 - Gray M (2018b). Geodiversity, geoheritage, geoconservation and their relationship to geotourism. In Dowling R & Newsome D (eds) Handbook of Geotourism (pp 48-60). Edward Elgar Publishing.
 - Guiver J & McGrath P (2015). Slow tourism. Exploring the discources. Dos Algarves: A Multidisciplinary e-Journal. 27: 11-34.
 - Guo W & Chung S (2019) Using Tourism Carrying Capacity to Strengthen UNESCO Global Geopark Management in Hong Kong. Geoheritage 11 (1): 193-205.
 - Harmon B & Viles H (2013) Beyond geomorphosites: trade-offs, optimization, and networking in heritage landscapes. Environment Systems & Decisions. 33 (2): 272-285.
 - Hayes S & Lovelock B (2017). ‘Demystifying’ worldmaking: exploring New Zealand’s clean and green imaginary through the lens of angling tourists. Tourism Recreation Research. 42 (3): 380-391.
 - Hedley A, Laing P & Laing N (2010). Castlepoint: Coastal Station and Settlement (p 272). Hedley’s Bookshop.
 - Henriques M.H & Brilha J (2017). UNESCO Global Geoparks: a strategy towards global understanding and sustainability. Episodes. 40(4): 349-355.
 - Henriques M.H, dos Reis R.P, Brilha J & Mota T (2011) . Geoconservation as an Emerging Geoscience. Geoheritage. 3(2): 117-128.
 - Henriques M.H, Luisa Canales M, Garcia-Frank A & Gomez-Heras M (2019) Accessible Geoparks in Iberia: a Challenge to Promote Geotourism and Education for Sustainable Development. Geoheritage. 11 (2): 471-484.
 - Hose T.A (1997) Geotourism - Selling the earth to Europe. In Marinos PG, Koukis GC, Tsiambaos GC & Stournaras GC (eds) Engineering Geology and the Environment (Vols 1-3, pp 2955-2960). Balkema, Rotterdam.
 - Hose T.A (2006) Geotourism and interpretation. In Dowling R & Newsome D (eds) Geotourism. Oxford: Elsevier Butterworth-Heinemann.
 - Hose T.A (2012a). 3G’s for Modern Geotourism. Geoheritage.4: 1-18.
 - Hose T.A (2012b). Geotourism and Geoconservation. Geoheritage. 4(1-2): 1-5.
 - Hose T.A, Markovic SB, Komac B & Zorn M (2011) . Geotourism - a short introduction. Acta Geographica Slovenica-Geografski Zbornik. 51(2): 339-342.
 - Ielenicz M (2009). Geotope, geosite, geomorphosite. The Annals of Valahia Universtiy of Targoviste, Geographical Series 9/2009: 7-22.
 - Ilies D & Josan N (2009). Geosites-geomorphosites and relief. GeoJournal of Tourism and Geosites. 3 (1): 78-85
 - Johnston M.R (1973) Geology of Castlepoint headland and reef, Wairarapa, New Zealand. New Zealand Journal of Geology and Geophysics. 16 (4): 909-916.
 - Joyce B (2010a). Volcano tourism in the new Kanawinka Global Geopark of Victoria and SE South Australia. In Erfurt-Cooper P & Cooper M (eds) Volcano and Geothermal Tourism, Sustainable Geo-Resources for Leisure and Recreation (pp 302-311). London, UK: Earthscan.
 - Joyce E (2010b). Australia’s Geoheritage: History of Study, A New Inventory of Geosites and Applications to Geotourism and Geoparks. Geoheritage. 2 (1-2): 39-56.
 - Justice S.C (2018). UNESCO Global Geoparks, Geotourism and Communication of the Earth Sciences: A Case Study in the Chablais UNESCO Global Geopark, France. Geosciences. 8 (5):149.
 - Keever P.J.M & Zouros N (2005). Geoparks: Celebrating Earth heritage, sustaining local communities. Episodes 28 (4): 274-278.
 - King D.N, Goff J.R, Chague-Goff C, McFadgen B, Jacobsen GE, Gadd P & Horrocks M (2017) Reciting the layers: Evidence for past tsunamis at Mataora-Wairau Lagoon, Aotearoa-New Zealand. Marine Geology 389: 1-16.
 - Kovac M, Marton E, Oszczypko N, Vojtko R, Hok J, Kralikova S, Plasienka D, Kluciar T, Hudackova N & Oszczypko-Clowes M (2017.) Neogene palaeogeography and basin evolution of the Western Carpathians, Northern Pannonian domain and adjoining areas. Global and Planetary Change 155: 133-154.
 - Kurz W, Neubauer F, Genser J & Dachs E (1998). Alpine geodynamic evolution of passive and active continental margin sequences in the Tauern Window (eastern Alps, Austria, Italy): a review. Geologische Rundschau 87 (2): 225-242
 - Lagally U (2018). Earth sciences for everybody as promotion of geotourism: the ‘Day of Geosites’ in Germany In Dowling R & Newsome D (eds) Handbook of Geotourism (pp 254-263). Edward Elgar Publishing.
 - Laird M.G, Bassett K.N, Schioler P, Morgans H.E.G, Bradshaw J.D & Weaver S.D (2003) Paleoenvironmental and tectonic changes across the Cretaceous/Tertiary boundary at Tora, southeast Wairarapa, New Zealand: a link between Marlborough and Hawke’s Bay. New Zealand Journal of Geology and Geophysics 46 (2): 275-293.
 - Lapcik P (2019). Facies anatomy of a progradational submarine channelized lobe complex: semi-quantitative analysis of the Ropianka Formation (Campanian-Paleocene) in Hucisko Jawornickie section, Skole Nappe, Polish Carpathians. Acta Geologica Polonica. 69 (1): 111-141.
 - Leach B.F (1981). The prehistroy of the southern Wairarapa. Journal of the Royal Society of New Zealand. 11 (1): 11-33.
 - Lee J.K, Woo KS & Lee H.K (2019). Present status and future improvement directions of educational signboards at geoheritage sites in protected areas of Korea. Journal of the Geological Society of Korea. 55 (1): 37-57.
 - Lee J.M & Begg J.G (2002). Geology of the Wairarapa area: scale 1:250,000. Lower Hutt: Institute of Geological & Nuclear Sciences Limited. Institute of Geological & Nuclear Sciences 1:250,000 geological map 11. 66 p. + 1 folded map
 - Letos D & Muntele I (2011). The capitalization of the touristic potential of Cozla mountain. GeoJournal of Tourism and Geosites. 7 (1): 143-153.
 - Lim K (2014). A study of Geotourism Growth through Recognition of Geoeducation and Geoconservation for the Geoheritage. Journal of Tourism and Leisure Research. 26(6): 43-59.
 - Litchfield NJ & Rieser U (2005). Optically stimulated luminescence age constraints for fluvial aggradation terraces and loess in the eastern North Island, New Zealand. New Zealand Journal of Geology and Geophysics. 48 (4): 581-589.
 
- Maghsoudi M, Moradi A, Moradipour F & Nezammahalleh MA (2019). Geotourism Development in World Heritage of the Lut Desert. Geoheritage. 11 (2): 501-516
 - McFadgen B.G (1980) A stone row system at Okoropunga on the southeast Wairarapa coast and inferences about coastal stone rows elsewhere in Central New Zealand. New Zealand Journal of Science. 23 (2): 189-197.
 - McFadgen B.G & Goff J.R (2007) Tsunamis in the New Zealand archaeological record. Sedimentary Geology 200 (3-4): 263-274.
 - Melendez G (2018) Palaeontological Heritage in Spain: from Conservation to Education. Geoheritage. 10 (3): 341-342.
 - Meschede M (2014). Accretionary Wedge. In Harff J, Meschede M, Petersen S & Thiede J (eds) Encyclopedia of Marine Geosciences. Dordrecht: Springer.
 - Micallef A, Mountjoy J.J, Barnes P.M, Canals M & Lastras G (2014). Geomorphic response of submarine canyons to tectonic activity: Insights from the Cook Strait canyon system, New Zealand. Geosphere. 10 (5): 905-929.
 - Miccadei E, Piacentini T & Esposito G (2011). Geomorphosites and Geotourism in the Parks of the Abruzzo Region (Central Italy). Geoheritage. 3(3): 233-251.
 - Migoń P (2018). Geo-interpretation: how and for whom? In Dowling R & Newsome D (eds) Handbook of Geotourism (pp 224-233). Edward Elgar Publishing.
 - Migon P & Pijet-Migon E (2016). Geoconservation and tourism at geothermal sites - lessons learnt from the Taupo Volcanic Zone, New Zealand. Proceedings of the Geologists Association. 127(3): 413-421.
 - Migon P & Pijet-Migon E (2019). Natural Disasters, Geotourism, and Geo-interpretation. Geoheritage. 11(2): 629-640.
 - Mikhailenko A.V & Ruban D.A (2019) Geo-Heritage Specific Visibility as an Important Parameter in Geo-Tourism Resource Evaluation. Geosciences. 9 (4): 146.
 - Molina E, Alegret L, Arenillas I, Arz J.A, Gallala N, Manuel Grajales-Nishimura J, Murillo-Muneton G & Zaghbib-Turki D (2009). The Global Boundary Stratotype Section and Point for the base of the Danian Stage (Paleocene, Paleogene, “Tertiary”, Cenozoic): auxiliary sections and correlation. Episodes. 32 (2): 84-95.
 
- Moore G.F, Boston B.B, Strasser M, Underwood M.B & Ratliff RA (2015). Evolution of tectono-sedimentary systems in the Kumano Basin, Nankai Trough forearc. Marine and Petroleum Geology. 67: 604-616.
 - Morgans H.E.G, Scott G.H, Beu A.B, Graham I.J, Mumme T.C, St George W & Strong C.P (1996). New Zealand Cenozoic time scale. Insitute of Geological and Nuclear Sciences science report 96/38.
 - Morrow J & Mowatt S (2015) The Implementation of Authentic Sustainable Strategies: i-SITE Middle Managers, Employees and the Delivery of 100% Pure New Zealand. Business Strategy and the Environment. 24(7): 656-666.
 - Mortimer N, Campbell HJ, Stagpoole M, Wood RA, Rattenbury MS, Sutherland R & Seton M (2017) Zealandia: Earth’s hidden continent. GSA Today. 27(3): 27-35.
 - Mortimer N, Rattenbury M.S, King PR, Bland KJ, Barrell D.J.A, Bache F, Begg J.G, Campbell H.J, Cox SC, Crampton JS, Edbrooke SW, Forsyth P.J, Johnston MR, Jongens R, Lee J.M, Leonard GS, Raine J.I, Skinner D.N.B, Timm C, Townsend D.B, Tulloch A.J, Turnbull I.M & Turnbull R.E (2014) High-level stratigraphic scheme for New Zealand rocks. New Zealand Journal of Geology and Geophysics. 57(4): 402-419.
 
- Nair V, Mohamed B & Thomas T.K (2016) How can the tourism industry respond to the global challenges arising from climate change and environmental degradation? Worldwide Hospitality and Tourism Themes. 8(5): 606-610.
 - Neef G (1992) Turbidite deposition in 5 Miocene, bathyal formations along an active plate margin, North Island, New Zealand - with notes on styles of deposition at the margins of east-coast bathyal basins. Sedimentary Geology. 78 (1-2): 111-136.
 - Neef G (1995). Cretaceous and Cenozoic geology east of the Tinui Fault Complex in northeastern Wairarapa, New Zealand. New Zealand Journal of Geology and Geophysics. 38(3): 375-394.
 - Neef G (1997). Stratigraphy, structural evolution, and tectonics of the northern part of the Tawhero basin and adjacent areas, northern Wairarapa, North Island, New Zealand. New Zealand Journal of Geology and Geophysics. 40 (3): 335-358.
 - Neef G (2000). Late Miocene marine fossil-rich, rock-fall, avalanche, mud-flow and debris-flow deposits adjoining and near the western margin of the Tawhero Basin, outer forearc North Island, New Zealand. Transactions of the Royal Society of Edinburgh-Earth Sciences 90: 189-201.
 - Nemcok M, Pospisil L, Lexa J & Donelick RA (1998). Tertiary subduction and slab break-off model of the Carpathian-Pannonian region. Tectonophysics. 295 (3-4): 307-340.
 - Németh K & Moufti M.R (2017). Geoheritage Values of a Mature Monogenetic Volcanic Field in Intra-continental Settings: Harrat Khaybar, Kingdom of Saudi Arabia. Geoheritage. 9 (3): 311-328.
 - Németh K, W.u J, Sun C & Liu J (2017) Update on the Volcanic Geoheritage Values of the Pliocene to Quaternary Arxan-Chaihe Volcanic Field, Inner Mongolia, China. Geoheritage. 9 (3): 279-297.
 - Newsome D & Dowling R (2006) The scope and nature of geotourism. In Dowling R & Newsome D (eds) Geotourism (pp 3-25). Oxford, England: Elsevier.
 - Newsome D & Dowling R (2018) Ch. Geoheritage and Geotourism. In Reynard E & Brilha J (eds) Geoheritage (pp 305-321). Elsevier.
 - Newsome D & Dowling R.K (2010a). Setting an agenda for geotourism. In Newsome D & Dowling RK (eds) Geotourism: the tourism of geology and landscape. Oxford: Goodfellow Publishers .
 - Newsome D & Dowling RK (2010b). The future of geotourism where to from here. In Newsome D & Dowling RK (eds) Geotourism: the tourism of geology and landscape. Oxford: Goodfellow Publishers.
 - Nicol A, Mazengarb C, Chanier F, Rait G, Uruski C & Wallace L (2007) Tectonic evolution of the active Hikurangi subduction margin, New Zealand, since the Oligocene. Tectonics. 26(4):1-24.
 - O’Connor P.J (2008). The role of geotourism in supporting regeneration in disadvantaged rural communities in Ireland. In Pineda FD & Brebbia CA (eds) Sustainable Tourism Iii ( pp 267-275)
 - Orchiston C & Higham J.E.S (2016). Knowledge management and tourism recovery (de)marketing: the Christchurch earthquakes 2010-2011. Current Issues in Tourism. 19 (1): 64-84.
 - Palmer J & Németh K (2018) Geoheritage values and geotourism potential of the Wairarapa “mudstone country” and Manawatu “terrace-land”. Geoscience Society of New Zealand Miscellaneous Publication. 151A: 209
 - Panizza M (2001). Geomorphosites: Concepts, methods and examples of geomorphological survey. Chinese Science Bulletin 46: 4-6.
 - Parkner T, Page M.J, Marutani T & Trustrum N.A (2006). Development and controlling factors of gullies and gully complexes, East Coast, New Zealand. Earth Surface Processes And Landforms. 31 (2): 187-199.
 - Pearce D.G (2000). Tourism plan reviews: methodological considerations and issues from Samoa. Tourism Management. 21 (2): 191-203.
 - Pellitero R, González-Amuchastegui M, Ruiz-Flaño P & Serrano E (2011). Geodiversity and Geomorphosite Assessment Applied to a Natural Protected Area: the Ebro and Rudron Gorges Natural Park (Spain). Geoheritage. 3 (3): 163-174.
 - Pereira P & Pereira D (2010). Methodological guidelines for geomorphosite assessment. Geomorphologie-Relief Processus Environnement. (2): 215-222.
 - Petrovic M.D, Vasiljevic D.A, Vujicic M.D, Hose T.A, Markovic SB & Lukic T (2013). Global geopark and candidate - Comparative analysis of Papuk Mountain Geopark (Croatia) and Fruska Gora Mountain (Serbia) by using GAM model. Carpathian Journal of Earth and Environmental Sciences. 8 (1): 105-116.
 - Pike S (2017). Destination positioning and temporality: Tracking relative strengths and weaknesses over time. Journal of Hospitality and Tourism Management 31: 126-133.
 - Pillans B, McGlone M, Palmer A, Mildenhall D, Alloway B & Berger G (1993). The Last Glacial Maximum In Central And Southern North Island, New-Zealand - A Paleoenvironmental Reconstruction Using The Kawakawa Tephra Formation As A Chronostratigraphic Marker. Palaeogeography Palaeoclimatology Palaeoecology 101 (3-4): 283-304.
 - Pirmez C, Beaubouef R.T, Friedmann S.J, Mohrig D.C & Weimer P (2000) Equilibrium Profile and Baselevel in Submarine Channels: Examples from Late Pleistocene Systems and Implications for the Architecture of Deepwater Reservoirs. In  Deep-Water Reservoirs of the World (pp 782-805 ). SEPM Society for Sedimentary Geology.
 - Pirmez C & Imran J (2003). Reconstruction of turbidity currents in Amazon Channel. Marine and Petroleum Geology. 20 (6-8): 823-849
 - Poch J & Llordes J.P (2011). The Basque Coast Geopark: support for good practices in geotourism. GeoJournal of Tourism and Geosites. 8(2): 272-280.
 - Pouderoux H, Proust J-N & Lamarche G (2014). Submarine paleoseismology of the northern Hikurangi subduction margin of New Zealand as deduced from Turbidite record since 16 ka. Quaternary Science Reviews. 84: 116-131.
 - Prosser C.D (2018). Geoconservation, Quarrying and Mining: Opportunities and Challenges Illustrated Through Working in Partnership with the Mineral Extraction Industry in England. Geoheritage. 10(2): 259-270.
 - Rapprich V, Valenta J, Broz M, Kadlecova E, de Vries BvW, Petronis MS & Rojik P (2019) A Crucial Site in the Argument Between Neptunists and Plutonists: Reopening of the Historical Adit in the Komorni hurka (Kammerbuhl) Volcano After 180 Years. Geoheritage. 11 (2): 347-358.
 - Rees C, Palmer J & Palmer A (2018). Plio-Pleistocene geology of the Lower Pohangina Valley, New Zealand. New Zealand Journal of Geology and Geophysics. 61 (1): 44-63.
 - Reynard E & Brilha J (2018). Geoheritage: A Multidisciplinary and Applied Research Topic. In Reynard E & Brilha J (eds) Geoheritage (pp 3-9). Elsevier.
 - Reynard E, Coratz P & Regolini-Bissing G (2009) Geomorphosites. In  Verlag Dr. Friedrich Pfeil, München, Germany ( p 240).
 - Reynard E & Coratza P (2013). Scientific research on geomorphosites. A review of the activities of the IAG Working Group on Geomorphosites over the last twelve years. Geografia Fisica E Dinamica Quaternaria. 36 (1): 159-168.
 - Reynard E, Coratza P, Cayla N, Clivaz M, Comanescu L, Darbellay L, Giusti C, Grecu F, Hoblea F & Pereira P (2018). InterGEO: a Digital Platform for University Education on Geomorphological Heritage. Geoheritage 10 (4): 645-657.
 
- Reynard E, Coratza P & Giusti C (2011a) Geomorphosites and geotourism. Geoheritage. 3 (3): 129-130.
 
- Reynard E, Fontana G, Kozlik L & Scapozza C (2007) A method for assessing «scientific» and «additional values» of geomorphosites. Geographica Helvetica. 62(3): 148-158
 
- Reynard E, Hoblea F, Cayla N & Gauchon C (2011b). Iconic sites for Alpine geology and geomorphology Rediscovering Heritage? Revue De Geographie Alpine-Journal of Alpine Research. 99 (1-4): 188-202
 - Reyners M (1998). Plate coupling and the hazard of large subduction thrust earthquakes at the Hikurangi subduction zone, New Zealand. New Zealand Journal of Geology and Geophysics 41 (4): 343-354
 - Robinson AG, Banks C.J, Rutherford M.M & Hirst J.P.P (1995) Stratigraphic and Structural Development of the Eastern Pontides, Turkey. Journal of the Geological Society. 152: 861-872
 - Rozycka M & Migon P (2018). Customer-Oriented Evaluation of Geoheritage-on the Example of Volcanic Geosites in the West Sudetes, SW Poland. Geoheritage 10 (1): 23-37.
  Sa dos Santos W.F, Carvalho IdS, Brilha J.B & Leonardi G (2016). Inventory and assessment of palaeontological sites in the Sousa Basin (Paraiba, Brazil): preliminary study to evaluate the potential of the area to become a geopark. Geoheritage. 8 (4): 315-332.
 - Schmitz B, Molina E & von Salis K (1998) The Zumaya section in Spain: A possible global stratotype section for the Selandian and Thanetian stages. Newsletters on Stratigraphy 36 (1): 35-42.
 - Serrano E, Jose J & Trueba G (2011). Environmental education and landscape leisure. Geotourist map and geomorphosites in the Picos de Europa National Park. GeoJournal of Tourism and Geosites. 8 (2): 295-308.
 - Shane P, Froggatt P, Black T & Westgate J (1995). Chronology of Pliocene and Quaternary bioevents and climatic events from fission-track ages on tephra beds, Wairarapa, New Zealand. Earth and Planetary Science Letters 130 (1-4): 141-154.
 - Stelting C.E, Bouma A.H & Stone CG (2000). Fine-grained turbidite systems: Overview. In Bouma AH & Stone CG (eds) Fine-Grained Turbidite Systems( pp 1-7).
 - Stewart E.J, Wilson J, Espiner S, Purdie H, Lemieux C & Dawson J (2016). Implications of climate change for glacier tourism. Tourism Geographies 18 (4): 377-398
 - Sutherland R, Stagpoole V, Uruski C, Kennedy C, Bassett D, Henrys S, Scherwath M, Kopp H, Field B, Toulmin S, Barker D, Bannister S, Davey F, Stern T & Flueh ER (2009) Reactivation of tectonics, crustal underplating, and uplift after 60 Myr of passive subsidence, Raukumara Basin, Hikurangi-Kermadec fore arc, New Zealand: Implications for global growth and recycling of continents. Tectonics 28.
 
- Suzuki D.A & Takagi H (2018). Evaluation of Geosite for Sustainable Planning and Management in Geotourism. Geoheritage. 10 (1): 123-135
 - Sylvester Z & Lowe D.R (2004). Textural trends in turbidites and slurry beds from the Oligocene flysch of the East Carpathians, Romania. Sedimentology 51 (5): 945-972.
 - Talling P.J, Masson D.G, Sumner E.J & Malgesini G (2012) Subaqueous sediment density flows: Depositional processes and deposit types. Sedimentology 59 (7): 1937-2003.
 - Tenkate W & Sprenger A (1993) Orbital cyclisites above and below the Cretaceous Paleogene boundary at Zumaya (N SPain), Agost and Relleu (SE Spain). Sedimentary Geology 87 (1-2): 69-101.
 - Tet’ak F, Pivko D & Kovacik M (2019). Depositional systems and paleogeography of Upper Cretaceous-Paleogene deep-sea flysch deposits of the Magura Basin (Western Carpathians). Palaeogeography Palaeoclimatology Palaeoecology 533: 109250
 - Tucker H, Shelton EJ & Bae H (2017). Post-disaster tourism: Towards a tourism of transition. Tourist Studies .17 (3): 306-327.
 
- Turner S (2013). Geoheritage and Geoparks: One (Australian) Woman’s Point of View. Geoheritage. 5 (4): 249-264
 - UNESCO (2018) UNESCO Global Geoparks. Retrieved from: http://www.unesco.org/new/en/natural-sciences/environment/earth-sciences/unesco-global-geoparks/
 - Vasiljevic D.A, Vujicic M.D, Bozic S, Jovanovic T, Markovic S.B, Basarin B, Lukic T & Carkadzic J (2018) Trying to underline geotourist profile of National park visitors: Case study of NP Fruska Gora, Serbia (Typology of potential geotourists at NP Fruska Gora). Open Geosciences. 10 (1): 222-233.
 - Vikneswaran N, Badaruddin M & Thomas T.K (2016). How can the tourism industry respond to the global challenges arising from climate change and environmental degradation? Worldwide Hospitality and Tourism Themes. 8 (5): 606-610.
 - Wang H.J, Rahn M & Zhou J (2018) Tectonothermal evolution of the Triassic flysch in the Bayan Har Orogen, Tibetan plateau. Tectonophysics 723: 277-287.
 - Warowna J, Zglobicki W, Kolodynska-Gawrysiak R, Gajek G, Gawrysiak L & Telecka M (2016). Geotourist values of loess geoheritage within the planned Geopark Malopolska Vistula River Gap, E Poland. Quaternary International. 399: 46-57.
 - Weaver D.B (2016) Beyond convention: reimagining indigenous tourism. Tourism Geographies. 18 (1): 28-37
- Wilson C.J.N (2001) The 26.5 ka Oruanui eruption, New Zealand: an introduction and overview. Journal of Volcanology and Geothermal Research. 112(1-4): 133-174
 - Yuhora K, Yamashita S & Takahashi S (2016). Geopark Education Adopts a Geographical Viewpoint at Muroto Senior High School. Journal of Geography-Chigaku Zasshi. 125 (6): 813-829.
 
- Zecha S & Regelous A (2018). Promoting Geodiversity Education by Using EarthCaching in National Geoparks. Geoheritage. 10 (4): 637-643.
 - Zeng F (2014). An evaluation of residents’ perceptions of the creation of a geopark: a case study on the geopark in Mt. Huaying Grand Canyon, Sichuan Province, China. Environmental Earth Sciences. 71 (3): 1453-1463.
 - Zglobicki W & Baran-Zglobicka B (2013). Geomorphological Heritage as a Tourist Attraction. A Case Study in Lubelskie Province, SE Poland. Geoheritage. 5 (2): 137-149.
 - Zili L, Zaghbib-Turki D, Alegret L, Arenillas I & Molina E (2009). Foraminiferal turnover across the Paleocene/Eocene boundary at the Zumaya section, Spain: record of a bathyal gradual mass extinction. Revista Mexicana De Ciencias Geologicas. 26 (3): 729-744.