Inventory of Onshore Hydrocarbon Seeps in Romania (HYSED-RO Database)

Seeps are the expression of the migration of hydrocarbons from subsurface accumulations to the surface in sedimentary basins. They may represent an important indication of the presence of petroleum (gas and oil) reservoirs and faults, and are a natural source of greenhouse gas (methane) and atmospheric pollutants (ethane, propane) to the atmosphere. Romania is one of the countries with the largest number of seeps in the world, due to the high petroleum potential and active tectonics. Based on a review of the available literature, and on the field surveys performed by the authors during the last 17 years, we report the first comprehensive GIS-based inventory of 470 seeps in Romania (HYSED-RO), including gas seeps (10.4% of the total), oil seeps (11.7%), mud volcanoes (50.4%), gas-rich springs (12.6%), asphalt (solid) seeps (4.3%), unclassified manifestations (4.0%), and uncertain seeps (6.6%). Seeps are typically located in correspondence with major faults and vertical and fractured stratigraphic contacts associated to petroleum reservoirs (anticlines) in low heat flow areas, and their gas-geochemistry reflects that of the subsurface reservoirs. The largest and most active seeps occur in the Carpathian Foredeep, where they release thermogenic gas, and subordinately in the Transylvanian Basin, where gas is mainly microbial. HYSED-RO may represent a key reference for baseline characterization prior to subsurface petroleum extraction, for environmental studies, and atmospheric greenhouse gas emission estimates in Romania.


Introduction
Subsurface hydrocarbon accumulations and reservoirs in petroleum fields are not always perfectly sealed. Communication pathways may exist in many cases at faults crossing or at the boundaries of cap rocks, facilitating the upward migration of fluids from the deep reservoir to the surface. In petroleum geology, the natural release of hydrocarbon-dominant fluids is named "seepage" [1][2][3]. The overall geological system leading to seepage is termed "petroleum seepage system"-PSS [4], which describes the complex relationship between the elements and processes that favor the migration of hydrocarbons along faults, fractures, and highly permeable volumes of rock. Historically, seeps have provided important clues in the global exploration for oil and gas deposits. Nowadays they still can serve as indicators of gas and/or oil subsurface accumulations, acting as a natural borehole from the surface to the reservoir. However, some seeps may be related to small accumulations, without a commercial value. The assessment of the origin of seeping gas is nevertheless a key task for understanding, without drilling, the subsurface hydrocarbon potential, genesis (e.g., source rock maturity), quality (secondary alteration of gas), and other aspects related to the hydrocarbon accumulations, e.g., the presence distribution and characteristics of the seeps. We will show, in particular, the relationships between seeps, petroleum fields, tectonic structures and heat flow.

History of Seep Exploration in Romania
On the territory of Romania, surface gas manifestations were known well before the discovery of methane (the burning gas from marshes) by Alessandro Volta in 1776. The first gas manifestation was documented in Transylvania in 1671 near Bazna, where an "everlasting fire" was observed by local shepherds [18]. Scholars of the 17th and 18th centuries successively wrote about gas manifestations in Transylvania in local newspapers, poems and scientific reports. The first notice on the gas seep from Bazna is from the diary of Johann Irthell, who visited Bazna on the 11th of April 1672 to see personally the "burning waters". The news spread rapidly in Transylvania, so that after a few years the mentioned phenomenon was described in a poem of Valentin Frank von Frankenstein, while the first scientific description was given by Georgius Vette and Henric Vollgnad in 1688, in the work De aquis ardendibus [19].
The first detailed scientific report was written by Luigi Fernando Marsigli in 1727, who described and mapped the local geology of the Bazna seep, outlining the burning gas [18,19]. Early scientific reports on Transylvanian gas manifestations were written in the mid/late-19th century, mainly regarding the area of Sibiu [20][21][22][23]. In the 20th century many gas manifestations became known in the Transylvanian region, and they were described as mud volcanoes, eternal flames, "bubbling pools" and oils seeps. Posewitz [24] catalogued the oil-, and bitumen manifestations in Transylvania, mentioning also some eternal fires, and gas-bearing springs.
In the Moldavian region, some gas manifestations were documented in the first decade of the 20th century [25,26]. In his PhD thesis, Vancea [27] described in detail the distribution of seeps in the Transylvanian Plain and, as supplementary material, he described a few seeps in Romania, and the importance of connate water in petroleum exploration [28]. Banyai [29] described the morphological features of some mud volcanoes from the Eastern Transylvanian region, indicating the areas where these features are more frequent. In the following decades, seeps were mentioned in geological works [30,31]. The first catalogue of gas and petroleum manifestations in Romania, containing about 1000 seeping sites was produced as an unpublished report by Tonescu, in 1953 [15]. In this catalogue, 73% of seeps were indicated as "oil springs", and the rest of 27% were gas emanations including mud volcanoes. Unfortunately, it seems that the original report is no longer retrievable. Paraschiv [15] recalled however the seeps distribution based on Tonescu's inventory, and provided a synthesis map ( Figure 1). Peaha [32] developed the first catalogue of Romanian mud volcanoes, including a set of chemical analyses of water and gas. Chromatographic analyses of a few Carpathian seeps were then included in a remarkable collection of natural gas data, mostly referring to reservoirs, by Filipescu and Huma [33].
A new approach based on gas flux measurements in seeping areas, combined with molecular and isotopic analyses has developed gradually after the year 2000 [34][35][36]. The first isotopic composition of methane from a Romanian seep was published in 2008 [37]. Geochemical and flux data were then reported for seeps from Transylvania [38][39][40], from the Carpathian Foredeep [41], and from the Moldavian Platform [40,42]. Rowland et al. [43] investigated a large number of springs in the Pannonian region, and measured the dissolved hydrocarbon concentrations, and also the δ 13 C-CH 4 in some samples. A more complete survey of the geochemistry of Romanian seeps was published by Baciu et al. [44]. seeps were indicated as "oil springs", and the rest of 27% were gas emanations including mud volcanoes. Unfortunately, it seems that the original report is no longer retrievable. Paraschiv [15] recalled however the seeps distribution based on Tonescu's inventory, and provided a synthesis map (Figure 1). Peaha [32] developed the first catalogue of Romanian mud volcanoes, including a set of chemical analyses of water and gas. Chromatographic analyses of a few Carpathian seeps were then included in a remarkable collection of natural gas data, mostly referring to reservoirs, by Filipescu and Huma [33].

Development of the HYSED-RO Database
HYSED-RO developed from an accurate and wide literature survey and direct field investigations. It contains only emissions from sedimentary hydrocarbon-prone areas, respecting the proper sense of "petroleum seep" notion. CO 2 -rich manifestations related to geothermal or volcanic areas will be included in a separate database, although they may contain certain amounts of hydrocarbons.
The version presented in this paper (HYSED-RO v. 1.0, 2017) reports 470 seeps, out of which 98 are active, 23 are inactive (extinct) and the activity of the remaining 349 is unknown, as explained below. The database will be progressively updated as new field observations will be acquired. HYSED-RO is a GIS-based database, and as such it can be examined on the basis of geological, geo-structural and seismic maps. This allows the assessment of relationships between seepage and a wide set of geological factors, like tectonic lineaments, petroleum field locations, and heat flow (relationships described in more detail in this paper).
Field work, including identification and description of several seepage sites, was conducted to check the information reported in the literature and also to look for new, undocumented seeps. In some cases, seeps are ephemeral features, and seeps that were reported as active in the past may not show any activity today, or cannot be even identified in the field. Based on the activity of the seeps we categorize them as "active", "inactive" and "with unknown activity". Active seeps are defined as those that were mentioned in the recent literature as being active, or those that periodically are active due to different factors (rainfall, increase of pressure in the subsurface, or seismicity), or their activity has been directly observed in the field. Inactive seeps are those that had no visible activity during the last decades or are mentioned as "fossil" or "extinct seeps" in the literature. The manifestations that were mentioned in the literature only in very generic terms, that we could not identify in the field, and it is not sure whether they are real seeps or not, are classified as "uncertain". When reliable information on the existence of a seep is available, but the type of seep is not clear (gas seep, mud volcano, etc.) it is categorized as "unclassified manifestation".
The database has a simple tabular structure (Table 1), and it is divided in two sections. The first section includes data related to the type and geographic location of the seeps. Each column or field of the data table contains an attribute. Each row or record uniquely corresponds to a seepage site. Each field representing an attribute is treated as one variable, which can be defined as: "Char (n)": a text string of n characters, or "Decimal (n, m)": a real number of n digits (including decimal separator) and m decimal places. For the most representative sites, this information is complemented by gas flux estimates resulted from direct measurements, and geochemical data, including gas and water molecular and isotopic analyses, included in Section 2 of the database. When unpublished, this kind of data can be delivered to interested entities or researchers on request.
Items (1) to (8) are available for all seeps. If one item is missing, the entry in the database is not valid. The ID (identification number) is a unique code given to each individual seep in the database, and it is composed of 11 characters, in the following format: where CC is a two-character ISO 3166-1 code for names of officially recognized countries (https://www.iso.org/iso-3166-country-codes.html). PPPP is the code of the geologic provinces defined by the USGS [45] that cover the territory of the country: 4061-Carpathian-Balkanian Basin; 4048-Pannonian Basin; 4057-Transylvania; 1013-Ukrainian Shield; 4064-West Black Sea Basin; 4047-North Carpathian Basin; 4063-Dobrogean Orogen; 1103-Dobrogea Foreland. TTT is the type of manifestation: gas seep-GAS; mud volcano-MVO; oil seep-OIL; and gas-bearing spring-SPR. ### is an alphanumerical character ranging from 001 to 999.

Data Gathering
Published literature has been used as primary information source for identifying the seeps in the field. Together with the literature data, any other information from newspapers, TV, internet, oral communications, etc. regarding the potential existence of seeps, has been considered. The general geological structure was also an important premise in the assessment of the potential seeping areas. The geographic coordinates of each seep were determined by a GPS handheld device. Coordinates of the closest locality have been considered for the uncertain sites that were not identified in the field. Very often, more vents may occur on a limited area, suggesting they belong to a unique seepage system. In some cases, this is very intuitive, as in the case of Paclele Mari, where dozens of vents occur on a convex plateau, and the whole structure can be defined as one mud volcano, considering there is a unique channel that feed all the vents.

Seep Distribution and Mapping
The majority of the seeps are mud volcanoes (50.4%), followed by gas-bearing springs (12.6%), oil seeps (11.7%), gas seeps (10.4%), and solid seeps (4.3%). The percentage of unclassified manifestations is 4.0%, and 6.6% of the seeps are uncertain ( Table 2). The distribution of seeps is in agreement with the occurrence of petroleum plays, and depends on the tectonic features of the areas of interest. In order to better understand the distribution of the seeps, HYSED-RO was implemented in a GIS interface, using MapInfo Professional V. 11 (Pitney Bowes Software, Troy, NY, USA). The layers used in the evaluation of the seep occurrences are: the geological map of Romania 1:1,000,000 [45], a Digital Elevation Model (DEM) of Romania, the map of the geologic provinces in Romania [46] (Figures 2 and 3), and the geothermal map of Romania [47]. The maps were digitized and geo-referenced in WGS84 Latitude and Longitude projection. The seeps have been categorized according to their type and status. Their geographic coordinates and features have then been stored in a GIS environment, able to manage and analyze all geographical data. In particular, the following base layers are included: tectonics [48,49]; geological provinces of Europe [46]; heat flow [47,[50][51][52].

Seeps vs. Tectonics
A close link has been observed between the occurrence of seeps and the tectonic characteristics. The mud volcanoes, in particular, are sensitive indicators of the stress field in folded areas. The anticlinal crests in hydrocarbon-prone areas meet favorable conditions for the appearance of migration pathways for the fluids. Conduits may take the form of pipes or dikes, reflected at the surface by the shape and distribution of the mud-releasing features [53].
Seeps are usually located along tectonic lineaments, as illustrated by several cases in the Transylvanian Basin and in the Eastern Carpathian Flysch and Foredeep zones (Figures 4 and 5).
In the inner part of the Transylvanian Basin, the main structural elements are the gas-bearing domes and anticlines. Mud volcanoes and dry seeps may occur on these gas-bearing structures, presumably controlled by the presence of local faults (e.g., Deleni, Tauni, Bazna seeps). On the south-western border of the basin, a large number of seeps are located along the sinuous Rusi-Cenade reverse fault system, cutting Pliocene deposits ( Figure 6). Faults associated to the numerous salt diapirs, located along the borders of the basin, provide pathways for the upward migration of gas (e.g., Praid and Homorod seeps-eastern border, Ocna Sibiului-southern border, Ocnisoara-western border).

Seeps vs. Tectonics
A close link has been observed between the occurrence of seeps and the tectonic characteristics. The mud volcanoes, in particular, are sensitive indicators of the stress field in folded areas. The anticlinal crests in hydrocarbon-prone areas meet favorable conditions for the appearance of migration pathways for the fluids. Conduits may take the form of pipes or dikes, reflected at the surface by the shape and distribution of the mud-releasing features [53].
Seeps are usually located along tectonic lineaments, as illustrated by several cases in the Transylvanian Basin and in the Eastern Carpathian Flysch and Foredeep zones (Figures 4 and 5). In the inner part of the Transylvanian Basin, the main structural elements are the gas-bearing domes and anticlines. Mud volcanoes and dry seeps may occur on these gas-bearing structures, presumably controlled by the presence of local faults (e.g., Deleni, Tauni, Bazna seeps). On the south-western border of the basin, a large number of seeps are located along the sinuous Rusi-Cenade reverse fault system, cutting Pliocene deposits ( Figure 6). Faults associated to the numerous salt diapirs, located along the borders of the basin, provide pathways for the upward migration of gas (e.g., Praid and Homorod seeps-eastern border, Ocna Sibiului-southern border, Ocnisoara-western border).

Seeps vs. Tectonics
A close link has been observed between the occurrence of seeps and the tectonic characteristics. The mud volcanoes, in particular, are sensitive indicators of the stress field in folded areas. The anticlinal crests in hydrocarbon-prone areas meet favorable conditions for the appearance of migration pathways for the fluids. Conduits may take the form of pipes or dikes, reflected at the surface by the shape and distribution of the mud-releasing features [53].
Seeps are usually located along tectonic lineaments, as illustrated by several cases in the Transylvanian Basin and in the Eastern Carpathian Flysch and Foredeep zones (Figures 4 and 5). In the inner part of the Transylvanian Basin, the main structural elements are the gas-bearing domes and anticlines. Mud volcanoes and dry seeps may occur on these gas-bearing structures, presumably controlled by the presence of local faults (e.g., Deleni, Tauni, Bazna seeps). On the south-western border of the basin, a large number of seeps are located along the sinuous Rusi-Cenade reverse fault system, cutting Pliocene deposits ( Figure 6). Faults associated to the numerous salt diapirs, located along the borders of the basin, provide pathways for the upward migration of gas (e.g., Praid and Homorod seeps-eastern border, Ocna Sibiului-southern border, Ocnisoara-western border).   In the Eastern Carpathian Flysch and Foredeep, the main tectonic features are represented by N-S trending anticlines and synclines and related thrust and reverse faults. Consequently, the seeps are located on lineaments following the same direction. A relevant example is the Berca-Arbanasi hydrocarbon-bearing structure (Figure 7), where the largest Romanian mud volcanoes are located [35]. Four mud volcanoes (from North to South: Beciu, Paclele Mici, Paclele Mari and Fierbatori) are positioned on the axis of the faulted anticline, generally at the intersection with transversal faults. Paroxistic activity of the mud volcanoes has been observed in relation to the important seismic events in the nearby Vrancea seismic area [55].
The strata in the flysch zone are strongly folded and tilted. Figure  Although the Neogene of the Moldavian Platform (Romanian segment of the East European Platform) shows some potential of gas generation, there are few commercial deposits due to the lack of proper traps and reservoirs [15]. In most of the cases, a few seeps are found on river terraces or Figure 6. Distribution of seeps in the Transylvanian Basin in relation to the gas commercial deposits (red pattern), major faults (black dashed lines), anticlines (red dashed lines), synclines (green dashed lines), after [38,54]. Symbols of the seeps are the same as in Figure 2.
In the Eastern Carpathian Flysch and Foredeep, the main tectonic features are represented by N-S trending anticlines and synclines and related thrust and reverse faults. Consequently, the seeps are located on lineaments following the same direction. A relevant example is the Berca-Arbanasi hydrocarbon-bearing structure (Figure 7), where the largest Romanian mud volcanoes are located [35]. Four mud volcanoes (from North to South: Beciu, Paclele Mici, Paclele Mari and Fierbatori) are positioned on the axis of the faulted anticline, generally at the intersection with transversal faults. Paroxistic activity of the mud volcanoes has been observed in relation to the important seismic events in the nearby Vrancea seismic area [55].
The strata in the flysch zone are strongly folded and tilted. Figure 8 depicts the main features and tectonic style along a WNW-ESE section through the most external flysch units and the western margin of the Foredeep [56], relevant for the gas seeps of Lepsa (about 8 km north of the section) and Andreiasu (about 15 km south of the section). Lepsa seep is located on a north-south oriented overthrust fold in the Paleogene flysch, dominated by a succession of thin layers of sandstone and clay. At Andreiasu, close to the Casin-Bisoca Fault, the strata are in a sub-vertical position and faulted, thus facilitating the gas migration through the Miocene sandstones and volcanic tuffs occurring in the area.
Although the Neogene of the Moldavian Platform (Romanian segment of the East European Platform) shows some potential of gas generation, there are few commercial deposits due to the lack of proper traps and reservoirs [15]. In most of the cases, a few seeps are found on river terraces or close to the river bed and probably related to shallow gas accumulations. The monocline structure of the Neogene sediments is intersected by the erosion along the valleys, thus providing connections to the surface of the shallow gas accumulations [42]. close to the river bed and probably related to shallow gas accumulations. The monocline structure of the Neogene sediments is intersected by the erosion along the valleys, thus providing connections to the surface of the shallow gas accumulations [42].  close to the river bed and probably related to shallow gas accumulations. The monocline structure of the Neogene sediments is intersected by the erosion along the valleys, thus providing connections to the surface of the shallow gas accumulations [42].

Seeps vs. Petroleum Systems
Approximately 30% of the Romanian territory is covered by 18 petroleum systems that range from Palaeozoic to Pliocene in age [14]. We observed that most of the seeps are located in correspondence with productive reservoirs of these petroleum systems [44]. The inventory mapping also shows that the type of seepage reflects the characteristics and source rock maturity of the petroleum systems. Microbial gas seeps (whose geochemical data are described in [44]) dominate the low temperature diagenetic system of the Transylvanian Basin, included in the USGS geologic province denomination 4057 (Figure 3). Thermogenic gas seeps ( [44,57]) are widespread in the catagenetic systems of Carpathian, Moesian and Moldavian geologic units belonging to the provinces 4061 and 1013.
Liquid (oil) and solid (asphalt) seeps are obviously located in correspondence with thermogenic petroleum systems (Eastern Carpathians, Apuseni Mountains and NW of the Transylvanian Basin) where source rock maturity entered the oil window.
Overall, seep type and geochemistry reflect the different geological and maturity conditions of the basins where seeps are located. A vertical sequence for petroleum systems has also been suggested in some basins by seeps displaying different maturity and secondary alterations ( [44]). Accordingly, the seeps can be used as exploration tool to assess the subsurface petroleum potential and quality.

Seep vs. Heat Flow
Heat flow determines the temperatures at depth that, together with the duration of the thermal event, influence the maturity of source rocks and the consequent gas and oil production in sedimentary basins. Too high temperatures (>250 • C, at metagenesis after catagenesis [58]), destroy hydrocarbons. Heat flow data of Romania ( [47] and references therein) have been compared with the location and type of seeps. All the investigated seeps are located in low and medium heat flow areas (Figure 9), consistent with typical thermal gradients of sedimentary petroleum provinces [58,59]. The central part of the Transylvanian Basin is a low heat flux region, the values increasing towards the borders, from 30 to 60-70 mW m −2 , more accentuated on the eastern side of the basin close to the volcanic range of the Eastern Carpathians [47]. The relatively cold and thick pile of sediments of the Transylvanian Basin has produced microbial methane. Higher temperatures, approaching catagenesis and oil window, may justify thermogenic gas seeps at the western border of the Transylvanian Basin [44]. The heat flow further increases towards the volcanic areas of the Eastern Carpathians, exceeding 120 mW m −2 , in the northern part and at the southern border of the volcanic range. Here there are no documented seeps anymore ( Figure 9) and geothermal, CO 2 -rich manifestations prevail [60]. Some of the uncertain seeps described in the old literature, however, occur within the higher thermal flux zones (Figure 10), as the northern and southern segments of the Eastern Carpathians. Such cases deserve specific investigation because at least some of the "presumed" seeps, may be actually manifestations of geothermal CO 2 -rich fluids.
The current heat flow in the Carpathian Flysch and Foredeep areas is around 40-60 mW m −2 , with some exceptions at the limit between the Southern Carpathian Foredeep and the Moesian Platform, where it may drop down to 20-30 mW m −2 . In the Moldavian Platform the heat flow is about 40-50 mW m −2 , with a slight increasing tendency from West to East [47].   Figure 9).

Gas Origin and Output to the Atmosphere
HYSED-RO contains molecular and isotopic data of the natural gas released by 17 major seeps and methane flux data from almost 100 seeps distributed in the several petroleum basins of Romania. Detailed discussion on gas origin and emission to the atmosphere is reported in Baciu et al. [44]. Here we report key concepts that explain the main meaning of the geochemical data of the various seeps reported in the inventory.
Stable C and H isotope composition of methane (coupled with molecular hydrocarbon composition and other geological data) revealed that seeps in the external area of the Carpathians (in the Flysch and Foredeep zones) and on the western margin of the Moldavian Platform, release gas with thermogenic origin. In some cases the gas shows signals of secondary methanogenesis after biodegradation [44]. The gases from the Transylvanian Basin are mainly microbial; however, thermogenic gas occurs at the eastern edge of the basin, close to the Neogene volcanic belt. In the central part of the Transylvanian Basin the gases contain measurable amounts of C2+ alkanes, which may imply that the gas is not totally microbial [44].   Figure 9).

Gas Origin and Output to the Atmosphere
HYSED-RO contains molecular and isotopic data of the natural gas released by 17 major seeps and methane flux data from almost 100 seeps distributed in the several petroleum basins of Romania. Detailed discussion on gas origin and emission to the atmosphere is reported in Baciu et al. [44]. Here we report key concepts that explain the main meaning of the geochemical data of the various seeps reported in the inventory.
Stable C and H isotope composition of methane (coupled with molecular hydrocarbon composition and other geological data) revealed that seeps in the external area of the Carpathians (in the Flysch and Foredeep zones) and on the western margin of the Moldavian Platform, release gas with thermogenic origin. In some cases the gas shows signals of secondary methanogenesis after biodegradation [44]. The gases from the Transylvanian Basin are mainly microbial; however, thermogenic gas occurs at the eastern edge of the basin, close to the Neogene volcanic belt. In the central part of the Transylvanian Basin the gases contain measurable amounts of C2+ alkanes, which may imply that the gas is not totally microbial [44].  Figure 9).

Gas Origin and Output to the Atmosphere
HYSED-RO contains molecular and isotopic data of the natural gas released by 17 major seeps and methane flux data from almost 100 seeps distributed in the several petroleum basins of Romania. Detailed discussion on gas origin and emission to the atmosphere is reported in Baciu et al. [44]. Here we report key concepts that explain the main meaning of the geochemical data of the various seeps reported in the inventory.
Stable C and H isotope composition of methane (coupled with molecular hydrocarbon composition and other geological data) revealed that seeps in the external area of the Carpathians (in the Flysch and Foredeep zones) and on the western margin of the Moldavian Platform, release gas with thermogenic origin. In some cases the gas shows signals of secondary methanogenesis after biodegradation [44]. The gases from the Transylvanian Basin are mainly microbial; however, thermogenic gas occurs at the eastern edge of the basin, close to the Neogene volcanic belt. In the central part of the Transylvanian Basin the gases contain measurable amounts of C 2+ alkanes, which may imply that the gas is not totally microbial [44].
The comparison of reservoir vs. seep geochemistry (CH 4 isotopes and molecular composition) confirms that most seeps are actually linked to reservoirs. For intense, high flux seeps, the gas has also the same hydrocarbon molecular composition of the reservoir (same C 1 /C 2+ ratio) [44]. In general, gas seep geochemistry reflects the different geological and maturity conditions of the basins where the seeps are located, and highlights that in the same basin, different gas sources with different maturity and secondary alterations may exist. This may be due to source rocks or reservoirs at different depths, suggesting a vertical sequence of petroleum systems.
HYSED-RO flux data show that the mud volcanoes on the Berca-Arbanasi hydrocarbon-bearing structure and the Sarmasel and Tauni seeps in Transylvania are the biggest gas emitters, with an estimated output in the order of hundreds of tons of methane per year. Total estimated methane emission from seeps in Romania is around 3000 t y −1 [44].
Available data of the stable C composition of methane can be used in top-down atmospheric models for the estimation of atmospheric methane emissions [61]. The average δ 13 C-CH 4 value for the Carpathian Foredeep is −41.1 (Alimpesti, Andreiasu, Bacau, Lopatari, Raiuti, Paclele mud volcanoes), while for the central-western Transylvanian Basin is −64.6 (Sarmasel, Deleni, Tauni seeps), and for the eastern Transylvanian Basin is −30.9 (Praid, Corund, Homorod seeps). Considering the relative CH 4 emission factors of the several seeps, whereby Carpathian Foredeep seeps are the most intense and large [44], a weighted δ 13 C-CH 4 value of about −44 is estimated.

Concluding Remarks
The present paper describes the first comprehensive GIS-based inventory of hydrocarbon seeps in Romania (HYSED-RO), based on (a) investigations performed in the field by the authors during the last 17 years and (b) information gathered in the literature covering more than three centuries. The database currently comprises a total number of 470 seeps, out of which 98 are active, 23 are inactive, and the activity of 349 is unknown. Most of the seeps are mud volcanoes (50.4%). In the other categories, there are gas-bearing springs (12.6%), oil seeps (11.7%), gas seeps (10.4%), solid seeps (4.3%), unclassified manifestations (4.0%), and uncertain seeps (6.6%).
The occurrence of the seeps is largely dependent on the spatial distribution of the hydrocarbon bearing structures and on the regional and local tectonics (faults and vertical, fractured stratigraphic contacts), which determine communication pathways for the fluids from the reservoirs to the surface. HYSED-RO shows, in particular, that the most important gas emitting structures are distributed over commercial reservoirs in the Carpathian Foredeep and Transylvanian Basin. Geochemical studies of the several seeps [44] showed that the emitted gas reflects the nature and quality of the reservoirs. As examples, the seeps in the Transylvanian Basin reveal the quality of the microbial gas in the reservoirs, extremely rich in methane (>90 vol %), without undesirable gases (e.g., CO 2 , H 2 S), while some seeps from the Eastern Carpathians release gas with indications of biodegradation, which may reflect shallow reservoirs of low quality petroleum. The intensity of seepage is also suggestive for the fluid pressure in the reservoirs, thus indicating the potential for their economic recovery. Basically, seeps are "natural sampling valves" of subsurface hydrocarbon accumulations. In this respect, HYSED-RO can be a valuable tool for petroleum exploration. The seep inventory shows the areas where natural hydrocarbon degassing is occurring, especially for small or ephemeral seeps (e.g., many small mud volcanoes in Transylvania or in the Moldavian Platform), often not known by national authorities for environmental protection and oil industry. Knowing the pre-existence of natural gas seepage-independent of man-made activity-is a critical prerequisite for baseline characterization needed to assess the environmental impacts due to petroleum exploitation, including shale fracking.
Finally, the inventory can also be useful for civil protection purposes, as some gas seeps may represent hazards for the population, infrastructure, and the built environment. The availability of relevant data on the occurrence of seeps may support the establishment of proper protective measures. This inventory also provides first data for estimating the greenhouse gas emission from geological sources at a national level [44].
The interactive map of HYSED-RO inventory is available at http://www.hysedro.blogspot.com. The database files including data points and metadata are available for scientific use upon request to the corresponding author. This way of dissemination should be perceived as an open invitation to researchers from different fields of the natural sciences to contribute with information on the gas-emitting features they may observe during their field work, but the same invitation is extended to hikers, local communities, and the general public. This kind of voluntary contribution may help us to continuously update the database, and to monitor the evolution of some of the seeps. It will also support the better knowledge and protection of these natural features. A similar database is being developed for seeps in Italy (HYSED-IT; [62]), another European country with widespread hydrocarbon seepage.