Front. Astron. Space Sci. Frontiers in Astronomy and Space Sciences Front. Astron. Space Sci. 2296-987X Frontiers Media S.A. 1198689 10.3389/fspas.2023.1198689 Astronomy and Space Sciences Original Research Stoichiometric model of a fully closed bioregenerative life support system for autonomous long-duration space missions Vermeulen et al. 10.3389/fspas.2023.1198689 Vermeulen Angelo C. J. 1 2 * Papic Alvaro 2 Nikolic Igor 1 Brazier Frances 1 1 Systems Engineering and Simulation, Faculty of Technology, Policy and Management, Delft University of Technology, Delft, Netherlands 2 SEADS (Space Ecologies Art and Design), Temse, Belgium

Edited by: Cyprien Verseux, University of Bremen, Germany

Reviewed by: Lucie Poulet, UMR6602 Institut Pascal (IP), France

Gabriela Soreanu, Gheorghe Asachi Technical University of Iași, Romania

*Correspondence: Angelo C. J. Vermeulen, a.c.j.vermeulen@tudelft.nl
16 08 2023 2023 10 1198689 01 04 2023 07 07 2023 Copyright © 2023 Vermeulen, Papic, Nikolic and Brazier. 2023 Vermeulen, Papic, Nikolic and Brazier

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Bioregenerative life support systems (BLSS) are vital for long-duration and remote space missions to increase mission sustainability. These systems break down human waste materials into nutrients and CO2 for plants and other edible organisms, which in turn provide food, fresh water, and oxygen for astronauts. The central idea is to create a materially closed loop, which can significantly reduce mission mass and volume by cutting down or even eliminating disposable waste. In most BLSS studies only a fraction of the resources, such as food, are provided by the system itself, with the rest taken on board at departure or provided through resupply missions. However, for autonomous long-duration space missions without any possibility of resupply, a BLSS that generates all resources with minimal or no material loss, is essential. The goal of this study is to develop a stoichiometric model of a conceptually fully closed BLSS that provides all the metabolic needs of the crew and organisms. The MELiSSA concept of the European Space Agency is used as reference system, consisting of five interconnected compartments, each inhabited by different types of organisms. A detailed review of publicly available MELiSSA literature from 1989 to 2022 revealed that no existing stoichiometric model met the study’s requirements. Therefore, a new stoichiometric model was developed to describe the cycling of the elements C, H, O, and N through all five MELiSSA compartments and one auxiliary compartment. A compact set of chemical equations with fixed coefficients was established for this purpose. A spreadsheet model simulates the flow of all relevant compounds for a crew of six. By balancing the dimensions of the different compartments, a high degree of closure is attained at steady state, with 12 out of 14 compounds exhibiting zero loss, and oxygen and CO2 displaying only minor losses between iterations. This is the first stoichiometric model of a MELiSSA-inspired BLSS that describes a continuous provision of 100% of the food and oxygen needs of the crew. The stoichiometry serves as the foundation of an agent-based model of the MELiSSA loop, as part of the Evolving Asteroid Starships (E|A|S) research project.

space exploration human spaceflight bioregenerative life support waste processing food production ecosystem modeling simulation MELiSSA section-at-acceptance Astrobiology

香京julia种子在线播放

    1. <form id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></form>
      <address id=HxFbUHhlv><nobr id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></nobr></address>

      1 Introduction

      Bioregenerative life support systems (BLSS) will be a key component of future long-duration space exploration, as they can reduce mission mass and volume (Imhof et al, 2017; Audas et al, 2022). A BLSS recycles human waste by feeding it to an artificial ecosystem consisting of a range of different organisms, such as higher plants, microalgae, and bacteria. These organisms gradually break down the waste and in the process provide fresh food, generate oxygen, and clean water for the crew (Nelson et al, 2010; Escobar and Nabity, 2017; Pannico et al, 2022). As a result, there is less reliance on constant resupply from Earth, which greatly reduces the material footprint needed for extended missions in space (Macelroy and Averner, 1978; Häder, 2020; Audas et al, 2022). Because of its constant fresh food provision, a BLSS also solves the problem that essential nutrients in processed and prepackaged food are deficient and degrade over time (Cooper et al, 2017; Carillo et al, 2020). Biological systems are also capable of self-repair, unlike mechanical systems that need separate systems to be repaired. Surviving organisms and tissues in a damaged BLSS can potentially recover to their original size on their own, even after significant damage (Bartsev et al, 1996).

      MELiSSA is a BLSS concept developed by ESA with a consortium of 15 international partners and the involvement of approximately 50 organizations (MELiSSA Foundation, 2023). It is an artificial ecosystem consisting of five interconnected compartments inhabited by higher plants, microalgae, microorganisms, and humans (abbreviated as C1 to C5, Figure 1; Table 1). Each of the compartments has its own specific metabolic function within the entire MELiSSA loop. First, all human waste is broken down through a sequence of three bioreactor types: a thermophilic anaerobic compartment (C1), photoheterotrophic compartment (C2) and nitrifying compartment (C3). This results in nutrients and CO2 for microalgae and plants (C4a and C4b) that, in turn, provide food, oxygen, and fresh water for the crew (C5), thus closing the loop (Hendrickx and Mergeay, 2007; Clauwaert et al, 2017; Vermeulen et al, 2020a). The operational pilot plant at Universitat Autònoma de Barcelona consists of a number of connected MELiSSA compartments and demonstrates a part of the system’s entire metabolic loop (Gòdia et al, 2004; Garcia-Gragera et al, 2021). Mathematical models of the MELiSSA loop have been developed and improved since the beginning of the MELiSSA project in 1989 (MELiSSA, 1989; Dussap and Gros, 1991). Some of these models are used for predictive control, to be run alongside the actual physical system. The MELiSSA models help in understanding the dynamics of mass flow fluctuations and achieving long-term reliability of the system (Cornet et al, 2001; Ciurans Molist et al, 2020; Poulet et al, 2020). Such ecosystem models can also be used to simulate the impact of unforeseen perturbations, deliberate interventions, or design changes in the system (Volk and Rummel, 1987; Pilo Teniente, 2015; Vermeulen et al, 2019).

      Overview of the MELiSSA loop with its five main compartments and a schematic overview of the different material flows. Some of these material flows have been adjusted in the stoichiometric model. See Section 3 Stoichiometric model assumptions for more details. Image by the MELiSSA Foundation.

      Systematic review of the stoichiometric descriptions of mass flows in the various compartments of the MELiSSA loop, as documented in publicly available literature. Only modeling studies containing at least one chemical equation are listed. “Full closure” indicates that all material outputs are utilized as inputs without the need for any external supplies. “Chemical equations” indicates the number of involved chemical equations that are factually listed in the publication. “Stoichiometrically balanced” indicates that the publication includes the stoichiometric coefficients for all fixed equations. Stoichiometric coefficients are described as dynamic when coefficients within at least one chemical equation fluctuate throughout a simulation run. The macromolecule composition is described as dynamic when the proportion of CHON elements in at least one chemical formula varies over time. C, compartment; N/A, not applicable.

      MELiSSA compartments Full closure Chemical equations Stoichiometrically balanced Stoichiometric coefficients Macromolecule composition References
      C1-C5 a No Partly Yes Fixed Fixed MELiSSA (1989)
      C1-C5 a No Complete Yes Fixed Fixed Dussap and Gros (1991)
      C1-C5 a No Complete Yes Fixed Fixed Dussap et al (1993)
      C1-C5 a No Complete Yes Fixed + dynamic Fixed + dynamic Poughon et al (1994)
      C1-C5 a No Complete Yes Fixed Fixed Fulget (1996)
      C1-C5 a No Examples Yes Fixed + dynamic Fixed + dynamic Fulget et al (1999)
      C1-C5 No Complete Incomplete Fixed + dynamic Fixed + dynamic Poughon et al (2000)
      C1-C5 No Complete Yes Fixed + dynamic Fixed + dynamic Poughon (2007b)
      C1-C5 No Complete Yes Fixed + dynamic Fixed + dynamic Guirado and Podhajsky (2008)
      C1-C5 No Examples Yes Fixed + dynamic Fixed + dynamic Poughon et al (2009)
      C1-C5 No Complete Yes Fixed + dynamic Fixed + dynamic Thiron (2020)
      C1-C4a No info Complete No None Fixed Hendrickx et al (2006)
      C2+C5 No Complete Yes Fixed Fixed Poughon (1994b)
      C4a-C5 No Partly Yes Fixed + dynamic Fixed + dynamic Poughon (2007a)
      C4a-C5 No Partly Yes Fixed Fixed Pilo Teniente (2015)
      C4a-C5 No Partly Yes Fixed Fixed Alemany et al (2019)
      C2 N/A Complete Yes Fixed Fixed Cornet and Dussap (2000)
      C3 N/A Complete Yes Fixed Fixed Cruvellier et al (2016)
      C4a N/A Complete Yes Fixed Fixed Poughon (1994a)
      C4a N/A Complete Yes Fixed Fixed Cornet et al (1998)
      C4a N/A Complete Yes Fixed Fixed Cornet et al (2003)
      C4a N/A Complete Yes Dynamic Dynamic Poughon et al (2020)
      C4b N/A Complete No None No info Maclean et al (2010)
      C4b N/A Complete Yes Fixed Fixed Poulet et al (2020)

      Only C4a (Limnospira compartment), no C4b (higher plant compartment).

      To describe the material flows in an ecosystem model, the stoichiometric relations governing these flows need to be developed (Volk and Rummel, 1987; Begon and Townsend, 2021). Several authors have described mass flow models of biological life support systems (Volk and Rummel, 1987; Garland, 1989; Loader et al, 1997; Finn, 1998), but not all studies contain the underlying chemical stoichiometric equations. The stoichiometric equations describing mass flows in the MELiSSA BLSS have been published in publicly-available MELiSSA literature, both in peer-reviewed papers and research reports disseminated by the MELiSSA Foundation (Table 1). In less than half of the listed studies, all five MELiSSA compartments were modeled. The other studies describe either a smaller number of compartments, or just one. The loop is never fully closed in any of the listed studies. There are a number of reasons for this. First of all, it may indicate that certain compounds are not part of the regenerative logic of the system and are not generated within the loop, and therefore need to be fully supplied from the outside (Dussap et al, 1993; Poughon et al, 2009). Secondly, it may indicate an output that does not get regenerated, such as a residual indigestible material (Fulget et al, 1999; Poughon et al, 2000). Because of the resulting accumulation of that material, external supply is required to keep the system running. And thirdly, it may be a deliberate system design choice that only a fraction of a particular compound or resource is generated by the BLSS. Several MELiSSA studies, for example, explicitly specify that only a part of the needed food is generated by the loop (Poughon et al, 2000; Guirado and Podhajsky, 2008; Thiron, 2020). The rest is supplied externally. One of the main motivations for such a choice is the limited space that is available in common spacecraft and habitat designs to integrate bioreactors and plant growth chambers (Johnson et al, 2021; Gorce et al, 2022).

      Almost three-quarters of the researched MELiSSA studies list all the stoichiometric equations of their models. The other studies limit themselves to a few selected equations as examples. The stoichiometric coefficients in the chemical equations can either be fixed beforehand, or dynamically calculated on the fly, during the simulation run. For example, in simulating the production of Limnospira indica in C4a, the elementary composition of its biomass is determined by the variable light irradiation inside the photobioreactor (Guirado and Podhajsky, 2008; Poughon et al, 2009). This then necessitates the recalculation of all subsequent stoichiometric equations that contain this dynamic empirical biomass formula.

      Whenever creating a stoichiometry to describe mass flows in an ecosystem, choices need to be made about which compounds will be considered, and which level of detail will be used to describe the ecosystem’s biochemical processes (Dubitzky et al, 2011; Begon and Townsend, 2021). This choice is influenced by the level of knowledge about the individual metabolism of the different species in the ecosystem, but ultimately depends on the objective and scope of the study. In light of this, the reviewed MELiSSA studies contain varying assumptions concerning the role of different compounds such as lipids, VFAs, and CO2. This is partially due to the fact that over time the MELiSSA models became more detailed, incorporating more processes. Early MELiSSA papers assume that all lipids consumed by the astronauts are oxidized and no traces of lipids are found in the feces (Dussap and Gros, 1991; Dussap et al, 1993). This is later adjusted, and lipids become an explicit part of the feces input in C1 (Poughon et al, 2000). The mix of VFAs that is generated as output of C1 varies quite significantly between the different studies. In the earliest studies only acetic acid and butyric acid are considered (Dussap and Gros, 1991; Dussap et al, 1993), but in later studies, a broader range of VFAs is included in the stoichiometry (Poughon, 2007b).

      The stoichiometric model presented in this paper provides the basis for an agent-based model (ABM) of the MELiSSA loop to study the interactions and emergent behavior of the MELiSSA system from a theoretical perspective. MELiSSA is fundamentally a complex system, and properties such as nonlinearity, stochasticity and chaos may play an important role in the overall behavior of the system, just like in any ecosystem (Bjørnstad, 2015). Because all functions in the MELiSSA loop are tightly coupled and since there are no large buffers in space (e.g., no large atmosphere or water bodies), slight variations in any of the loop’s pathways can have dramatic consequences for the entire loop (Macelroy and Averner, 1978; Poulet et al, 2018). It is precisely such tightly coupled systems that are potentially susceptible to irregular fluctuations (Bjørnstad, 2015). ABMs offer a modeling approach that helps to gain insight in such dynamics and helps to explore different design strategies for improved system robustness. ABMs are used to understand those aspects of the behavior of a system, that are not easily represented in population-level differential equations, such as variation among individuals or variation of individuals during their life cycle (Deangelis and Grimm, 2014). Such exploratory models of complex systems lead to insights about essential mechanisms and principles, a use of modeling that is different from modeling aimed at making accurate predictions (Forrest and Mitchell, 2016). MELiSSA modeling studies to date focus on predictive modeling because this is a critical step to develop a durable and reliable BLSS. The MELiSSA modelling approach is mechanistic, allowing for better understanding and improved predictions. With its exploratory approach, this study can be considered as being complementary to the existing MELiSSA modeling efforts. In order to build an ABM of MELiSSA, a stoichiometric model is needed that describes essential mass flows between all involved agents (bacteria, microalgae, higher plants and humans). The key metabolic processes in all MELiSSA compartments need to be described with enough detail to enable studying the interactions. At the same time, they need to remain abstracted enough to tlead to an efficient model and retain relevance for the objectives of the ABM which are centered around investigating the impact of heterogeneity in the system and capturing nonlinear dynamics.

      The research objective of this paper is to develop a stoichiometric model of a hypothetical BLSS that provides a crew with fresh food and oxygen during long-duration space missions, without the need for any resupply. The envisioned stoichiometric model should conform to the following requirements:

      • The system provides 100% of all the food and oxygen needed by the crew, and 100% of the resources needed by every other organism in the loop.

      • The stoichiometric model is fully closed. Every single output is used as input, and no additional external supplies are needed during steady state.

      • No kinetic limitations of the stoichiometric reactions are used. Macromolecules do not have a dynamic composition (e.g., biomass composition does not depend on light), and all stoichiometric coefficients are fixed. Such a static approach is sufficient for the current goals of the ABM of MELiSSA that is developed on top of this stoichiometry.

      The current literature (listed in Table 1) has, to the authors’ knowledge, not included such a stoichiometric model of the MELiSSA loop. As described above, the loop is always considered to be partially open, primarily for practical reasons. The research reported in this paper is more theoretical in nature, focusing on the internal dynamics of hypothetically fully closed systems.

      2 Materials and methods

      This section describes the empirical formulas, stoichiometric equations, and the setup of a static spreadsheet model.

      2.1 Empirical formulas

      Table 2 provides an overview of the empirical formulas of the large biomolecules. More detailed descriptions of the origin of each formula are then discussed. The empirical formulas of biomass and feces are provided for comparative purposes only, since both biomass and feces are represented as proportions of carbohydrates, proteins and lipids.

      Empirical formulas of all large biomolecules used in the stoichiometric model.

      Compound Empirical formula References and notes
      Carbohydrates CH1.6667O0.8333 General polysaccharides
      Proteins CH1.5900O0.3100N0.2500 Hu et al, 2010, Fu et al, 2016
      Lipids CH1,9216O0,1177 Tripalmitin
      Biomass Rhodospirillum rubrum CH1.6472O0.3598N0.1788 18% carbohydrates, 72% proteins, 10% lipids, literature
      Biomass Limnospira sp. CH1.6472O0.3598N0.1788 18% carbohydrates, 72% proteins, 10% lipids, literature
      Biomass higher plants (edible) CH1.6889O0.6090N0.05487 70% carbohydrates, 20% proteins, 10% lipids, literature
      Biomass higher plants (non-edible) CH1.6667O0.8333 100% carbohydrates, 0% proteins, 0% lipids, literature
      Feces CH1.6676O0.6028N0.07715 66% carbohydrates, 28% proteins, 6% lipids, calculated
      2.1.1 Carbohydrates, proteins, and lipids

      The formula for carbohydrates is the general formula for polysaccharides, used in many of the MELiSSA modeling studies (e.g., in (Dussap et al, 1993) and (Thiron, 2020)). There are, however, a range of different general formulas used for proteins in the MELiSSA literature (Fulget, 1996; Duatis et al, 2008; Thiron, 2020). Our approach was to first calculate the average protein composition of the different organisms generating biomass in our model, and then compare this result to general protein formulas published in multiple BLSS studies. In our stoichiometric model, the higher plants are represented by a generalized “ideal plant”, in line with previous MELiSSA studies (Poughon, 2007b; Guirado and Podhajsky, 2008). It is conceived as an average of the nutritional properties of bread wheat, durum wheat, potato and soybean. The CHON composition of the protein of this ideal plant is then averaged with that of Limnospira and Rhodobacter, all using values listed in MELiSSA literature (Dussap et al, 1993; Duatis et al, 2008). This results in CH1.5799O0.3134N0.2558. This composition comes very close to CH1.5900O0.3100N0.2500 used as general protein formula in studies of the Lunar Palace BLSS (Hu et al, 2010; Fu et al, 2016), and hence this latter formula was selected. Palmitic acid is often used as lipid in stoichiometries of ecosystems (Volk and Rummel, 1987; Guirado and Podhajsky, 2008). However, vegetable oils mainly consist of more complex triacylglycerols (Zambelli et al, 2015). Therefore in this stoichiometric model, lipids are represented by tripalmitin, consisting of one glycerol molecule bonded to three palmitic acid molecules. Just like with palmitic acid, the fermentation of glycerol also leads to VFAs, further contributing to the overall VFA output of C1.

      2.1.2 Biomass

      Throughout the MELiSSA literature the empirical formulas of the biomass of specific organisms differ to varying degrees. The empirical formula that describes the biomass of the purple sulphur bacteria Rhodobacter varies, for example, between (Dussap et al, 1993), (Hendrickx et al, 2006), and (Thiron, 2020). This is to be expected since empirical formulas are always approximations. Moreover, in several MELiSSA modeling studies, biomass composition is dependent on environmental parameters such as light input, and consequently the biomass composition varies throughout the course of a simulation (Table 1). In the current stoichiometric model, biomass is consistently represented as a combination of carbohydrates, proteins, and lipids, and not as one single molecule. This makes it possible to trace the levels of each individual compound throughout the loop. To establish the composition of the different types of biomass, in this study an estimate was made based on averaged values from literature. As validation, the corresponding empirical formula was then compared with formulas published in MELiSSA literature.

      In correspondence with the MELiSSA study of (Hendrickx et al, 2006), the same biomass composition was chosen for both Rhodospirillum rubrum and Limnospira sp. Based on literature (e.g., Dussap et al, 1993; Kobayashi and Kobayashi, 1995; Teuling et al, 2017), the following biomass composition was established (dry weight, or DW): 18% carbohydrates, 72% proteins, and 10% lipids. Both organisms contain high amounts of protein, hence their interesting potential as food sources (Vrati, 1984; Clauwaert et al, 2017; Muys et al, 2019). Using the previously established CHON compositions of carbohydrates, proteins and lipids, the overall biomass composition corresponds to the empirical formula CH1.6472O0.3598N0.1788. This is in the same range as the biomass formulas used in the MELiSSA modeling studies of (Hendrickx et al, 2006) and (Thiron, 2020).

      Carbohydrates are typically the most abundant compound in plants, encompassing both structural (e.g., cellulose) and non-structural (e.g., starch and glucose) carbohydrates (Moore and Hatfield, 1994). As explained before, an ‘ideal plant’ is used in the current model. Four energy- and/or protein-rich crops were selected and combined based on (Paradiso et al, 2013): bread and durum wheat, potato and soybean. Averaging the carbohydrate, protein and lipid fractions of these four crops (Molders et al, 2012; Stasiak et al, 2012; Page and Feller, 2013; Paradiso et al, 2013), leads to the following edible plant biomass composition (DW): 70% carbohydrates, 20% proteins, and 10% lipids. Taking into account the above CHON compositions of carbohydrates, proteins and lipids, this corresponds to the empirical macroformula CH1.6889O0.6090N0.05487. This is in line with the “food” descriptions in other MELiSSA studies (Fulget, 1996; Thiron, 2020). The inedible part of the plant is considered to be 100% carbohydrates, assuming that it is only composed of cellulose.

      2.1.3 Feces

      Just like with biomass, feces are expressed in the model as a combination of carbohydrates, proteins and lipids. The composition of feces was calculated by stoichiometrically balancing the equation for human metabolism (see 5.5). The proportions of carbohydrates, proteins and lipids as reported in literature [review in (Rose et al, 2015)] were used as a guideline to solve the stoichiometry. The resulting empirical formula is CH1.6676O0.6028N0.07715.

      2.2 Establishing the stoichiometric equations

      As a general method to balance the equations that represent the reactions occurring in each compartment, each one is stated as a linear programming problem–or non-linear when the constraints of the particular reaction so require–where the goal is to minimize the total amount of matter used. The variables are defined as the number of moles of each compound on either side of the equation. The objective function is subject to the linear constraint of keeping the number of atoms of a given element equal on both sides of the equation. Additionally for each problem ad hoc constraints are described e.g., to keep proportions of VFA’s produced similar to experimental data on C2. The solving engine for these problems is MS Excel’s built-in optimization tool. For linear problems Simplex is used and for non-linear problems GRN non-linear, both on their default configurations settings. In reporting throughout this paper 4 significant decimals are used, but during calculations in Excel 14 decimals were used.

      This stoichiometry is a full revision of a version of that was concisely published in two IAC conference papers (Vermeulen et al, 2018b; 2019) and presented at two MELiSSA conferences (Vermeulen et al, 2018a; 2020b). The revisions were based on a detailed validation session and on further literature study. The main changes were the following:

      • All biomass compositions were re-evaluated (Limnospira, Rhodospirillum and higher plants) to more accurately represent their carbohydrate, protein and lipid content.

      • No more distinction between fecal and bacterial protein, in line with the singular use of carbohydrates and lipids in the stoichiometry.

      • The VFA spectrum was extended beyond acetic acid and butyric acid.

      • Instead of representing urine as NH3, it is now represented as urea, in order to bring carbon into the urine cycling. Hence, a new subsystem with its own stoichiometry representing the breakdown of urea was added.

      2.3 Static spreadsheet model to test the closure level

      To calculate the steady state flow values for each stream, each compartment’s equation is linked, thus building an equation based model that represents the whole cycle. To solve it, the one a priori known value is used: the consumption rate of the crew compartment, which has been set at 3,000 kcal/day/person. Knowing these requirements for C5, the required food outputs of C4a and C4b are calculated. By continuing to work backwards through the entire stoichiometry, the flow rates through the entire loop could be established over a period of 1 day for a crew of six.

      Since there is a subcycle between C1 and C2 the exchange flows between them need to be calculated (Figure 2). Having the required steady state flows into C1, and out of C2 (except the subcycle flows into each other), the process of exchange between C1 and C2 can be examined. This requires an iterative approach in which, given an initial output of VFAs from C1 to C2, a subsequent C2 to C1 flow of biomass is calculated. This, in turn, updates the flow from C1 to C2 again. This eventually converges into a stable exchange between C1 and C2, with complete processing of all VFAs. This then results in cycle-wide steady state flow rates.

      Diagram representing the subcycle between C1 and C2, and the process of converging to a stable exchange of VFAs and biomass between the two compartments.

      3 Stoichiometric model assumptions

      This section describes the general assumptions behind the model and those per compartment.

      3.1 General assumptions

      This stoichiometry is based on the four CHON elements; sulfur and phosphor are not yet included. Throughout the entire loop, food and feces are described as a proportion of carbohydrates, proteins and lipids. The water in this stoichiometric model is what is used in the chemical reactions. The actual water needs of a human crew (drinking and hygiene water) have not been included. Net biomass production only occurs in C2, C4a, and C4b because this biomass is a material input for other compartments. For the other compartments it is assumed there is no growth, and hence no net biomass production occurs. This has an impact on the stoichiometric approach for each compartment, which differs when describing a growth metabolism vs. a maintenance metabolism without growth (Cruvellier et al, 2016). It is further assumed that all waste materials are fully digested in C1, and no undigestible residues are left. And there is no limitation in terms of available space and hence the number of bioreactors and growth chambers that can be included.

      3.2 MELiSSA loop adjustments

      Several adjustments have been made to the official concept of the MELiSSA loop, as illustrated in Figure 1. They can be summarized as follows (more details further down in this section):

      • The Rhodospirillum biomass that is produced in C2 is not used for human consumption, but is sent to C1 where it is anaerobically digested, together with the other organic waste.

      • C2 does not produce any net CO2.

      • Urea is included as an individual waste product and is processed by C3.

      • C3 consists of two subcompartments, C3a and C3b, to respectively process ammonia and urea.

      • An auxiliary compartment has been added to burn surplus hydrogen and generate water.

      3.3 Assumptions for Compartment 1: anaerobic fermentation

      In C1, the thermophilic anaerobic bacteria can convert up to 70% of the waste under experimental conditions (Lasseur et al, 2010). This is because cellulose and hemicellulose can only be hydrolyzed slowly, further hampered by the fact that both are often intermeshed with lignin in lignocellulose, which is even more difficult to break down (Liebetrau et al, 2017). In this model it is assumed that all organic waste is fully converted in C1. In the MELiSSA system up to 71% of the organic nitrogen in the waste is converted into ammonium (Clauwaert et al, 2017). In our stoichiometric model it is assumed that this is 100%. The VFA spectrum produced by C1 consists of four VFAs: acetic acid, propionic acid, butyric acid and valeric acid. These VFAs are characteristically produced during acidogenesis (Jiang et al, 2013; Khan et al, 2016).

      3.4 Assumptions for Compartment 2: <italic>Rhodospirillum</italic> culture

      In the original conceptualization of MELiSSA, it was assumed that in addition to the microalgae and crops from C4, the edible purple sulfur bacteria from C2 would also serve as a food supply for the crew (Lasseur et al, 1996; Albiol et al, 2000; Figure 1). However, because of close proximity, there is always a chance that pathogens associated with fecal matter in C1 pass through membrane filtration and end up with the edible bacteria in C2 (Clauwaert et al, 2017). Therefore, the produced biomass from this compartment is sent back to C1 for anaerobic degradation. This is in correspondence with other MELiSSA modeling studies (Poughon et al, 2000; Waters et al, 2004; Thiron, 2020). The overall diet of the crew should only contain a limited amount of microbial single-cell protein, because of increased nucleic acid contents and lower palatability (Tusé and Miller, 1984; Tranquille and Emeis, 1996). The crew is already using another microbial food source (Limnospira from C4a, Figure 1), and this might be an additional argument to avoid consuming the bacteria from C2. Theoretically, the assimilation of the VFAs by purple sulfur bacteria in C2 could generate a net CO2 output (Hendrickx and Mergeay, 2007). However, in practice, this has not been observed yet in MELiSSA bioreactor experiments with C1 effluent as VFA source (Mastroleo et al, 2020).

      3.5 Assumptions for Compartment 3: processing ammonia and urea

      Analogous to C4, this compartment also consists of two subcompartments, C3a and C3b. C3a processes ammonia and produces nitrate. C3b is a urine treatment unit based on C3a, and also produces nitrate. It is assumed that all incoming ammonia and urea are converted to nitrate.

      3.6 Assumptions for Compartment 5: human metabolism

      The assumptions about the human metabolism are summarized in Table 3. The reference caloric value of the higher plants is estimated to be 4,000 kcal/kg dry weight (DW). This is based on the reported caloric values of the four crops that constitute the ideal plant (adjusted for their moisture content): bread wheat and durum wheat (Kaleta and Górnicki, 2013; Knowledge4 Policy, 2021; WebMD, 2021), potato (Decker and Ferruzzi, 2013; U.S. Food & Drug Administration, 2017), and soybean (Food Data Central, 2019b). The reference caloric value of Limnospira (commercially available as Spirulina) is estimated to be 3,000 kcal/kg DW (Food Data Central, 2019a). Solid output of the human metabolism is described as feces and urea. All other body wastes such as dead skin cells, sweat solids, hair, nails, etc. (Hu et al, 2010) are assumed to be part of the feces. H2O is only integrated in the stoichiometry as a metabolite, and as such the amounts do not represent actual water consumption. It is assumed that the anabolism and catabolism of the crew are in dynamic equilibrium, and that the crew neither gains or loses weight.

      Assumptions about the human metabolism. Listed values are for one crew member.

      Parameter Value References and notes
      Daily caloric intake 3,000 kcal Poughon et al, 1994; Ewert et al, 2022
      Daily caloric intake from higher plants 2,700 kcal 90%
      Daily caloric intake from Limnospira 300 kcal 10%
      Oxygen consumption 0.84 kg/day or 52.50 mol/day Jones, 2003
      CO2 production 1 kg/day or 22.72 mol/day Jones, 2003
      Urea production 0.1239–0.7758 mol/day Based on Putnam, 1971; MedlinePlus, 2023.
      Feces production Stoichiometrically determined
      4 Results

      This section describes the results of the stoichiometry development, and the simulation results.

      4.1 Stoichiometric equations

      For each compartment a set of chemical equations with fixed stoichiometric coefficients has been developed, making sure that all types of output are also used as input in the loop, thus enabling full closure. The equations are partially based on existing equations from the MELiSSA literature, but were tailored to serve the needs of this study based on additional literature.

      4.1.1 Compartment 1: fermentation with thermophilic anaerobic bacteria

      The first compartment of the MELiSSA loop is an anaerobic digester with a consortium of thermophilic bacteria that receives fecal material and other wastes, such as inedible plant material (Michel et al, 2005; Mastroleo, 2009; Poughon et al, 2013). Anaerobic digestion is a natural process that typically consists of four subsequent stages: hydrolysis, acidogenesis, acetogenesis and methanogenesis. Volatile fatty acids are intermediate products in this process, while methane is the main end product (Hendrickx and Mergeay, 2007; Khan et al, 2016). In the MELiSSA system a partial anaerobic transformation is used with a maximal yield of volatile fatty acids (VFAs), all while avoiding the production of methane. The VFAs are useful biomolecules that can be utilized as nutrients for other microorganisms downstream in the MELiSSA loop. Such partial anaerobic digestion is achieved by operating the bioreactor under slightly acidified conditions (pH 5.3) which inhibits methanogenic activity (Mastroleo, 2009; Clauwaert et al, 2017).

      Organic waste typically consists of polymers (carbohydrates and proteins), as well as fats (Liebetrau et al, 2017). These cannot be taken up by the microbial cells and have to be broken down into smaller compounds first. During hydrolysis, bacteria excrete hydrolytic enzymes which split these substances into sugars, amino acids, fatty acids, and glycerol. These water-soluble hydrolysis products are then converted into volatile fatty acids, CO2 and hydrogen during anaerobic fermentation (or acidogenesis) (Lauwers et al, 2013; Liebetrau et al, 2017). C1 receives input from three different sources: feces from the crew, inedible plant material, and biomass from C2. The next paragraphs describe the stoichiometries of the fermentation of, respectively, carbohydrates, proteins, and lipids.

      4.1.1.1 Carbohydrate fermentation

      Carbohydrate fermentation in the MELiSSA loop has been described in multiple studies, and results in the production of VFAs, CO2 and hydrogen. In Hendrickx et al (2006) only one generic formula representing all VFAs is used, while in all the other consulted MELiSSA studies just two VFAs are produced (e.g., Dussap et al, 1993; Poughon et al, 2009). In some studies this is complemented with the production of bacterial biomass (Guirado & Podhajsky, 2008; Thiron, 2020). As stated above, in the current stoichiometry, four VFAs are used and no biomass production is needed in the equation of this compartment. As a consequence, a new stoichiometry had to be calculated. Proportions of acetic acid, propionic acid, butyric acid and valeric acid were used from (Yin et al, 2016), describing the anaerobic fermentation of glucose as model carbohydrate at pH 6. In below equation, the proportions of these four VFAs were used since the MELiSSA C1 bioreactor also operates at a low pH of 5.3 (Clauwaert et al, 2017). CO2 and hydrogen are produced in equal proportions, in line with other published balanced carbohydrate fermentation reactions (Dussap et al, 1993; Yin et al, 2016). C H 1.6667 O 0.8333 + 0.1666 H 2 O 0.04651 C 2 H 4 O 2 a c e t i c a c i d + 0.0711 C 3 H 6 O 2 p r o p i o n i c a c i d + 0.09535 C 4 H 8 O 2 b u t y r i c a c i d + 0.006313 C 5 H 10 O 2 v a l e r i c a c i d + 0.2807 C O 2 + 0.2807 H 2

      4.1.1.2 Protein fermentation

      Just like with carbohydrates, there are several stoichiometric descriptions of the fermentation of proteins in the MELiSSA literature. Older studies use acetic acid and butyric acid as VFA output (e.g., Dussap & Gros, 1991), subsequent studies use up to 5 VFAs (e.g., Guirado & Podhajsky, 2008; Thiron, 2020). In these latter studies, microbial biomass is also included as an output. As biomass output is not necessary in this part of our stoichiometric model, and as a slightly differing empirical formula for protein is used, a new stoichiometric balance needs to be calculated. Proteins are polymers consisting of amino acids. During anaerobic fermentation, proteins are first broken down into their individual amino acids. The subsequent fermentation of these amino acids can be described using a single reaction step, and results in the production of VFAs, CO2, hydrogen and ammonia (Ramsay and Pullammanappallil, 2001; Poughon et al, 2009). In stoichiometric equations, CO2, hydrogen and ammonia are generally being produced in a similar order of magnitude (Poughon, 2007b; Thiron, 2020). Tepari (2019) investigated protein fermentation during anaerobic wastewater treatment, and used bovine serum albumin (BSA) as model protein. At pH 5, only acetic acid (37.4%), propionic acid (13.7%) and butyric acid (48.9%) were produced, and no valeric acid. These proportions were used to calculate the following stoichiometric equation: C H 1.59 O 0.31 N 0.25 + 0.4962 H 2 O 0.1056 C 2 H 4 O 2 a c e t i c a c i d + 0.03867 C 3 H 6 O 2 p r o p i o n i c a c i d + 0.138 C 4 H 8 O 2 b u t y r i c a c i d + 0.1209 C O 2 + 0.037 H 2 + 0.25 N H 3

      4.1.1.3 Lipid fermentation

      As described above, tripalmitin is used as model lipid in this model. Since in previous MELiSSA modeling studies, palmitic acid is used, a new stoichiometric equation has to be calculated. When lipids such as triacylglycerols are anaerobically digested, they are first hydrolyzed to glycerol and long-chain fatty acids (LCFAs) (Mackie et al, 1991; Li et al, 2002; Cirne et al, 2007). Subsequently, the glycerol is degraded into VFAs, while the LCFAs are broken down into acetic acid and hydrogen through beta-oxidation (Weng and Jeris, 1976; Ceron-Chafla et al, 2021; Holohan et al, 2022). Firstly, the glycerol and LCFA fermentations will be individually described, after which they will be integrated into one equation.

      The hydrolysis of tripalmitin can be written as follows: C 51 H 98 O 6 t r i p a l m i t i n + 3 H 2 O C 3 H 8 O 3 g l y c e r o l + 3 C 16 H 32 O 2 p a l m i t i c a c i d

      In experiments by Yin et al, 2016, glycerol was subjected to (non-strict) anaerobic fermentation at pH 6. The study presents three individual stoichiometric reactions, each describing the degradation of glycerol into acetic acid, propionic acid, and butyric acid. The experiment resulted in a production of 32% acetic acid, 52% propionic acid, and 16% butyric acid. Integrating the three individual stoichiometric reactions, taking into account the above mentioned proportion, leads to the following single equation: C 3 H 8 O 3 g l y c e r o l + 0.2727 H 2 O 0.2727 C 2 H 4 O 2 a c e t i c a c i d + 0.4481 C 3 H 6 O 2 p r o p i o n i c a c i d + 0.1396 C 4 H 8 O 2 b u t y r i c a c i d + 0.5519 C O 2 + 1.3766 H 2 + 0.4481 H 2 O

      The stoichiometric equation for palmitic beta-oxidation is described by Li et al (2018): C 16 H 32 O 2 p a l m i t i c a c i d + 14 H 2 O 8 C 2 H 4 O 2 a c e t i c a c i d + 14 H 2

      Integrating Eqs 35 leads to: C H 1.9216 O 0.1177 + 0.8789 H 2 O 0.4759 C 2 H 4 O 2 a c e t i c a c i d + 0.008786 C 3 H 6 O 2 p r o p i o n i c a c i d + 0.002737 C 4 H 8 O 2 b u t y r i c a c i d + 0.01082 C O 2 + 0.8505 H 2

      4.1.2 Compartment 2: VFA processing by photoheterotrophic bacteria

      The purpose of MELiSSA C2 is to recycle the VFAs produced in C1. The purple non-sulphur phototrophic bacteria Rhodospirillum rubrum is utilized for this. Under anaerobic conditions, it uses light as an energy source to convert VFAs into biomass (Mastroleo, 2009; De Meur et al, 2020). Because we use a custom composition of Rhodospirillum biomass in this model, a new stoichiometry had to be calculated. Anaerobic digestion of organic matter may lead to fluctuating proportions of different VFAs (Alloul et al, 2018). Because of the variable VFA spectrum coming out of C1, it was decided to develop a individual equation for each VFA. The assimilation process consumes CO2 and NH3, and for acetic acid also hydrogen (Dussap et al, 1993; Fulget, 1996; Thiron, 2020). As explained above, no net CO2 production occurs in the following equations. C 2 H 4 O 2 a c e t i c a c i d + H 2 + 0.2774 C O 2 + 0.4073 N H 3 0.3327 C H 1.6667 O 0.8333 c a r b o h y d r a t e s + 1.6291 C H 1.5900 O 0.3100 N 0.2500 p r o t e i n s + 0.3155 C H 1,9216 O 0,1177 l i p i d s R h o d o s p i r i l l u m + 1.7353 H 2 O 0.8556 C 3 H 6 O 2 p r o p i o n i c a c i d + 0.1611 C O 2 + 0.4879 N H 3 0.3986 C H 1.6667 O 0.8333 c a r b o h y d r a t e s + 1.9515 C H 1.5900 O 0.3100 p r o t e i n s N 0.2500 + 0.3780 C H 1,9216 O 0,1177 l i p i d s R h o d o s p i r i l l u m + 1.052 H 2 O 0.5563 C 4 H 8 O 2 b u t y r i c a c i d + 0.3086 C O 2 + 0.4532 N H 3 0.3702 C H 1.6667 O 0.8333 c a r b o h y d r a t e s 1.8127 C H 1.5900 O 0.3100 p r o t e i n s N 0.2500 + 0.3511 C H 1,9216 O 0,1177 l i p i d s R h o d o s p i r i l l u m + 0.8182 H 2 O 0.4084 C 5 H 10 O 2 v a l e r i c a c i d + 0.3762 C O 2 + 0.4325 N H 3 0.3533 C H 1.6667 O 0.8333 c a r b o h y d r a t e s + 1.7299 C H 1.5900 O 0.3100 N 0.2500 p r o t e i n s + 0.335 C H 1,9216 O 0,1177 l i p i d s R h o d o s p i r i l l u m + 0.6991 H 2 O

      4.1.3 Compartment 3: nitrification of ammonia and ureolysis

      In some cases, ammonia and urea can be taken up directly by plants or microorganisms grown for food production. However, in a BLSS it is desirable to convert ammonia and urea to nitrate because it is a less volatile and less reactive molecule (Clauwaert et al, 2017; De Paepe et al, 2018; Sachdeva et al, 2021). MELiSSA C3 receives ammonia from C2 and urine from C5, and transforms both into nitrate for the plants and microalgae in C4 (Duatis and Moreno, 2009; 2010). Ammonia is processed in a nitrifying bioreactor, while urine is processed separately in a modified version of the nitrifying bioreactor, the so-called Urine Treatment Unit (UTU) (Guirado and Podhajsky, 2008). As both bioreactors process different compounds, two separate compartments are defined in the stoichiometric model: 3a and 3b.

      4.1.3.1 Compartment 3a: nitrification of ammonia by nitrifying bacteria

      The nitrifying bioreactor as defined in the MELiSSA project is a cylindrical packed-bed bioreactor with polystyrene beads that are colonized by a biofilm of Nitrosomonas europaea and Nitrobacter winogradsky (Cruvellier et al, 2016; Garcia-Gragera et al, 2021). The nitrification of ammonia happens in two steps, each step carried out by one of the bacterial species (Cruvellier et al, 2016): 2 N H 3 + 3 O 2 2 H N O 2 + 2 H 2 O N i t r o s o m o n a s 2 H N O 2 + O 2 2 H N O 3 N i t r o b a c t e r

      Taken together this results in: N H 3 + 2 O 2 H N O 3 + H 2 O

      4.1.3.2 Compartment 3b: breakdown of urea by ureolytic heterotrophic bacteria

      In a recent MELiSSA study of a UTU concept (Christiaens et al, 2019), showed how human urine can be nitrified with the use of a synthetic microbial community of ureolytic heterotrophic bacteria and nitrifying bacteria cultivated in a continuous stirred tank reactor. First, during ammonification, the ureolytic heterotrophic bacteria break down urea into ammonia and CO2, a process catalyzed by the enzyme urease (Nicolau et al, 2010; De Paepe et al, 2018). Subsequently, the ammonia is further converted into nitrate by the nitrifying bacteria as described above. The stoichiometric equation describing the ammonification of urea can be found in (Guirado and Podhajsky, 2008): C O N H 2 2 u r e a + H 2 O 2 N H 3 + C O 2

      Taken together with (13), this results in: C O N H 2 2 u r e a + 4 O 2 2 H N O 3 + H 2 O + C O 2

      4.1.4 Compartment 4: food production 4.1.4.1 Compartment 4a: photoautotrophic cyanobacteria growth

      The edible microalgae Limnospira platensis and indica (previously Arthrospira platensis and indica) are cyanobacteria that consume nitrate and CO2 that are produced in other compartments of the loop, and generate oxygen and food for the crew (Clauwaert et al, 2017; Alemany et al, 2019; Sachdeva et al, 2021). Limnospira is grown inside a gas lift photobioreactor with illumination and a gas injection system (Garcia-Gragera et al, 2021; Poughon et al, 2021). Since a custom biomass composition is used in our model, a new stoichiometric balance had to be calculated (generated biomass is dry mass): C O 2 + 0.1788 H N O 3 + 0.7342 H 2 O 0.1461 C H 1.6667 O 0.8333 c a r b o h y d r a t e s + 0.7154 C H 1.5900 O 0.3100 N 0.2500 p r o t e i n s + 0.1385 C H 1,9216 O 0,1177 l i p i d s L i m n o s p i r a + 1.4554 O 2

      4.1.4.2 Compartment 4b: higher plant growth

      Edible crops are grown in C4b of the MELiSSA loop (Pannico et al, 2022). The goal of this compartment is in line with that of C4a: recycling nitrate and CO2, and producing oxygen and food (Lasseur et al, 2010; Peiro et al, 2020). In the current model, a harvest index of 50% is used, which means that half of the plant biomass is edible, while the other half is not. The non-edible part is considered as waste and is used as feed for the thermophilic anaerobic bacteria in C1. This can be summarized in the general equation: C O 2 + H 2 O + H N O 3 b i o m a s s h i g h e r p l a n t s e d i b l e + 0.8965 b i o m a s s h i g h e r p l a n t s n o n e d i b l e + O 2

      This corresponds to the following stoichiometric equation (generated biomass is dry mass): 2.9101 C O 2 + 0.0842 H N O 3 + 2.4 H 2 O 0.9629 C H 1.6667 O 0.8333 c a r b o h y d r a t e s + 0.3368 C H 1.5900 O 0.3100 N 0.2500 p r o t e i n s + 0.2348 C H 1,9216 O 0,1177 l i p i d s h i g h e r p l a n t s e d i b l e + 1.3756 C H 1.6667 O 0.8333 c a r b o h y d r a t e s h i g h e r p l a n t s n o n - e d i b l e + 3.1961 O 2

      4.1.5 Compartment 5: human metabolism

      The general equation for the human metabolism can be summarized as follows (Guirado and Podhajsky, 2008): b i o m a s s h i g h e r p l a n t s e d i b l e + b i o m a s s L i m n o s p i r a + O 2 f e c e s + C O N H 2 2 + H 2 O + C O 2

      To determine the final stoichiometric equation, a few preparatory calculations are needed. Firstly, the right proportion between edible higher plant biomass and Limnospira biomass needs to be established. In the model, 90% of the caloric content of the food comes from the higher plants, and 10% comes from Limnospira. Higher plants and Limnospira are considered to have a caloric content of respectively 4,000 and 3000 kcal/kg. This amounts to a target daily food intake of 675 g DW (27.86 mol) of higher plant biomass and 100 g DW (4.56 mol) of Limnospira biomass per crew member. Using the biomass compositions, this can then be translated into corresponding values for carbohydrates, proteins, and lipids. For oxygen, CO2, and urea, reported daily values were transformed into moles (Table 2), and this was then used for calculating a stoichiometry with realistic proportions between these compounds and biomass uptake. The amount of feces is calculated by solving the stoichiometric equation. Below equation describes the metabolism of one crew member: 17.6580 C H 1.6667 O 0.8333 c a r b o h y d r a t e s + 6.1760 C H 1.5900 O 0.3100 N 0.2500 p r o t e i n s + 4.3059 C H 1,9216 O 0,1177 l i p i d s h i g h e r p l a n t s e d i b l e + 0.6594 C H 1.6667 O 0.8333 c a r b o h y d r a t e s + 3.2286 C H 1.5900 O 0.3100 N 0.2500 p r o t e i n s + 0.6253 C H 1,9216 O 0,1177 l i p i d s L i m n o s p i r a + 24.2893 O 2 6.1661 C H 1.6667 O 0.8333 c a r b o h y d r a t e s + 3.1982 C H 1.5900 O 0.3100 N 0.2500 p r o t e i n s + 0.9989 C H 1,9216 O 0,1177 l i p i d s f e c e s + 0.7758 C O N H 2 2 u r e a + 17.2869 H 2 O + 21.5142 C O 2

      4.1.6 Auxiliary processes 4.1.6.1 Excess hydrogen burn-off

      One auxiliary process has been added to the loop with the function to burn any excess hydrogen. This is in correspondence with other MELiSSA modeling studies such as (Poughon et al, 2000). 2 H 2 + O 2 2 H 2 O

      4.2 Simulation results

      The results of the stoichiometric spreadsheet model running in steady state conditions can be found in Table 4. The overall input and output values of all compounds in the loop are listed, for a crew of six over a period of 1 day (the time unit used in the ABM for which this stoichiometry was developed). All different stoichiometric equations were successfully balanced with each other, which resulted in a very high closure level with 12 out of 14 values reaching full closure, with no material loss. Oxygen displays an almost neglectable loss of 0.00021%, and CO2 0.00001%. This is assumed to be caused by the numerical limitations of the software tools (e.g., number of decimals and rounding). The values correspond to respectively 2.13 mg O2 and 1.47 mg CO2 material loss per person per day. Over the course of a year, this would amount to 4.67 g O2 and 3.21 g CO2 for the entire crew. The VFA spectrum (weight-based) at steady state consists of 57.20% acetic acid, 12.59% propionic acid, 29.98% butyric acid and 0.23% valeric acid.

      Input and output values of the different compounds, after reaching steady state. Values for a crew of six, over a period of 1 day.

      Compound Input (mol) Output (mol) Material loss(mol) Material loss(%)
      Carbohydrates 1832.2836 1832.2836 0.0000 0.00000%
      Protein 7587.0675 7587.0675 0.0000 0.00000%
      Lipids 1490.3754 1490.3754 0.0000 0.00000%
      Acetic acid 1570.2506 1570.2506 0.0000 0.00000%
      Propionic acid 426.4730 426.4730 0.0000 0.00000%
      Butyric acid 1207.5497 1207.5497 0.0000 0.00000%
      Valeric acid 10.8735 10.8735 0.0000 0.00000%
      Urea 4.6548 4.6548 0.0000 0.00000%
      NH3 1882.6600 1882.6601 0.0001 0.00000%
      HNO3 14.1068 14.1068 0.0000 0.00000%
      H2O 5591.1996 5591.1994 -0.0001 0.00000%
      H2 2004.4967 2004.4967 0.0000 0.00000%
      O2 391.0724 391.0716 -0.0008 -0.00021%
      CO2 1543.0613 1543.0611 -0.0002 -0.00001%
      5 Discussion

      This study presents a stoichiometric model of the MELiSSA BLSS, which describes the main metabolic processes in all compartments using a series of static chemical equations. In steady state, 100% of the necessary food and oxygen for the crew, as well as 100% of the compounds required for all other organisms in the loop, are generated. Using the assumption that all waste compounds in the loop can be broken down again, full closure was achieved in the static model by balancing the size of each compartment relative to all the other compartments. This is the first time a fully closed static MELiSSA model is described that can provide all the food and oxygen for astronauts. The model is a useful tool for theoretical research on autonomous BLSSs for long-duration space exploration, allowing tracing of the elements over the entire cycle.

      However, in practice, no system can be 100% closed. There are always potential outfluxes from the system, for example, through a buildup of recalcitrant organic matter (Hendrickx and Mergeay, 2007; Lasseur et al, 2010; Zhang et al, 2018) or the accumulation of different precipitates such as carbonates and phosphates (De Paepe et al, 2018; Christiaens et al, 2019). Recalcitrant organic matter consists mainly of plant fibers that do not rapidly biodegrade. Plant cell walls are composed of cellulose, hemicellulose, pectin, and lignin and are among the least degradable polymers in a BLSS (Hendrickx and Mergeay, 2007; Zhang et al, 2018). As a result, up to about 70% of the organic waste could be digested in MELiSSA lab experiments (Farges et al, 2008; Lasseur et al, 2010). The MELiSSA community has investigated additional methods to break down the remaining recalcitrant matter using hydrothermal and chemical oxidation, and anaerobic and hyperthermophilic cellulose-degrading bacteria (Hendrickx & Mergeay, 2007; Lasseur et al, 2010). The most efficient results were obtained by using supercritical water oxidation (SCWO) in which up to 98% of all organic matter could be broken down (Hendrickx and Mergeay, 2007; Zhang et al, 2018). Hydrogen peroxide was used as an oxidizer in these tests, which is a compound that can be generated within a BLSS loop (Tikhomirov et al, 2011; Vijapur et al, 2017; Nelson et al, 2020). Nevertheless, it is safe to assume that, in reality, there will always be some level of material loss in any BLSS loop. It would therefore be interesting to investigate what accumulates over a long period of time and then calculate the dimensions of the reserve storage required to mitigate the resulting material shortages.

      The interaction between compartments C1 and C2 represents a subcycle in this model. The Rhodospirillum biomass generated in C2 is fed back into C1 as organic waste. The goal of this subcycle is to produce a surplus of CO2, H2, and NH3 that is needed in the other compartments of the MELiSSA loop. Steady state of this subcycle is reached gradually until enough biomass from C2 is sent back to C1 to ensure that the exact amount of VFAs is produced to sustain a stable Rhodospirillum culture. If too few VFAs are provided, C2 will underperform due to a lack of Rhodospirillum growth. If too many VFAs are provided, the Rhodospirillum culture will not be able to process all of them, leading to accumulation and negative impacts on the mass flows and efficiency of the entire loop. At steady state, this subcycle generates the following surplus (per day, for a crew of six): 213 mol CO2, 434 mol H2, and 126 mol NH3. The subcycle also explains the dominant amount of protein present in the loop, despite the food provided to the crew containing much less protein (Table 4). Due to the need for a significant mass of Rhodospirillum to process the VFAs and their high protein content (72%), a lot of protein is present in the loop in general. Bio-electrochemical oxidation is a potentially more efficient approach to convert the VFAs into CO2. In a study by Luther et al, 2018, using bio-electrochemical oxidation within a microbial electrolysis cell, 80%–100% of the VFAs in the effluent of MELiSSA’s C1 were converted into CO2 in 7 days. With further process development and optimization, bio-electrochemical oxidation can potentially entirely replace the use of Rhodospirillum in C2.

      (Jiang et al, 2013) conducted a series of experiments on anaerobic digestion of food waste. At a pH of 5.0, the resulting VFA spectrum consisted of 60.40% acetic acid, 8.32% propionic acid, 31.13% butyric acid, and 0.15% valeric acid, which falls within the same range as the results achieved in our model at steady state. Another study by Lim et al, 2008 also involved anaerobic digestion of food waste. At a pH of 5.5, the resulting VFA spectrum had a more variable composition, but acetic acid was generally the dominant VFA. These reported results are consistent with the results of our static model, which suggests that the theoretical VFA outcome is a potentially realistic scenario. However, it must be noted that VFA spectra generated during anaerobic fermentation are highly variable and depend on various factors such as substrate composition and inoculum (Poughon et al, 2013; Khan et al, 2016), pH and temperature (Shin et al, 2004; Jiang et al, 2013; Khan et al, 2016), and bioreactor design and operation (Lim et al, 2008; Bharathiraja et al, 2016). For instance, in experiments using food waste by Shin et al (2004), and in MELiSSA experiments using a synthetically composed organic waste (Luther et al, 2018; Zhang et al, 2019), both acetic acid and butyric acid were the main VFAs, with butyric acid being dominant.

      Growing enough plants to meet the food requirements of a human crew typically results in the production of more oxygen than the crew needs. This is because the harvest index of the crops is often less than 100%, meaning that only a portion of the grown biomass is used for food (Jones, 2003). However, the unused biomass still generates oxygen. In a closed loop system, this could lead to an accumulation of excess oxygen (as e.g., observed in Poughon et al, 2000; Gros et al, 2003) and, eventually, the collapse of the system. In our model, the excess oxygen is used to hydrolyze solid waste in C1, in the form of H2O. The hydrogen-burning auxiliary helps to transform any remaining oxygen gas into H2O. Since oxygen originates from splitting H2O during photosynthesis, the amount of surplus oxygen matches the amount of surplus hydrogen in the system.

      In our model, the human crew consumes 21.51 mol/person/day of CO2, which is very close to the target value of 22.72 mol/person/day used for calculating the stoichiometry. However, the human crew only uses 24.29 mol/person/day of oxygen, much lower than the target value of 52.50 mol/person/day. The proportion between consumed oxygen and produced CO2 should be 2.31, but in our model, this proportion is 1.13, precisely because of the lower oxygen consumption. Remarkably, other BLSS modeling studies show the same tendency. In MELiSSA studies, the proportion is 1.20 (Fulget, 1996; Thiron, 2020) and in Lunar Palace studies, it ranges from 1.13 to 1.38 (Hu et al, 2010; Fu et al, 2016). This suggests that the description of the O2 consumption in all of these studies is not entirely complete. This could be due to the fact that in reality, more compounds are oxidized than represented in the simplified stoichiometric equation for human metabolism. The total mass of consumed food per person per day corresponds to 781 g DW/day. This is slightly higher than the range of baseline values used in MELiSSA studies and reference documents. Adjusted to a daily food intake of 3,000 kcal/person, the values in these studies and documents range from 509 to 748 g DW/day (e.g. (Dussap et al, 1993; Poughon et al, 2000; Thiron, 2020)).

      The composition of feces was calculated by solving the stoichiometric equation for human metabolism, using the general characterization of feces reported in the literature as a starting point (Rose et al, 2015). After solving the equation for human metabolism, the resulting feces composition was found to be 66% carbohydrates, 28% proteins, and 6% lipids, corresponding to the empirical formula CH1.6676O0.6028N0.07715. This formula is similar to the one used in BLSS studies by Hu et al, 2010 and Fu et al, 2016, which is CH1.70O0.60N0.05. Because the MELiSSA diet is a vegan diet, it is expected that there is a high fiber intake and a resulting high amount of fiber in the feces (Brodribb et al, 1980; Kay, 1982; Forsum et al, 1990). As a result, the overall carbohydrate fraction could be quite significant. It should be noted that the fraction of protein in the feces could be an underestimation, as bacteria can compose up to half of the fecal solids and may have a significant protein content (Cummings, 2001; Rose et al, 2015). Intestinal bacterial growth is not included in the human stoichiometry, which is justifiable as it means more digestion is done by the same groups of microbes in C1. The amount of feces produced is 253 g DW/person/day. Jenkins et al, 2001 measured 160 ± 24 g DW/person/day for a vegetarian diet. It should be noted that no other solid body wastes have been included in the stoichiometry, and as such the ‘feces’ in the equation represent a more general representation of all wastes combined.

      Further development of the stoichiometric model could focus on increasing its granularity. This could involve expanding the model to include more elements beyond CHON, incorporating a greater diversity of crops, and using a more detailed description of human metabolism. Nitrogen, phosphorus, and potassium are the three major macro-elements for plants (Usherwood and Segars, 2001; Kanazawa et al, 2008). Therefore, including the latter two elements in the next version of the model seems like a logical step. Fertilizers typically contain these three elements, but it is crucial that they are applied in the correct ratio (Usherwood and Segars, 2001; Chun et al, 2017), which must also be factored into the model. Diverse plant crops could be used instead of the concept of an ‘ideal plant’, allowing for a more detailed description of the crew’s daily diet. Integrating water usage beyond it's role as a metabolite will result in a more accurate description of the entire water budget. For example, water transpiration, water excretion, and free water inside the plant could be factored in. By integrating solid body wastes other than feces, the stoichiometry of human metabolism could also be further refined.

      Data availability statement

      The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

      Author contributions

      AV and AP conceived and designed the study and performed the MELiSSA literature review. AV did the background research for the stoichiometry, while AP developed the tool in Excel to balance the chemical equations. The stoichiometric spreadsheet model was created in Excel by AP. AV wrote the first draft of the manuscript, and AP contributed sections. All authors contributed to the article and approved the submitted version.

      The first, and shorter, version of the stoichiometry was compiled in collaboration with Jason Kiem from SmartCrops BV. Prof. Siegfried Vlaeminck from the University of Antwerp provided validation of this initial version and offered suggestions for further development. The authors would like to express their gratitude to Christophe Lasseur and Christel Paille from MELiSSA at ESTEC, as well as Prof. Claude-Gilles Dussap from Polytech Clermont, for providing input and feedback throughout the research. Ralph Lindeboom from Delft University of Technology also provided valuable feedback on a previous version of the manuscript.

      Conflict of interest

      The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

      Publisher’s note

      All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

      References Albiol J. Gòdia F. Montesinos J. Perez J. Vernerey A. Cabello F. (2000). Biological life support system demostration facility: the melissa pilot plant. SAE Tech. Pap. 10.4271/2000-01-2379 Alemany L. Peiro E. Arnau C. Garcia D. Poughon L. Cornet J.-F. (2019). Continuous controlled long-term operation and modeling of a closed loop connecting an air-lift photobioreactor and an animal compartment for the development of a life support system. Biochem. Eng. J. 151, 107323. 10.1016/j.bej.2019.107323 Alloul A. Ganigué R. Spiller M. Meerburg F. Cagnetta C. Rabaey K. (2018). Capture–ferment–upgrade: a three-step approach for the valorization of sewage organics as commodities. Environ. Sci. Technol. 52, 67296742. 10.1021/acs.est.7b05712 Audas C. Ortega Ugalde S. Paille C. Lamaze B. Lasseur C. (2022). Life support systems beyond low Earth orbit advocates for an improved resources management approach. Ecol. Eng. Environ. Prot. 2022, 513. 10.32006/eeep.2022.1.0513 Bartsev S. I. Gitelson J. I. Lisovsky G. M. Mezhevikin V. V. Okhonin V. A. (1996). Perspectives of different type biological life support systems (BLSS) usage in space missions. Acta Astronaut. 39, 617622. 10.1016/S0094-5765(97)00012-X Begon M. Townsend C. R. (2021). Ecology: from individuals to ecosystems. New Jersey, United States: John Wiley & Sons. Bharathiraja B. Sudharsanaa T. Bharghavi A. Jayamuthunagai J. Praveenkumar R. (2016). Biohydrogen and Biogas – an overview on feedstocks and enhancement process. Fuel 185, 810828. 10.1016/j.fuel.2016.08.030 Bjørnstad O. N. (2015). Nonlinearity and chaos in ecological dynamics revisited. Proc. Natl. Acad. Sci. 112, 62526253. 10.1073/pnas.1507708112 Brodribb J. Condon R. E. Cowles V. DeCosse J. J. (1980). Influence of dietary fiber on transit time, fecal composition, and myoelectrical activity of the primate right colon. Dig. Dis. Sci. 25, 260266. 10.1007/BF01308515 Carillo P. Morrone B. Fusco G. M. De Pascale S. Rouphael Y. (2020). Challenges for a sustainable food production system on board of the international space station: a technical review. Agronomy 10, 687. 10.3390/agronomy10050687 Ceron-Chafla P. Chang Y. Rabaey K. van Lier J. B. Lindeboom R. E. F. (2021). Directional selection of microbial community reduces propionate accumulation in glycerol and glucose anaerobic bioconversion under elevated pCO2 . Front. Microbiol. 12, 675763. 10.3389/fmicb.2021.675763 Christiaens M. E. R. De Paepe J. Ilgrande C. De Vrieze J. Barys J. Teirlinck P. (2019). Urine nitrification with a synthetic microbial community. Syst. Appl. Microbiol. 42, 126021. 10.1016/j.syapm.2019.126021 Chun J.-H. Kim S. Arasu M. V. Al-Dhabi N. A. Chung D. Y. Kim S.-J. (2017). Combined effect of Nitrogen, Phosphorus and Potassium fertilizers on the contents of glucosinolates in rocket salad (Eruca sativa Mill). Saudi J. Biol. Sci. 24, 436443. 10.1016/j.sjbs.2015.08.012 Cirne D. G. Paloumet X. Björnsson L. Alves M. M. Mattiasson B. (2007). Anaerobic digestion of lipid-rich waste—effects of lipid concentration. Renew. Energy 32, 965975. 10.1016/j.renene.2006.04.003 Ciurans Molist C. Bazmohammadi N. Vasquez J. C. Dussap C.-G. Guerrero J. Gòdia F. (2020). Hierarchical control of space closed ecosystems - expanding microgrid concepts to bioastronautics. IEEE Ind. Electron. Mag. 15, 1627. 10.1109/MIE.2020.3026828 Clauwaert P. Muys M. Alloul A. De Paepe J. Luther A. Sun X. (2017). Nitrogen cycling in Bioregenerative Life Support Systems: challenges for waste refinery and food production processes. Prog. Aerosp. Sci. 91, 8798. 10.1016/j.paerosci.2017.04.002 Cooper M. Perchonok M. Douglas G. L. (2017). Initial assessment of the nutritional quality of the space food system over three years of ambient storage. npj Microgravity 3, 17. 10.1038/s41526-017-0022-z Cornet J.-F. Dussap C. G. Gros J.-B. (1998). “Kinetics and energetics of photosynthetic micro-organisms in photobioreactors,” in Bioprocess and algae reactor Technology, apoptosis advances in biochemical engineering biotechnology (Berlin, Heidelberg: Springer), 153224. 10.1007/BFb0102299 Cornet J.-F. Dussap C. G. Leclercq J.-J. (2001). “Simulation, design and model based predictive control of photobioreactors,” in Engineering and Manufacturing for biotechnology focus on biotechnology. Editors Hofman M. Thonart P. (Dordrecht: Springer Netherlands), 227238. 10.1007/0-306-46889-1_15 Cornet J.-F. Favier L. Dussap C.-G. (2003). Modeling stability of photoheterotrophic continuous cultures in photobioreactors. Biotechnol. Prog. 19, 12161227. 10.1021/bp034041l Cornet J. F. Dussap C. G. (2000). Kinetic and stoichiometric analysis of Rhodospirillum rubrum growth under carbon substrate limitations in rectangular photobioreactors. MELiSSA Technical Note 45.5. Laboratoire de Génie Chimique et Biochimique Available at: https://www.melissafoundation.org/download/469 (Accessed March 23, 2023). Cruvellier N. Poughon L. Creuly C. Dussap C.-G. Lasseur C. (2016). Growth modelling of Nitrosomonas europaea ATCC® 19718 and nitrobacter winogradskyi ATCC® 25391: a new online indicator of the partial nitrification. Bioresour. Technol. 220, 369377. 10.1016/j.biortech.2016.08.090 Cummings J. (2001). “The effect of dietary fiber on fecal weight and composition,” in CRC handbook of dietary fiber in human nutrition. Editor Spiller G. (Florida, United States: CRC Press), 183252. 10.1201/9781420038514.ch4.4 De Meur Q. Deutschbauer A. Koch M. Bayon-Vicente G. Cabecas P. Ruddy W. (2020). New perspectives on butyrate assimilation in Rhodospirillum rubrum S1H under photoheterotrophic conditions. BMC Microbiol. 20, 126. 10.1186/s12866-020-01814-7 De Paepe J. Lindeboom R. E. F. Vanoppen M. De Paepe K. Demey D. Coessens W. (2018). Refinery and concentration of nutrients from urine with electrodialysis enabled by upstream precipitation and nitrification. Water Res. 144, 7686. 10.1016/j.watres.2018.07.016 Deangelis D. Grimm V. (2014). Individual-based models in ecology after four decades. F1000prime Rep. 6, 39. 10.12703/P6-39 Decker E. A. Ferruzzi M. G. (2013). Innovations in food chemistry and processing to enhance the nutrient profile of the white potato in all forms. Adv. Nutr. 4, 345S350S. 10.3945/an.112.003574 Duatis J. Guirado V. Podhajsky S. (2008). MELiSSA adaptation for space. Phase II. Preliminary life support system design. MELiSSA Technical Note 88.3. NTE SA (Werfen Group) Available at: https://www.melissafoundation.org/download/174. Duatis J. Moreno D. (2009). MELiSSA adaptation for space. Phase II. LIfe support system recommended design for a demonstration in a moon base. Barcelona, Spain: NTE-SENER. MELiSSA Technical Note 88.5. Duatis J. Moreno D. (2010). MELiSSA adaptation for space. Phase II. Recommendations and future work. PreliminaryDevelopment plan. MELiSSA technical note 88.6. NTE-SENER Available at: https://www.melissafoundation.org/download/238 (Accessed March 23, 2023). Dubitzky W. Southgate J. Fuß H. (Editors) (2011). “Adaptation and self-organizing systems,” Understanding the dynamics of biological systems (New York, NY: Springer). 10.1007/978-1-4419-7964-3 Dussap C.-G. Cornet J.-F. Gros J.-B. (1993). Simulation of mass fluxes in the MELISSA microorganism based ecosystem. Warrendale, PA SAE Int., 932125. 10.4271/932125 Dussap C. G. Gros J. B. (1991). Simulation des flux de matières dans l’écosystème MELiSSA. MELiSSA Technical Note 14.1. Laboratoire de Génie Chimique et Biochimique Available at: https://www.melissafoundation.org/download/600 (Accessed March 23, 2023). Escobar C. Nabity J. (2017). “Past, present, and future of closed human life support ecosystems-a review,” in Proceedings of the 47th International Conference on Environmental Systems, Charleston, South Carolina, July 2017, 1620. Ewert M. K. Chen T. T. Powell C. D. (2022). Life support baseline values and assumptions document. Available at: https://ntrs.nasa.gov/citations/20210024855 (Accessed April 1, 2023). Farges B. Poughon L. Creuly C. Cornet J.-F. Dussap C.-G. Lasseur C. (2008). Dynamic aspects and controllability of the MELiSSA project: a bioregenerative system to provide life support in space. Appl. Biochem. Biotechnol. 151, 686699. 10.1007/s12010-008-8292-2 Finn C. K. (1998). Steady-state system mass balance for the BIO-plex. Warrendale, PA: SAE International. 10.4271/981747 FoodData Central (2019a). Seaweed, spirulina, dried. Available at: https://fdc.nal.usda.gov/fdc-app.html#/food-details/170495/nutrients (Accessed March 23, 2023). FoodData Central (2019b). Soybeans, mature cooked, boiled, without salt. Available at: https://fdc.nal.usda.gov/fdc-app.html#/food-details/174271/nutrients (Accessed March 23, 2023). Forrest S. Mitchell M. (2016). Adaptive computation: the multidisciplinary legacy of John H. Holland. Commun. ACM 59, 5863. 10.1145/2964342 Forsum E. Eriksson C. Göranzon H. Sohlström A. (1990). Composition of faeces from human subjects consuming diets based on conventional foods containing different kinds and amounts of dietary fibre. Br. J. Nutr. 64, 171186. 10.1079/BJN19900019 Fu Y. Li L. Xie B. Dong C. Wang M. Jia B. (2016). How to establish a bioregenerative life support system for long-term crewed missions to the moon or mars. Astrobiology 16, 925936. 10.1089/ast.2016.1477 Fulget N. (1996). Complete loop control: first study. MELiSSA technical note 28.3. ADERSA. Available at: https://www.melissafoundation.org/download/552 (Accessed March 23, 2023). Fulget N. Poughon L. Richalet J. Lasseur Ch. (1999). Melissa: global control strategy of the artificial ecosystem by using first principles models of the compartments. Adv. Space Res. 24, 397405. 10.1016/S0273-1177(99)00490-1 Garcia-Gragera D. Arnau C. Peiro E. Dussap C.-G. Poughon L. Gerbi O. (2021). Integration of nitrifying, photosynthetic and animal compartments at the MELiSSA pilot plant. Front. Astron. Space Sci. 8, 750616. 10.3389/fspas.2021.750616 Garland J. L. (1989). A simple, mass balance model of carbon flow in a controlled ecological life support system. Available at: https://ntrs.nasa.gov/citations/19900001255 (Accessed March 23, 2023). Gòdia F. Albiol J. Pérez J. Creus N. Cabello F. Montràs A. (2004). The MELISSA pilot plant facility as an integration test-bed for advanced life support systems. Adv. Space Res. 34, 14831493. 10.1016/j.asr.2003.08.038 Gorce B. Garnier L. Ssi Yan Kai H. Schini P.-Y. Lefevre A. Afuera E. (2022). “Conceptual design of an environment control and life support system for a mars transit mission,” in MELiSSA Conference 2022, Toulouse, France, November 2022, 810. Available at: https://www.melissafoundation.org/download/932 (Accessed March 23, 2023). Gros J. B. Poughon L. Lasseur C. Tikhomirov A. A. (2003). Recycling efficiencies of C,H,O,N,S, and P elements in a biological life support system based on micro-organisms and higher plants. Adv. Space Res. 31 (1), 195199. 10.1016/S0273-1177(02)00739-1 Guirado V. Podhajsky S. (2008). MELiSSA adaptation for space. Phase II. Summary of the European life support systems technologies. MELiSSA Technical Note 88.2. NTE SA (Werfen Group) Available at: https://www.melissafoundation.org/download/173 (Accessed March 23, 2023). Häder D. (2020). On the way to mars—flagellated algae in bioregenerative life support systems under microgravity conditions. Front. Plant Sci. 10, 1621. 10.3389/fpls.2019.01621 Hendrickx L. De Wever H. Hermans V. Mastroleo F. Morin N. Wilmotte A. (2006). Microbial ecology of the closed artificial ecosystem MELiSSA (Micro-Ecological Life Support System Alternative): reinventing and compartmentalizing the Earth’s food and oxygen regeneration system for long-haul space exploration missions. Res. Microbiol. 157, 7786. 10.1016/j.resmic.2005.06.014 Hendrickx L. Mergeay M. (2007). From the deep sea to the stars: human life support through minimal communities. Curr. Opin. Microbiol. 10, 231237. 10.1016/j.mib.2007.05.007 Holohan B. C. Duarte M. S. Szabo-Corbacho M. A. Cavaleiro A. J. Salvador A. F. Pereira M. A. (2022). Principles, advances, and perspectives of anaerobic digestion of lipids. Environ. Sci. Technol. Wash. 56, 47494775. 10.1021/acs.est.1c08722 Hu E. Bartsev S. Liu H. (2010). Conceptual design of a bioregenerative life support system containing crops and silkworms. Adv. Space Res. 45, 929939. 10.1016/j.asr.2009.11.022 Imhof B. Weiss P. Vermeulen A. Flynn E. Hyams R. Kerrigan C. (2017). “Space architectures,” in Star Ark: a Living, Self-Sustaining Spaceship . Editor Armstrong R. (Cham: Springer International Publishing), 287340. 10.1007/978-3-319-31042-8_12 Jenkins D. Kendall C. Popovich D. Vidgen E. Mehling C. Vuksan V. (2001). Effect of a very-high-fiber vegetable, fruit, and nut diet on serum lipids and colonic function. Metabolism Clin. Exp. 50, 494503. 10.1053/meta.2001.21037 Jiang J. Zhang Y. Li K. Wang Q. Gong C. Li M. (2013). Volatile fatty acids production from food waste: effects of pH, temperature, and organic loading rate. Bioresour. Technol. 143, 525530. 10.1016/j.biortech.2013.06.025 Johnson C. Boles H. Spencer L. Poulet L. Romeyn M. Bunchek J. (2021). Supplemental food production with plants: a review of nasa research. Front. Astronomy Space Sci. 8. 10.3389/fspas.2021.734343 Jones H. (2003). “Design rules for life support systems,” in 33rd International Conference on Environmental Systems (ICES), Vancouver, B.C., Canada, 7-10 July 2003. Available at: https://ntrs.nasa.gov/citations/20040012725 (Accessed March 23, 2023). Kaleta A. Górnicki K. (2013). Criteria of determination of safe grain storage time - a review. Adv. Agrophysical Res., 295318. Kanazawa S. Ishikawa Y. Tomita-Yokotani K. Hashimoto H. Kitaya Y. Yamashita M. (2008). Space agriculture for habitation on Mars with hyper-thermophilic aerobic composting bacteria. Adv. Space Res. 41, 696700. 10.1016/j.asr.2007.09.040 Kay R. (1982). Dietary fiber. J. Lipid Res. 23, 221242. 10.1016/S0022-2275(20)38151-7 Khan M. A. Ngo H. H. Guo W. S. Liu Y. Nghiem L. D. Hai F. I. (2016). Optimization of process parameters for production of volatile fatty acid, biohydrogen and methane from anaerobic digestion. Bioresour. Technol. 219, 738748. 10.1016/j.biortech.2016.08.073 Knowledge4Policy (2021). Nutritional value of whole grains. Available at: https://knowledge4policy.ec.europa.eu/health-promotion-knowledge-gateway/whole-grain-nutritional-value-whole-2_en (Accessed March 23, 2023). Kobayashi M. Kobayashi M. (1995). “Waste remediation and treatment using anoxygenic phototrophic bacteria,” in Anoxygenic photosynthetic bacteria advances in photosynthesis and respiration. Editors Blankenship R. E. Madigan M. T. Bauer C. E. (Dordrecht: Springer Netherlands), 12691282. 10.1007/0-306-47954-0_62 Lasseur C. Brunet J. De Weever H. Dixon M. Dussap G. Godia F. (2010). MELiSSA: the European project of closed life support system. Gravitational Space Res. 23. http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.464.3630andrep=rep1andtype=pdf Lasseur Ch. Verstraete W. Gros J. B. Dubertret G. Rogalla F. (1996). Melissa: a potential experiment for a precursor mission to the moon. Adv. Space Res. 18, 111117. 10.1016/0273-1177(96)00097-X Lauwers J. Appels L. Thompson I. P. Degrève J. Van Impe J. F. Dewil R. (2013). Mathematical modelling of anaerobic digestion of biomass and waste: power and limitations. Prog. Energy Combust. Sci. 39, 383402. 10.1016/j.pecs.2013.03.003 Li Y. Y. Sasaki H. Yamashita K. Seki K. Kamigochi I. (2002). High-rate methane fermentation of lipid-rich food wastes by a high-solids co-digestion process. Water Sci. Technol. 45, 143150. 10.2166/wst.2002.0420 Liebetrau J. Weinrich S. Sträuber H. Kretzschmar J. (2017). “Anaerobic fermentation of organic material: biological processes and their control parameters,” in Encyclopedia of sustainability science and Technology. Editor Meyers R. A. (New York, NY: Springer), 130. 10.1007/978-1-4939-2493-6_962-1 Lim S.-J. Kim B. J. Jeong C.-M. Choi J. Ahn Y. H. Chang H. N. (2008). Anaerobic organic acid production of food waste in once-a-day feeding and drawing-off bioreactor. Bioresour. Technol. 99, 78667874. 10.1016/j.biortech.2007.06.028 Loader C. Garland J. Raychaudhuri S. Wheeler R. (1997). A simple mass balance model of nitrogen flow in a bioregenerative life support system. Life support & biosphere Sci. Int. J. earth space 4, 3141. Luther A. Beyaert A. Brutsaert M. Lasseur C. Rebeyre P. Clauwaert P. (2018). “Bio-electrochemical oxidation for CO2 recovery in regenerative life support systems,” in 42nd COSPAR Scientific Assembly Abstracts (COSPAR (Committee on Space Research)), California, USA, July 2018. Available at: http://hdl.handle.net/1854/LU-8573095 (Accessed March 23, 2023). Macelroy R. D. Averner M. M. (1978). Space ecosynthesis: an approach to the design of closed ecosystems for use in space. Available at: https://ntrs.nasa.gov/citations/19780018797 (Accessed March 23, 2023). Mackie R. I. White B. A. Bryant M. P. (1991). Lipid metabolism in anaerobic ecosystems. Crit. Rev. Microbiol. 17, 449479. 10.3109/10408419109115208 Maclean H. Dochain D. Waters G. Dixon M. Chaerle L. Van Der Straeten D. (2010). “Identification of simple mass balance models for plant growth - towards food production on manned space missions,” in Proceedings IFAC-CAB 2010 Conference, Leuven, Belgium, July 2010 (International Federation of Automatic Control), 79. 335–340. Available at: https://dial.uclouvain.be/pr/boreal/object/boreal:87876 (Accessed March 23, 2023). Mastroleo F. (2009). Molecular characterization of the life support bacterium Rhodospirillum rubrum S1H cultivated under space related environmental conditions. Université de Mons-Hainaut, Mons, Belgium. Available at: https://researchportal.sckcen.be/files/4517866/Molecular_characterization_of_the_life_support_bacterium_Rhodospirillum_rubrum_S1H_cultivated_under_space_related_environmental_conditions_thesis.Felice.20.01.09.BIS.without.pdf (Accessed March 23, 2023). Mastroleo F. Moussalli C. Raeymaekers L. Smolders C. Leysen L. Coninx I. (2020). “Investigating volatile fatty acids conversion to CO2 by the MELiSSA bacterium Rhodospirillum rubrum,” in MELiSSA Conference 2020, France, November 2020, 35. online. MedlinePlus (2023). Urine 24-hour volume. Available at: https://medlineplus.gov/ency/article/003425.htm (Accessed March 23, 2023). MELiSSA Foundation (2023). Melissa space research program. Available at: https://www.melissafoundation.org/page/melissa-project (Accessed March 23, 2023). MELiSSA (1989). Modelling. MELiSSA technical note 4. Available at: https://www.melissafoundation.org/download/609 (Accessed March 23, 2023). Michel N. De Wever H. Dotremont C. Van Hoof V. (2005). Engineering of the waste compartment. MELiSSA technical note 71.9.4. ECo process assistance. VITO Available at: https://www.melissafoundation.org/download/327 (Accessed March 23, 2023). Molders K. Quinet M. Decat J. Secco B. Dulière E. Pieters S. (2012). Selection and hydroponic growth of potato cultivars for bioregenerative life support systems. Adv. Space Res. 50, 156165. 10.1016/j.asr.2012.03.025 Moore K. J. Hatfield R. D. (1994). “Carbohydrates and forage quality,” in Forage quality, evaluation, and utilization (New Jersey, United States: John Wiley & Sons, Ltd), 229280. 10.2134/1994.foragequality.c6 Muys M. Sui Y. Schwaiger B. Lesueur C. Vandenheuvel D. Vermeir P. (2019). High variability in nutritional value and safety of commercially available Chlorella and Spirulina biomass indicates the need for smart production strategies. Bioresour. Technol. 275, 247257. 10.1016/j.biortech.2018.12.059 Nelson G. J. Vijapur S. H. Hall T. D. Brown B. R. Peña-Duarte A. Cabrera C. R. (2020). Electrochemistry for space life support. Electrochem. Soc. Interface 29, 4752. 10.1149/2.F06201IF Nelson M. Pechurkin N. S. Allen J. P. Somova L. A. Gitelson J. I. (2010). Closed ecological systems, space life support and biospherics. Environ. Biotechnol., 517565. 10.1007/978-1-60327-140-0_11 Nicolau E. Rodríguez-Martínez J. A. Fonseca J. Richardson T.-M. J. Flynn M. Griebenow K. (2010). Bioelectrochemical oxidation of urea with urease and platinized boron doped diamond electrodes for water recycling in space applications. ECS Trans. 33, 18531859. 10.1149/1.3484676 Page V. Feller U. (2013). Selection and hydroponic growth of bread wheat cultivars for bioregenerative life support systems. Adv. Space Res. 52, 536546. 10.1016/j.asr.2013.03.027 Pannico A. Cimini G. Quadri C. Paradiso R. Bucchieri L. Rouphael Y. (2022). A plant characterization unit for closed life support: hardware and control design for atmospheric systems. Front. Astron. Space Sci. 9, 820752. 10.3389/fspas.2022.820752 Paradiso R. Micco V. Buonomo R. Aronne G. Barbieri G. De Pascale S. (2013). Soilless cultivation of soybean for bioregenerative life-support systems: a literature review and the experience of the MELiSSA project - food characterisation phase I. Plant Biol. Stuttg. Ger. 16, 6978. 10.1111/plb.12056 Peiro E. Pannico A. Colleoni S. Bucchieri L. Rouphael Y. De Pascale S. (2020). Air distribution in a fully-closed higher plant growth chamber impacts crop performance of hydroponically-grown lettuce. Front. Plant Sci. 11, 537. 10.3389/fpls.2020.00537 Pilo Teniente S. (2015). Simulation of the gas phase integration between compartments CIVa and CV of the MELiSSA Pilot Plant. Barcelona, Spain: Universitat Politècnica de Catalunya. Available at: https://upcommons.upc.edu/handle/2117/86047 (Accessed March 23, 2023). Poughon L. Creuly C. Farges B. Dussap C.-G. Schiettecatte W. Jovetic S. (2013). Test of an anaerobic prototype reactor coupled with a filtration unit for production of VFAs. Bioresour. Technol. 145, 240247. 10.1016/j.biortech.2012.12.052 Poughon L. Creuly C. Godia F. Leys N. Dussap C.-G. (2021). Photobioreactor Limnospira indica growth model: application from the MELiSSA plant pilot scale to ISS flight experiment. Front. Astronomy Space Sci. 8. 10.3389/fspas.2021.700277 Poughon L. Dussap C. G. Cornet J. F. Gros J. B. (1994). Melissa: behavior of the ecosystem under different light radiant energy inputs. Warrendale, PA SAE Int., 941347. 10.4271/941347 Poughon L. (2007a). Dynamic modelling of a coupled MELiSSA crew - compartment C4a with Matlab-Simulink. MELiSSA Technical Note 83.1. Laboratoire de Génie Chimique et Biochimique Available at: https://www.melissafoundation.org/download/199 (Accessed March 23, 2023). Poughon L. Farges B. Dussap C. G. Godia F. Lasseur C. (2009). Simulation of the MELiSSA closed loop system as a tool to define its integration strategy. Adv. Space Res. 44, 13921403. 10.1016/j.asr.2009.07.021 Poughon L. Gros J. B. Dussap C. G. (2000). “MELISSA loop: first estimate of flow rates and concentrations through the loop,” in 30th International Conference on Environmental Systems (ICES), Toulouse, France, 10-13 July 2000 (SAE Technical Paper Series) (Warrendale, PA: SAE International). Poughon L. Laroche C. Creuly C. Dussap C.-G. Paille C. Lasseur C. (2020). Limnospira indica PCC8005 growth in photobioreactor: model and simulation of the ISS and ground experiments. Life Sci. Space Res. 25, 5365. 10.1016/j.lssr.2020.03.002 Poughon L. (2007b). MELiSSA loop mass balance modelling with MatLab/simulink. MELiSSA Technical Note 79.2. Laboratoire de Génie Chimique et Biochimique Available at: https://www.melissafoundation.org/download/197 (Accessed March 23, 2023). Poughon L. (1994a). MELiSSA simulation and modelling. Spirulina modelling. Variable global stoichiometric equation of Spirulina platensis in different light conditions. Effects of light on diet composition. Best conditions to insure maximal mass recycling. MELiSSA Technical Note 17.3. Laboratoire de Génie Chimique et Biochimique Available at: https://www.melissafoundation.org/download/572 (Accessed March 23, 2023). Poughon L. (1994b). MELiSSA simulation and modelling. Study of partial conversion. The stoichiometric assumptions: origin and effects. The main consequences of the partial conversion of components. MELiSSA Technical Note 17.2. Laboratoire de Génie Chimique et Biochimique Available at: https://www.melissafoundation.org/download/571 (Accessed March 23, 2023). Poulet L. Dussap C.-G. Fontaine J.-P. (2018). A physical modeling approach for higher plant growth in reduced gravity environments. Astrobiology 18, 10931100. 10.1089/ast.2017.1804 Poulet L. Dussap C.-G. Fontaine J-P. (2020). Development of a mechanistic model of leaf surface gas exchange coupling mass and energy balances for life-support systems applications. Acta Astronaut. 175, 517530. 10.1016/j.actaastro.2020.03.048 Putnam D. F. (1971). Composition and concentrative properties of human urine. Available at: https://ntrs.nasa.gov/citations/19710023044 (Accessed March 23, 2023). Ramsay I. R. Pullammanappallil P. C. (2001). Protein degradation during anaerobic wastewater treatment: derivation of stoichiometry. Biodegradation 12, 247257. 10.1023/A:1013116728817 Rose C. Parker A. Jefferson B. Cartmell E. (2015). The characterization of feces and urine: a review of the literature to inform advanced treatment Technology. Crit. Rev. Environ. Sci. Technol. 45, 18271879. 10.1080/10643389.2014.1000761 Sachdeva N. Poughon L. Gerbi O. Dussap C.-G. Lasseur C. Leroy B. (2021). Ground demonstration of the use of Limnospira indica for air revitalization in a bioregenerative life-support system setup: effect of non-nitrified urine–derived nitrogen sources. Front. Astronomy Space Sci. 8. 10.3389/fspas.2021.700270 Shin H.-S. Youn J.-H. Kim S.-H. (2004). Hydrogen production from food waste in anaerobic mesophilic and thermophilic acidogenesis. Int. J. Hydrogen Energy 29, 13551363. 10.1016/j.ijhydene.2003.09.011 Stasiak M. Gidzinski D. Jordan M. Dixon M. (2012). Crop selection for advanced life support systems in the ESA MELiSSA program: durum wheat (Triticum turgidum var durum). Adv. Space Res. 49, 16841690. 10.1016/j.asr.2012.03.001 Tepari E. A. (2019). Evaluation of biohydrogen production from Co-fermentation of carbohydrates and proteins. Available at: https://ir.lib.uwo.ca/etd/6566. Teuling E. Wierenga P. A. Schrama J. W. Gruppen H. (2017). Comparison of protein extracts from various unicellular green sources. J. Agric. Food Chem. 65, 79898002. 10.1021/acs.jafc.7b01788 Thiron B. (2020). Global control of a life support system: flow control and optimisation. Stockholm, Sweden: KTH School of Electrical Engineering and Computer Science. Available at: http://urn.kb.se/resolve?urn=urn:nbn:se:kth:diva-278509 (Accessed March 23, 2023). Tikhomirov A. A. Ushakova S. A. Velichko V. V. Tikhomirova N. A. Kudenko Yu. A. Gribovskaya I. V. (2011). Assessment of the possibility of establishing material cycling in an experimental model of the bio-technical life support system with plant and human wastes included in mass exchange. Acta Astronaut. 68, 15481554. 10.1016/j.actaastro.2010.10.005 Tranquille N. Emeis J. J. (1996). Rhodospirillum rubrum food acceptability study for MELiSSA. MELiSSA Technical Note 30.1. TNO Prevention and Health Available at: https://www.melissafoundation.org/download/554 (Accessed March 23, 2023). Tusé D. Miller M. W. (1984). Single‐cell protein: current status and future prospects. C R C Crit. Rev. Food Sci. Nutr. 19, 273325. 10.1080/10408398409527379 U.S. Food & Drug Administration (2017). Nutrition information for raw vegetables. Available at: https://www.fda.gov/food/food-labeling-nutrition/nutrition-information-raw-vegetables (Accessed March 23, 2023). Usherwood N. Segars W. (2001). Nitrogen interactions with phosphorus and potassium for optimum crop yield, nitrogen use effectiveness, and environmental stewardship. TheScientificWorldJournal 1, 5760. Suppl 2. 10.1100/tsw.2001.97 Vermeulen A. C. J. Papic A. Brazier F. (2018a). “Modeling and simulating the MELiSSA loop to understand the effects of system interaction on survivability during long-duration interstellar missions: an agent-based approach,” in AgroSpace-MELiSSA Conference, Rome, Italy, May 16–18, 2018. 10.5281/zenodo.7791088 Vermeulen A. C. J. Sirenko M. Papic A. Kiem J. Theys A. Brazier F. (2018b). “Evolving asteroid Starships: a bio-inspired approach for interstellar space systems,” in Proceedings of the 69th International Astronautical Congress (IAC), Bremen, Germany, October 1–5, 2018. Vermeulen A. C. Papic A. Kiem J. Hallak D. Brazier F. (2019). “Modeling and simulating a regenerative life support system to understand the effects of system interaction on survivability during deep space missions: an agent-based approach,” in Proceedings of the 70th International Astronautical Congress (IAC), Washington DC, October 21–25, 2019. Vermeulen A. C. J. Hubers C. de Vries L. Brazier F. (2020a). What horticulture and space exploration can learn from each other: the Mission to Mars initiative in the Netherlands. Acta Astronaut. 177, 421424. 10.1016/j.actaastro.2020.05.015 Vermeulen A. C. J. Papic A. Kiem J. Hallak D. Brazier F. (2020b). “Exploring the impact of irregular metabolic efficiencies and the space environment on the survivability of a regenerative life support system through agent-based modeling,” in MELiSSA Conference 2020, November 3–5, 2020,online. 10.5281/zenodo.7791110 Vijapur S. H. Hall T. D. Snyder S. Inman M. Taylor E. J. Skinn B. (2017). Electrochemical peroxide generation. ECS Trans. 77, 947962. 10.1149/07711.0947ecst Volk T. Rummel J. D. (1987). Mass balances for a biological life support system simulation model. Adv. Space Res. 7, 141148. 10.1016/0273-1177(87)90045-7 Vrati S. (1984). Single cell protein production by photosynthetic bacteria grown on the clarified effluents of biogas plant. Appl. Microbiol. Biotechnol. 19, 199202. 10.1007/BF00256454 Waters G. C. R. Dixon M. A. Masot A. Albiol J. Godia F. (2004). “Static mass balance studies of the MELiSSA pilot plant: integration of a higher plant chamber,” in 34th International Conference on Environmental Systems (ICES), Colorado Springs, CO, 19-22 July 2004 (SAE Technical Paper Series) (Warrendale, PA: SAE International). WebMD (2021). Are there health benefits of durum wheat? Available at: https://www.webmd.com/diet/are-there-health-benefits-of-durum-wheat (Accessed March 23, 2023). Weng C. Jeris J. S. (1976). Biochemical mechanisms in the methane fermentation of glutamic and oleic acids. Water Res. 10, 918. 10.1016/0043-1354(76)90151-2 Yin J. Yu X. Wang K. Shen D. (2016). Acidogenic fermentation of the main substrates of food waste to produce volatile fatty acids. Int. J. Hydrogen Energy 41, 2171321720. 10.1016/j.ijhydene.2016.07.094 Zambelli A. León A. Garcés R. (2015). “Mutagenesis in sunflower,” in Sunflower. Editors Sunflower Martínez-Force E. Dunford N. T. Salas J. J. (Urbana, Illinois: AOCS Press), 2752. 10.1016/B978-1-893997-94-3.50008-8 Zhang D. Clauwaert P. Luther A. López Barreiro D. Prins W. Brilman D. W. F. (2018). Sub- and supercritical water oxidation of anaerobic fermentation sludge for carbon and nitrogen recovery in a regenerative life support system. Waste Manag. 77, 268275. 10.1016/j.wasman.2018.04.008 Zhang D. Luther A. Clauwaert P. Ronsse F. (2019). Mild temperature hydrothermal oxidation of anaerobic fermentation filtrate for carbon and nitrogen recovery in a regenerative life support system. J. Supercrit. Fluids 145, 3947. 10.1016/j.supflu.2018.11.022
      ‘Oh, my dear Thomas, you haven’t heard the terrible news then?’ she said. ‘I thought you would be sure to have seen it placarded somewhere. Alice went straight to her room, and I haven’t seen her since, though I repeatedly knocked at the door, which she has locked on the inside, and I’m sure it’s most unnatural of her not to let her own mother comfort her. It all happened in a moment: I have always said those great motor-cars shouldn’t be allowed to career about the streets, especially when they are all paved with cobbles as they are at Easton Haven, which are{331} so slippery when it’s wet. He slipped, and it went over him in a moment.’ My thanks were few and awkward, for there still hung to the missive a basting thread, and it was as warm as a nestling bird. I bent low--everybody was emotional in those days--kissed the fragrant thing, thrust it into my bosom, and blushed worse than Camille. "What, the Corner House victim? Is that really a fact?" "My dear child, I don't look upon it in that light at all. The child gave our picturesque friend a certain distinction--'My husband is dead, and this is my only child,' and all that sort of thing. It pays in society." leave them on the steps of a foundling asylum in order to insure [See larger version] Interoffice guff says you're planning definite moves on your own, J. O., and against some opposition. Is the Colonel so poor or so grasping—or what? Albert could not speak, for he felt as if his brains and teeth were rattling about inside his head. The rest of[Pg 188] the family hunched together by the door, the boys gaping idiotically, the girls in tears. "Now you're married." The host was called in, and unlocked a drawer in which they were deposited. The galleyman, with visible reluctance, arrayed himself in the garments, and he was observed to shudder more than once during the investiture of the dead man's apparel. HoME香京julia种子在线播放 ENTER NUMBET 0016www.jxxxlwpq.org.cn
      mimalm.org.cn
      www.fushipifa.net.cn
      www.lpchain.com.cn
      www.gdvnet.com.cn
      gpmygj.com.cn
      ntjfc.com.cn
      www.usmpbo.com.cn
      szlddq.com.cn
      www.szcourt.org.cn
      处女被大鸡巴操 强奸乱伦小说图片 俄罗斯美女爱爱图 调教强奸学生 亚洲女的穴 夜来香图片大全 美女性强奸电影 手机版色中阁 男性人体艺术素描图 16p成人 欧美性爱360 电影区 亚洲电影 欧美电影 经典三级 偷拍自拍 动漫电影 乱伦电影 变态另类 全部电 类似狠狠鲁的网站 黑吊操白逼图片 韩国黄片种子下载 操逼逼逼逼逼 人妻 小说 p 偷拍10幼女自慰 极品淫水很多 黄色做i爱 日本女人人体电影快播看 大福国小 我爱肏屄美女 mmcrwcom 欧美多人性交图片 肥臀乱伦老头舔阴帝 d09a4343000019c5 西欧人体艺术b xxoo激情短片 未成年人的 插泰国人夭图片 第770弾み1 24p 日本美女性 交动态 eee色播 yantasythunder 操无毛少女屄 亚洲图片你懂的女人 鸡巴插姨娘 特级黄 色大片播 左耳影音先锋 冢本友希全集 日本人体艺术绿色 我爱被舔逼 内射 幼 美阴图 喷水妹子高潮迭起 和后妈 操逼 美女吞鸡巴 鸭个自慰 中国女裸名单 操逼肥臀出水换妻 色站裸体义术 中国行上的漏毛美女叫什么 亚洲妹性交图 欧美美女人裸体人艺照 成人色妹妹直播 WWW_JXCT_COM r日本女人性淫乱 大胆人艺体艺图片 女同接吻av 碰碰哥免费自拍打炮 艳舞写真duppid1 88电影街拍视频 日本自拍做爱qvod 实拍美女性爱组图 少女高清av 浙江真实乱伦迅雷 台湾luanlunxiaoshuo 洛克王国宠物排行榜 皇瑟电影yy频道大全 红孩儿连连看 阴毛摄影 大胆美女写真人体艺术摄影 和风骚三个媳妇在家做爱 性爱办公室高清 18p2p木耳 大波撸影音 大鸡巴插嫩穴小说 一剧不超两个黑人 阿姨诱惑我快播 幼香阁千叶县小学生 少女妇女被狗强奸 曰人体妹妹 十二岁性感幼女 超级乱伦qvod 97爱蜜桃ccc336 日本淫妇阴液 av海量资源999 凤凰影视成仁 辰溪四中艳照门照片 先锋模特裸体展示影片 成人片免费看 自拍百度云 肥白老妇女 女爱人体图片 妈妈一女穴 星野美夏 日本少女dachidu 妹子私处人体图片 yinmindahuitang 舔无毛逼影片快播 田莹疑的裸体照片 三级电影影音先锋02222 妻子被外国老头操 观月雏乃泥鳅 韩国成人偷拍自拍图片 强奸5一9岁幼女小说 汤姆影院av图片 妹妹人艺体图 美女大驱 和女友做爱图片自拍p 绫川まどか在线先锋 那么嫩的逼很少见了 小女孩做爱 处女好逼连连看图图 性感美女在家做爱 近距离抽插骚逼逼 黑屌肏金毛屄 日韩av美少女 看喝尿尿小姐日逼色色色网图片 欧美肛交新视频 美女吃逼逼 av30线上免费 伊人在线三级经典 新视觉影院t6090影院 最新淫色电影网址 天龙影院远古手机版 搞老太影院 插进美女的大屁股里 私人影院加盟费用 www258dd 求一部电影里面有一个二猛哥 深肛交 日本萌妹子人体艺术写真图片 插入屄眼 美女的木奶 中文字幕黄色网址影视先锋 九号女神裸 和骚人妻偷情 和潘晓婷做爱 国模大尺度蜜桃 欧美大逼50p 西西人体成人 李宗瑞继母做爱原图物处理 nianhuawang 男鸡巴的视屏 � 97免费色伦电影 好色网成人 大姨子先锋 淫荡巨乳美女教师妈妈 性nuexiaoshuo WWW36YYYCOM 长春继续给力进屋就操小女儿套干破内射对白淫荡 农夫激情社区 日韩无码bt 欧美美女手掰嫩穴图片 日本援交偷拍自拍 入侵者日本在线播放 亚洲白虎偷拍自拍 常州高见泽日屄 寂寞少妇自卫视频 人体露逼图片 多毛外国老太 变态乱轮手机在线 淫荡妈妈和儿子操逼 伦理片大奶少女 看片神器最新登入地址sqvheqi345com账号群 麻美学姐无头 圣诞老人射小妞和强奸小妞动话片 亚洲AV女老师 先锋影音欧美成人资源 33344iucoom zV天堂电影网 宾馆美女打炮视频 色五月丁香五月magnet 嫂子淫乱小说 张歆艺的老公 吃奶男人视频在线播放 欧美色图男女乱伦 avtt2014ccvom 性插色欲香影院 青青草撸死你青青草 99热久久第一时间 激情套图卡通动漫 幼女裸聊做爱口交 日本女人被强奸乱伦 草榴社区快播 2kkk正在播放兽骑 啊不要人家小穴都湿了 www猎奇影视 A片www245vvcomwwwchnrwhmhzcn 搜索宜春院av wwwsee78co 逼奶鸡巴插 好吊日AV在线视频19gancom 熟女伦乱图片小说 日本免费av无码片在线开苞 鲁大妈撸到爆 裸聊官网 德国熟女xxx 新不夜城论坛首页手机 女虐男网址 男女做爱视频华为网盘 激情午夜天亚洲色图 内裤哥mangent 吉沢明歩制服丝袜WWWHHH710COM 屌逼在线试看 人体艺体阿娇艳照 推荐一个可以免费看片的网站如果被QQ拦截请复制链接在其它浏览器打开xxxyyy5comintr2a2cb551573a2b2e 欧美360精品粉红鲍鱼 教师调教第一页 聚美屋精品图 中韩淫乱群交 俄罗斯撸撸片 把鸡巴插进小姨子的阴道 干干AV成人网 aolasoohpnbcn www84ytom 高清大量潮喷www27dyycom 宝贝开心成人 freefronvideos人母 嫩穴成人网gggg29com 逼着舅妈给我口交肛交彩漫画 欧美色色aV88wwwgangguanscom 老太太操逼自拍视频 777亚洲手机在线播放 有没有夫妻3p小说 色列漫画淫女 午间色站导航 欧美成人处女色大图 童颜巨乳亚洲综合 桃色性欲草 色眯眯射逼 无码中文字幕塞外青楼这是一个 狂日美女老师人妻 爱碰网官网 亚洲图片雅蠛蝶 快播35怎么搜片 2000XXXX电影 新谷露性家庭影院 深深候dvd播放 幼齿用英语怎么说 不雅伦理无需播放器 国外淫荡图片 国外网站幼幼嫩网址 成年人就去色色视频快播 我鲁日日鲁老老老我爱 caoshaonvbi 人体艺术avav 性感性色导航 韩国黄色哥来嫖网站 成人网站美逼 淫荡熟妇自拍 欧美色惰图片 北京空姐透明照 狼堡免费av视频 www776eom 亚洲无码av欧美天堂网男人天堂 欧美激情爆操 a片kk266co 色尼姑成人极速在线视频 国语家庭系列 蒋雯雯 越南伦理 色CC伦理影院手机版 99jbbcom 大鸡巴舅妈 国产偷拍自拍淫荡对话视频 少妇春梦射精 开心激动网 自拍偷牌成人 色桃隐 撸狗网性交视频 淫荡的三位老师 伦理电影wwwqiuxia6commqiuxia6com 怡春院分站 丝袜超短裙露脸迅雷下载 色制服电影院 97超碰好吊色男人 yy6080理论在线宅男日韩福利大全 大嫂丝袜 500人群交手机在线 5sav 偷拍熟女吧 口述我和妹妹的欲望 50p电脑版 wwwavtttcon 3p3com 伦理无码片在线看 欧美成人电影图片岛国性爱伦理电影 先锋影音AV成人欧美 我爱好色 淫电影网 WWW19MMCOM 玛丽罗斯3d同人动画h在线看 动漫女孩裸体 超级丝袜美腿乱伦 1919gogo欣赏 大色逼淫色 www就是撸 激情文学网好骚 A级黄片免费 xedd5com 国内的b是黑的 快播美国成年人片黄 av高跟丝袜视频 上原保奈美巨乳女教师在线观看 校园春色都市激情fefegancom 偷窥自拍XXOO 搜索看马操美女 人本女优视频 日日吧淫淫 人妻巨乳影院 美国女子性爱学校 大肥屁股重口味 啪啪啪啊啊啊不要 操碰 japanfreevideoshome国产 亚州淫荡老熟女人体 伦奸毛片免费在线看 天天影视se 樱桃做爱视频 亚卅av在线视频 x奸小说下载 亚洲色图图片在线 217av天堂网 东方在线撸撸-百度 幼幼丝袜集 灰姑娘的姐姐 青青草在线视频观看对华 86papa路con 亚洲1AV 综合图片2区亚洲 美国美女大逼电影 010插插av成人网站 www色comwww821kxwcom 播乐子成人网免费视频在线观看 大炮撸在线影院 ,www4KkKcom 野花鲁最近30部 wwwCC213wapwww2233ww2download 三客优最新地址 母亲让儿子爽的无码视频 全国黄色片子 欧美色图美国十次 超碰在线直播 性感妖娆操 亚洲肉感熟女色图 a片A毛片管看视频 8vaa褋芯屑 333kk 川岛和津实视频 在线母子乱伦对白 妹妹肥逼五月 亚洲美女自拍 老婆在我面前小说 韩国空姐堪比情趣内衣 干小姐综合 淫妻色五月 添骚穴 WM62COM 23456影视播放器 成人午夜剧场 尼姑福利网 AV区亚洲AV欧美AV512qucomwwwc5508com 经典欧美骚妇 震动棒露出 日韩丝袜美臀巨乳在线 av无限吧看 就去干少妇 色艺无间正面是哪集 校园春色我和老师做爱 漫画夜色 天海丽白色吊带 黄色淫荡性虐小说 午夜高清播放器 文20岁女性荫道口图片 热国产热无码热有码 2015小明发布看看算你色 百度云播影视 美女肏屄屄乱轮小说 家族舔阴AV影片 邪恶在线av有码 父女之交 关于处女破处的三级片 极品护士91在线 欧美虐待女人视频的网站 享受老太太的丝袜 aaazhibuo 8dfvodcom成人 真实自拍足交 群交男女猛插逼 妓女爱爱动态 lin35com是什么网站 abp159 亚洲色图偷拍自拍乱伦熟女抠逼自慰 朝国三级篇 淫三国幻想 免费的av小电影网站 日本阿v视频免费按摩师 av750c0m 黄色片操一下 巨乳少女车震在线观看 操逼 免费 囗述情感一乱伦岳母和女婿 WWW_FAMITSU_COM 偷拍中国少妇在公车被操视频 花也真衣论理电影 大鸡鸡插p洞 新片欧美十八岁美少 进击的巨人神thunderftp 西方美女15p 深圳哪里易找到老女人玩视频 在线成人有声小说 365rrr 女尿图片 我和淫荡的小姨做爱 � 做爱技术体照 淫妇性爱 大学生私拍b 第四射狠狠射小说 色中色成人av社区 和小姨子乱伦肛交 wwwppp62com 俄罗斯巨乳人体艺术 骚逼阿娇 汤芳人体图片大胆 大胆人体艺术bb私处 性感大胸骚货 哪个网站幼女的片多 日本美女本子把 色 五月天 婷婷 快播 美女 美穴艺术 色百合电影导航 大鸡巴用力 孙悟空操美少女战士 狠狠撸美女手掰穴图片 古代女子与兽类交 沙耶香套图 激情成人网区 暴风影音av播放 动漫女孩怎么插第3个 mmmpp44 黑木麻衣无码ed2k 淫荡学姐少妇 乱伦操少女屄 高中性爱故事 骚妹妹爱爱图网 韩国模特剪长发 大鸡巴把我逼日了 中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片 大胆女人下体艺术图片 789sss 影音先锋在线国内情侣野外性事自拍普通话对白 群撸图库 闪现君打阿乐 ady 小说 插入表妹嫩穴小说 推荐成人资源 网络播放器 成人台 149大胆人体艺术 大屌图片 骚美女成人av 春暖花开春色性吧 女亭婷五月 我上了同桌的姐姐 恋夜秀场主播自慰视频 yzppp 屄茎 操屄女图 美女鲍鱼大特写 淫乱的日本人妻山口玲子 偷拍射精图 性感美女人体艺木图片 种马小说完本 免费电影院 骑士福利导航导航网站 骚老婆足交 国产性爱一级电影 欧美免费成人花花性都 欧美大肥妞性爱视频 家庭乱伦网站快播 偷拍自拍国产毛片 金发美女也用大吊来开包 缔D杏那 yentiyishu人体艺术ytys WWWUUKKMCOM 女人露奶 � 苍井空露逼 老荡妇高跟丝袜足交 偷偷和女友的朋友做爱迅雷 做爱七十二尺 朱丹人体合成 麻腾由纪妃 帅哥撸播种子图 鸡巴插逼动态图片 羙国十次啦中文 WWW137AVCOM 神斗片欧美版华语 有气质女人人休艺术 由美老师放屁电影 欧美女人肉肏图片 白虎种子快播 国产自拍90后女孩 美女在床上疯狂嫩b 饭岛爱最后之作 幼幼强奸摸奶 色97成人动漫 两性性爱打鸡巴插逼 新视觉影院4080青苹果影院 嗯好爽插死我了 阴口艺术照 李宗瑞电影qvod38 爆操舅母 亚洲色图七七影院 被大鸡巴操菊花 怡红院肿么了 成人极品影院删除 欧美性爱大图色图强奸乱 欧美女子与狗随便性交 苍井空的bt种子无码 熟女乱伦长篇小说 大色虫 兽交幼女影音先锋播放 44aad be0ca93900121f9b 先锋天耗ばさ无码 欧毛毛女三级黄色片图 干女人黑木耳照 日本美女少妇嫩逼人体艺术 sesechangchang 色屄屄网 久久撸app下载 色图色噜 美女鸡巴大奶 好吊日在线视频在线观看 透明丝袜脚偷拍自拍 中山怡红院菜单 wcwwwcom下载 骑嫂子 亚洲大色妣 成人故事365ahnet 丝袜家庭教mp4 幼交肛交 妹妹撸撸大妈 日本毛爽 caoprom超碰在email 关于中国古代偷窥的黄片 第一会所老熟女下载 wwwhuangsecome 狼人干综合新地址HD播放 变态儿子强奸乱伦图 强奸电影名字 2wwwer37com 日本毛片基地一亚洲AVmzddcxcn 暗黑圣经仙桃影院 37tpcocn 持月真由xfplay 好吊日在线视频三级网 我爱背入李丽珍 电影师傅床戏在线观看 96插妹妹sexsex88com 豪放家庭在线播放 桃花宝典极夜著豆瓜网 安卓系统播放神器 美美网丝袜诱惑 人人干全免费视频xulawyercn av无插件一本道 全国色五月 操逼电影小说网 good在线wwwyuyuelvcom www18avmmd 撸波波影视无插件 伊人幼女成人电影 会看射的图片 小明插看看 全裸美女扒开粉嫩b 国人自拍性交网站 萝莉白丝足交本子 七草ちとせ巨乳视频 摇摇晃晃的成人电影 兰桂坊成社人区小说www68kqcom 舔阴论坛 久撸客一撸客色国内外成人激情在线 明星门 欧美大胆嫩肉穴爽大片 www牛逼插 性吧星云 少妇性奴的屁眼 人体艺术大胆mscbaidu1imgcn 最新久久色色成人版 l女同在线 小泽玛利亚高潮图片搜索 女性裸b图 肛交bt种子 最热门有声小说 人间添春色 春色猜谜字 樱井莉亚钢管舞视频 小泽玛利亚直美6p 能用的h网 还能看的h网 bl动漫h网 开心五月激 东京热401 男色女色第四色酒色网 怎么下载黄色小说 黄色小说小栽 和谐图城 乐乐影院 色哥导航 特色导航 依依社区 爱窝窝在线 色狼谷成人 91porn 包要你射电影 色色3A丝袜 丝袜妹妹淫网 爱色导航(荐) 好男人激情影院 坏哥哥 第七色 色久久 人格分裂 急先锋 撸撸射中文网 第一会所综合社区 91影院老师机 东方成人激情 怼莪影院吹潮 老鸭窝伊人无码不卡无码一本道 av女柳晶电影 91天生爱风流作品 深爱激情小说私房婷婷网 擼奶av 567pao 里番3d一家人野外 上原在线电影 水岛津实透明丝袜 1314酒色 网旧网俺也去 0855影院 在线无码私人影院 搜索 国产自拍 神马dy888午夜伦理达达兔 农民工黄晓婷 日韩裸体黑丝御姐 屈臣氏的燕窝面膜怎么样つぼみ晶エリーの早漏チ○ポ强化合宿 老熟女人性视频 影音先锋 三上悠亚ol 妹妹影院福利片 hhhhhhhhsxo 午夜天堂热的国产 强奸剧场 全裸香蕉视频无码 亚欧伦理视频 秋霞为什么给封了 日本在线视频空天使 日韩成人aⅴ在线 日本日屌日屄导航视频 在线福利视频 日本推油无码av magnet 在线免费视频 樱井梨吮东 日本一本道在线无码DVD 日本性感诱惑美女做爱阴道流水视频 日本一级av 汤姆avtom在线视频 台湾佬中文娱乐线20 阿v播播下载 橙色影院 奴隶少女护士cg视频 汤姆在线影院无码 偷拍宾馆 业面紧急生级访问 色和尚有线 厕所偷拍一族 av女l 公交色狼优酷视频 裸体视频AV 人与兽肉肉网 董美香ol 花井美纱链接 magnet 西瓜影音 亚洲 自拍 日韩女优欧美激情偷拍自拍 亚洲成年人免费视频 荷兰免费成人电影 深喉呕吐XXⅩX 操石榴在线视频 天天色成人免费视频 314hu四虎 涩久免费视频在线观看 成人电影迅雷下载 能看见整个奶子的香蕉影院 水菜丽百度影音 gwaz079百度云 噜死你们资源站 主播走光视频合集迅雷下载 thumbzilla jappen 精品Av 古川伊织star598在线 假面女皇vip在线视频播放 国产自拍迷情校园 啪啪啪公寓漫画 日本阿AV 黄色手机电影 欧美在线Av影院 华裔电击女神91在线 亚洲欧美专区 1日本1000部免费视频 开放90后 波多野结衣 东方 影院av 页面升级紧急访问每天正常更新 4438Xchengeren 老炮色 a k福利电影 色欲影视色天天视频 高老庄aV 259LUXU-683 magnet 手机在线电影 国产区 欧美激情人人操网 国产 偷拍 直播 日韩 国内外激情在线视频网给 站长统计一本道人妻 光棍影院被封 紫竹铃取汁 ftp 狂插空姐嫩 xfplay 丈夫面前 穿靴子伪街 XXOO视频在线免费 大香蕉道久在线播放 电棒漏电嗨过头 充气娃能看下毛和洞吗 夫妻牲交 福利云点墦 yukun瑟妃 疯狂交换女友 国产自拍26页 腐女资源 百度云 日本DVD高清无码视频 偷拍,自拍AV伦理电影 A片小视频福利站。 大奶肥婆自拍偷拍图片 交配伊甸园 超碰在线视频自拍偷拍国产 小热巴91大神 rctd 045 类似于A片 超美大奶大学生美女直播被男友操 男友问 你的衣服怎么脱掉的 亚洲女与黑人群交视频一 在线黄涩 木内美保步兵番号 鸡巴插入欧美美女的b舒服 激情在线国产自拍日韩欧美 国语福利小视频在线观看 作爱小视颍 潮喷合集丝袜无码mp4 做爱的无码高清视频 牛牛精品 伊aⅤ在线观看 savk12 哥哥搞在线播放 在线电一本道影 一级谍片 250pp亚洲情艺中心,88 欧美一本道九色在线一 wwwseavbacom色av吧 cos美女在线 欧美17,18ⅹⅹⅹ视频 自拍嫩逼 小电影在线观看网站 筱田优 贼 水电工 5358x视频 日本69式视频有码 b雪福利导航 韩国女主播19tvclub在线 操逼清晰视频 丝袜美女国产视频网址导航 水菜丽颜射房间 台湾妹中文娱乐网 风吟岛视频 口交 伦理 日本熟妇色五十路免费视频 A级片互舔 川村真矢Av在线观看 亚洲日韩av 色和尚国产自拍 sea8 mp4 aV天堂2018手机在线 免费版国产偷拍a在线播放 狠狠 婷婷 丁香 小视频福利在线观看平台 思妍白衣小仙女被邻居强上 萝莉自拍有水 4484新视觉 永久发布页 977成人影视在线观看 小清新影院在线观 小鸟酱后丝后入百度云 旋风魅影四级 香蕉影院小黄片免费看 性爱直播磁力链接 小骚逼第一色影院 性交流的视频 小雪小视频bd 小视频TV禁看视频 迷奸AV在线看 nba直播 任你在干线 汤姆影院在线视频国产 624u在线播放 成人 一级a做爰片就在线看狐狸视频 小香蕉AV视频 www182、com 腿模简小育 学生做爱视频 秘密搜查官 快播 成人福利网午夜 一级黄色夫妻录像片 直接看的gav久久播放器 国产自拍400首页 sm老爹影院 谁知道隔壁老王网址在线 综合网 123西瓜影音 米奇丁香 人人澡人人漠大学生 色久悠 夜色视频你今天寂寞了吗? 菲菲影视城美国 被抄的影院 变态另类 欧美 成人 国产偷拍自拍在线小说 不用下载安装就能看的吃男人鸡巴视频 插屄视频 大贯杏里播放 wwwhhh50 233若菜奈央 伦理片天海翼秘密搜查官 大香蕉在线万色屋视频 那种漫画小说你懂的 祥仔电影合集一区 那里可以看澳门皇冠酒店a片 色自啪 亚洲aV电影天堂 谷露影院ar toupaizaixian sexbj。com 毕业生 zaixian mianfei 朝桐光视频 成人短视频在线直接观看 陈美霖 沈阳音乐学院 导航女 www26yjjcom 1大尺度视频 开平虐女视频 菅野雪松协和影视在线视频 华人play在线视频bbb 鸡吧操屄视频 多啪啪免费视频 悠草影院 金兰策划网 (969) 橘佑金短视频 国内一极刺激自拍片 日本制服番号大全magnet 成人动漫母系 电脑怎么清理内存 黄色福利1000 dy88午夜 偷拍中学生洗澡磁力链接 花椒相机福利美女视频 站长推荐磁力下载 mp4 三洞轮流插视频 玉兔miki热舞视频 夜生活小视频 爆乳人妖小视频 国内网红主播自拍福利迅雷下载 不用app的裸裸体美女操逼视频 变态SM影片在线观看 草溜影院元气吧 - 百度 - 百度 波推全套视频 国产双飞集合ftp 日本在线AV网 笔国毛片 神马影院女主播是我的邻居 影音资源 激情乱伦电影 799pao 亚洲第一色第一影院 av视频大香蕉 老梁故事汇希斯莱杰 水中人体磁力链接 下载 大香蕉黄片免费看 济南谭崔 避开屏蔽的岛a片 草破福利 要看大鸡巴操小骚逼的人的视频 黑丝少妇影音先锋 欧美巨乳熟女磁力链接 美国黄网站色大全 伦蕉在线久播 极品女厕沟 激情五月bd韩国电影 混血美女自摸和男友激情啪啪自拍诱人呻吟福利视频 人人摸人人妻做人人看 44kknn 娸娸原网 伊人欧美 恋夜影院视频列表安卓青青 57k影院 如果电话亭 avi 插爆骚女精品自拍 青青草在线免费视频1769TV 令人惹火的邻家美眉 影音先锋 真人妹子被捅动态图 男人女人做完爱视频15 表姐合租两人共处一室晚上她竟爬上了我的床 性爱教学视频 北条麻妃bd在线播放版 国产老师和师生 magnet wwwcctv1024 女神自慰 ftp 女同性恋做激情视频 欧美大胆露阴视频 欧美无码影视 好女色在线观看 后入肥臀18p 百度影视屏福利 厕所超碰视频 强奸mp magnet 欧美妹aⅴ免费线上看 2016年妞干网视频 5手机在线福利 超在线最视频 800av:cOm magnet 欧美性爱免播放器在线播放 91大款肥汤的性感美乳90后邻家美眉趴着窗台后入啪啪 秋霞日本毛片网站 cheng ren 在线视频 上原亚衣肛门无码解禁影音先锋 美脚家庭教师在线播放 尤酷伦理片 熟女性生活视频在线观看 欧美av在线播放喷潮 194avav 凤凰AV成人 - 百度 kbb9999 AV片AV在线AV无码 爱爱视频高清免费观看 黄色男女操b视频 观看 18AV清纯视频在线播放平台 成人性爱视频久久操 女性真人生殖系统双性人视频 下身插入b射精视频 明星潜规测视频 mp4 免賛a片直播绪 国内 自己 偷拍 在线 国内真实偷拍 手机在线 国产主播户外勾在线 三桥杏奈高清无码迅雷下载 2五福电影院凸凹频频 男主拿鱼打女主,高宝宝 色哥午夜影院 川村まや痴汉 草溜影院费全过程免费 淫小弟影院在线视频 laohantuiche 啪啪啪喷潮XXOO视频 青娱乐成人国产 蓝沢润 一本道 亚洲青涩中文欧美 神马影院线理论 米娅卡莉法的av 在线福利65535 欧美粉色在线 欧美性受群交视频1在线播放 极品喷奶熟妇在线播放 变态另类无码福利影院92 天津小姐被偷拍 磁力下载 台湾三级电髟全部 丝袜美腿偷拍自拍 偷拍女生性行为图 妻子的乱伦 白虎少妇 肏婶骚屄 外国大妈会阴照片 美少女操屄图片 妹妹自慰11p 操老熟女的b 361美女人体 360电影院樱桃 爱色妹妹亚洲色图 性交卖淫姿势高清图片一级 欧美一黑对二白 大色网无毛一线天 射小妹网站 寂寞穴 西西人体模特苍井空 操的大白逼吧 骚穴让我操 拉好友干女朋友3p