Multiscale Porosity Microfluidics to Study Bacterial Transport in Heterogeneous Chemical Landscapes.

Microfluidic models are proving to be powerful systems to study fundamental processes in porous media, due to their ability to replicate topologically complex environments while allowing detailed, quantitative observations at the pore scale. Yet, while porous media such as living tissues, geological substrates, or industrial systems typically display a porosity that spans multiple scales, most microfluidic models to date are limited to a single porosity or a small range of pore sizes. Here, a novel microfluidic system with multiscale porosity is presented. By embedding polyacrylamide (PAAm) hydrogel structures through in-situ photopolymerization in a landscape of microfabricated polydimethylsiloxane (PDMS) pillars with varying spacing, micromodels with porosity spanning several orders of magnitude, from nanometers to millimeters are created. Experiments conducted at different porosity patterns demonstrate the potential of this approach to characterize fundamental and ubiquitous biological and geochemical transport processes in porous media. Accounting for multiscale porosity allows studies of the resulting heterogeneous fluid flow and concentration fields of transported chemicals, as well as the biological behaviors associated with this heterogeneity, such as bacterial chemotaxis. This approach brings laboratory studies of transport in porous media a step closer to their natural counterparts in the environment, industry, and medicine.

From the Mental to the Conceptual Model: The Challenge of Teaching Hydrogeology in the Field.

Field-based learning in hydrogeology enables students to develop their understanding and application of practical methodologies, and to enhance many of the generic skills (e.g., teamwork, problem-solving). However, teaching and learning hydrogeology in general, and especially in the field, presents cognitive difficulties, such as the diversity in student education and experience, the hidden nature of water movement and transport of chemicals, and the preexisting students' mental models of the subsurface, in particular. At any given experimental or teaching site there is only one reality for which lecturers can have an approximate conceptual model, including aquifer(s) geometry and functioning (e.g., flow direction). However, students' preconceptions (i.e., mental model), in some cases misconceptions, influence not only their outcome from the learning strategy designed, but also the conceptual model expression (i.e., flow chart, block diagram, or similar) for the study area or site. In practice, two general "teaching challenges" are identified to enable students' transition from the mental to the conceptual model: (1) identify and dispel any prior misconceptions and (2) show how to go from the partial information to the integration of new information for the development of the conceptual model. The inclusion of specific prior-to-field lessons in the classroom is recommended and in general, done. However, introducing a prior-to-field survey to learn about students' backgrounds, and methodologies for the development and expression of hydrogeological conceptual models and for testing multiple plausible conceptual models will help students transition from the mental to the conceptual model.

Global warming accelerates soil heterotrophic respiration.

Carbon efflux from soils is the largest terrestrial carbon source to the atmosphere, yet it is still one of the most uncertain fluxes in the Earth's carbon budget. A dominant component of this flux is heterotrophic respiration, influenced by several environmental factors, most notably soil temperature and moisture. Here, we develop a mechanistic model from micro to global scale to explore how changes in soil water content and temperature affect soil heterotrophic respiration. Simulations, laboratory measurements, and field observations validate the new approach. Estimates from the model show that heterotrophic respiration has been increasing since the 1980s at a rate of about 2% per decade globally. Using future projections of surface temperature and soil moisture, the model predicts a global increase of about 40% in heterotrophic respiration by the end of the century under the worst-case emission scenario, where the Arctic region is expected to experience a more than two-fold increase, driven primarily by declining soil moisture rather than temperature increase.

Internal Biofilm Heterogeneities Enhance Solute Mixing and Chemical Reactions in Porous Media.

Bacterial biofilms can form in porous media that are of interest in industrial applications ranging from medical implants to biofilters as well as in environmental applications such as in situ groundwater remediation, where they can be critical locations for biogeochemical reactions. The presence of biofilms modifies porous media topology and hydrodynamics by clogging pores and consequently solutes transport and reactions kinetics. The interplay between highly heterogeneous flow fields found in porous media and microbial behavior, including biofilm growth, results in a spatially heterogeneous biofilm distribution in the porous media as well as internal heterogeneity across the thickness of the biofilm. Our study leverages highly resolved three-dimensional X-ray computed microtomography images of bacterial biofilms in a tubular reactor to numerically compute pore-scale fluid flow and solute transport by considering multiple equivalent stochastically generated internal permeability fields for the biofilm. We show that the internal heterogeneous permeability mainly impacts intermediate velocities when compared with homogeneous biofilm permeability. While the equivalent internal permeability fields of the biofilm do not impact fluid-fluid mixing, they significantly control a fast reaction. For biologically driven reactions such as nutrient or contaminant uptake by the biofilm, its internal permeability field controls the efficiency of the process. This study highlights the importance of considering the internal heterogeneity of biofilms to better predict reactivity in industrial and environmental bioclogged porous systems.

Morphogenesis of Biofilms in Porous Media and Control on Hydrodynamics.

The functioning of natural and engineered porous media, like soils and filters, depends in many cases on the interplay between biochemical processes and hydrodynamics. In such complex environments, microorganisms often form surface-attached communities known as biofilms. Biofilms can take the shape of clusters, which alter the distribution of fluid flow velocities within the porous medium, subsequently influencing biofilm growth. Despite numerous experimental and numerical efforts, the control of the biofilm clustering process and the resulting heterogeneity in biofilm permeability is not well understood, limiting our predictive abilities for biofilm-porous medium systems. Here, we use a quasi-2D experimental model of a porous medium to characterize biofilm growth dynamics for different pore sizes and flow rates. We present a method to obtain the time-resolved biofilm permeability field from experimental images and use the obtained permeability field to compute the flow field through a numerical model. We observe a biofilm cluster size distribution characterized by a spectrum slope evolving in time between -2 and -1, a fundamental measure that can be used to create spatio-temporal distributions of biofilm clusters for upscaled models. We find a previously undescribed biofilm permeability distribution, which can be used to stochastically generate permeability fields within biofilms. An increase in velocity variance for a decrease in physical heterogeneity shows that the bioclogged porous medium behaves differently than expected from studies on heterogeneity in abiotic porous media.

Upscaling Mixing-Controlled Reactions in Unsaturated Porous Media.

We study mixing-controlled chemical reactions in unsaturated porous media from a pore-scale perspective. The spatial heterogeneity induced by the presence of two immiscible phases, here water and air, in the pore space generates complex flow patterns that dominate reactive mixing across scales. To assess the impact of different macroscopic saturation states (the fraction of pore volume occupied by water) on mixing-controlled chemical reactions, we consider a fast irreversible reaction between two initially segregated dissolved species that mix as one solution displaces the other in the heterogeneous flow field of the water phase. We use the pore-scale geometry and water distributions from the laboratory experiments reported by Jiménez-Martínez et al. (Geophys. Res. Lett. 42: 5316-5324, 2015). We analyze reactive mixing in three complementary ways. Firstly, we post-process experimentally observed spatially distributed concentration data; secondly, we perform numerical simulations of flow and reactive transport in the heterogeneous water phase, and thirdly, we use an upscaled mixing model. The first approach relies on an exact algebraic map between conservative and reactive species for an instantaneous irreversible bimolecular reaction that allows to estimate reactive mixing based on experimental conservative transport data. The second approach is based on reactive random walk particle tracking simulations in the numerically determined flow field in the water phase. The third approach uses a dispersive lamella approach that accounts for the impact of flow heterogeneity on mixing in terms of effective dispersion coefficients, which are estimated from both experimental data and numerical random walk particle tracking simulations. We observe a significant increase in reactive mixing for decreasing saturation, which is caused by the stronger heterogeneity of the water phase and thus of the flow field. This is consistently observed in the experimental data and the direct numerical simulations. The dispersive lamella model, parameterized by the effective interface width, provides robust estimates of the evolution of the product mass obtained from the experimental and numerical data.

Mixing Controlled Adsorption at the Liquid-Solid Interfaces in Unsaturated Porous Media.

The unsaturated zone, located between the soil surface and the phreatic level, plays an important role in defining the fate of any substance entering the subsoil. In addition to the processes of flow and transport taking place in the liquid phase, surface reactions such as adsorption to the solid phase may occur and increase the residence time of the substance entering the system. In this study, we aim to understand the pore-scale mechanisms that control adsorption in unsaturated systems. We combine 2D pore-scale experimental images with numerical simulations to analyze flow, transport, and adsorption under different liquid saturation degrees. We demonstrate the role of mixing on adsorption at the liquid-solid interfaces by analyzing the deformation in time of a pulse-injected surfactant. We also analyze the impact of the isotherm functional shape and the inclusion of the liquid-gas interfaces as adsorption sites on this surface reaction. The enhancement of mixing as saturation decreases is accompanied by a reduction in the amount of adsorbed mass, located mainly along preferential flow paths, where the solute is primarily transported. For the same isotherm, a nonlinear behavior of adsorption as a function of liquid saturation has been observed. This is explained by the nonlinear variation of the volume fraction of liquid behaving as preferential path or stagnation zone as liquid saturation decreases, despite the linear decrease in the surface area of solids accessible for adsorption. Supplementary Information: The online version contains supplementary material available at 10.1007/s11242-022-01747-x.

The role of biofilm matrix composition in controlling colony expansion and morphology.

Biofilms are biological viscoelastic gels composed of bacterial cells embedded in a self-secreted polymeric extracellular matrix (ECM). In environmental settings, such as in the rhizosphere and phyllosphere, biofilm colonization occurs at the solid-air interface. The biofilms' ability to colonize and expand over these surfaces depends on the formation of osmotic gradients and ECM viscoelastic properties. In this work, we study the influence of biofilm ECM components on its viscoelasticity and expansion, using the model organism Bacillus subtilis and deletion mutants of its three major ECM components, TasA, EPS and BslA. Using a multi-scale approach, we quantified macro-scale viscoelasticity and expansion dynamics. Furthermore, we used a microsphere assay to visualize the micro-scale expansion patterns. We find that the viscoelastic phase angle Φ is likely the best viscoelastic parameter correlating to biofilm expansion dynamics. Moreover, we quantify the sensitivity of the biofilm to changes in substrate water potential as a function of ECM composition. Finally, we find that the deletion of ECM components significantly increases the coherence of micro-scale colony expansion patterns. These results demonstrate the influence of ECM viscoelasticity and substrate water potential on the expansion of biofilm colonies on wet surfaces at the air-solid interface, commonly found in natural environments.

A Microfluidic Platform to Study Bioclogging in Porous Media.

Bacterial biofilms are found in several environmental and industrial porous media, including soils and filtration membranes. Biofilms grow under certain flow conditions and can clog pores, thereby redirecting the local fluid flow. The ability of biofilms to clog pores, the so-called bioclogging, can have a tremendous effect on the local permeability of the porous medium, creating a pressure buildup in the system, and impacting the mass flow through it. To understand the interplay between biofilm growth and fluid flow under different physical conditions (e.g., at different flow velocities and pore sizes), in the present study, a microfluidic platform is developed to visualize biofilm development using a microscope under externally-imposed, controlled physical conditions. The biofilm-induced pressure buildup in the porous medium can be measured simultaneously using pressure sensors and, later, correlated with the surface coverage of the biofilm. The presented platform provides a baseline for a systematic approach to investigate bioclogging caused by biofilms in porous media under flow conditions and can be adapted to studying environmental isolates or multispecies biofilms.

Spatiotemporal Mineral Phase Evolution and Arsenic Retention in Microfluidic Models of Zerovalent Iron-Based Water Treatment.

Arsenic (As) is a toxic element, and elevated levels of geogenic As in drinking water pose a threat to the health of several hundred million people worldwide. In this study, we used microfluidics in combination with optical microscopy and X-ray spectroscopy to investigate zerovalent iron (ZVI) corrosion, secondary iron (Fe) phase formation, and As retention processes at the pore scale in ZVI-based water treatment filters. Two 250 µm thick microchannels filled with single ZVI and quartz grain layers were operated intermittently (12 h flow/12 h no-flow) with synthetic groundwater (pH 7.5; 570 µg/L As(III)) over 13 and 49 days. Initially, lepidocrocite (Lp) and carbonate green rust (GRC) were the dominant secondary Fe-phases and underwent cyclic transformation. During no-flow, lepidocrocite partially transformed into GRC and small fractions of magnetite, kinetically limited by Fe(II) diffusion or by decreasing corrosion rates. When flow resumed, GRC rapidly and nearly completely transformed back into lepidocrocite. Longer filter operation combined with a prolonged no-flow period accelerated magnetite formation. Phosphate adsorption onto Fe-phases allowed for downstream calcium carbonate precipitation and, consequently, accelerated anoxic ZVI corrosion. Arsenic was retained on Fe-coated quartz grains and in zones of cyclic Lp-GRC transformation. Our results suggest that intermittent filter operation leads to denser secondary Fe-solids and thereby ensures prolonged filter performance.

Competition between growth and shear stress drives intermittency in preferential flow paths in porous medium biofilms.

Bacteria in porous media, such as soils, aquifers, and filters, often form surface-attached communities known as biofilms. Biofilms are affected by fluid flow through the porous medium, for example, for nutrient supply, and they, in turn, affect the flow. A striking example of this interplay is the strong intermittency in flow that can occur when biofilms nearly clog the porous medium. Intermittency manifests itself as the rapid opening and slow closing of individual preferential flow paths (PFPs) through the biofilm-porous medium structure, leading to continual spatiotemporal rearrangement. The drastic changes to the flow and mass transport induced by intermittency can affect the functioning and efficiency of natural and industrial systems. Yet, the mechanistic origin of intermittency remains unexplained. Here, we show that the mechanism driving PFP intermittency is the competition between microbial growth and shear stress. We combined microfluidic experiments quantifying Bacillus subtilis biofilm formation and behavior in synthetic porous media for different pore sizes and flow rates with a mathematical model accounting for flow through the biofilm and biofilm poroelasticity to reveal the underlying mechanisms. We show that the closing of PFPs is driven by microbial growth, controlled by nutrient mass flow. Opposing this, we find that the opening of PFPs is driven by flow-induced shear stress, which increases as a PFP becomes narrower due to microbial growth, causing biofilm compression and rupture. Our results demonstrate that microbial growth and its competition with shear stresses can lead to strong temporal variability in flow and transport conditions in bioclogged porous media.

Impact of a transformation from flood to drip irrigation on groundwater recharge and nitrogen leaching under variable climatic conditions.

The sustainability of agriculture in the Mediterranean climate is challenged by high irrigation water demands and nitrogen fertilizer losses to the environment, causing significant pressure on groundwater resources and groundwater-dependent ecosystems. Advanced irrigation technologies and improved fertilizer management have been promoted as key solutions to reduce the agricultural impact on aquatic systems. However, it remains unclear how different irrigation-fertilizer practices perform on the long-term under a highly variable climate, such as the Mediterranean one. Here, we conduct hydrological simulations over a fifty-year period to quantify the magnitude and dynamics of groundwater recharge and nitrogen leaching under five real-case irrigation-fertilizer practices observed in Valencia (eastern Spain). The Valencian Region is the largest citrus-producing region of Europe and current irrigation-fertilizer practices reflect the ongoing transformation of irrigation systems from flood to drip irrigation. Our simulations highlight three major implications of the irrigation transformation for groundwater resources. First, the transformation from flood to drip irrigation reduces the recharge fraction (19% vs. 16%) and especially the nitrogen leaching fraction (33% vs. 18%) on the long term. Second, the long-term performance of the two irrigation practices is subject to substantial inter-annual differences controlled by precipitation variability. The sensitivity of recharge and nitrogen leaching to annual meteorological conditions is stronger in drip irrigation, which eventually leads to a similar performance of flood and drip irrigation in wet years if fertilizer inputs are similar. Third, we identify a pronounced year-to-year nitrogen memory in the soil, whereby an enhanced (decreased) nitrogen leaching is observed after anomalously dry (wet) years, affecting the performance of irrigation-fertilizer practices. Overall, the study demonstrates the highly variable nature of the performance of irrigation-fertilizer practices, and the major findings can guide future efforts in designing sustainable water management strategies for agricultural areas with a Mediterranean climate.

Controlling pore-scale processes to tame subsurface biomineralization.

Microorganisms capable of biomineralization can catalyze mineral precipitation by modifying local physical and chemical conditions. In porous media, such as soil and rock, these microorganisms live and function in highly heterogeneous physical, chemical and ecological microenvironments, with strong local gradients created by both microbial activity and the pore-scale structure of the subsurface. Here, we focus on extracellular bacterial biomineralization, which is sensitive to external heterogeneity, and review the pore-scale processes controlling microbial biomineralization in natural and engineered porous media. We discuss how individual physical, chemical and ecological factors integrate to affect the spatial and temporal control of biomineralization, and how each of these factors contributes to a quantitative understanding of biomineralization in porous media. We find that an improved understanding of microbial behavior in heterogeneous microenvironments would promote understanding of natural systems and output in diverse technological applications, including improved representation and control of fluid mixing from pore to field scales. We suggest a range of directions by which future work can build from existing tools to advance each of these areas to improve understanding and predictability of biomineralization science and technology.

Phase Saturation Control on Mixing-Driven Reactions in 3D Porous Media.

Transported chemical reactions in unsaturated porous media are relevant to environmental and industrial applications. Continuum scale models are based on equivalent parameters derived from analogy with saturated conditions and cannot appropriately account for incomplete mixing. It is also unclear how the third dimension controls mixing and reactions. We obtain three-dimensional (3D) images by magnetic resonance imaging using an immiscible nonwetting liquid as a second phase and a fast irreversible bimolecular reaction. We study the impact of phase saturation on the dynamics of mixing and the reaction front. We quantify the temporally resolved effective reaction rate and describe it using the lamellar theory of mixing, which explains faster than Fickian (t0.5) rate of product formation by accounting for the deformation of the mixing interface between the two reacting fluids. For a given Péclet, although stretching and folding of the reactive front enhance as saturation decreases, enhancing the product formation, the product formation is larger as saturation increases. After breakthrough, the extinction of the reaction takes longer as saturation decreases because of the larger nonmixed volume behind the front. These results are the basis for a general model to better predict reactive transport in unsaturated porous media not achievable by the current continuum paradigm.

Crustal fingering facilitates free-gas methane migration through the hydrate stability zone.

Widespread seafloor methane venting has been reported in many regions of the world oceans in the past decade. Identifying and quantifying where and how much methane is being released into the ocean remains a major challenge and a critical gap in assessing the global carbon budget and predicting future climate [C. Ruppel, J. D. Kessler. Rev. Geophys. 55, 126-168 (2017)]. Methane hydrate ([Formula: see text]) is an ice-like solid that forms from methane-water mixture under elevated-pressure and low-temperature conditions typical of the deep marine settings (>600-m depth), often referred to as the hydrate stability zone (HSZ). Wide-ranging field evidence indicates that methane seepage often coexists with hydrate-bearing sediments within the HSZ, suggesting that hydrate formation may play an important role during the gas-migration process. At a depth that is too shallow for hydrate formation, existing theories suggest that gas migration occurs via capillary invasion and/or initiation and propagation of fractures (Fig. 1). Within the HSZ, however, a theoretical mechanism that addresses the way in which hydrate formation participates in the gas-percolation process is missing. Here, we study, experimentally and computationally, the mechanics of gas percolation under hydrate-forming conditions. We uncover a phenomenon-crustal fingering-and demonstrate how it may control methane-gas migration in ocean sediments within the HSZ.

Reduced gravity promotes bacterially mediated anoxic hotspots in unsaturated porous media.

Human endeavours into deep space exploration and the prospects of establishing colonies on nearby planets would invariably involve components of bioregenerative life support for food production, cabin atmosphere renewal, and waste recycling. Growing plants and their microbiomes in porous media under different gravitational fields may present new challenges due to effects of liquid distribution on gaseous exchange with roots and microorganisms. We provide the first direct evidence that capillary driven liquid reconfiguration in porous media under reduced gravity conditions reduces oxygen diffusion pathways and enhances anoxic conditions within bacterial hotspots. Parabolic flight experiments using model porous media inoculated with aerobic and facultative anaerobic bacteria reveal the systematic enhancement of anoxic conditions during the reduced gravity periods in the presence but not in the absence of bacterial activity. The promotion of anoxic conditions under reduced gravity may lead to higher nitrous oxide and methane emissions relative to Earth conditions, on the other hand, anoxic conditions could be beneficial for perchlorate bioremediation of Martian soil. The results highlight changes in soil bacterial microhabitats under reduced gravity and the challenges of managing bioregenerative life support systems in space.

Transport of Nano- and Microplastic through Unsaturated Porous Media from Sewage Sludge Application.

Wastewater treatment plants have been identified as important hubs for small particulate plastic, down to the nanometer scale, from urban areas to the environment. The reuse of sludge as fertilizer in agricultural practices can lead to accumulation of plastic in the soil. In this study, nanoplastic particles and microplastic fibers were synthesized with a passive inorganic tracer to aid in faster and more quantitative analysis using inductively coupled plasma mass spectrometry (ICP-MS). Using the anaerobic digestate of a pilot wastewater treatment plant spiked with metal-doped plastic, the excess sludge was dewatered, ensuring realistic associations between sludge and plastic. The resulting sludge cake was affixed atop an unsaturated porous-medium column of glass beads to assess: (i) the release of particulate plastic from the sludge, and (ii) the accumulation and mobility of plastic and organic matter through the column (analogous to a soil). A total of three particulate plastic treatments were assessed, in triplicate, where the plastic and mobile organic fractions were monitored for 14 pore water volumes. Due to size-limited transport, low deattachment from the sludge and reduced mobility through the column were found for microplastic fibers (>95% retention). However, cotransport between the mobile organic fraction and nanoplastic particles was observed, with 50% of both retained in the column. These results contribute to the understanding of the fate of particulate plastics and to assessing the associated environmental risks of particle mobility and percolation, particularly for nanoplastics.

Hydrogeological modelling for the watershed management of the Mar Menor coastal lagoon (Spain).

The Mar Menor is the largest lagoon along the Spanish Mediterranean coast. It suffers from eutrophication and algal blooms associated with intensive agricultural activities and urban pressure in the surrounding Campo de Cartagena plain. A balanced discharge of groundwater, carrier of algal nutrients such as nitrate, is essential to ensure the integrity of the coastal lagoon and the availability of groundwater resources inland. We here present a 3D hydrogeological model of the unconfined Quaternary aquifer that discharges into the lagoon. The model couples both surface water balance and groundwater dynamics and has been calibrated to available data in the period 2000-2016. The calibrated model allows understanding of the current state of the aquifer and its link to the lagoon. The potential discharge has been quantified in both space and time and falls between 69.5 and 84.9 hm3/yr during dry and wet periods, respectively (with values of nitrate discharge of 11.4-11.8 Mkg/yr in the absence of aquifer sink terms, e.g., leakage to deeper aquifers and pumping from groundwater wells). The predictive capabilities of the calibrated model can be used to test the impact of different integrated management scenarios on the surface-groundwater dynamics of the catchment. Three plausible management scenarios are proposed that include localized and distributed groundwater pumping (drains and groundwater wells, respectively). Results show the effectiveness of the scenarios in reducing the groundwater and nitrate discharge into the lagoon. The disadvantages of the proposed scenarios, including potential seawater intrusion, need to be balanced with their relative merits for the sustainable development of the region and the survival of the Mar Menor ecosystem. The modelling approach proposed provides a valuable tool for the integrated and holistic management of the Campo de Cartagena-Mar Menor catchment and should be of great interest to similar hydrological systems with high ecological value.

Eigenvector centrality for geometric and topological characterization of porous media.

Solving flow and transport through complex geometries such as porous media is computationally difficult. Such calculations usually involve the solution of a system of discretized differential equations, which could lead to extreme computational cost depending on the size of the domain and the accuracy of the model. Geometric simplifications like pore networks, where the pores are represented by nodes and the pore throats by edges connecting pores, have been proposed. These models, despite their ability to preserve the connectivity of the medium, have difficulties capturing preferential paths (high velocity) and stagnation zones (low velocity), as they do not consider the specific relations between nodes. Nonetheless, network theory approaches, where a complex network is a graph, can help to simplify and better understand fluid dynamics and transport in porous media. Here we present an alternative method to address these issues based on eigenvector centrality, which has been corrected to overcome the centralization problem and modified to introduce a bias in the centrality distribution along a particular direction to address the flow and transport anisotropy in porous media. We compare the model predictions with millifluidic transport experiments, which shows that, albeit simple, this technique is computationally efficient and has potential for predicting preferential paths and stagnation zones for flow and transport in porous media. We propose to use the eigenvector centrality probability distribution to compute the entropy as an indicator of the "mixing capacity" of the system.

Geo-material microfluidics at reservoir conditions for subsurface energy resource applications.

Microfluidic investigations of flow and transport in porous and fractured media have the potential to play a significant role in the development of future subsurface energy resource technologies. However, the majority of experimental systems to date are limited in applicability due to operating conditions and/or the use of engineered material micromodels. We have developed a high pressure and temperature microfluidic experimental system that allows for direct observations of flow and transport within geo-material micromodels (e.g. rock, cement) at reservoir conditions. In this manuscript, we describe the experimental system, including our novel micromodel fabrication method that works in both geo- and engineered materials and utilizes 3-D tomography images of real fractures as micromodel templates to better represent the pore space and fracture geometries expected in subsurface formations. We present experimental results that highlight the advantages of using real-rock micromodels and discuss potential areas of research that could benefit from geo-material microfluidic investigations. The experiments include fracture-matrix interaction in which water imbibes into the shale rock matrix from etched fractures, supercritical CO2 (scCO2) displacing brine in idealized and realistic fracture patterns, and three-phase flow involving scCO2-brine-oil.