Microorganism production and remediation processes are of critical importance to the next generation of sustainable industries. Undertaking mathematical treatment of dynamic biosystems operating at any spatial or temporal scale is essential to guarantee their performance and safety. However, constructing physical models remains a challenge due to the extreme complexity of process biological mechanisms. Data-driven models also encounter severe limitations because datasets from large-scale bioprocesses are often scarce without complete information and on a restricted operational space. To fill this gap, the current research compares the performance of advanced physical and data-driven models for dynamic bioprocess simulations subject to incomplete and scarce datasets, which to the best of our knowledge has never been addressed before. In specific, kinetic models were constructed by integrating different classic models, and state-of-the-art hyperparameter selection frameworks were developed to design artificial neural networks and Gaussian process regression models. An algae-bacteria consortium wastewater treatment process was selected to test the accuracy of these modeling strategies, as it is one of the most sophisticated biosystems due to the intricate mutualistic and competitive interactions. Based on the current results and available data, a heuristic model selection procedure is provided. This study paves the way to facilitate future bioprocess modeling.