- Home
- A-Z Publications
- Basin Research
- Previous Issues
- Volume 31, Issue 3, 2019
Basin Research - Volume 31, Issue 3, 2019
Volume 31, Issue 3, 2019
-
-
Tectonic and oceanographic process interactions archived in Late Cretaceous to Present deep‐marine stratigraphy on the Exmouth Plateau, offshore NW Australia
AbstractDeep‐marine deposits provide a valuable archive of process interactions between sediment gravity flows, pelagic sedimentation and thermohaline bottom‐currents. Stratigraphic successions can also record plate‐scale tectonic processes (e.g. continental breakup and shortening) that impact long‐term ocean circulation patterns, including changes in climate and biodiversity. One such setting is the Exmouth Plateau, offshore NW Australia, which has been a relatively stable, fine‐grained carbonate‐dominated continental margin from the Late Cretaceous to Present. We combine extensive 2D (~40,000 km) and 3D (3,627 km2) seismic reflection data with lithologic and biostratigraphic information from wells to reconstruct the tectonic and oceanographic evolution of this margin. We identified three large‐scale seismic units (SUs): (a) SU‐1 (Late Cretaceous)—500 m‐thick, and characterised by NE‐SW‐trending, slope‐normal elongate depocentres (c. 200 km long and 70 km wide), with erosional surfaces at their bases and tops, which are interpreted as the result of contour‐parallel bottom‐currents, coeval with the onset of opening of the Southern Ocean; (b) SU‐2 (Palaeocene—Late Miocene)—800 m‐thick and characterised by: (a) very large (amplitude, c. 40 m and wavelength, c. 3 km), SW‐migrating, NW‐SE‐trending sediment waves, (b) large (4 km‐wide, 100 m‐deep), NE‐trending scours that flank the sediment waves and (c) NW‐trending, 4 km‐wide and 80 m‐deep turbidite channel, infilled by NE‐dipping reflectors, which together may reflect an intensification of NE‐flowing bottom currents during a relative sea‐level fall following the establishment of circumpolar‐ocean current around Antarctica; and (c) SU‐3 (Late Miocene—Present)—1,000 m‐thick and is dominated by large (up to 100 km3) mass‐transport complexes (MTCs) derived from the continental margin (to the east) and the Exmouth Plateau Arch (to the west), and accumulated mainly in the adjacent Kangaroo Syncline. This change in depositional style may be linked to tectonically‐induced seabed tilting and folding caused by collision and subduction along the northern margin of the Australian plate. Hence, the stratigraphic record of the Exmouth Plateau provides a rich archive of plate‐scale regional geological events occurring along the distant southern (2,000 km away) and northern (1,500 km away) margins of the Australian plate.
-
-
-
Along‐strike variability in shelf‐margin morphology and accretion pattern: An example from the northern margin of the South China Sea
Authors Haiteng Zhuo, Yingmin Wang, Zhen Sun, Ying Wang, Qiang Xu, Pengfei Hou, Xingxing Wang, Zhongxian Zhao, Wei Zhou and Shang XuAbstractInterplays among diachronous tectonism, uneven sediment supply, and local marine hydraulic processes make the northern margin of the South China Sea (SCS) an ideal location to investigate the complexity of along‐strike variability in shelf margins. This study examines shelf‐margin morphology, stratigraphy, and sedimentation from the northern SCS using multichannel seismic reflection profiles complemented with the data from commercial and ocean drilling sites. Analysis of seismic reflection profiles reveals three categories of shelf‐margin cross‐sectional profiles, the concave‐up, linear, and sigmoidal, according to which five margin sectors were recognized. Results show that these margin segments differ in relief, shelf‐edge trajectory, submarine canyon development, and long‐term accretion pattern. The westernmost margin sector, or the Yinggehai (YGH)‐western Qiongdongnan (QDN) margin, has appeared to be supply dominated since its commencement at ca. 10.5 Ma, which is characterized by well‐developed prograding clinoforms, low‐angle shelf‐edge trajectories, and an absence of canyons. Presence of concave‐up profiles is also suggestive of high sediment influx. In contrast, the eastern QDN margin was primarily regulated by local subsidence and faulting, leading to a stationary shelf‐edge migrating pattern and linear upper‐slope morphology. Densely distributed slope‐confined gullies indicate the margin’s disequilibrium and erosive nature. Further east, the Pearl River Mouth (PRM) margin formed much earlier (ca. 30 Ma) and experienced a more complicated accretion history, including three phases which were dominated by sequential marginal faulting (before ca. 30 Ma), basement structure (ca. 30–23 Ma), and sediment supply (ca. 23 Ma to the present). The overall sigmoidal morphology and truncated stratigraphy of this margin probably resulted from the sculpting of local marine processes, especially ocean currents and internal waves. The exception of the central PRM margin where concave‐up profiles develop is mainly related to canyon erosion. Overall, this study highlights the vital role of local forcing factors in controlling along‐margin variations and determining the final fates of different margin segments. A comparison between the northern SCS and other well‐established examples reveals that concave‐upward shelf‐margin shapes, which are usually associated with high sediment supply, little influence from hydraulic regimes, or sometimes, high degree of canyon development, may be an indicator of good reservoir potential beyond the shelf edge.
-
-
-
Paleogene evolution and demise of the proto‐Paratethys Sea in Central Asia (Tarim and Tajik basins): Role of intensified tectonic activity at ca. 41 Ma
AbstractThe proto‐Paratethys Sea covered a vast area extending from the Mediterranean Tethys to the Tarim Basin in western China during Cretaceous and early Paleogene. Climate modelling and proxy studies suggest that Asian aridification has been governed by westerly moisture modulated by fluctuations of the proto‐Paratethys Sea. Transgressive and regressive episodes of the proto‐Paratethys Sea have been previously recognized but their timing, extent and depositional environments remain poorly constrained. This hampers understanding of their driving mechanisms (tectonic and/or eustatic) and their contribution to Asian aridification. Here, we present a new chronostratigraphic framework based on biostratigraphy and magnetostratigraphy as well as a detailed palaeoenvironmental analysis for the Paleogene proto‐Paratethys Sea incursions in the Tajik and Tarim basins. This enables us to identify the major drivers of marine fluctuations and their potential consequences on Asian aridification. A major regional restriction event, marked by the exceptionally thick (≤ 400 m) shelf evaporites is assigned a Danian‐Selandian age (ca. 63–59 Ma) in the Aertashi Formation. This is followed by the largest recorded proto‐Paratethys Sea incursion with a transgression estimated as early Thanetian (ca. 59–57 Ma) and a regression within the Ypresian (ca. 53–52 Ma), both within the Qimugen Formation. The transgression of the next incursion in the Kalatar and Wulagen formations is now constrained as early Lutetian (ca. 47–46 Ma), whereas its regression in the Bashibulake Formation is constrained as late Lutetian (ca. 41 Ma) and is associated with a drastic increase in both tectonic subsidence and basin infilling. The age of the final and least pronounced sea incursion restricted to the westernmost margin of the Tarim Basin is assigned as Bartonian–Priabonian (ca. 39.7–36.7 Ma). We interpret the long‐term westward retreat of the proto‐Paratethys Sea starting at ca. 41 Ma to be associated with far‐field tectonic effects of the Indo‐Asia collision and Pamir/Tibetan plateau uplift. Short‐term eustatic sea level transgressions are superimposed on this long‐term regression and seem coeval with the transgression events in the other northern Peri‐Tethyan sedimentary provinces for the 1st and 2nd sea incursions. However, the 3rd sea incursion is interpreted as related to tectonism. The transgressive and regressive intervals of the proto‐Paratethys Sea correlate well with the reported humid and arid phases, respectively in the Qaidam and Xining basins, thus demonstrating the role of the proto‐Paratethys Sea as an important moisture source for the Asian interior and its regression as a contributor to Asian aridification.
-
-
-
Fluid flow and pore pressure development throughout the evolution of a trough mouth fan, western Barents Sea
[Conceptual model of TMF evolution during a full glacial/interglacial cycle. (a) Interglacial stage with hemipelagic sedimentation. During winter, dense shelf water flows due to sea ice formation and brine release maintain free of sediment some of the upper and middle slope gullies excavated during the deglaciation. (b) GM, the material transported by ice streams is dumped over the shelf edge as debris flows which can erode the underlying sediments. (c) Deglaciation: the turbid meltwater plumes leave a bed of plumite/turbidite sediments covering the shelf and TMF area, while the most energetic flows excavate gullies on the upper slope. The thickness of this unit increases towards the south. d: Submarine landslides triggered by earthquakes from isostatic rebound induced by ice sheet retreat. (1) hyperpycnal flow; (2) hemipelagic (interglacial) sediments; (3) gullies; (4) contour currents, (5) subglacial (diamicton) till; (6) debris flows; (7) meltwater plumes; (8) gully erosion and plumite/turbidite sedimentation; (9) iceberg rafting; (10) earthquake; (11) landslides; (12) glacial trough. In (b) to (d), sea‐ice is not shown for a better visualisation of the slope processes. Overpressure shading and fluid flow vectors are depicted.
Using a combination of geophysical and geotechnical data from Storfjorden Trough Mouth Fan off southern Svalbard, we investigate the hydrogeology of the continental margin and how this is affected by Quaternary glacial advances and retreats over the continental shelf. The geotechnical results show that plumites, deposited during the deglaciation, have high porosities, permeabilities and compressibilities with respect to glacigenic debris flows and tills. These results together with margin stratigraphic models obtained from seismic reflection data were used as input for numerical finite element models to understand focusing of interstitial fluids on glaciated continental margins. The modelled evolution of the Storfjorden TMF shows that tills formed on the shelf following the onset of glacial sedimentation (ca. 1.5 Ma) acted as aquitards and therefore played a significant role in decreasing the vertical fluid flow towards the sea floor and diverting it towards the slope. The model shows that high overpressure ratios (up to λ ca. 0.6) developed below the shelf edge and on the middle slope. A more detailed model for the last 220 kyrs accounting for ice loading during glacial maxima shows that the formation of these aquitards on the shelf focused fluid flow towards the most permeable plumite sediments on the slope. The less permeable glacigenic debris flows that were deposited during glacial maxima on the slope hinder fluid evacuation from plumites allowing high overpressure ratios (up to λ ca. 0.7) to develop in the shallowest plumite layers. These high overpressures likely persist to the Present and are a critical precondition for submarine slope failure.
]
-
-
-
Salt thickness and composition influence rift structural style, northern North Sea, offshore Norway
Abstract“Salt” giants are typically halite‐dominated, although they invariably contain other evaporite (e.g. anhydrite, bittern salts) and non‐evaporite (e.g. carbonate, clastic) rocks. Rheological differences between these rocks mean they impact or respond to rift‐related, upper crustal deformation in different ways. Our understanding of basin‐scale lithology variations in ancient salt giants, what controls this and how this impacts later rift‐related deformation, is poor, principally due to a lack of subsurface datasets of sufficiently regional extent. Here we use 2D seismic reflection and borehole data from offshore Norway to map compositional variations within the Zechstein Supergroup (ZSG) (Lopingian), relating this to the structural styles developed during Middle Jurassic‐to‐Early Cretaceous rifting. Based on the proportion of halite, we identify and map four intrasalt depositional zones (sensu Clark et al., Journal of the Geological Society, 1998, 155, 663) offshore Norway. We show that, at the basin margins, the ZSG is carbonate‐dominated, whereas towards the basin centre, it becomes increasingly halite‐dominated, a trend observed in the UK sector of the North Sea Basin and in other ancient salt giants. However, we also document abrupt, large magnitude compositional and thickness variations adjacent to large, intra‐basin normal faults; for example, thin, carbonate‐dominated successions occur on fault‐bounded footwall highs, whereas thick, halite‐dominated successions occur only a few kilometres away in adjacent depocentres. It is presently unclear if this variability reflects variations in syn‐depositional relief related to flooding of an underfilled presalt (Early Permian) rift or syn‐depositional (Lopingian) rift‐related faulting. Irrespective of the underlying controls, variations in salt composition and thickness influenced the Middle Jurassic‐to‐Early Cretaceous rift structural style, with diapirism characterising hangingwall basins where autochthonous salt was thick and halite‐rich and salt‐detached normal faulting occurring on the basin margins and on intra‐basin structural highs where the salt was too thin and/or halite‐poor to undergo diapirism. This variability is currently not captured by existing tectono‐stratigraphic models largely based on observations from salt‐free rifts and, we argue, mapping of suprasalt structural styles may provide insights into salt composition and thickness in areas where boreholes are lacking or seismic imaging is poor.
-
-
-
Provenance analysis of detrital monazite, zircon and Cr‐spinel in the northern Tibetan Plateau: Implications for the Paleozoic tectonothermal history of the Altyn Tagh and Qimen Tagh Ranges
Authors Wen Zhu, Chaodong Wu, Jialin Wang and Tianqi Zhou[AbstractSediment provenance studies have proven to be an effective method to extract the sediment provenance and tectonic process information recorded by detrital minerals. In this contribution, we conducted detrital monazite and zircon U‐Pb geochronology and detrital Cr‐spinel major element chemistry analyses on samples from the Qaidam Basin to reconstruct the spatial and temporal evolution of the Altyn Tagh Range and the Qimen Tagh Range in the northern Tibetan Plateau. Based on the significant variation in [Th/U]N, [Gd/Lu]N and [Eu/Eu*]N and the U‐Pb ages of the monazite and zircon, the South Altyn Tagh subduction‐collision belt and the North Qimen Tagh Range were, respectively, the main provenances of the Ganchaigou section and the Dongchaishan‐Weitai section in the Qaidam Basin in the Cenozoic. Paleozoic peak metamorphism, retrograde granulite‐facies metamorphism and amphibolite‐facies metamorphism in the South Altyn Tagh subduction‐collision belt were well recorded by the detrital monazite. In comparison, the detrital zircon is a better indicator of igneous events than detrital monazite. Synthesizing the detrital monazite, zircon and Cr‐spinel data, we concluded that the South Altyn Tagh Ocean and Qimen Tagh Ocean existed in the early Paleozoic and that the Altyn Tagh terrane and Qimen Tagh terrane experienced different Paleozoic tectonothermal histories. The collision between the Qaidam terrane and the Azhong terrane occurred at ca. 500 Ma. The Middle Ordovician was the key period of transformation from the collision‐induced compressional environment to an extensional environment in the area of the South Altyn Tagh Range. In the early Paleozoic, the Qimen Tagh area was characterized by the subduction of oceanic crust.
,The tectonic evolution model of the South Altyn Tagh Range in the Palaeozoic. The Altyn Tagh Range experienced three major tectonic stages in the Palaeozoic. The collision between the Qaidam terrane and the Azhong terrane occurred at ~500 Ma. The Middle Ordovician was the key period of transformation from the collision‐induced compressional environment to an extensional environment. In the late Silurian‐Middle Devonian, the South Altyn Tagh was in a post‐collisional crustal extension setting.
-
-
-
Mid Campanian‐Lower Maastrichtian magnetostratigraphy of the James Ross Basin, Antarctica: Chronostratigraphical implications
AbstractThe James Ross Basin, in the northern Antarctic Peninsula, exposes which is probably the world thickest and most complete Late Cretaceous sedimentary succession of southern high latitudes. Despite its very good exposures and varied and abundant fossil fauna, precise chronological determination of its infill is still lacking. We report results from a magnetostratigraphic study on shelfal sedimentary rocks of the Marambio Group, southeastern James Ross Basin, Antarctica. The succession studied covers a ~1,200 m‐thick stratigraphic interval within the Hamilton Point, Sanctuary Cliffs and Karlsen Cliffs Members of the Snow Hill Island Formation, the Haslum Crag Formation, and the lower López de Bertodano Formation. The basic chronological reference framework is given by ammonite assemblages, which indicate a Late Campanian – Early Maastrichtian age for the studied units. Magnetostratigraphic samples were obtained from five partial sections located on James Ross and Snow Hill islands, the results from which agree partially with this previous biostratigraphical framework. Seven geomagnetic polarity reversals are identified in this work, allowing to identify the Chron C32/C33 boundary in Ammonite Assemblage 8‐1, confirming the Late Campanian age of the Hamilton Point Member. However, the identification of the Chron C32/C31 boundary in Ammonite Assemblage 8‐2 assigns the base of the Sanctuary Cliffs Member to the early Maastrichtian, which differs from the Late Campanian age previously assigned by ammonite biostratigraphy. This magnetostratigraphy spans ~14 Ma of sedimentary succession and together with previous partial magnetostratigraphies on Early‐Mid Campanian and Middle Maastrichtian to Danian columns permits a complete and continuous record of the Late Cretaceous distal deposits of the James Ross Basin. This provides the required chronological resolution to solve the intra‐basin and global correlation problems of the Late Cretaceous in the Southern Hemisphere in general and in the Weddellian province in particular, given by endemism and diachronic extinctions on invertebrate fossils, including ammonites. The new chronostratigraphic scheme allowed us to calculate sediment accumulation rates for almost the entire Late Cretaceous infill of the distal James Ross Basin (the Marambio Group), showing a monotonous accumulation for more than 8 Myr during the upper Campanian and a dramatic increase during the early Maastrichtian, controlled by tectonic and/or eustatic causes.
-
-
-
Discrete element modelling of extensional, growth, fault‐propagation folds
By Stuart HardyAbstractExtensional fault‐propagation folds are now recognised as being an important part of basin structure and development. They have a very distinctive expression, often presenting an upward‐widening monocline, which is subsequently breached by an underlying, propagating fault. Growth strata, if present, are thought to provide a crucial insight into the manner in which such structures grow in space and time. However, interpreting their stratigraphic signal is neither straightforward nor unique. Both analogue and numerical models can provide some insight into fold growth. In particular, the trishear kinematic model has been widely adopted to explain many aspects of the evolution and geometry of such fault‐propagation folds. However, in some cases the materials/rheologies used to represent the cover do not reproduce the key geometric/stratigraphic features of such folds seen in nature. This appears to arise from such studies not addressing adequately the very heterogenous mechanical stratigraphy seen in many sedimentary covers. In particular, flexural slip between beds/layers is often not explicitly modelled but, paradoxically, it appears to be an important deformation mechanism operative in such settings. Here, I present a 2D discrete element model of extensional fault‐propagation folding which explicitly includes flexural slip between predefined sedimentary units or layers in the cover. The model also includes growth strata and shows how they may reflect the various evolutionary stages of fold and fault growth. When flexural slip is included in the modelling scheme, the resultant breached monoclines and their growth strata are strikingly similar to some of those seen in nature. Results are also compared with those obtained using simple, homogeneous, frictional‐cohesive and elastic cover materials. Both un‐lithified and lithified growth strata are considered and clearly show that, rather than just being passive recorders of structural evolution, growth strata can themselves have an important effect on fault‐related fold growth. Implications for the evolution of and strain within, the resultant growth structures are discussed. A final focus of this study is the relationship that trishear might have with the upward‐widening zone of flexural slip activation away from a fault tip singularity.
-
-
-
Strain analysis of a seismically imaged mass‐transport complex, offshore Uruguay
More LessAbstractStrain style, magnitude and distribution within mass‐transport complexes (MTCs) are important for understanding the process evolution of submarine mass flows and for estimating their runout distances. Structural restoration and quantification of strain in gravitationally driven passive margins have been shown to approximately balance between updip extensional and downdip contractional domains; such an exercise has not yet been attempted for MTCs. We here interpret and structurally restore a shallowly buried (c. 1,500 mbsf) and well‐imaged MTC, offshore Uruguay using a high‐resolution (12.5 m vertical and 15 × 12.5 m horizontal resolution) three‐dimensional seismic‐reflection survey. This allows us to characterise and quantify vertical and lateral strain distribution within the deposit. Detailed seismic mapping and attribute analysis shows that the MTC is characterised by a complicated array of kinematic indicators, which vary spatially in style and concentration. Seismic‐attribute extractions reveal several previously undocumented fabrics preserved in the MTC, including internal shearing in the form of sub‐orthogonal shear zones, and fold‐thrust systems within the basal shear zone beneath rafted‐blocks. These features suggest multiple transport directions and phases of flow during emplacement. The MTC is characterised by a broadly tripartite strain distribution, with extensional (e.g. normal faults), translational and contractional (e.g. folds and thrusts) domains, along with a radial frontally emergent zone. We also show how strain is preferentially concentrated around intra‐MTC rafted‐blocks due to their kinematic interactions with the underlying basal shear zone. Overall, and even when volume loss within the frontally emergent zone is included, a strain difference between extension (1.6–1.9 km) and contraction (6.7–7.3 km) is calculated. We attribute this to a combination of distributed, sub‐seismic, ‘cryptic’ strain, likely related to de‐watering, grain‐scale deformation and related changes in bulk sediment volume. This work has implications for assessing MTCs strain distribution and provides a practical approach for evaluating structural interpretations within such deposits.
-
-
-
Peat layer accumulation and post‐burial deformation during the mid‐late Holocene in the Po coastal plain (Northern Italy)
Authors Luigi Bruno, Bruno Campo, Andrea Di Martino, Wan Hong and Alessandro AmorosiAbstractPeat horizons are characteristic features of delta plains worldwide. In this study, we tested the use of peat‐based correlations to assess the deformation of Holocene strata in the Po coastal plain (Northern Italy). The Holocene stratigraphy, about 30 km inland from the modern coastline consists of a peat‐bearing, estuarine and deltaic succession, up to 23 m thick. Through the analysis of 31 core data and 100 piezocone penetration tests, we identified and mapped three 10–40 cm‐thick peat layers (T1–T3) dated to 6.6–5.8, 5.5–5.0 and 3.3–2.7 cal kyr BP respectively. These peat horizons were found to be suitable stratigraphic markers within the Holocene succession over an area of about 200 km2. The mid‐late Holocene palaeogeography, reconstructed through high‐resolution peat correlation, supported by 72 radiocarbon dates, highlights a typical upper delta plain environment, with ribbon‐shaped distributary channels and swamp interdistributary areas. Peat layers are inclined towards E‐NE with gradients that increase downsection from ~0.016% (T3) to 0.021% (T1). The gradient of the oldest peat horizon is one order of magnitude larger than the slope of the modern delta plain (~0.0025%). We infer that peat horizons accumulated during periods of low sediment supply mainly controlled by autogenic processes and were deformed after deposition. Differential compaction of underlying sedimentary strata and recent tectonic activity of the buried Apenninic thrust systems are the most likely drivers of strata deformation. Based on isochore maps, we document that higher sedimentation rates in topographically depressed areas compensated, in part at least, the ongoing deformation, keeping unaltered the topographic gradient and the depositional environment. This study demonstrates that peat‐based correlation and mapping can shed lights on the mechanisms of strata accumulation and deformation in deltaic settings, constituting a robust basis for reconstructing delta evolution.
-
-
-
Spatial distribution of micrometre‐scale porosity and permeability across the damage zone of a reverse‐reactivated normal fault in a tight sandstone: Insights from the Otway Basin, SE Australia
More LessAbstractKnowledge of the permeability structure of fault‐bearing reservoir rocks is fundamental for developing robust hydrocarbon exploration and fluid monitoring strategies. Studies often describe the permeability structure of low porosity host rocks that have experienced simple tectonic histories, while investigations of the influence of faults with multiple‐slip histories on the permeability structure of porous clastic rocks are limited. We present results from an integrated petrophysical, microstructural, and mineralogical investigation of the Eumeralla Formation (a tight volcanogenic sandstone) within the hanging wall of the Castle Cove Fault which strikes 30 km NE–SW in the Otway Basin, southeast Australia. This late Jurassic to Cenozoic‐age basin has experienced multiple phases of extension and compression. Core plugs and thin sections oriented relative to the fault plane were sampled from the hanging wall at distances of up to 225 m from the Castle Cove Fault plane. As the fault plane is approached, connected porosities increase by ca. 10% (17% at 225 m to 24% at 0.5 m) and permeabilities increase by two orders of magnitude (from 0.04 mD at 225 m to 1.26 mD at 0.5 m). Backscattered Scanning Electron Microscope analysis shows that microstructural changes due to faulting have enhanced the micrometre‐scale permeability structure of the Eumeralla Formation. These microstructural changes have been attributed to the formation of microfractures and destruction of original pore‐lining chlorite morphology as a result of fault deformation. Complex deformation, that is, formation of macrofractures, variably oriented microfractures, and a hanging wall anticline, associated with normal faulting and subsequent reverse faulting, has significantly influenced the off‐fault fluid flow properties of the protolith. However, despite enhancement of the host rock permeability structure, the Eumeralla Formation at Castle Cove is still considered a tight sandstone. Our study shows that high‐resolution integrated analyses of the host rock are critical for describing the micrometre‐scale permeability structure of reservoir rocks with high porosities, low permeabilities, and abundant clays that have experienced complex deformation.
-
Volumes & issues
-
Volume 36 (2024)
-
Volume 35 (2023)
-
Volume 34 (2022)
-
Volume 33 (2021)
-
Volume 32 (2020)
-
Volume 31 (2019)
-
Volume 30 (2018)
-
Volume 29 (2017)
-
Volume 28 (2016)
-
Volume 27 (2015)
-
Volume 26 (2014)
-
Volume 25 (2013)
-
Volume 24 (2012)
-
Volume 23 (2011)
-
Volume 22 (2010)
-
Volume 21 (2009)
-
Volume 20 (2008)
-
Volume 19 (2007)
-
Volume 18 (2006)
-
Volume 17 (2005)
-
Volume 16 (2004)
-
Volume 15 (2003)
-
Volume 14 (2002)
-
Volume 13 (2001)
-
Volume 12 (2000)
-
Volume 11 (1999)
-
Volume 10 (1998)
-
Volume 9 (1997)
-
Volume 8 (1996)
-
Volume 7 (1994)
-
Volume 6 (1994)
-
Volume 5 (1993)
-
Volume 4 (1992)
-
Volume 3 (1991)
-
Volume 2 (1989)
-
Volume 1 (1988)