There are a growing number of large-scale, complex hydrologic models that are capable of simulating integrated surface and subsurface flow. Many are coupled to land-surface energy balance models, biogeochemical and ecological process models, and atmospheric models. Although they are being increasingly applied for hydrologic prediction and environmental understanding, very little formal verification and/or benchmarking of these models has been performed. Here we present the results of an intercomparison study of seven coupled surface-subsurface models based on a series of benchmark problems. All the models simultaneously solve adapted forms of the Richards and shallow water equations, based on fully 3-D or mixed (1-D vadose zone and 2-D groundwater) formulations for subsurface flow and 1-D (rill flow) or 2-D (sheet flow) conceptualizations for surface routing. A range of approaches is used for the solution of the coupled equations, including global implicit, sequential iterative, and asynchronous linking, and various strategies are used to enforce flux and pressure continuity at the surface-subsurface interface. The simulation results show good agreement for the simpler test cases, while the more complicated test cases bring out some of the differences in physical process representations and numerical solution approaches between the models. Benchmarks with more traditional runoff generating mechanisms, such as excess infiltration and saturation, demonstrate more agreement between models, while benchmarks with heterogeneity and complex water table dynamics highlight differences in model formulation. In general, all the models demonstrate the same qualitative behavior, thus building confidence in their use for hydrologic applications.

Emphasizing the physical intricacies of integrated hydrology and feedbacks in simulating connected, variably saturated groundwater‐surface water systems, the Integrated Hydrologic Model Intercomparison Project initiated a second phase (IH‐MIP2), increasing the complexity of the benchmarks of the first phase. The models that took part in the intercomparison were ATS, Cast3M, CATHY, GEOtop, HydroGeoSphere, MIKE‐SHE, and ParFlow. IH‐MIP2 benchmarks included a tilted v‐catchment with 3‐D subsurface; a superslab case expanding the slab case of the first phase with an additional horizontal subsurface heterogeneity; and the Borden field rainfall‐runoff experiment. The analyses encompassed time series of saturated, unsaturated, and ponded storages, as well as discharge. Vertical cross sections and profiles were also inspected in the superslab and Borden benchmarks. An analysis of agreement was performed including systematic and unsystematic deviations between the different models. Results show generally good agreement between the different models, which lends confidence in the fundamental physical and numerical implementation of the governing equations in the different models. Differences can be attributed to the varying level of detail in the mathematical and numerical representation or in the parameterization of physical processes, in particular with regard to ponded storage and friction slope in the calculation of overland flow. These differences may become important for specific applications such as detailed inundation modeling or when strong inhomogeneities are present in the simulation domain.

[1] The accurate characterization of fractured geologic medium, imaging of fracture patterns and their connectivity have been a challenge for decades. Recently, hydraulic tomography has been proposed as a new method for imaging the hydraulic conductivity (K) and specific storage (S s ) distributions of fractured geologic media. While encouraging results have been obtained in the field, the method has not been rigorously assessed in a controlled laboratory setting. In this study, we assess the performance of transient hydraulic tomography (THT) in a fractured dolomitic rock block. The block is characterized through flow-through tests and multiple pumping tests. The pumping test data were then analyzed with the THT code of Zhu and Yeh (2005) to image the fracture patterns and their connectivity through the delineation of K and S s distributions (or tomograms). Results show that the THT analysis of pumping tests yields high-K and low-S s zones that capture the fracture pattern and their connectivity quite well and those patterns become more vivid as additional pumping test data are added to the inverse model. The performance of the estimated K and S s tomograms are then assessed by: (1) comparing the tomograms obtained from synthetic to real data; (2) comparing the tomograms from two different pumping configurations; (3) comparing the estimated geometric mean of the hydraulic conductivity (K G ) from the K tomogram to the effective hydraulic conductivity (K eff ) estimated from the flow-through tests; and (4) predicting five independent pumping tests not used in the construction of the K and S s tomograms. The performance assessment of the K and S s tomograms reveals that THT is able to image high-K and low-S s zones that correspond to fracture locations in the fractured rock block and that the tomograms can be used to predict drawdowns from pumping tests not used in the construction of the tomograms with reasonable fidelity.

Crosstalk between cancer cells and carcinoma-associated fibroblasts (CAFs) has earned recognition as an interaction that plays a pivotal role in carcinogenesis. Thus, we attempted to clarify whether increase in the level of CAFs promotes cancer progression by proportionally enhancing the interaction between cancer cells and CAFs. We first analyzed clinical correlation between the levels of fibroblasts and cancer progression and found that the level of CAFs made a noticeable difference on the prognosis of patients with oral squamous cell carcinoma (OSCC). In vivo animal study also demonstrated that tumor volume depended on the dose of CAFs that was co-injected with OSCC cells. The same tendency was observed in an in vitro study. We also found that interleukin-1α (IL-1α) secreted from OSCC cells had dual effects on CAFs: IL-1α not only promoted the proliferation of CAFs but also upregulated the secretion of cytokines in CAFs such as CCL7, CXCL1, and IL-8. The induction activity of cytokine secretion by IL-1α surpassed that of proliferation in OSCC cells. In summary, we unraveled an important interactive mechanism of carcinogenesis: IL-1α released from carcinoma stimulates the proliferation of CAFs and the simultaneous increase in cytokine secretion from CAFs promotes cancer progression in human OSCC. On the basis of these findings, we propose that the level of CAFs is eligible for being selected as a prognostic factor that will be useful in routine diagnosis. We also propose that blockage of reciprocal interaction between cancer cells and CAFs will provide an insight for developing novel chemotherapeutic strategy.

No abstract