A Method for Assessing Regional Bioenergy Potentials Based on GIS Data and a Dynamic Yield Simulation Model

The assessment of regional bioenergy potentials from different types of natural land cover is an integral part of simulation tools that aim to assess local renewable energy systems. This work introduces a new workflow, which evaluates regional bioenergy potentials and its impact on water demand based on geographical information system (GIS)-based land use data, satellite maps on local crop types and soil types, and conversion factors from biomass to bioenergy. The actual annual biomass yield of crops is assessed through an automated process considering the factors of local climate, crop type, soil, and irrigation. The crop biomass yields are validated with historic statistical data, with deviation less than 7% in most cases. Additionally, the resulting bioenergy potentials yield between 10.7 and 12.0 GWh/ha compared with 13.3 GWh/ha from other studies. The potential contribution from bioenergy on the energy demand were investigated in the two case studies, representing the agricultural-dominant rural area in North Germany and suburban region in South Germany: Simulation of the future bioenergy potential for 2050 shows only smaller effects from climate change (less than 4%) and irrigation (below 3%), but the potential to cover up to 21% of the transport fuels demand in scenario supporting biodiesel and bioethanol for transportation.


Introduction
Although the metabolism of industrial societies strongly relies on minerals and fossilized biomass, annually harvested biomass from vegetation contributes about 10% to primary energy use in the European Union [1]. Biomass can be derived from different resources, e.g., agricultural land or forest, and transferred into different forms of bioenergy, e.g., biogas, liquid, and solid fuel. In decentralized and renewable energy systems, bioenergy can play important roles in fueling the parts of the transport sector that cannot be easily electrified, or in securing controllable electricity supply that can counteract fluctuations of intermittent wind and photovoltaic power sources [2]. It is therefore important to assess regional biomass potentials and to understand the possible variables that might influence the potential in order to help local governments and planning authorities to make informed choices regarding Table 1. Land use and crop category used in this study. Data adapted from [15,16]. Conflicts between maps from different sources are common, since they were derived with different methods and were based on different primary sources. Generally, the DLM data has a high level of accuracy and reliability in terms of overall land use type, e.g., farming land, vineyard, or built-up area, when compared with satellite data, e.g., from Google Maps. Therefore, the polygons from this source serve as the basic unit when merging the two sources [15]. In the case of a conflict regarding overall land use, DLM data is prioritized. Crop information [16] is then attached to each DLM polygon as an additional attribute. In case multiple crop types from [16] exist on the same DLM polygon, which for agricultural land mostly refers to individual fields, the land use type with the largest area share is assigned to this polygon. Generally, only areas classified as agricultural in the DLM map are attached with additional crop type information from [16].

Crop Type Only Specified in
Plant-soil relationships in the surface soil layer affect crop productivity [17]. For example, yields of Miscanthus range from 5 to 13 t/(ha a) on poor soil or marginal land, while from 7 to 44 t/(ha a) on arable land with higher-quality soils [18]. Therefore, local soil types should be considered to achieve a more accurate biomass or bioenergy yield simulation. For this, a map showing the distribution of typical soil types (soil texture) in the top soils of Germany (resolution 1:100,000) from the Federal Institute for Geosciences and Natural Resources [19] was overlaid with the enriched DLM map. Considering the large heterogeneity in the data and the resulting uncertainty in the precision for a given site, the depiction of the obtained soil texture is presented at the level of the soil types group, according to the German soil classification system (KA5) [19] and as shown in Appendix A. Figure 1 shows the original DLM, crop data, soil data, and the superimposed data at the example of the city of Marbach, Ludwigsburg county, in the south-western state of Baden-Württemberg. It is The accuracy of the thus created CityGML data set was validated by comparing it to statistical data for Ludwigsburg county. The county's total land area, classified into the main forms of land use, was compared with the total land area from the state's 2018 land use report [22]. As Table 2 shows, the total area dedicated to agriculture and forest differed by only 3.2%, and the farming area by less than 7%. Grassland and garden areas show larger differences between, with possible reasons being (i) that the DLM dataset only counted polygons with areas of more than one hectare, with smaller fields not included in the dataset in the first place and (ii) that the DLM dataset contained more categories than the state land use report. For instance, orchard meadows, tree nurseries, and fruit plantations are part of the created GML map, but not the land use report. Since orchard meadow could for example be regarded as grassland with fruit trees, combining both categories yield a sum (8210 ha) that is closer to the grass land area from the land use report (7967 ha). Similarly, adding nurseries and fruit plantations to the garden category reduced this initial difference.  The accuracy of the thus created CityGML data set was validated by comparing it to statistical data for Ludwigsburg county. The county's total land area, classified into the main forms of land use, was compared with the total land area from the state's 2018 land use report [22]. As Table 2 shows, the total area dedicated to agriculture and forest differed by only 3.2%, and the farming area by less than 7%. Grassland and garden areas show larger differences between, with possible reasons being (i) that the DLM dataset only counted polygons with areas of more than one hectare, with smaller fields not included in the dataset in the first place and (ii) that the DLM dataset contained more categories than the state land use report. For instance, orchard meadows, tree nurseries, and fruit plantations are part of the created GML map, but not the land use report. Since orchard meadow could for example be regarded as grassland with fruit trees, combining both categories yield a sum (8210 ha) that is closer to the grass land area from the land use report (7967 ha). Similarly, adding nurseries and fruit plantations to the garden category reduced this initial difference.

Assessment Method for Local Biomass Potential
Biomass potentials can be distinguished between theoretical, technical, economic, exploitable, and sustainable potentials.
The theoretical potential describes the potential that exists in a given region within a certain time period of physically usable energy supply, e.g., the energy stored in the entire plant mass. It is determined by physical limits and marks the upper limit of bioenergy's theoretically realizable contribution to energy supply. Due to insurmountable technical, ecological, economic, and administrative barriers, this potential can generally only be tapped to a limited extent. It therefore has no practical relevance for assessing the actual usability of biomass [23].
The technical potential describes that part of the theoretical potential that can be used after taking into account technical restrictions, e.g., salvage rate, storage losses, and conversion losses. In addition, existing structural and ecological restrictions and other legal requirements and possible social restrictions are taken into account, as they represent barriers to the use of bioenergy similar to technically induced restrictions [23]. In the following, the technical bioenergy potential is the potential after processing (e.g., pelletizing, pyrolysis, and methanization) in the form of secondary energy carriers (biomethane, biodiesel, bioethanol, and burning fuel) or primary energy carriers (energy wood and residues) destined for combustion. In order to calculate the biomass yield considering climate, soil situation, land management, and irrigation patterns, an external crop yield and water demand model, AquaCrop, was integrated into our energy simulation environment (see also the next section). The amount of biomass produced over the course of a growing season on a given area is given as one of the outputs of AquaCrop, which can be considered as the yield of the crop and starting point for assessing bioenergy potentials [24]. By applying conversion factors to theoretical and technical energy yields, the energy potential at these two levels can be calculated and given as the output of the energy simulation model.
It has to be noted that the dynamic yield modeling only applies to field crops, grass, and short rotation coppice, while not to forests and fruit plantations, as the AquaCrop's model does not apply to tree-based biomass creation. Moreover, woody biomass as residue is collected only under the situations of the clearing of fields and of irregular intervals, when trees are felled because of excessive age, diseases, and/or storm damage [23]. The woody biomass from these sources is accumulated through several years. In this paper harvesting of energy crops were simulated at the annual base. Therefore, in this simulation the yield of vineyard, bushes, fruit planation, and orchard meadow have the static raw biomass yield from the literature [23]. The potentials calculated in this paper as well as their calculation methods are summarized in Table 3.

Dynamic Yield Model
AquaCrop is a model that describes the interactions between the plant, atmospheric factors, and the soil. From the root zone, the plant extracts water and nutrients. Calculation scheme of AquaCrop is affected by water stress and temperature stress [24]. The accumulation of biomass and transpiration of crop are influenced by crop type, atmospheric temperature, rainfall, irradiation, CO 2 concentration, irrigation, ground water level, and fertilization. In this study the ground water level and fertilization are not included in the yield model. With the help of the dynamic yield model a more accurate biomass yield can be retrieved for the further simulation. The local climate, soil texture, or even irrigation can result in very different yield values than the typical statistical values.
The aim of the new bioenergy potential workflow is to extract the amount of above-ground biomass indicated, which is the starting point for further bioenergy conversion processes. The biomass Energies 2020, 13, 6488 6 of 24 produced is proportional to the cumulative amount of crop transpiration (ΣTr), with the biomass water productivity (WP) in dry mass as the proportionality factor. Normalization for climatic conditions yields normalized biomass water productivity (WP*), valid for diverse locations, seasons, and CO 2 concentrations [24]. Water in the soil that is easily extracted by the plant is called readily available water (RAW). Hereby, water stress is defined as the percentage of RAW in the root zone, with 0% meaning the soil water is at field capacity, which is ideal for crop growing, and 100% representing the threshold of stomatal closure.

Calculation of Bioenergy Potentials
The calculation of theoretical potential is required for the calculation of technical potentials in a next step. The theoretical energy potential of a certain form of biomass can be described by its heating value, multiplied with its dry mass production yield and its water content percentage, with production yields calculated by taking environmental factors including temperature, precipitation, irradiation, soil texture, and irrigation into account (as discussed in Section 3.3).
As explained in Section 3.2, the technical potential of crops and grass can be derived from the theoretical potential by applying conversion coefficients including conversion losses, and the percentage that are used for energetic purposes. As a first step, the distribution between energy and food is defined for specific crops. Secondly, for the part that is used energetically, in most cases more than one type of secondary (bio-)energy carrier can be obtained per crop, so the respective shares need to be determined. Thus, multiplying bioenergy conversion coefficients with the amount of available biomass allows one to calculate the technical potential of a certain crop in a certain form of bioenergy. Here, the solid bioenergy potential is set equal to the theoretical potential of any crop type. This process is pursued for each land field, each crop type, and each possible form of bioenergy. It has to be noted that the thus calculated technical potential is the potential of the crop used for energetic purposes only. Regardless of whether the crop is used for food purpose or energy purpose, its byproducts, e.g., grain straw or leaf mass, are always considered to have technical energy potential of residue.
Of the harvested forest wood, around 26% are utilized for an energetic purpose in Germany [25], making it an important contributor to local bioenergy potentials, especially in rural and mountainous regions. However, given that the growth of forest biomass is also dependent on other factors than soil, climate, and temperature, notably groundwater levels [26], it is much more difficult to model dynamically. Therefore, a static value for the annual technical potential of forest-based bioenergy (in the form of energy wood) is assumed in the following. Similar to forest, we used literature values of the biomass yield potential of orchard, vineyard, and bushes/hedges. A limited share of the biomass from annual pruning in landscape conservation areas is used as solid residue, since its collection is often complex and expensive.
Finally, all formulas and related parameters for the calculation of biomass and bioenergy potentials for the most widely spread crops in Germany are given in Appendix B.

Simulation Environment and Interface
As mentioned earlier, the assessment of biomass potentials is to be included in an existing modeling environment in order to compare different RES potentials and contrast these with energetic demands in a given region (Figure 2). The modeling environment SimStadt developed at HFT Stuttgart allows one to assess electricity, cooling and heat demands, and renewable energy potentials (photovoltaic) on a single-building level using 3D city models (in the CityGML format [21]). SimStadt provides a modular workflow management for various, primarily energetic, system analysis purposes. Each workflow serves a specific purpose, e.g., heating demand of buildings or photovoltaic potential, while certain modules are shared between workflows, e.g., importing data or data preprocessing [27].

Simulation Environment and Interface
As mentioned earlier, the assessment of biomass potentials is to be included in an existing modeling environment in order to compare different RES potentials and contrast these with energetic demands in a given region (Figure 2). The modeling environment SimStadt developed at HFT Stuttgart allows one to assess electricity, cooling and heat demands, and renewable energy potentials (photovoltaic) on a single-building level using 3D city models (in the CityGML format [21]). SimStadt provides a modular workflow management for various, primarily energetic, system analysis purposes. Each workflow serves a specific purpose, e.g., heating demand of buildings or photovoltaic potential, while certain modules are shared between workflows, e.g., importing data or data preprocessing [27]. For the newly established workflow on regional bioenergy potentials, most of the predefined modules are not applicable due to the fact that the input data is land use polygons instead of building geometries, the exception being the import module that can read CityGML files regardless of the type of objects (building or land use polygon) and the weather model that imports the meteorological data in TMY3 format generated by Meteonorm for the specific region in hourly or monthly resolution. The meteorological data are stored in SimStadt and can be called in later steps.
To model bioenergy potentials more accurately than by using static values for all crops, a new module "YieldGenerator" was developed. Firstly, climate data, including daily precipitation and temperature data from the step "WeatherProcessor" is read, converted, and saved in a form that AquaCrop can process. Secondly, the crop reference evapotranspiration (ETo) is determined from meteorological data, including latitude, longitude, altitude, temperature, relative humidity, wind speed, and cloud shading. The FAO Penman-Monteith method [28] was used for the computation of ETo from meteorological data. These meteorological data were called and calculated according to the Penman-Monteith method on a daily basis over a whole year. Thirdly, irrigation patterns are assigned to crops. In this paper we were interested in the net water demand of a crop to determine the amount of potential irrigation required, on top of precipitation, in a given environment. Net irrigation requirement was calculated by adding a certain amount of water to the soil profile each day in case the local soil water content itself was not sufficient for a given crop type. Fourthly, project files for each crop on all possible soil types in a given region were generated by the workflow, which serve as the main input to the AquaCrop model; all required crop parameters are listed in Appendix C. To limit the complexity of the model, we assumed the soil profile of each soil type was deep (>1.5 For the newly established workflow on regional bioenergy potentials, most of the predefined modules are not applicable due to the fact that the input data is land use polygons instead of building geometries, the exception being the import module that can read CityGML files regardless of the type of objects (building or land use polygon) and the weather model that imports the meteorological data in TMY3 format generated by Meteonorm for the specific region in hourly or monthly resolution. The meteorological data are stored in SimStadt and can be called in later steps.
To model bioenergy potentials more accurately than by using static values for all crops, a new module "YieldGenerator" was developed. Firstly, climate data, including daily precipitation and temperature data from the step "WeatherProcessor" is read, converted, and saved in a form that AquaCrop can process. Secondly, the crop reference evapotranspiration (ETo) is determined from meteorological data, including latitude, longitude, altitude, temperature, relative humidity, wind speed, and cloud shading. The FAO Penman-Monteith method [28] was used for the computation of ETo from meteorological data. These meteorological data were called and calculated according to the Penman-Monteith method on a daily basis over a whole year. Thirdly, irrigation patterns are assigned to crops. In this paper we were interested in the net water demand of a crop to determine the amount of potential irrigation required, on top of precipitation, in a given environment. Net irrigation requirement was calculated by adding a certain amount of water to the soil profile each day in case the local soil water content itself was not sufficient for a given crop type. Fourthly, project files for each crop on all Energies 2020, 13, 6488 8 of 24 possible soil types in a given region were generated by the workflow, which serve as the main input to the AquaCrop model; all required crop parameters are listed in Appendix C. To limit the complexity of the model, we assumed the soil profile of each soil type was deep (>1.5 m) and uniform, regardless of the possibility that soils could have several layers with different characteristics. Fifthly, the AquaCrop Plug-In was called from SimStadt; as output, the annual biomass incremental yield in t/(ha a) and water demand, including potential irrigation demands, in millimeters of all crops on all soil types were stored in an XML configuration file extending other input settings, e.g., the conversion pathways to different bioenergy forms per relevant crop, static default biomass yields, or bioenergy conversion factors from raw biomass.
Another module, "BiomassProcessor" then processes all land use polygons. Users can modify parameters, such as the annual forest wood energetic use rate, the share of energy crops such as corn and rapeseed that are actually used for energetic purposes, or the grass land energy usage rate. The default variables values are shown in Appendix D where data are collected and adapted from [23,25,[29][30][31]. Further input parameters can also be imported from an XML configuration file step. The module analyses each land field polygon, tagged with a certain type of vegetation and soil. Therefore, the module was able to find the corresponding biomass yield of the crop on the soil, the possible bioenergy usages and bioenergy conversion coefficient from the XML configuration file. It then calculated the corresponding technical bioenergy energy potential, with the output being exported to a CSV file.

Approach to Data Validation
The biomass yield results of AquaCrop for various crop and soil types were validated by comparison to statistical data. The biomass yield given by AquaCrop is the above-ground value in dry mass (DM), whereas statistical values are typically given in fresh mass, i.e., including average water contents [23,32] and thus requiring conversion into dry mass values. In SimStadt, yields were simulated for three German counties with different climate patterns, under no specific field management, no irrigation, and no shallow ground water available.
The county of Ludwigsburg is located in the centre of Baden-Württemberg, in Germany's south-western corner. It covers an area of 687 km 2 , with a population of about 550,000 inhabitants. 55% of the county's land area is agricultural land, and 18% forest [33]. The county of Dithmarschen in the state of Schleswig-Holstein is bordering the North Sea. Its population of 133,000 is spread over an area of 1428 km 2 , with 78% of the county's land area being agricultural land, and 3% forest [34]. The county of Ilm-Kreis has the land area of 805 km 2 , among which 43% is agricultural land and 42% is forest [35]. The total population in year 2019 is about 106,000 [35].
Thus, Ludwigsburg represents a suburban and densely populated county in Germany's south, whereas Dithmarschen represents a rural county in Germany's north with low forest cover and high agriculture cover, and Ilm-Kreis represents a neither very rural nor (sub)urban county in Germany's hilly and relatively densely forested center. However, both counties use a sizeable share of their land for agriculture and should thus have meaningful bioenergy potentials Table 4 shows the comparison of yield simulation results with actual yields for selected crops in and the three German counties of Ludwigsburg (south-west), Dithmarschen (north), and Ilm-Kreis (mid-eastern Germany). Only yields on soil types that actually exist in the respective county are compared. As yields vary between location due to a difference in climate, a location-dependent actual yield was introduced to better validate yields in different locations. Actual crop yield data were obtained from the global yield gap atlas (GYGA) [36]. GYGA provides robust estimates of untapped crop production potential on existing farmland, based on current climate data and available soil and water resources globally [37]. In Germany, crop yield data is only available for wheat, maize, and barley. Since barley did not feature in the crop category in Table 1, only yields of rainfed wheat and maize could be compared with modeling results. The actual yield only considered the yield of the harvestable organ, e.g., grain of wheat and maize. Harvest index values, i.e., the ratio between weight of Energies 2020, 13, 6488 9 of 24 harvestable organ and the weight of the whole plant, of 52% and 47% for maize and wheat respectively were retrieved from [38,39]. Furthermore, it has to be noted that GYGA did not differentiate between spring and winter cereals. The climate data used for validation was the average climate between year 2000 and 2010.
As Table 4 shows, all modeling results lay within the yield range given by GYGA. The difference between yields on different soils of the same crop was highest for winter cereals, with 5.2 t DM /(ha a), while for most crops the yield difference between soils was in the order of 1 t DM /(ha a). In Ludwigsburg county, the actual yield of spring cereal was 15.3 t DM /(ha a), compared with 15.5 t DM /(ha a) based on our model, i.e., a deviation of 1.3%. The modeled yield of spring cereals was even closer to the actual yield for Ilm-Kreis with 16.0 and 16.1 t DM /(ha a), respectively, i.e., a difference of 0.1 t DM /(ha a) or 0.6%. Larger deviations occurred for winter cereal in Ludwigsburg and in Ilm-Kreis, and maize in Ilm-Kreis. The modeled yield of maize in Ilm-Kreis had the largest deviation of 7.4 t DM /(ha a), or 35%, compared with statistical yield of 21 t DM /(ha a). Deviation can be explained by the deviation of crop harvest rates and low spatial resolution. No yield data of maize in Dithmarschen was available based on GYGA data, even though [15] shows maize being grown in Dithmarschen. Differences between simulated and actual yields can furthermore be explained by the biomass modeling process: while simulated yields were the accumulated yields over the growing periods, including all the parts above ground and without land management and fertilization, any form of land management to remove, e.g., unwanted parts of plants during the growing period would change the actual yield.
A validation of modeled bioenergy potentials on a subnational and substrate level with actual bioenergy usage levels was inherently limited, since for example information on in/exports of bioenergy as primary or secondary energy carriers into and out of a (in our case) county is typically not available, and technical potentials are rarely fully exploited for a variety of reasons. On a national level in Germany, few studies have assessed aggregated bioenergy potentials in long-term scenarios [30,40,41]. Generally, bioenergy potentials vary depending on natural protection regulation, which impacts the available land area for energy crops, yielding for example values from 543, without consideration of energy crops, to 1425 PJ/a in a "highest probability" scenario, to 1900 PJ/a under an "optimal land use" scenario [41]. As the model presented here considered the total natural land use area, the total natural vegetation cover area of Germany, 298,065 ha, would be applied to derive energy yields on a national level. Given that, the average technical bioenergy yield should be between 5.1 and 17.7 MWh/(ha a) based on above-mentioned potentials [41], with the energy yield of 13.3 MWh/(ha a) as the most realistic considering the constrains. As Table 5 shows, the model in this paper yielded the total bioenergy technical potentials between 10.7 and 12.0 GWh/ha in all three counties, which fit the result of other studies quite well.

Scenarios Setting
For each of the counties of Ludwigsburg and Dithmarschen, one base case and three sensitivity scenarios were defined to quantify the influence of (i) climate change, (ii) priority setting on the forms of bioenergy produced from the available biomass resources, and (iii) use of irrigation on energy yields. Hereby, the aim is not to come up with comprehensive and realistic future scenarios, but rather to test whether the proposed modeling approach reacts reasonably to parameter changes. Table 6 gives an overview of the four scenarios. Here, only changes with regards to the base case are listed, with all other parameters remaining the same as in the base case (Table A5). The values from Table A5 were applied in the base case while Table A6 gives the share of various secondary energy carriers produced from the most important energy crops in the base case scenario (Germany, 2018 data) [29]. Table 7 gives the key parameters of climate at current state (average between 2000 and 2010) and in forecasted case 2050. It has to be noted that beyond the crops listed in Table A6, it is possible to for example produce bioethanol from plants rich in starch, e.g., potato. However almost all the potato for non-food use goes to material production for industrial starch [29]. Furthermore, residue of energy crops is mostly directly used as solid fuel, without any further conversion, while about 40% of maize silage residue is used for biogas production [29,31] (the other maize residue products are assumed to become solid fuel). Table 6. Base case and sensitivity scenarios.

Scenario Name Explanation
Base case Values of Tables A5 and A6 in Appendix D applied [23] Climate 2050 Climate forecast data in 2050 including temperature, precipitation, and CO 2 concentration change. The key parameters of climate situation in both counties are listed in Table 7.

Optimization for fuel consumption
If an energy crop can be a source for biodiesel and bioethanol, all of its yield will be used to this end. If the crop cannot be used for the production of this biofuel carrier, it would follow the same distribution as given in Table A6 Water-energy nexus The impact of different irrigation levels on bioenergy potential. Water stress is set at different levels in percentage to simulate water demand under different irrigation conditions. The irrigation water demand is the minimum amount of water that has to remain in the root zone throughout the growing cycle, and as such the water stress that is allowed in the season.

Ludwigsburg and Dithmarschen Test Cases
Land coverages and populations data for Ludwigsburg and Dithmarschen were already introduced in Section 2.6.
As for Ludwigsburg, it has set itself a goal (in 2013) of reducing greenhouse gas emissions by 90% until 2050 through increased energy efficiency and increasing the share renewable energy sources, e.g., through using biomass for local heat and power generation [43]. The extrapolated 2018 end energy demand (electricity, heating, gasoline, and diesel) in Ludwigsburg is 8506 GWh, with diesel and gasoline accounting for 106 GWh and 58 GWh, respectively [43,44]. The extrapolated end energy Energies 2020, 13, 6488 11 of 24 demand in Dithmarschen in 2018 was 5714 GWh. Among this, end energy demand for transport was 971 GWh, or 11% of the total for diesel and 6% for petrol [45,46]. Table 7. Temperature, precipitation, and CO 2 concentration levels in the three counties, average climate between 2000 and 2010 and 2050. The climate data is generated by Meteonorm [42].

The Impact of Climate Change
Changes in ambient temperature, precipitation, and atmospheric CO 2 concentration (see Table 7), allow one to assess the impact of climate change on crop and by this bioenergy yields. Climate change will increase average temperatures in all three regions, by about 0.5 • C until 2050 compared to the average value between 2000 and 2010, while the change in precipitation patterns was mixed: yearly precipitation levels were expected to increase from 794 to 839 mm in Dithmarschen, and from 570 to 590 mm in Ilm-Kreis, while they were expected to fall from 729 to 716 mm in Ludwigsburg. Figure 3 shows how various crops react to changes in these parameters as indicated in Table 6 for Ludwigsburg and Dithmarschen. The spans of black lines indicate different yields on different soil types, while the points show the median value of crops yield under the average climate between 2000 and 2010 and climate in 2050. These developments would imply changes of ±0.2 t DM /(ha a) by 2050 for most crops compared with current yields, except for maize, with around 1.5 t DM /(ha a) higher yields in Dithmarschen. In fact, only the change in climate until 2050 might make maize cultivation economically attractive in Dithmarschen. Taking these temperature and precipitation changes into account, our model still gives plausible results. In both counties, the yield of rapeseed increased between 11% and 13%. The specific yield of maize in Dithmarschen shows an increase of 35%, as the region moved into a climatic zone favorable to maize cultivation; however, in Ludwigsburg, maize yields decreased by 12%. Opposite to that, grass would produce 21% less biomass in Dithmarschen but 8% more in Ludwigsburg. This can be explained by grass tending to be in favor of moderate cool and humid climate. In both counties, yields of spring cereal, winter cereal, sugar beet, and potato changed by less than ±3% until 2050.
As Appendix E shows, the total bioenergy potential increased by only 0.2%, from 646 to 648 GWh/a, due to climate change in Ludwigsburg. Therein, bioethanol and solid fuel potentials do not change at all, while the biodiesel potential increased by 13%, from 1.5 to 1.7 GWh/a, which could be explained by the fact that the yield of rapeseed, being the main source for biodiesel (see Table A6), was positively affected by a temperature increase from 10.1to 10.8 • C as expected for Ludwigsburg until 2050. In Dithmarschen, the technical energy potential changed by 4%, from 1327 to 1381 GWh/a by 4%. Half of the increase could be explained by the increase in the yield of maize, leading to higher biogas yields. Additionally, higher yields of SRC resulted in higher solid fuel potentials. In both regions, the potentials of energy wood from forest are generally limited as (i) less than 4% of forest biomass can be harvested annually (see Table A5), (ii) only about 25% (see Table A5) of harvested wood biomass was directed to energetic uses, and (iii) forests were not the dominant type of land cover type in both regions, with shares of 33% and 3% in Ludwigsburg and Dithmarschen, respectively. region moved into a climatic zone favorable to maize cultivation; however, in Ludwigsburg, maize yields decreased by 12%. Opposite to that, grass would produce 21% less biomass in Dithmarschen but 8% more in Ludwigsburg. This can be explained by grass tending to be in favor of moderate cool and humid climate. In both counties, yields of spring cereal, winter cereal, sugar beet, and potato changed by less than ±3% until 2050. As Appendix E shows, the total bioenergy potential increased by only 0.2%, from 646 to 648 GWh/a, due to climate change in Ludwigsburg. Therein, bioethanol and solid fuel potentials do not change at all, while the biodiesel potential increased by 13%, from 1.5 to 1.7 GWh/a, which could be explained by the fact that the yield of rapeseed, being the main source for biodiesel (see Table A6), was positively affected by a temperature increase from 10.1to 10.8 °C as expected for Ludwigsburg

Optimizing Biofuel for Tranportation Sector
Biofuels have been required to account for at least 2% of the total transportation fuels used in EU member states since 2005. With the Renewable Energy Directive 2018/2001 (RED II), adopted in December 2018, the EU is continuing the political framework for the use of renewable energy sources in the transport sector for the period from 2021 to 2030 [47], while first-generation bioethanol, i.e., ethanol from agricultural crops, will be phased out until 2030. According to [48], sugarcane ethanol will no longer appear in the mix from 2025, as straw ethanol is considered the most attractive unrestricted non-food and non-feed-based option for gasoline. In Germany, the additional introduction of a biofuel quota in 2007 required oil companies to ensure that 4.4% of diesel sales are made of biodiesel. In addition, they needed to ensure that 1.2% (from 2008, 2%; from 2009, 2.8%; and from 2010, 3.6%) of the sales of motor fuel originate from biofuels [49]. In a scenario that optimizes biofuels for transportation, all energy crops that can be converted into biofuels (biodiesel and bioethanol) were used to 100% for that purpose. In this scenario, 1st generation bioethanol is still considered as a source for transportation in the short term and can be produced from grains, maize, sugar beet, potato, and SRC, while biodiesel is produced from rapeseed. The climate data is unchanged as in the base case with average climate between 2000 and 2010.
As Figure 4 shows, bioethanol potentials increased from 48.1 to 98.5 GWh/a, i.e., by 105%, in Ludwigsburg. Since crops used for bioethanol could alternatively be used for biogas production, biogas potentials decreased by 62%. Similarly, in Dithmarschen the amount of bioethanol potential increased by 100%. However, the technical potential of biodiesel in both counties stayed unchanged.
Differences in total technical bioenergy potentials stem from different conversion efficiencies from biomass to biogas and bioethanol, respectively, as indicated in Table 8: for instance, for cereals, maize and SRC, a conversion to biogas yields 4-18 percentage points higher energy in the secondary energy carrier than in the case of bioethanol. The actual use of maize as an energy crop (see Table A6) also prioritizes for bioethanol production, regardless of the fact that biogas production would be more effective. In case of being used as food feedstock, however, maize silage is used for biogas production [29].
sugar beet, potato, and SRC, while biodiesel is produced from rapeseed. The climate data is unchanged as in the base case with average climate between 2000 and 2010.
As Figure 4 shows, bioethanol potentials increased from 48.1 to 98.5 GWh/a, i.e., by 105%, in Ludwigsburg. Since crops used for bioethanol could alternatively be used for biogas production, biogas potentials decreased by 62%. Similarly, in Dithmarschen the amount of bioethanol potential increased by 100%. However, the technical potential of biodiesel in both counties stayed unchanged. Differences in total technical bioenergy potentials stem from different conversion efficiencies from biomass to biogas and bioethanol, respectively, as indicated in Table 8: for instance, for cereals,  Next to conversion efficiencies, a decisive factor is the demand for certain types of biofuel. In 2018, Ludwigsburg county used 1722 GWh of diesel and 997 GWh of gasoline for transportation [43]. In the base case, if all bioethanol and biodiesel had been used for transportation, biofuels would have covered 3% of the total fuel demand, while biofuels could have covered 6% of the fuel consumption in the scenario prioritizing biodiesel and bioethanol production. In Dithmarschen 106 GWh of diesel and 58 GWh of gasoline were consumed in 2017. Being more rural than Ludwigsburg, bioethanol and biodiesel could have covered 11% of the total fuel demand in the base case and 21% in the scenario prioritizing biodiesel and bioethanol production.

The Impact of Irrigation
In a fourth scenario, external irrigation is applied in addition to natural precipitation. A crop's irrigation demand is determined by the minimal amount of external water that has to remain in the root zone throughout the growing cycle, so that the given water stress is maintained in the growing season. Water stress levels from 10% to 90%, in 10% steps, and from 92% to 98% in 2% steps, were set to simulate water demand under different irrigation conditions. The higher the water stress level, the less water is allowed to stay in the soil. Adequate or more water in the soil (lower actual water stress level than set value) brought by natural rainfall is allowed. In this case no external irrigation is needed. Figure 5 shows the relation between relative technical bioenergy yield and the amount of irrigation. The relative technical bioenergy yield is defined as the total technical bioenergy yield under current water stress level with the total technical bioenergy yield at a water stress level of 0% (adequate water during growing period). As Figure 4 shows, irrigation led to a technical bioenergy potential increased The relative indifference of energy potentials to irrigation can be explained by the relative biomass. Relative biomass is the ratio between the actual biomass amount and the ideal reference case (no water stress, no soil fertility, no soil salinity stress, and no weed infestation). Here, we discarded salinity and weed infestation in the first place, so the reference only presented the ideal yield with enough water available. Table 9 shows relative biomass without irrigation in the base case climate for both Ludwigsburg and Dithmarschen. It can be seen that some crops (spring cereal, maize, sugar beet, and potato) already yield (almost) 100% of biomass without irrigation, while other crops, e.g., grass and SRC, required additional water supply to reach a 100% yield. Here, soil types also influenced biomass yields. For most crops, and regardless of soil types, applying irrigation can only increase biomass yields by less than 7%, thus contributing little to the overall technical bioenergy potential. To be more detailed, from a bioenergy perspective, irrigation should for example only by applied to crops on silty clay soil and SRC on clayish silt in Ludwigsburg.

Discussion
This paper proposed a new workflow that allows one to assess technical bioenergy potentials at the regional level in high resolution, considering local climates and soil conditions, and crop The relative indifference of energy potentials to irrigation can be explained by the relative biomass. Relative biomass is the ratio between the actual biomass amount and the ideal reference case (no water stress, no soil fertility, no soil salinity stress, and no weed infestation). Here, we discarded salinity and weed infestation in the first place, so the reference only presented the ideal yield with enough water available. Table 9 shows relative biomass without irrigation in the base case climate for both Ludwigsburg and Dithmarschen. It can be seen that some crops (spring cereal, maize, sugar beet, and potato) already yield (almost) 100% of biomass without irrigation, while other crops, e.g., grass and SRC, required additional water supply to reach a 100% yield. Here, soil types also influenced biomass yields. For most crops, and regardless of soil types, applying irrigation can only increase biomass yields by less than 7%, thus contributing little to the overall technical bioenergy potential. To be more detailed, from a bioenergy perspective, irrigation should for example only by applied to crops on silty clay soil and SRC on clayish silt in Ludwigsburg.

Discussion
This paper proposed a new workflow that allows one to assess technical bioenergy potentials at the regional level in high resolution, considering local climates and soil conditions, and crop distributions. It builds on commonly available input data such as digital landscape models. In the newly established workflow the advantages of two tools are combined: SimStadt assesses regional energy potentials and demands primarily based on geoinformatic data, and AquaCrop simulates biomass yields and water demands based on soil types, crop types, and climate situations [50]. Integrating AquaCrop enables more accurate bioenergy potential calculations than merely applying average statistical yields that disregard local environmental parameters such as climate and soil.
Yield modeling results were compared with actual yield from GYGA, with deviations of usually less than 1.3%. On the aggregated level, bioenergy potential studies for Germany have shown average biomass technical energy yields of between 5.1 and 17.7 MWh/(ha a), with most realistic yields of 13.3 MWh/(ha a) considering existing constrains. The newly established workflow yields values between 10.7 and 12.0 MWh/(ha a), thus in line with [20,40,41]. For the case studies of Ludwigsburg and Dithmarschen, technical bioenergy potentials were calculated as 465 GWh/a and 1327 GWh/a respectively, which would be 5% and 23% of the total extrapolated 2018 and 2017 end energy demands. As a comparison around 3.9% of electricity is covered by biomass in county Ludwigsburg in 2014 [43]; while 6.6% of the thermal energy demand is covered by biomass in Dithmarschen in 2010 [45]. It has to be noticed that the biomass in the reports [43,45] does not necessarily come from local sources. The results show that at least for the case of Dithmarschen, the more rural county of the two, there might still be a substantial untapped potential in local biomass utilization.
Climate change will lead to a 0.2% increase in the technical bioenergy potential in Ludwigsburg until 2050, and a 4% increase in Dithmarschen according to the result of this paper. For Dithmarschen, the increment of temperature from climate condition favors the maize; in Ludwigsburg a very slight increase of biomass potential is driven by two opposing effects: the temperature increment is favorable however reduced precipitation is detrimental to biomass yields, consistent with results in [51].
By directing bioenergy production to bioethanol and biodiesel for transportation, bioethanol productions could almost double in both regions covering 6% and 21% of the total transportation energy consumption in Ludwigsburg and in Dithmarschen respectively. Under the current application case, solid fuel or biogas could be used in a combined heat and power unit (CHP), while bioethanol and biodiesel play a more important role in transportation until 2030. However, with the introduction of act RED II, bioethanol and biodiesel production from crop is not allowed beyond 2030. While biogas is expected to be the dominate fuel for non-electric transportation until 2050 [52]. Thus, the local energy planner should also consider the policy restrains to make policies to direct biomass use.
With regards to irrigation, our model shows that yield gains were at best modest for the most important relevant crops in the studied regions: in Ludwigsburg 128 m 3 /ha of irrigation water only brought the total technical energy potential increment of less than 1% (See Figure 5). The benefits seemed more significant in Dithmarschen-a 4% total technical energy potential increment with 441 m 3 /ha irrigation (see Figure 5). Different amounts of external irrigation to keep the water stress level between 10% and 98% did not bring significant gain among them. Most crops had close to 100% relative biomass rate in both regions, which means crop yields could not be improved through additional irrigation. This is in accordance to the presented simulation results showing limited benefits from irrigation. Irrigation should thus only be applied selectively, e.g., on SRC on silty clay in Ludwigsburg, where the relative biomass rate could be increased by 52%.
There are limitations to the method in this paper: (i) The availability of reasonably detailed soil and crop maps determined the applicability, even though the method itself is generic and globally applicable. So far, the bioenergy potential in any region only in Germany can be assessed. (ii) For forests and orchards, static biomass yield values from the literature were assigned, rather than the in-depth modeling applied to crops such as cereals or maize. As forest growing cycles span multiple years and yields also depend on groundwater levels, rather than merely on precipitation. A model that assesses forest similar to crops would be much more complex. Therefore, the influence from precipitation and irrigation is only restricted to crop. (iii) The secondary energy carriers that might be produced do not (yet) include hydrogen from biomass. Even though it is now widely recognized that considerable amounts of hydrogen can be produced from renewable sources through biological processes such as bacterial hydrogen fermentation, at the present stage of development of hydrogen fermentation technology has not been applied at the industrial-scale with economic advantage [53]. Nevertheless, the presented method can and should be extended to cover hydrogen potential analysis in order to confront with the biofuel development trend.

Conclusions
The methodology and workflow presented in this paper present a novel, generic method to evaluate bioenergy potentials for any region with good accuracy, by simultaneously considering constraints of water availability, climate, and soil. Furthermore, these bioenergy potentials can compare other RES to potentials and local heating, cooling, or electricity demands within one simulation platform based on a consistent set of input data that is moreover relatively accessible for many regions. The new method does not, however, aim to provide a highly accurate assessment of regional bioenergy potentials based on customized data for a specific region.
In the context of the water-energy nexus, our research extends an energy-centric modeling environment to water-related aspects, namely crop transpiration and potential irrigation benefits or requirements. A combination with the urban residential and non-residential water demand simulation workflow presented in [54], which also bases on similar input data, allows one to study the nexus between bioenergy, crop irrigation, and urban water consumption in the regional context. With the newly established biomass, roof PV, electricity, and heating demand workflow [27] in place, the modeling environment SimStadt can now assess regional renewable energy balances with the constraint of water and climate in regions without rich wind, ground-based PV, or hydro potentials. As the next steps, (i) a workflow similar to the presented bioenergy one will be developed for wind power. (ii) The nexus will be extended to food-energy nexus, by, e.g., (a) assessing local balances between food-use and energy-use of crops to fulfill local demands and to meet strategic goals of emission reduction by regionalization of food supply [55]. (b) Investigating the impact of diet change, e.g., more plant-based diets, on local food and energy potentials will be carried out.
P theoretical is the theoretical energy potential of a specific land field in GJ/ha·a. i is the crop type. A vegetation is the area in ha of the vegetation type. Y vegetation is the dry matter production yield of a specific vegetation type in t/(ha·a). H u is the calorific value in gigajoules per tonne [GJ/t lutro ].
where P biogas, GL OM is the annual biogas fuel potential [GJ/a]. i is the crop type. A i is the area of the crop that can be as the source of biogas production in hectares (ha), E crop is the crop harvest yield of fresh mass of each crops [kg/ha·a], TS and oTS are its dry mass rate and organic dry mass rate, and E CH4,i is its methane yield [l/kg oTS]. n energy,i is the share of actual energetic use per crop, and H u is the calorific value of methane in gigajoules per tonne [GJ/l CH 4 ].
where P VO is the fuel potential from rapeseed cultivation for vegetable oil production [GJ/a]. A biodiesel is the area of the area of biodiesel production crops (rapeseed). P Primary is the theoretical yield of the crop in GJ/(ha a), C biodiesel is the conversion efficiency from dry raw mass to biodiesel.
where P ethanol is the fuel potential from energy crops for ethanol production [GJ/a], A ethanol is the area of the area of bioethanol production crops (sugar beet, grain maize and cereals), and C energy is the ethanol conversion factor from dry raw mass. P Primary is the theoretical yield of the crop in GJ/(ha a).
where P residue is the fuel potential of residues as solid fuels in gigajoules per year in GJ/a, A is the area of individual object types in ha, Y is the dry matter production yield of a specific vegetation type in t/(ha·a), C residue is the residue yield in wet mass per total yield in wet mass, TS is the dry mass rate, n energy is the percentage of energy usage of residue, and H u is the heat value in gigajoule per kilogram.

Parameter Winter Cereal Spring Cereal
Soil water depletion factor for canopy expansion, upper limit 0 0.2 Soil water depletion factor for canopy expansion, lower limit 0.35 0.55 Shape factor for water stress coefficient for canopy expansion 2.5 3.5 Soil water depletion factor for pollination (p-pol), upper threshold 0.9 0.9 Shape factor for water stress coefficient for stomatal closure 2 5 Shape factor for water stress coefficient for canopy senescence 2 3 Appendix D Table A5. Default parameter values used in modelling local bioenergy potentials.

Parameter Default Value Explanation
Conifer trees harvest rate 4.5% [23] The percentage in volume of conifer trees harvested annually out of all conifer trees Deciduous trees harvest rate 3.0% [23] The percentage in volume of deciduous trees harvested annually out of all deciduous trees Forest energy usage rate 25.6% [25] The percentage in volume of solid forest wood with diameters > 7 cm that is used for energy purposes Energy crop rate 14.0% [29] The percentage of farmland area used for energy crop cultivation (e.g., rapeseed, maize). Energy crops are used exclusively for energetic purposes. Since no data source gives information on the end product of a crop (energy or food) per field, we assume, in line with statistical data, that 14% of each field's area is used for energetic purposes. Residue energy usage rate 62.0% [30] The percentage of residue by-products which are used for energetic purposes. Rate of maize residue for Biogas production 39.4% [29,31] The percentage of maize residue (silage) for biogas production. The rest of maize residue of maize is used as solid fuel.

Appendix E
Deciduous trees harvest rate 3.0% [23] The percentage in volume of deciduous trees harvested annually out of all deciduous trees Forest energy usage rate 25.6% [25] The percentage in volume of solid forest wood with diameters > 7 cm that is used for energy purposes Energy crop rate 14.0% [29] The percentage of farmland area used for energy crop cultivation (e.g., rapeseed, maize). Energy crops are used exclusively for energetic purposes.
Since no data source gives information on the end product of a crop (energy or food) per field, we assume, in line with statistical data, that 14% of each field's area is used for energetic purposes. Residue energy usage rate 62.0% [30] The percentage of residue by-products which are used for energetic purposes.
Rate of maize residue for Biogas production 39.4% [29,31] The percentage of maize residue (silage) for biogas production. The rest of maize residue of maize is used as solid fuel.