Organic petrology and geochemistry of Triassic and Jurassic coals of the Tabas Basin, Northeastern/Central Iran

Jurassic and Triassic coals and organic matter-rich shales of the Tabas Basin, Iran, were investigated with respect to their thermal maturation and petrographic composition. The Triassic coals of the Parvadeh coal field range in vitrinite reflectance between 1.17% and 1.37% and the Jurassic coals of the Mazino coal field between 2.08% and 2.29% VRr. Maceral analysis revealed a predominance of vitrinite in all samples, with slightly higher percentages of inertinite in the Jurassic samples. Rock–Eval analyses confirm the presence of type III kerogen in the Triassic coals, while the Jurassic coals contain strongly carbonized residual kerogen. The range of sulfur contents for coals from the Parvadeh coal field (0.37% to 4.64%) and Mazino coal field (0.45% to 2.92%) is related to the effect of marine water in peat. The studied samples are characterized by the predominance of short- over long-chained n-alkanes. The relatively high Pr/Ph ratios indicate predominance of terrestrial organic matter whereas Pr/n-C17 and Ph/n-C18 ratios prove oxic conditions during deposition. In addition, the DBT/Phen ratio shows that these coals formed in fluvial/deltaic environments. Molecular geochemical parameters such as carbon preference index (CPI, 0.99–1.04), methyl phenanthrene index (MPI, 1.20–1.60), methyl naphthalene ratio (MNR, 1.61–3.45), and ethyl naphthalene ratio (ENR, 4–6.78) confirm the high maturity of the samples. Burial and thermal history reconstruction indicates necessity of an erosional thickness of about 4000 m in the Parvadeh area. Towards the Mazino area, a higher basal heat flow up to 80 mW m−2 is assumed for the Paleogene leading to higher maturities.

a Geological map of Tabas block with sampling locations marked, b Geological cross sections along lines AA 0 and BB 0 . Modified after Konon et al. (2016) Organic petrology and geochemistry of Triassic and Jurassic coals of the Tabas Basin… 355 logs, the Mazino coal deposit consists of more than 120 coal layers. In this study, we focus on the three layers exploitable within the Chahrekhneh mine (K 1 , K 2 , and K 3 ). As compared to Permian and Carboniferous coals on the one hand and Cenozoic coals and lignites on the other hand, Mesozoic coals have been much less studied in the past. Furthermore, only little information is available on coals of the Middle East, of which those of Turkey are best studied. Although the mineralogy, inorganic geochemistry, and stratigraphy of the Mesozoic coals of Iran have been investigated by several authors (e.g. Pazand 2015; Goodarzi et al. 2006;Ghassabiyan et al. 2010; Moore and Esmaeili 2012;Ardebili et al. 2012;Moosavirad et al. 2013;Ghavidel-Syooki et al. 2015;Rajabzadeh et al. 2016), only little information is available on the characteristics of organic material that formed these coals (Solaymani and Taghipour 2012; Stasiuk et al. 2006;Alizadeh et al. 2011Alizadeh et al. , 2015Ghavidel-Syooki et al. 2015). In order to investigate depositional conditions, plant input and thermal maturity of the aforementioned Iranian coal deposits, a wide range of organic petrological and organic geochemical methods was applied. This paper will provide first detailed data on the organic matter in the Triassic and Jurassic coals of the Tabas Basin. In addition, a first numerical model on burial and uplift of the Tabas Basin is presented based on the maturity data and geological information.

Geological setting
The Tabas Basin is part of the Central Iranian Micro Continent (Takin 1972;Fürsich et al. 2009;Wilmsen et al. 2009) and located at the central part of the Tabas Block, which is separated by the Nayband Fault to the east and the Kalmard-Kouhbanan Fault to the west from the Yazd and the Lut Blocks, respectively (Alavi 1996;Ghassabiyan et al. 2010, see Fig. 1). The basin is characterized by E-Wtrending faults and fault-related folds in the east and N-Strending faults in the west (Fig. 1). The basin formed as a result of the closure of the Paleo-Tethys Ocean during the Late Permian to the Late Triassic ), when the Iranian plate moved northwards, from an equatorial position to about 25°northern latitude during the Jurassic (Fig. 2, Mattei et al. 2015;Wilmsen et al. 2009;Bordenave and Hegre 2005), always at the northern border of the Neo-Thetys Ocean (Stampfli and Borel 2002). The Tabas Basin holds five major coal deposits, namely the Parvadeh, Nayband, Mazino, Abdooghi, and Kalshaneh coalfields (Ahangaran et al. 2011 rad et al. 2013;Goodarzi et al. 2006;Moore and Esmaeili 2012;Ghassabiyan et al. 2010). This group deposited in a continental lagoonal setting and consists, from old to young, of the Nayband, the Ab-e-Haji, the Badamu, and the Hojedk Formations (Fig. 3). Ghassabiyan et al. (2010) suggested that the Shemshak Group megasequence probably formed in a foreland basin in front of tectonically active blocks. The Shemshak Group within the Tabas Basin overlies the Shotory Formation and is unconformably overlain by the Parvadeh Formation of the Magu Group, which formed in a shallow marine environment (Stöcklin et al. 1965;Ghassabiyan et al. 2010;Fig. 2). The Cretaceous sequence in the Tabas Basin includes shallow marine deposits, e.g. red coarse-grained sandstone, pebbly sandstone to limestone, marlstone, and conglomerate. The Cenozoic sediments in the Tabas Block consist of marl, gypsiferrous marl, and thick-bedded Neogene/Quaternary sediments (Ghassabiyan et al. 2010). In the Parvadeh coalfield, coal layers occur within the Upper Triassic Nayband Formation that has a thickness of about 1800 m and is separated into four members which are, from bottom to top: the Bidestan Member, build up by 120 m of shale, siltstone, and sandstone with sandy limestone and marl; the Howze Sheikh Member, consisting of 240 m of shale, followed by siltstone and finally sandstone; the Howze Khan Member, with 230 m of shale and sandstone; and the Ghadir Member with a thickness of more than 1000 m consisting of sandstone, shale, limestone, and coal layers (Rajabzadeh et al. 2016;Ghavidel-Syooki et al. 2015).
Coal-bearing layers of the Mazino coalfield are hosted by the Middle Jurassic Hojedk Formation, with a thickness of about 60 m, consisting of sandstone, siltstone, shale, and coal seams (Ghassabiyan et al. 2010).

Samples
Fifteen samples of coal and organic matter-rich clay stone were collected from different mines within the Parvadeh and Mazino deposits after removing the weathered surface. The sample set covers several coal seams, which are accessible in the Chahrekhneh mine and several mines in the Parvadeh coalfield (Table 1). All samples were used for organic petrology, elemental analysis, ash content measurement, and Rock-Eval pyrolysis. In addition, molecular organic geochemistry was analyzed on 11 samples from the Parvadeh coalfield.

Petrography
For petrographic analysis, small pieces of rock and coal were embedded in a two-phase epoxy resin (Araldite Ò ) and, after hardening, ground and polished following to the procedures described in detail in Littke et al. (2012). Analyses include maceral-group counting and vitrinite reflectance measurements. Maceral analyses were carried out by counting 500 points per sample on an equidistal grid, following the Stopes-Heerlen system (Taylor et al. 1998), with an incident light microscope (Axio Imager, Zeiss), equipped with an automated stage at 500 9 magnification. Vitrinite reflectance (VRr) was measured on 100 random points per sample with an Axioplan microscope (Zeiss), calibrated with a YAG (Yttrium-Aluminum-Garnet) standard (0.89%). More information on microscopic device and measuring procedure can be found in Sachse et al. (2011).

Bulk and elemental analyses
Prior to analyses, the samples were crushed to fine powders using a rotation mill. The measurement of total inorganic (TIC) and organic carbon (TOC) contents was performed using a LiquiTOC II analyzer (elementar) and the total sulfur content (TS) was measured with a total evaporation analyzer (Leco200). Ash content was measured based on the DIN51719 standard which is in accordance to ISO 1171. The waterfree sample (1 g) was first heated within 60 min to 500°C.
Temperature is then increased until 815°C is reached. This temperature is held for 60 min. Each sample was further measured in triplets.
Rock-Eval pyrolysis was performed using a Rock-Eval VI apparatus (VINCI), using 20 mg pulverized sample of the bituminous coals and 100 mg of the semianthracites, respectively. Details of the analytical and device as well as the calculation of Rock-Eval parameters [S1, S2, S3, T max , Hydrogen Index (HI), Oxygen Index (OI), and Production Index (PI)] are described in Behar et al. (2001).

Molecular geochemistry
An amount of 3-5 g (based on TOC contents) of homogenized sample was extracted with a fully automatic DIO-NEX ASE 150 instrument (Thermo Scientific). Extraction was performed using dichloromethane under high temperature and pressure conditions (100°C, 100 bars). The solvent was then changed to n-hexane and the total amount was reduced to 0.5 ml using rotary evaporation. The samples were subsequently dried with anhydrous sodium sulfate and finally desulfurized. The extracts were separated into two fractions (aliphatic and aromatic) by liquid chromatography on a micro column using 2 g of activated silica gel. For activation, the column was kept in the oven at 200°C for 15 h prior to chromatography. The aliphatic fraction was eluted with five ml of n-pentane. 5 ml of a mixture of n-pentane and dichloromethane (40:60; v:v) were used to separate the aromatic hydrocarbons.
Both fractions were then concentrated to less than 0.2 ml and analyzed by GC-FID. The analyses were performed with a Fisons Instruments GC 8000 series equipped with a split/splitless injector, flame ionization detector (FID), and Zebron ZB-1 HT Inferno fused silica column (30 m 9 0.25 mm i.d., film thickness 0.25 lm, Phenomenex Ò ). Gas chromatography conditions were set to 270°C injector temperature, 1 ll splitless injection (60 s) at 80°C, 3 min isothermal hold, then programmed at 5°C min -1 to 310°C, and held there for 20 min. Hydrogen was used as carrier gas. GC-MS was performed on the aromatic fraction using a quadrupol mass spectrometer Trace MS (Thermoquest) linked to a Mega Series HRGC 5160 gas chromatograph (Carlo Erba, IT), equipped with a 30 m 9 0.25 mm 9 0.25 lm film Zebron ZB-5 fused silica capillary column. The chromatographic conditions were set to 270°C injector temperature, 1 ll split injection at 80°C, 3 min isothermal hold, then programmed at 3°C min -1 to 320°C, 20 min isothermal hold. The mass spectrometer was operated in single ion mode (EI?, 70 eV) with a source temperature of 200°C. Helium was used as carrier gas. The GC-FID and GC-MS results were evaluated using Atlas and X-calibur softwares, respectively.

Basin modeling
To better understand the burial and thermal history and examine different scenarios of erosion magnitude and basal heat flow values in the Tabas Basin, 1D modeling was performed in two different settings, i.e. the Parvadeh and Mazino coal fields. The present-day geometry of the coal fields including the horizon depths and lithologies as well as geological events in terms of deposition and erosion/ hiatus were used as the main input data. The models were reconstructed using the PetroMod software (Schlumberger; Version 2016). Vitrinite reflectance data were used as calibration tool (Senglaub et al. 2006;Hantschel and Kauerauf 2009;Hartkopf-Fröder et al. 2015). The basic concept of 1D basin modeling for outcrop areas has been described in detail in Nöth et al. (2001). To quantify burial and thermal evolution in the Parvadeh and Mazino areas, two pseudo wells have been reconstructed. The Silurian-Devonian succession acts as the lower boundary and the models consist of 12 events, including one erosion event, which took place after Alpine orogeny from 34 Ma to present-day. Due to the lack of sufficient information from deep wells some input data including formation thickness, lithologies and ages had to be estimated based on the general geological information on Central Iran (Aghanabati 1998;Darvishzadeh 1991;Shahrabi 1999). Also to avoid confusion caused by formation names and to keep the models simple, the geological events were named based on the age instead of local names (Table 2). SWIT has been calculated applying the latitude of 33°(Central Asia) (Wygrala 1989). A present-day heat flow of around 48 mW m -2 is assumed based on the crustal thickness variation in Central Iran (Tunini et al. 2015; see also Molinaro et al. 2005   McKenzie 1980). A single rift phase occurred in the late Paleozoic (Bordenave and Hegre 2005). The input data for the Parvadeh area is provided in Table 2 and petrophysical rock properties in Table 3. Calibration of the model results was performed by comparing the measured and calculated reflectance applying the EASY %R o algorithm of Sweeney and Burnham (1990).

Petrography
The microscopic analyses showed that the major maceral group in the samples from both, the Parvadeh and Mazino deposits is vitrinite (V) although the samples from the Mazino deposit show a slight enrichment in inertinite (I).
Liptinite was not detected in the samples. Microscopically counted mineral contents are not reported, because only small coal pieces were studied, which do not represent the entire samples. VR r of Triassic samples in the Parvadeh deposit ranges from 1.17% to 1.37%, corresponding to medium volatile bituminous coal. Reflectance values of the Jurassic coals from the Mazino coalfield are higher and vary between 2.08% and 2.29%, corresponding to semianthracite. Table 4 shows the relative contents of maceral groups in the analyzed Triassic and Jurassic samples. Pyrite and clay minerals are commonly present. The vitrinite of the Parvadeh samples shows small pores in low abundance, while vitrinite of the samples from the Mazino area shows large pores in high concentrations (Fig. 4).

Bulk and elemental analyses
TOC contents vary between 7.55% and 82.42% for the Triassic samples and between 25.85% and 82.00% for the Organic petrology and geochemistry of Triassic and Jurassic coals of the Tabas Basin… 361 Jurassic samples (Table 4), whereas TIC and, thus, carbonate contents are very low. The ash content of the Triassic samples varies between 7.00% and 84.00%, while that of the Jurassic samples ranges between 5.09% and 62.00%. TS contents of the Parvadeh and Mazino samples range from 0.37% to 4.64% and from 0.45% to 2.92%, respectively (Table 4). The sum of ash and TOC ranges from about 85% to 92% for the samples. This range can be expected due to the presence of about 8% to 15% (in total) organic nitrogen, sulfur, oxygen and hydrogen (Ragot 1977;Durand and Monin 1980;see Taylor et al. 1998). HI varies between 82 and 224 mg HC/g TOC for the Triassic and between 6 and 11 mg HC/g TOC for the Jurassic samples. OI varies between 0 and 10 mg CO 2 /g TOC for the Triassic and between 1 and 4 CO 2 /g TOC for the Jurassic samples (Table 4). Low OI values of coals have been commonly observed (Jasper et al. 2009) and do not reflect low O/C ratios; they are related to the release of carbon monoxide rather than carbon dioxide and of oxygen-bearing organic molecules during Rock-Eval pyrolysis (Zieger et al. 2018). T max values range between 461 and 480°C for the Triassic and between 592 to 611°C for the Jurassic samples (Table 4).

Molecular geochemistry
A typical GC-FID trace for the Triassic Parvadeh coals is shown in Fig. 5. For the Parvadeh samples the pristane/phytane ratio (Pr/Ph) varies between 2.51 and 5.46, whereas the Pr/n-C 17 and Ph/n-C 18 ratios vary between 0.11-0.50 and 0.04-0.95, respectively (Table 5). Carbon preference index (CPI) values range from 0.99 to 1.04. Waxiness index (WI) and n-C 17 /n-C 27 ratio, which compare the amount of short-to that of long-chained alkanes, vary between 0.22-0.51 and 4.52-17.06, respectively. The methylphenanthrene index (MPI) varies between 1.20 and 1.60 for the coal and shale samples (Table 5). Methylnaphthalene isomer ratio (MNR) and ethylnaphthalene isomer ratio (ENR) vary between 1.61 and 3.45 and 4.00 and 6.78, respectively. The dibenzothiophene/phenanthrene ratios (DBT/Phen) range from 0.15 to 0.22 (Table 5).  (Table 4). The overall moderate HI and low OI values (Table 3) confirm that the organic matter studied corresponds to kerogen type III. Jasper et al. (2009)  reported lower HI and, at the same time, higher OI values of organic-rich sediments as compared to the associated coals for a set of Carboniferous samples from Germany. The same observation can be made for the analyzed samples (Fig. 6a, b): the highest HI values occur in coal samples and the average OI values are higher in the organic matter-rich sediments and lower in the coals. This might be explained by stronger degradation of the kerogen in the presence of mineral matter and/or mineral matrix effects (Espitalié et al. 1984;Jasper et al. 2009).
Coal-rich sedimentary sequences are ideally suited to study the relationship between optical and geochemical maturity parameters. This is due to the fact that the optical parameters, especially vitrinite reflectance, can be determined with great certainty on coals, whereas there are many pitfalls in other sedimentary sequences, in particular marine sediments, e.g. to the scarcity of vitrinite or the presence of resedimented vitrinite particles (Hartkopf-Fröder et al. 2015). Accordingly, some molecular geochemical parameters, suitable for rather high levels of maturity as reached by the Triassic coals (about 1.2%VR r )  are discussed below. A well-established maturity-indicating ratio based on the distribution of n-alkanes is the CPI, introduced by Bray and Evans (1961). This index is usually high ([ 1) for low mature organic-rich sediments, and its value converges towards 1 for samples of high thermal maturity. CPI ratios of the Parvadeh samples show no predominance of odd-or even-numbered n-alkanes (Table 5) indicating high thermal maturity. There is no difference with respect to this parameter between coals and organic matter-rich shales. As indicated by Fig. 5, all samples from the Parvadeh coalfield show a dominance of short-chain (n-C 11 -n-C 23 ) over long chain n-alkanes ([ n-C 23 ), which can be explained by the high thermal maturity rather than the influence of aquatic organic matter input. MPI-1, MNR, and ENR are maturity indices calculated based on GC-MS results measured on the aromatic fraction extracted from coal and shale samples. The methylphenanthrene index (MPI-1) is a well-established indicator for maturity, in particular for type III kerogen. MPI-1 increases with increasing maturity until VR r values of about 1.35% are reached and then decreases with further thermal maturity (Radke et al. 1982a;Radke and Welte 1983). While values of the organic-rich sediments from the Parvadeh deposit correlate well with those of sediments from the Canada Basin (Fig. 7a), the MPI-1 ratio of the coals more or less constantly increase with VR r (Fig. 7b). The turning point from positive to negative correlations with vitrinite reflectance seems to be at higher thermal maturities for coals as compared to dispersed organic matter of the same kerogen type. Willsch and Radke (1995) reported a stagnation of MPI-1 values of coals around 1.65% VR r , and data published by Radke et al. (1982b) showed that MPI-1 values constantly increased for a sample-set of coals with maturities ranging from 0.63% to 1.70% VR r . From MPI-1 values, vitrinite reflectance (R c ) can be calculated ( Table 5). The calculated reflectance (R c ) shows values similar to the measured reflectance for the coal samples but is slightly lower for the organic matterrich shale samples (Table 5). A positive correlation with VR r is also observed for the MNR values (Fig. 7c); while the values are in the range expected for coals of medium volatile bituminous rank (Radke et al. 1994).

Depositional environment
The high vitrinite contents of the Jurassic and Triassic samples are typical for humic coals and associated organic matter-rich sediments. High Pr/Ph ratios of the Parvadeh samples give further evidence for a predominantly terrigenous plant input. The Pr/Ph ratio has also been used to determine the redox conditions during sedimentation (Didyk et al. 1978), although other factors such as biogenic precursor material and thermal alteration have to be considered as well (Ten Haven et al. 1987). The ratios of Pr/n-C 17 and Ph/n-C 18 point towards a deposition under oxic conditions, typical for coals and organic matter-rich sediments from coal-bearing strata rich in kerogen type III. Figure 8 shows a positive correlation of Pr/n-C 17 vs. Ph/n-C 18 ratios of the Parvadeh samples and samples of lower mature bituminous coals and shales from the Ruhr Basin (Germany). Both ratios are exceptionally low in two of the Parvadeh shales if compared to the coals, probably indicating stronger thermal degradation of pristane and phytane in the shales as compared to the coals rather than more reducing depositional environments. Information on the depositional setting of organic sediments can also be obtained by cross-plotting the DBT/Phen ratio versus the Pr/Ph ratio (Hughes et al. 1995). Because of high Pr/Ph and very low DBT/Phen values, the coal and shale samples from the Parvadeh coal field fall into the fluvial/deltaic zone ( Fig. 9) thus supporting the indications drawn from the parameters discussed above. Aside from maceral composition and molecular geochemical parameters, information on depositional environment can also be deduced from the elemental composition, in particular sulfur-organic carbon relationships (Berner 1984). Sulfurcontaining amino acids and other organic compounds of the peat-forming plants as well as some pyrite/marcasite are the most abundant types of sulfur present in freshwaterinfluenced coals (Casagrande et al. 1980;Chou 2012).
Most of the sulfur in coals is of syngenetic origin, formed during early diagenesis at the peat stage (Casagrande et al. 1980). The total sulfur (TS) content can be used as Fig. 8 Pristane/n-C17 versus Phytane n-C18 plot after Shanmugam (1985) showing the Parvadeh samples in comparison to less mature samples from the Ruhr Basin (Germany) Fig. 9 Organic-rich sediments and coal samples of the Parvadeh coalfield plotted in the DBT/Phen against Pr/Ph diagram of Hughes et al. (1995) showing the depositional setting indicative for marine influence on the coal forming peat, as seawater is enriched in sulfate and interactions of seawater and organic matter leads to the formation of abundant pyrite due to the activities of sulfate reducing bacteria (Casagrande et al. 1977;Given and Miller 1985). The amount of sulfur in coals not affected by marine waters is around 0.5%, and increases to 0.5%-1.5% and [ 1.5% in coals which became subject to slight or strong marine influence, respectively (Sykes et al. 2014;Stock et al. 2016). The interpretive TOC/TS plot developed by Jasper et al. (2010) reveals that the Triassic coals plot in the field of topogenic to transitional mires, but not in the field of ombrogenous mires. Some coals might have been influenced by marine flooding as indicated by quite high TS/ TOC ratios (Fig. 10). In addition, one of the organic matter-rich Parvadeh shales plots along the normal marine line (Berner 1984), also suggesting periodic marine flooding.

Temperature history
Maturity data clearly prove an advanced maturation of coal-bearing rocks in both the Parvadeh and Mazino area. Such high maturation can be due to either a deep burial of the rocks (at low-moderate basal heat flows) or a moderate burial at high heat flows or a combination of both (Senglaub et al. 2006). Considering the present-day geology and heat flow history of the region (Table 2) an erosion thickness of up to 4000 m is necessary to reach the measured maturity data and calibration in Parvadeh area. This thickness of now eroded sedimentary rocks is slightly higher than the reported total thickness of the upper Jurassic-Eocene succession of type sections in Central Iran (Aghanabati 1998;Shahrabi 1999). Erosion is assumed to have occurred during the Alpine Orogeny (Fig. 11a). The coal-bearing Triassic succession was accordingly at maximum burial depth (4500 m) during the Late Eocene/Early Oligocene (Fig. 11a). Following the Alpine Orogeny, uplift and erosion in the region the burial depth decreased to around 500 m at present-day. Due to the absence of major peaks in Jurassic to Eocene heat flow, the upper Triassic succession also experienced the highest temperature (150°C) during the Late Eocene/Early Oligocene with the onset of uplift and erosion. Temperature decreased to about 35°C at present-day. The maximum thermal maturity of 1.2% VR r was also reached during the Late Eocene/Early Oligocene remaining constant up to present-day.
Vitrinite reflectance values are even higher (Table 3) in the Mazino area, and thus maximum paleo-temperatures must have also exceeded those in the Parvadeh area (150°C). In this region, the thermal maturity of the Upper Jurassic succession reached more than 2% VR r . Such maturity variations can be explained by either greater burial depth, higher heat flows, thrust faulting and hydrothermal fluids or a combination of these factors (Ruppert et al. 2010;Hower et al. 2019). Assuming the same burial depth (about 4500 m), such maturities will be reached by applying basal heat flow values of 80 instead of 50 mW m -2 during maximum burial in the Late Eocene (Fig. 11b); or by assuming a much thicker eroded sequence of about 6000 m in this area (at the same heat flows in the Parvadeh area). Although there is no justified evidence, we cannot rule out the possibility of short-term thermal events caused by volcanic intrusions and hydrothermal fluids. Whereas the data for the Parvadeh area seem to fit well with geological data and the paleo-geothermal regime, more information is needed to understand the data for the Mazino region (e.g. apatite fission track data).

Conclusions
Triassic coals and organic matter-rich sedimentary rocks of the Parvadeh coalfield are characterized by abundant vitrinite and less inertinite and a vitrinite reflectance of about 1.2%-1.3% VR r . Ash contents of these medium volatile bituminous coals are usually high, whereas sulfur contents are variable, but mostly high. Jurassic coals and organic matter-rich sedimentary rocks of the Mazino coal field are more mature, with vitrinite reflectance values of about 2.0% to 2.2% VR r . Some of these rocks are rich in inertinite, especially the coals. Similar to the Parvadeh samples, ash contents are mostly high and sulfur contents variable, but often high. Based on the TOC-sulfur relationship, it can be deduced that most coals are derived from topogeneous or transitional mires, but not from ombrogeneous mires. Organic geochemical maturity parameters were tested on these samples and compared to vitrinite reflectance. Differences were observed for coals and other rocks with respect to the evolution of the methyl phenanthrene index at high levels of maturity, i.e. MPI-1 further increases for coals whereas in shales, it reaches the turnaround point at about 1.3% VR r and then decreases. In case of the Triassic coals, MNR values are higher for shales than for coals of the same maturity, whereas Pr/n-C 17 and Ph/n-C 18 ratios are lower. Modeling at the location of Parvadeh coalfield indicates the necessity of about 4000 m erosion which probably occurred during the Alpine Orogeny in order to reach the present-day thermal maturity of the Triassic successions applying basal heat flow of approx. 50 mW m -2 . At the location of Mazino coalfield the present-day thermal maturity which is higher than the Parvadeh area will reach either with higher erosion amounts of about 6000 m or heat flow of up to 80 mW m -2 . Magmatic intrusions can have caused the high heat flow in the area of the Mazino coals.