Enhancing methane production from the invasive macroalga Rugulopteryx okamurae through anaerobic co-digestion with olive mill solid waste: process performance and kinetic analysis

Vol.:(0123456789) 1 3 https://doi.org/10.1007/s10811-021-02548-3 Enhancing methane production fromtheinvasive macroalga Rugulopteryx okamurae throughanaerobic co‑digestion witholive mill solid waste: process performance andkinetic analysis D.delaLama‑Calvente1· M.J.Fernández‑Rodríguez1,2· J.Llanos1· J.M.Mancilla‑Leytón3· R.Borja1 Received: 28 March 2021 / Revised and accepted: 18 June 2021 © The Author(s) 2021 Abstract The biomass valorisation of the invasive brown alga Rugulopteryx okamurae (Dictyotales, Phaeophyceae) is key to curbing the expansion of this invasive macroalga which is generating tonnes of biomass on southern Spain beaches. As a feasible alternative for the biomass management, anaerobic co-digestion is proposed in this study. Although the anaerobic digestion of macroalgae barely produced 177mL of CH4 g−1 VS, the co-digestion with a C-rich substrate, such as the olive mill solid waste (OMSW, the main waste derived from the two-phase olive oil manufacturing process), improved the anaerobic digestion process. The mixture improved not only the methane yield, but also its biodegradability. The highest biodegradability was found in the mixture 1 R. okamurae—1 OMSW, which improved the biodegradability of the macroalgae by 12.9% and 38.1% for the OMSW. The highest methane yield was observed for the mixture 1 R. okamurae—3 OMSW, improving the methane production of macroalgae alone by 157% and the OMSW methane production by 8.6%. Two mathematical models were used to fit the experimental data of methane production time with the aim of assessing the processes and obtaining the kinetic constants of the anaerobic co-digestion of different combination of R. okamurae and OMSW and both substrates independently. First-order kinetic and the transference function models allowed for appropriately fitting the experimental results of methane production with digestion time. The specific rate constant, k (first-order model) for the mixture 1 R. okamurae- 1.5 OMSW, was 5.1 and 1.3 times higher than that obtained for the mono-digestion of single OMSW and the macroalga, respectively. In the same way, the transference function model revealed that the maximum methane production rate (Rmax) was also found for the mixture 1 R. okamurae—1.5 OMSW (30.4mL CH4 g−1 VS day−1), which was 1.6 and 2.2 times higher than the corresponding to the mono-digestions of the single OMSW and sole R. okamurae (18.9 and 13.6mL CH4 g−1 VS day−1), respectively. Keywords Invasive macroalgae· Ruguloterix okamurae· Dictyotales· Anaerobic co-digestion· Methane yield· Process kinetics· Transference function· First-order model Introduction The greatest biodiversity loss in the world is produced by the introduction of invasive alien species (IAS) (Convention on Biological Diversity 2017), with the seas the areas that suffer the most its consequences. A clear recent example of IAS in Mediterranean aquatic ecosystems is the dictotalean macroalga Rugulopteryx okamurae, a north-western Pacific native species (Huang 1994). Rugulopteryx okamurae specifically comes from temperate waters near China, Korea, the Philippines, and Japan. This macroalga has been recently included in the Spanish catalogue of invasive alien species according to the ministerial order TED/1126/2020 of the * M. J. Fernández-Rodríguez mjf[email protected] 1 Instituto de La Grasa (CSIC), Campus Universidad Pablo de Olavide, Edificio 46. Ctra. de Utrera Km. 1, 41013Sevilla, Spain 2 Departamento de Sistemas Físicos, Químicos y Naturales, Universidad Pablo de Olavide, Ctra. de Utrera, km 1, 41013Sevilla, Spain 3 Departmento de Biología Vegetal y Ecología, Facultad de Biología, Universidad de Sevilla, 41080Sevilla, Spain / Published online: 25 July 2021 Journal of Applied Phycology (2021) 33:4113–4124 1 3 Ministry of Ecological Transition published in the state official bulletin (BOE, number 314 of 1 December 2020). Rugulopteryx okamurae has proliferated on the Atlantic and Mediterranean coast in a short time displacing the local biota and generating an unprecedented ecological impact, causing a great concern for the sea-dependent local activities in the aforementioned geographical enclave (El Aamri etal. 2018; Navarro-Barranco etal 2019; García-Gómez etal. 2020; Sempere-Valverde etal. 2021). The high content of spatane type diterpenes, with herbivorous deterrent activity, of this species has allowed its high and rapid propagation and production. In fact, only in 1year after being observed for first time in the Straitof Gibraltar area in the autumn of 2015 in Ceuta, the beach cleaning machines of this city removed more than 5000 t of upstream macroalga biomass of the IAS (Ocaña etal. 2016). It took only 4years to occupy the entire Cadiz coastline and to expand along the coasts of Malaga and Granada, affecting very relevant spaces of great ecological value included in the Natura 2000 Network. This rapid expansion is unprecedented; the most similar case may be the spread of the seaweed Sargassum spp. in the Caribbean (van Tussenbroek etal. 2017). Rugulopteryx okamurae not only does have negative effects on the biodiversity of the affected areas, but it also has a great impact on human activities and economic affairs (García-Gómez etal. 2018). Economic impact associated to fishing activities and beach management in Southern Spain has been estimated to be higher than one million Euros per year (Altamirano Jeschke etal. 2016). Currently, tonnes of collected biomass of R. okamurae end up in landfills. However, this macroalga is known to contain polysaccharides (alginates and fucoidans) which are used in the food industry, nutraceutical and pharmaceutical industry, medical treatments, etc. (Hernández-Carmona etal. 2012). It is rich in secondary metabolites among which the terpene and polyphenolic compounds stand out presenting photoprotective properties, antioxidants, and immunestimulants (De Paula etal. 2011). Interest in the study of these natural products has been growing in recent years due to their ecological role and their utility for human beings (Hernández-Carmona etal. 2012). Different groups, associations, and local administrations affected by this invasion are proposing for this macroalgae to not be treated as a waste but rather as a valued raw material for cosmetic, dermatological, or composting use, for example. It is about finding sustainable uses and exploitation of the large available biomass that will finance the removal of arrivals from the coast: the so-called elimination for its valorisation. In addition to the uses or exploitations already mentioned, another alternative could be as a substrate in anaerobic digestion (AD). AD is an anoxic biological process with the participation of different microorganisms that transform organic matter until obtaining biogas as a final product (mainly methane, 60–70%) with a high calorific value, liable to be used energetically by combustion in engines, turbines, or boilers, either individually or mixed with other fuels (Abad etal. 2019). As a circular economy model, AD has been widely used for the degradation and stabilization of domestic and industrial wastes (Parra Huertas 2015). Recently, this biological process has been carried out in different species of red, green, and brown macroalgae obtaining methane yield coefficients of between 0.10 and 0.30 m3 CH4 kg−1 volatile solid (Morand and Briand 1996; Barbot etal. 2016; Saratale etal. 2018). Therefore, it could be a good alternative that would help eliminate and enhance the thousands of tonnes of biomass of R. okamurae deposited on the beaches, decreasing, in turn, the spread of this macroalga due to the removal of vegetative propagules. The co-digestion of different residues for improvement in the yields obtained from individual anaerobic digestion processes is a frequent practice (Saratale etal. 2018). In this way, the combination of different organic substrates generating a homogeneous mixture increases the performance and methane yield of the process. Moreover, anaerobic digestion of single macroalgae has been studied, although they are substrates with a low C/N ratio, which largely limits the proper functioning of the anaerobic process, generating low methane yields (Tabassum etal. 2016; 2017). Given these results, some authors have recently proposed co-digestion processes of macroalgae with nitrogen-poor but carbon-rich substrates to improve and compensate for the C/N ratio of anaerobic processes (Grosser 2017; Saratale etal. 2018). In this sense, R. okamurae, as nitrogen-rich biomass (C/N around 10:1–15:1, Geider and La Roche 2002), could potentially be used in AD as a co-substrate for C-rich by-products. One of these by-products could be the olive mill solid waste (OMSW), the main solid by-product generated in the process of producing olive oil by centrifugation in two phases, which is a mixture of olive pulp, skin, and pit (Fernández- Rodríguez etal. 2020). This by-product has a high organic content and its treatment and management represents an important environmental problem. Spain, as the world’s leading producer of olive oil, generates 2–4 million tonnes of OMSW per year (Moreno-Maroto etal. 2019). The C/N ratio of the OMSW is above the value indicated by Paul and Dutta (2018) as optimal for the process (C/N of OMSW 32, C/N optimal 25–30). Furthermore, the lignocellulosic content of the OMSW, along with the presence of certain AD inhibitor compounds, such as polyphenols, limits the AD process of this compound alone (Varnero etal. 2014). The use of macroalgae as a nitrogen-rich co-substrate would help balance the C/N ratio towards the values indicated as optimal and would help to dilute the concentration of toxic substances that the OMSW possesses (Ferreira etal. 2018; Li etal. 2018), thus, obtaining an improvement in the kinetics of the process, higher methane performance coefficients, 4114 Journal of Applied Phycology (2021) 33:4113–4124 1 3 and an increase in the biodegradability of the substrates (Fernández-Rodríguez etal. 2019a, b). Therefore, the aim of this study was to evaluate the methane potential and the stability of the anaerobic co-digestion of OMSW and the invasive macroalga R. okamurae in order to propose a possible alternative for waste management that contributes to avoid the spread and new colonization of this invasive alien species as well as to assess the viability of a higher methane production from this alga. The influence of the C/N ratio on the kinetics of the anaerobic processes and ultimate methane yields of both substrates individually and different mixtures of them was also evaluated in biochemical methane potential (BMP) tests. Materials andmethods Influents The experimental olive oil mill factory located in the Instituto de la Grasa was selected for the collection of the two-phase OMSW. The OMSW was collected and frozen to − 4°C until used. Then, OMSW was sieved (2mm) in order to remove any piece of olive pit trapped in the pulp prior use. The invasive macroalgaeR. okamuraewas provided by the Laboratory of Marine Biology of the University of Seville. The inoculum used during the anaerobic study was collected from an industrial up-flow anaerobic sludge blanket reactor treating brewery wastewater located at Heineken factory (Seville, southern Spain). The inoculum was chosen due to its high methanogenic activity, which was determined in previous anaerobic assays (Fernández-Rodríguez etal. 2021b). Both inoculum and substrates were characterized and the results are summarized in Table1. Experimental procedure The mesophilic batch experiments (BMP tests) were performed with an inoculum to substrate ratio fixed at 2 (as volatile solid (VS)) as described elsewhere (Fernández-Rodríguez etal. 2019a). A 10% micronutrient solution (Fernández-Rodríguez etal. 2019b) and nitrogen gas was added at the beginning of the experiment in order to keep anaerobic conditions. The anaerobic reactors of 250mL total volume were placed in a thermostatic water bath at a controlled temperature of 35 ± 2°C. Three replicates for each setting were carried out. Data are presented as means ± standard deviations of the means (n = 3). The agitation rate was constant by magnetic bars at 380rpm. For the purpose of the experiment, the produced biogas was passed through a 3N NaOH solution in order to trap the CO2, while the remaining gas was assumed to be methane. The anaerobic digestion experiments were stopped when the gas production remained essentially unchanged; this period was c.a. 28days across the board. Five different VS ratios of the substrates R. okamurae-OMSW were evaluated in triplicate. The different co-digestion mixtures studied were as follows: 1 R. okamurae-0 OMSW; 1 R. okamurae-3 OMSW; 1 R. okamurae-1.5 OMSW; 1 R. okamurae-1 OMSW; and 0 R. okamurae-1 OMSW. As controls, a triplicate with only inoculum and trace elements solution but without the addition of substrate was performed, in order to evaluate the endogenous methane production. Analytical methods The substrates and inoculum were analyzed before the beginning of the experiments. The digestates or effluents of each set of experiments were also analyzed at the end of each experiment. The performed analyses were as follows: total solids (TS) and volatile solids (VS) were performed according to the standard method 2540E (APHA 2012), soluble chemical oxygen demand (SCOD) was also determined following the standard method 5220D (APHA 2012), using the closed digestion and the colorimetric method. Total chemical oxygen demand (COD) was analyzed as described by Raposo etal. (2008). pH and total alkalinity (TA) were determined in fresh samples using a pH metermodel Crison 20 Basic. TA was measured by pH titration to 4.3 (APHA 2012). Total ammonia nitrogen (TAN) was measured by Table 1 Physicochemical properties of the anaerobic sludge, olive mill solid waste, and the invasive alien macroalgae (Rugulopteryx okamurae). Values represent mean ± standard error TS: total solid; VS: volatile solid; TCOD: total chemical oxygen demand Parameter Anaerobic sludge Olive mill solid waste Rugulopteryx okamurae TS (g kg−1) VS (g kg−1) VS/VT ratio TCOD (g O2 kg−1) pH C (%) 30.5 ± 0.7 22.7 ± 1.8 0.75 32.8 ± 1.1 7.24 ± 0.01 nd 245.2 ± 1.7 197.6 ± 2.5 0.81 324.1 ± 0.7 4.95 ± 0.2 49.70 ± 0.38 197.6 ± 7.3 102.1 ± 4.1 0.52 204.4 ± 1.8 nd 30.29 ± 0.50 N (%) C/N ratio nd nd 1.58 ± 0.06 31.4 ± 0.8 1.99 ± 0.03 15.2 ± 0.4 4115Journal of Applied Phycology (2021) 33:4113–4124 1 3 distillation and titration according to the standard method 4500-NH3 (APHA 2012). C and N were determined through an LECO CHNS-932 Elemental Analyzer (Leco Corporation, USA). Soluble parameters were determined after sample centrifugation (Eppendorf, 9000xg, 10min) and filtration (47mm glass fiber filter). Kinetic models First‑order kinetic model In order to study the process kinetics and estimate the process performance in the anaerobic digestion and co-digestion of the two substrates studied, the following first-order kinetic model was used: where G is the cumulative specific methane production (mL CH4 g−1 VSadded), Gm is the ultimate methane production (mL CH4 g−1 VSadded), k is the specific rate constant (day−1), and t is the digestion time (day). This kinetic model is generally applied to evaluate the kinetics of the batch anaerobic digestion processes for different types of biodegradable substrates (Li etal. 2012; Scarcelli etal. 2020). This model is based on the premise that methane generation is proportional to the amount of substrate and not limited by microbial cell mass (Wang etal. 2017). Transference function model The transference function (TF) model was also used to fit the experimental data of methane production during biochemical methane production (BMP) tests (Eq.(2)). The transference function (Reaction curve-type model) (RC), utilized mainly for control purposes, contemplates that any process might be analyzed as a system receiving inputs and producing outputs (Donoso-Bravo etal. 2010). The TF model has been successfully applied by several authors for the biomethanization of different organic wastes (Donoso-Bravo etal. 2010; Li etal. 2012; Pinto-Ibieta etal. 2016). The TF model is given by the following expression: where B (mL CH4 g−1 VSadded) is the cumulative specific methane production, Bmax (mL CH4 g−1 VSadded) is the ultimate methane production, Rmax is the maximum methane production rate (mL CH4 (g−1 VSadded·day−1)), t (day) is the digestion time, and γ (day) is the lag time. Determination coefficients (R2) and standard errors of estimates (SEE) were calculated to evaluate the goodness- of-fit and the accuracy of the results for both models. The (1) G −Gm [ 1−exp (−k⋅t) ] (2) B =Bmax ∗ ( 1−exp [ − R max (t− 𝛾 ) Bmax ]) kinetic parameters for each experiment and mathematical adjustment were determined numerically from the experimental data obtained by non-linear regression using the software Sigma-Plot (version 11). Results Geographical distribution andmain characteristics ofthesubstrates The invasive alien macroalga R. okamurae was detected in the south of Spain for the first time in the Strait of Gibraltar area in 2015 (Altamirano Jeschke etal. 2016; García-Gómez etal. 2020), and in just a few years, it has spread rapidly over a wide area of the southern coast of Spain. Rugulopteryx okamurae is becoming the dominant macroalga species in the area of the Strait of Gibraltar (Sempere-Valverde etal. 2021). On the west coast of southern Spain, the presence of this invasive alien species has been detected up to Chipiona, and in the east, its presence has been described up to the coast of Benalmadena (García-Gómez etal. 2020). The presence of this macroalga is accumulating tons of biomass on the Andalusian coasts of southern Spain. The management and valorisation of the macroalga biomass is vital to reduce the spread of R. okamurae by vegetative propagules. The TS content of R. okamurae biomass was 197.6 ± 7.3g kg−1 and the VS/TS ratio was 0.52. The lower VS content reduces its potential as a raw material for the AD process (Lee etal. 2015). The C content of the invasive macroalgae biomass was 30.29 ± 0.50% and the N content of 1.99 ± 0.03%, reaching a C/N ratio of 15.2 ± 0.4 (Table1). The OMSW presents a high VS/TS ratio content (0.81). The C content of the OMSW was 49.70 ± 0.38%, meanwhile its N content was 1.58 ± 0.06%, keeping a C/N ratio of 30.4. The main characteristics of the inoculum and both substrates (TS, VS, C/N, COD, pH) are given in Table1. Co‑digestion methane yield Accumulated volumes of methane for the different substrates and combinations of substrates tested are presented in Fig.1. The methane yield versus time curves of the evaluated experiments were characterized by having a typically exponential form up to 23days and then keeping the methane production about constant until day 28 (Fig.1). The methane yield for all the substrates and mixtures tested reached an asymptotic value until the experiment concluded (Fig.1). The highest methane production in the study was found in the 1 R. okamurae-3 OMSW co-digestion mixture, for which 454 ± 69mL CH4 g−1 VS was obtained. The co-digestion mixture 1 R. okamurae- 3 OMSW improved the methane yield of 0 R. okumurae-1 OMSW and 1 R. okamurae-0 OMSW by 8.6% and 157%, respectively (0 R. okamurae-1 4116 Journal of Applied Phycology (2021) 33:4113–4124 1 3 OMSW = 418 ± 56mL CH4 g−1 VS and 1 R. okamurae-0 OMSW = 177 ± 15mL CH4 g−1 VS). The co-digestion mixture 1 R. okamurae- 1.5 OMSW achieved a methane yield of 368 ± 42mL CH4 g−1 VS. Similar value was found in the codigestion mixture 1 R. okamurae- 1 OMSW (380 ± 15mL CH4 g−1 VS), not observing significant differences between these co-digestion mixtures. The lowest methane yield values were obtained for the AD of the invasive alien macroalgae alone (R. okamurae, 177 ± 15mL CH4 g−1 VS). Benefits ofanaerobic co‑digestion The inoculum is a microbial consortium that carries out the anaerobic digestion process. One of the major drawbacks of AD can be found in the susceptibility of the microbial community to changes in operational parameters and the presence of inhibitory substances, mainly in the final step of the process, methanogenesis. The accumulation of inhibitory substances during AD process is often a result from feeding the reactor with unbalanced C/N ratio substrates. There are numerous anaerobic codigestion studies where it has been seen that a feeding substrate C/N ratio of between 25 and 30 is the optimal one to avoid the formation of inhibitory substances in the process and, therefore, to facilitate the anaerobic process of methane production (Allen etal. 2013). Rugulopteryx okamurae biomass used in this study had a C/N ratio of 15.2 (Table1), being below the optimum for AD process. Co-digestion with a C-rich substrate, such as OMSW (C/N ratio 31.4), helped to balance the ratio, obtaining values of 27.4, 24.9, and 23.3 for co-digestions 1 R. okamurae- 3 OMSW; 1 R. okamurae- 1.5 OMSW; and 1 R. okamurae- 1 OMSW, respectively (Table2). Table2 shows the biodegradability values (percentage of organic matter eliminated after AD process) which were obtained for each co-digestion mixture and for the sole OMSW and macroalgae. The highest biodegradability was obtained with the 1 R. okamurae- 1 OMSW mixture (87%). Very similar values were obtained with the co-digestions of 1 R. okamurae-3 OMSW (85%) and 1 R. okamurae- 1.5 OMSW (81%). On the other hand, with the mono-digestions of R. okamurae and OMSW, biodegradability values of 77% and 63%, respectively, were reached. The AD of the sole R. okamurae presented a C/N ratio of 15.2 at the beginning of the experiment, being the final concentration of ammonium in its digestate 1645 ± 22mg L−1. At the end of the study, in OMSW reactors the final ammonium concentration was 970 ± 33mg L−1, indicating the lowest nitrogen concentration among the different substrates tested. The co-digestion mixture 1 R. okamurae- 3 OMSW, with a C/N ratio of 27.4, presented a final ammonium concentration of 1125 ± 12mg L−1. The 1 R. okamurae-1.5 OMSW mixture reached a C/N ratio of 24.9, obtaining a final ammonium concentration in the anaerobic reactors of 1352 ± 25mg L−1. And finally, the last mixture tested (1 R. okamurae- 1 OMSW) presented a C/N ratio at the beginning of the experiment of 23.3 and the final concentration of ammonium in the anaerobic reactor was 1338 ± 14mg L−1. Study ofpossible synergic effects The experimental methane yields observed for each codigestion mixture (Fig.1) were compared to the calculated or theoretical methane yields based on the OMSW and R. okamurae methane yields separately, according to Eq.(3): where 418 and 177 are the experimental methane yields (mL CH4 g−1 VSadded) obtained for OMSW and R. okamurae, respectively. X OMSW and Y R. okamurae are the proportions of OMSW and R. okamurae, respectively, in each codigestion mixture. (3) Calculated methane yield ( mL CH4g −1 VSadded ) =X ⋅ OMSW*(418)+Y ⋅R . okamurae *(177) Fig. 1 Methane yield versus time curves for Rugulopteryx okamurae (1 R.o.-0 OMSW), olive mill solid waste (0 R.o.-1 OMSW), and the different co-digestion mixtures tested. Values represent mean ± standard error Table 2 C/N ratio of Rugulopteryx okamurae (R. o.), olive mill solid waste (OMSW), and the different co-digestion mixtures. Biodegradability and ammonium concentration values at the end of the experiment. Values represent means ± standard error Biodegradability (%) Ammonium (mg L−1) C/N ratio 1 R. o.-0 OMSW 11 R. o.e-1 OMSW 1R. o.-1.5 OMSW 1 R. o.-3 OMSW 0 R. o.-1 OMSW 77 87 81 85 63 1645 ± 22 1338 ± 14 1352 ± 25 1125 ± 12 970 ± 33 15.2 ± 0.4 23.3 ± 0.6 24.9 ± 0.4 27.4 ± 0.2 31.4 ± 0.5 4117Journal of Applied Phycology (2021) 33:4113–4124 1 3 In all cases tested, experimental methane yield or BMP values were higher than the calculated methane yields for the three co-digestion mixtures tested (Table3). 26.8% for the co-digestion mixture 3 OMSW-1 R. okamurae, 13.9% for the co-digestion mixture 1.5 OMSW-1 R. okamurae, and 27.5% for the co-digestion mixture 1 OMSW-1 R. okamurae. Therefore, according to the increase in the BMP values, the biodegradability of the abovementioned co-digestion mixtures was also much higher than the biodegradability of the sole macroalga. Kinetic studies First‑order kinetic model Table4 shows the kinetic parameters obtained from the application of the first-order model to the experimental data of methane production-time corresponding to the different batch anaerobic digestion assays of all the mixtures and single substrates carried out. Values after ± represent the standard deviations of each parameter. Figure2 shows as an example the experimental points of methane production versus time in the same plot that the theoretical curve predicted by the first-order model for some of the different substrates (100% R. okamurae) and mixtures (1 R. okamurae—1.5 OMSW) tested. The good fit of the data indicated that the model adequately predicted the experimental data. Among the different co-digestion mixtures tested, the highest specific rate constants, k, was found for the mixtures 1 R. okamurae—1.5 OMSW and 1 R. okamurae- 3 OMSW with values of 0.066 and 0.042 days−1, respectively. These values were 5.1 and 3.2 times higher than that obtained for the single or mono anaerobic digestion of OMSW. On the other hand, the kinetic constant for the mixture 1 R. okamurae- 1.5 OMSW was 29% higher than that obtained for the anaerobic digestion of the sole macroalga. Table 3 Increase in the experimental methane yield compared to the calculated or theoretical value for the different mixtures of olive mill solid waste (OMSW) and Rugulopteryx okamurae (R. o.) Calculated (mL CH4 g−1 VS) Experimental (mL CH4 g−1 VS) Increase in methane production (%) 1 R. o.-0 OMSW 1 R. o.-1 OMSW 1 R. o.-1.5 OMSW 1 R. o.- 3 OMSW 0 R. o.- 1 OMSW 177 298 322 358 418 177 380 367 454 418 0 27.5 13.9 26.8 0 Table 4 Values of the kinetic constant obtained from the first-order model for the different substrates and co-digestion mixtures tested. Values represent mean ± standard error OMSW:olive mil solid waste; R. o.: Rugulopteryx okamurae; S.E.E.: standard error of estimate; Gmax: ultimate methane production (mL CH4 g−1 VS); k: specific rate constant or apparent kinetic constant (day−1); R2: determination coefficient Gmax (mL CH4 g−1 VS) K (day−1) R2S.E.E 1 R. o.-0 OMSW 1 R. o.-1 OMSW 1 R.o.-1.5 OMSW 1 R. o.- 3 OMSW 0 R. o.- 1 OMSW 190 ± 8 370 ± 11 460 ± 15 512 ± 22 486 ± 26 0.051 ± 0.006 0.017 ± 0.004 0.066 ± 0.004 0.042 ± 0.004 0.013 ± 0.004 0.998 0.994 0.995 0.995 0.992 7.42 14.03 12.22 15.38 16.68 Fig. 2 Variation in the experimental methane production with time (points) as well as the theoretical curve (solid lines) predicted by the first-order kinetic model for some of the cases tested (1 Rugulopteryx okamurae (R.o.)- 0 olive mill solid waste (OMSW) and 1 R.o.- 1.5 OMSW). Values represent mean ± standard error 4118 Journal of Applied Phycology (2021) 33:4113–4124 1 3 Transference function model Table5 summarizes the parameters obtained from the application of the transference function model to the experimental data of cumulative methane production time (Fig.1). The R2 values were higher than 0.988 in all cases. In addition, the low values of the standard deviations of the kinetic parameters and standard errors of estimates also pointed out a good fit of the experimental data to this suggested model for all cases assayed (Table5). Figure3 shows the experimental points of methane production versus digestion time as well as the theoretical curve predicted by the transference function model for some of the different mixtures (1 R. okamurae—3 OMSW and 1 R. okamurae—1.5 OMSW) tested as an example. The good adjustment of the data indicated that the model appropriately predicted the experimental data. Discussion The low VS content of the macroalgae is a major drawback for optimizing the anaerobic digestion process. Similar VS/TS ratio values of those obtaining in the experiment were previously described by Milledge and Harvey (2016), obtaining low methane yield values. The low C/N ratio value indicates the high nitrogen content of the biomass, placing the value of the C/N ratio below of the optimum value established for the biomethanization process (Zheng etal. 2021). In this work, co-digestion is proposed for optimizing the anaerobic process as a solution to valorise the R. okamurae biomass to prevent it from accumulating on Andalusian beaches and to be able to stop its expansion and for the need to adjust the C/N ratio to optimize the anaerobic digestion process. Due to its proximity and because it is one of the more polluting and main by-products in Andalusia (south Spain), the OMSW (main semi-solid by-product of the olive oil production process) is proposed in this study as C-rich co-substrate for anaerobic co-digestion together with the invasive alien macroalgae biomass. The high amount of VS present in OMSW makes it a highly potential substrate for AD process (Lee etal. 2015). This high C/N ratio value indicated the lignocellulosic nature and composition of the OMSW. Therefore, co-digestion was used in this study because the optimal C/N ratio for AD process is known to be between 25 and 30. All the C/N ratio values obtained with the co-digestion mixtures assayed were within the values established as optimal for the AD process (Allen etal. 2013). In previous biomethanization studies, the high nitrogen content of the macroalgae was indicated as one of the main factors that inhibit methane production (Costa etal. 2012). Co-digestion would not only help to balance the C/N ratio, but it also helps to dilute the presence of toxins and inhibitory substances for the anaerobic digestion process, such as salt, in the case of macroalgae and phenols and lignocellulosic material for the OMSW (Fernández-Rodríguez etal. 2020). The highest methane production in the study was found in the 1 R. okamurae-3 OMSW co-digestion mixture. The co-digestion mixture 1 R. okamurae- 3 OMSW improved the methane yield of 0 R. okumurae-1 OMSW and 1 R. okamurae-0 MSW by 8.6% and 157%, respectively. Sharply lower values were found by Elalami etal. (2020), who recorded a methane production value of 188 ± 19mL CH4 g−1 VS for the mixture 50% olive pomace (by-product of the traditional olive oil extraction mill) and 50% macroalgal residues of Gelidium sesquipedale. The lowest methane yield values were obtained for the AD of the invasive alien macroalgae alone. Previous studies carried out in batch mode using as substrate the macroalgae Saccharina latissima (brown alga) and Palmaria palmata (red alga) reported a methane yield of 209 ± 15mL CH4 g−1 VS and 257 ± 22mL CH4 g−1 VS, respectively (Jard etal. 2012). These values are somewhat higher than those obtained for R. okamurae (177 ± 15mL CH4 g−1 VS). Maiguizo- Diagne etal. (2019) used green macroalgae as AD substrate (Ulva lactuca and Codium tomentosusm) in BMP test, and achieved a final methane yield of 216mL CH4 g−1 VS. Recently, Thompson etal. (2021) studied different VS co-digestions mixtures of pelagic Sargassum (PS) and food waste (FW) in BMP tests. The mixture 75 PS:25 FW obtained a final methane yield of 97.46 ± 1.05mL CH4 g−1 VS. For the 50 PS:50 FW co-digestion mixture, the biomethane obtained was 182.33 ± 2.61mL CH4 g−1 VS. Table 5 Values of the parameters obtained from the transference function model for the different substrates studied. Values represent mean ± standard error OMSW: olive mil solid waste; R. o.: Rugulopteryx okamurae; Bmax: is the ultimate methane production; Rmax: is the maximum methane production rate; ʎ: is the lag time; S.E.E.: standard error of estimate; R2: determination coefficient Bmax (mL CH4 g−1 VS) Rmax (mL CH4 g−1 VS d−1)λ (d) R2S.E.E 1 R. o.-0 OMSW 1 R. o.-1 OMSW 1 R.o.-1.5 OMSW 1 R. o.- 3 OMSW 0 R. o.- 1 OMSW 179 ± 9 397 ± 17 461 ± 18 616 ± 13 496 ± 9 13.6 ± 0.7 18.4 ± 1.4 30.4 ± 1.3 29.6 ± 1.7 18.9 ± 1.5 8.3 × 10−10 1.5 × 10−8 3.8 × 10−10 0.0518 1.8 × 10−9 0.988 0.994 0.995 0.995 0.992 7.5 14.3 12.4 15.7 17.0 4119Journal of Applied Phycology (2021) 33:4113–4124 1 3 Meanwhile, the highest methane production was obtained in the mixture with the lowest PS concentration, 25 PS:75 FW; 201.67 ± 6.36mL CH4 g−1 VS. A similar pattern has been found in the present study, where the highest biomethane production was found in the mixture with the lowest content of the macroalgae R. okamurae (1 R. okamurae-3 OMSW; 454 ± 69mL CH4 g−1 VS). The low biodegradability values obtained (Table2) showed that the substrates used had a high complexity content. This result agrees with Sun etal. (2019). They reported a 75% of biodegradability of Laminaria digitata during anaerobic assays. Fernández-Rodríguez etal. (2019b) also found an OMSW biodegradability of 56%, very similar to that obtained in this experiment. In addition, the mentioned study also found an improvement in the biodegradability of different microalgae-OMSW mixtures compared with those obtained for the substrates separately. The biodegradability depends on the composition of the substrate used, being highly affected by the presence of organic compounds which are difficult to biodegrade, such as lignin (Sun etal. 2019). During AD process of the single invasive macroalgae (1 R. okamurae- 0 OMSW), a methane production of 177 ± 15mL CH4 g−1 VS was obtained, and at the end of the experiment, the biodegradability of the macroalgae was 77%. Previous studies found that the methane production during the anaerobic digestion of brown algae is limited by the inaccessibility of bacteria to complex sugars (Thompson etal. 2020). In addition, it should be considered that marine macroalgae tend to have a high content of salt, which can accumulate in the reactor, being able to achieve the inhibition concentration for the AD process (Ometto etal. 2018; Sun etal. 2019). According to Varnero Moreno (2011), for substrates with low C/N ratios, bacterial activity can be inhibited, due to the excessive amount of ammonia that can accumulate in the anaerobic reactor. It has been reported in the literature that ammonium values higher than 1500mg L−1 are toxic for the anaerobic digestion process (Yenigün and Demirel 2013). However, during the AD of sole OMSW, the C/N ratio is above the upper limit established in the bibliography as optimal for the anaerobic digestion process (C/N = 31.4). OMSW is a substrate with high carbon content which would lead to a slower degradation due to the low nitrogen content reducing the multiplication and development rate of bacteria (Varnero Moreno 2011). In the different co-digestion mixtures tested, the C/N ratio was balanced, and therefore, values below the limit established as ammonium toxic in the anaerobic reactors were obtained. In all cases tested, experimental methane yield or BMP values were higher than the calculated methane yields for the three co-digestion mixtures tested (Table3). The synergy effects of the OMSW and R. okamurae co-digestions with the above percentages were clearly shown with these results, being the co-digestion mixture 1 OMSW-1 R. okamurae the combination that achieved the highest synergic effect. In addition, the methane yields of the three mixtures tested were always higher than the methane yield of the sole R. okamurae (177mL CH4 g−1 VS), which is within the range of the methane yields reported for anaerobic digestion of other marine macroalgae such as Sargassum, Gracilaria, Laminaria, Ascophyllum, and Ulva, for which methane yields of between 140 and 280mL CH4 g−1 VS Fig. 3 Variation in the experimental methane production with time (points) as well as the theoretical curve (solid lines) predicted by the transference function model for some of the cases tested (1 Rugulopteryx okamurae (R.o.)-3 olive mill solid waste (OMSW) and 1 R.o.- 1.5 OMSW). Values represent mean ± standard error 4120 Journal of Applied Phycology (2021) 33:4113–4124 1 3 were achieved (Saratale etal. 2018). High contents in salts, phenolic compounds, heavy metals (in some cases), sulfides, and volatile acids stand out, along with low nitrogen content which has been reported as the characteristics that hinder and inhibit the anaerobic digestion process of macroalgae, and determine their low biodegradability (Saratale etal. 2018). Particularly, regarding R. okamurae, recent information (https:// www. elest recho digit al. com/ tag/ rugul opter ixokamu rae/ accessed 13 June 2021) indicates that the increase in sea water temperature, the lengthening of the summer months, and the changes in the wind regimes are promoting the proliferation of this macroalga, in places such as the Strait of Gibraltar, a phenomenon common in the hottest months of the year. This is causing the population of the area and tourists to have to live every summer with these seaweeds that are gaining more and more ground (García- Gómez etal. 2020). Added to the serious damage that the algae is causing to sectors such as fishing or tourism is the costs that the withdrawal of tonnes of this algal biomass means for municipalities such as Tarifa and Algeciras, which only in 2020 have already invested more than 100,000 euros each to carry out these tasks, given the flooding suffered by the beaches every time the Levante wind drags tons of these macroalgae towards the coast (https:// www. elest recho digit al. com/ tag/ rugul opter ixokamu rae/ accessed 13 June 2021). The first order kinetic model (Table4) shown the high R2 values as well as the low values of standard deviations of the kinetic constants and standard errors of estimates indicated that the experimental data correctly fit the proposed model. Beltran etal. (2016) showed that the co-digestion of waste activated sludge (WAS) and the microalga Chlorella sorokiniana (Ch.s.) with percentages of 75%WAS-25% Ch.s. and 50% WAS-50% Ch.s. were 12% and 42% higher than that obtained for the single microalga. The high kinetic constant values for the abovementioned mixtures (1 R. okamurae—1.5 OMSW and 1 R. okamurae- 3 OMSW) can be attributed to the appropriate C/N ratios of the both mixtures (24.9 and 27.4, respectively) compared with that observed for the single substrates (31.4 and 15.2 for 1 OMSW and 1 R. okamurae, respectively). Recent studies have demonstrated efficient anaerobic processes when the C/N ratio is within the range of 25 and 30 (Yan etal. 2015). Carbon helps to supply energy to the anaerobic microorganisms whereas nitrogen is crucial for the growth of the bacterial population. In addition, very recent studies (Fernández-Rodríguez etal. 2021a) have shown that the anaerobic co-digestion of 95% OMSW-5% Dunaliella salina achieved specific rate constant values of 0.19 days−1, a value 2.8 times higher than the maximum value reached in the present study for the mixture 1.5 OMSW-1 R. okamurae. The mentioned kinetic constant value for the mixture 95% OMSW-5% D. salina was also 12% higher than that obtained for the single OMSW, similar trend than that observed in the present research work. Dunaliella salina is a halophilic flagellate alga able to tolerate and acclimate to a wide range of salinities (Fernández-Rodríguez etal. 2021a). On the other hand, the anaerobic co-digestion of the macroalga Sargassum sp. with glycerol (Gly) (with 0.5% TS Sargassum sp. and 3.0g L−1 of Gly) and with waste fried oil (WFO) (1.31% TS Sargassum sp. and 0.88g L−1 of WFO) increased 38% and 19%, respectively, the specific rate constant when compared to the value obtained in the mono-digestion of the single macroalga (Oliveira etal. 2015). In addition, similar first order kinetic constant values to that obtained for 1 R. okamurae (0.051 day−1) were reported for anaerobic mono-digestion of fresh Ulva lactuca (0.06 day−1), commonly known as sea lettuce. However, the kinetic constant increased up to values of 0.12 and 0.15 day−1 when fresh Ulva was co-digested with dairy slurry at ratios of 75% fresh Ulva-25% slurry and 50% fresh Ulva-50% slurry, respectively (Allen etal. 2013). The good adjustment of the data indicated that the transference model appropriately predicted the experimental data. As can be seen, the highest value of the maximum methane production rate (Rmax) was found for the mixture 1 R. okamurae—1.5 OMSW (30.4mL CH4 g−1 VS day−1), which was 1.6 and 2.2 times higher than the corresponding to the mono-digestions of the single OMSW (18.9mL CH4 g−1 VS day−1) and sole R. okamurae (13.6mL CH4 g−1 VS day−1), respectively. In the same way, it was recently reported that the maximum methane production rate, Rmax, for the anaerobic co-digestion of the mixture 95% OMSW- 5% D. salina increased by 34.7% compared with the value obtained for the single OMSW (Fernández-Rodríguez etal. 2021a). The lowest values of Rmax and ultimate methane production (Bmax) were found for the single monodigestion of R. okamurae. Most macroalgae usually have high contents in salts, which can lead to an accumulation of them in anaerobic digesters, causing an inhibition of the microorganisms responsible for the process of biomethanization (Tabassum etal. 2017). Others macroalgae, such as U. latuca and Ascophyllum nodosum, may contain high levels of sulfide and polyphenols, respectively, compounds that from certain concentrations are inhibitors of the anaerobic digestion process (Tabassum etal. 2016). In general, this type of organisms also tends to have a high protein content that often limits the energy produced in anaerobic digestion processes due to their inhibition by the high concentration of ammoniacal nitrogen released (Tabassum etal. 2017), which have an inhibitory effect on anaerobic microorganisms, especially methanogenic archaea. On the other hand, the maximum methane production rate (Rmax) values found in anaerobic co-digestion of 4121Journal of Applied Phycology (2021) 33:4113–4124