Field-Calibrated Simplification of the F(x) Correction in Subsurface Drainage Design for Urban Heritage Sites | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Field-Calibrated Simplification of the F(x) Correction in Subsurface Drainage Design for Urban Heritage Sites Ahmed M. Noureldin This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-7906837/v1 This work is licensed under a CC BY 4.0 License Status: Under Review Version 1 posted 11 You are reading this latest preprint version Abstract Rising groundwater levels pose a significant threat to urban heritage sites, highlighting the need for efficient subsurface drainage systems. The Hooghoudt equation is a foundational tool for such designs, yet its practical use is hindered by the complexity of computing the correction term F(x), which is critical in urban environments with irregular geometries. This study presents a two-phase investigation. The first phase involves a field-calibrated analysis of three historic mosques in Cairo, applying a refined 10-term series for F(x). Findings reveal that omitting this term can lead to discharge errors exceeding 40%. In the second phase, a simplified empirical formula for F(x) is developed using data from over 22,000 simulation scenarios. A hybrid exponential decay model delivers outstanding accuracy (R² >0.9999; maximum error < 0.01%). The resulting three-parameter equation, F(x) ≈ 5.432 + 2.832⋅exp(− 6.414x), offers a practical alternative to complex series computations. Accompanying tools, such as nomographs and lookup tables enable rapid and reliable drainage planning to support the preservation of heritage sites vulnerable to rising groundwater. Subsurface Drainage Hooghoudt Equation F(x) Correction Empirical Modeling Heritage Site Conservation Field Calibration Monte Carlo Simulations and Design Nomograph Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 1. Introduction 1.1. The Heritage Conservation and Groundwater Challenge Historic urban centers, such as Cairo, are repositories of invaluable cultural heritage. However, these sites are increasingly under threat from environmental pressures, most notably rising groundwater levels [ 1 , 2 ]. The combination of urbanization, leaking infrastructure, and changes in regional hydrology has led to a persistent rise in the water table, which can cause catastrophic damage to the foundations and structures of ancient buildings. Effective management of the water table through engineered subsurface drainage is therefore not just a matter of geotechnical engineering, but a cornerstone of modern heritage conservation. 1.2. The Hooghoudt Equation and the F(x) Problem The Hooghoudt equation has long been a fundamental tool in drainage engineering, providing a steady-state approximation of the relationship between drain spacing, water table height, and soil hydraulic properties [ 3 ]. A critical, yet often overlooked, component of this equation is the F(x) correction term. This term accounts for the convergence of groundwater flow lines as they approach the drains, a phenomenon that is particularly pronounced in the constrained geometries typical of urban settings. The accurate determination of F(x) is essential for calculating the correct equivalent depth ( d ) of the impermeable layer, which in turn governs the entire drainage calculation. The challenge lies in the fact that the F(x) term is defined by a complex, slowly converging infinite series. Its direct calculation is cumbersome and not conducive to the rapid iteration required in practical design workflows. This computational barrier has led to a significant gap between drainage theory and common engineering practice [ 4 , 5 ]. 1.3. Current Practice and the Knowledge Gap In practice, designers often resort to simplified methods to deal with the F(x) term. These include using pre-computed nomographs [ 6 ], which may not be accurate for site-specific conditions, or, more alarmingly, neglecting the term altogether. While this simplifies the calculation, it introduces significant and often unquantified errors into the design, potentially leading to under-designed drainage systems that fail to protect the heritage structure. There is a clear knowledge gap in the literature regarding the magnitude of the error introduced by these simplifications, particularly in the context of urban heritage sites. Furthermore, there is a lack of practical, validated alternatives to the full series calculation that would allow engineers to maintain accuracy without sacrificing simplicity. 1.4. Research Objectives This study aims to bridge this gap through a comprehensive, two-part investigation. The primary objectives are: (1) To quantify the impact of the F(x) term on drainage design accuracy through a field-calibrated analysis of three historic mosques in Cairo, using a fully converged F(x) series calculation. (2) To develop and validate a simplified, robust, and highly accurate empirical formulation for the F(x) term that eliminates the need for complex series calculations. (3) To provide a set of practical design tools , including nomographs and lookup tables, to facilitate the immediate application of the empirical formulation in engineering practice. By pursuing these objectives, this paper presents a complete narrative, from problem identification and quantification in the field to the development and delivery of a practical, validated solution. 2. Methodology The methodological framework of this study integrates field-based calibration, analytical derivation, probabilistic modeling, and empirical regression to ensure both scientific rigor and practical applicability. The work was structured into four main components: (1) field data collection and site characterization, (2) theoretical computation of the converged F(x) correction term, (3) field-calibrated discharge and uncertainty analysis, and (4) development of a simplified empirical formulation. 2.1. Study Sites and Field Data Collection Three heritage mosques in Cairo—Amr Ibn al-As, Sayyida Ruqayya, and Hassan Pasha Taher were selected as representative case studies due to their historical and architectural significance, and their documented vulnerability to rising groundwater levels [1,2]. Field investigations included piezometric monitoring of the groundwater table and geotechnical characterization to determine soil stratigraphy and hydraulic conductivity. The procedures followed standard ILRI and FAO methodologies for subsurface drainage evaluation [3,4,10]. The resulting dataset was used to calibrate the hydraulic parameters (K, L, Di, Dd) for each site. This field-calibration process aligns with the approach adopted by Moriasi et al. [11] for validating drainage equations through observed field responses. 2.2. Theoretical Framework: Converged F(x) Computation The theoretical foundation is based on the Hooghoudt equation [5], which provides a steady-state relationship between drain spacing (L), hydraulic conductivity (K), and water table height (m). A critical component of this equation is the F(x) correction term , which accounts for the convergence of flow lines near the drain. The term is defined by an infinite series that converges slowly, as documented in early analytical works by Hooghoudt [5] and later improved by Van der Molen and Wesseling [6].In this study, a 10-term series expansion was adopted to ensure full numerical convergence of F(x), following the convergence criteria proposed by Van der Molen and Wesseling [6] and the design recommendations of ILRI [3]. This converged F(x) serves as the theoretical benchmark for calculating equivalent depth (d), in line with the comparative analyses by Wesseling [7] and the improved nomographic approaches of Ritzema [3]. 2.3. Field-Calibrated Discharge and Uncertainty Analysis Using the field-calibrated hydraulic parameters obtained from the three mosque sites, the required drainage discharge was computed for two conditions: (1) including the converged F(x) correction term, and (2) neglecting it entirely. The resulting difference was used to quantify the design error introduced by simplified practices that ignore F(x).A sensitivity analysis was performed to assess the effect of key design parameters (L, K, and d) on calculated discharge, consistent with the analytical procedures outlined by Smedema and Rycroft [8]. The approach was validated against empirical observations and prior numerical comparisons of Hooghoudt-type equations [11]. To incorporate the uncertainty associated with field-measured parameters, a Monte Carlo simulation with 10,000 iterations was conducted, applying the probabilistic framework introduced by Metropolis and Ulam [12] and refined through modern sensitivity analysis methods [13]. The resulting discharge distributions were summarized as P10, P50, and P90 percentiles , enabling risk-based design envelopes that define expected, median, and conservative design values. This risk-informed design concept aligns with the guidelines of the U.S. Army Corps of Engineers for hydraulic systems under uncertainty [15], thereby introducing a statistically robust framework for drainage design in heritage environments. 2.4. Empirical F(x) Formulation Development To avoid the gap between theoretical accuracy and practical usability, an extensive simulation dataset of 22,500 unique scenarios was generated. Key geometric and hydraulic parameters drain depth (Dd), spacing (L), and aquifer thickness above the drain (Di−Dd) were varied systematically across ranges representative of Cairo-type conditions, as summarized in Table 1. Table 1: Parameter Ranges for Simulation Dataset Generation Parameter Description Min Max Points/Values Dd Drain Depth 0.5 3.0 15 L Drain Spacing 20 100 20 Di - Dd Aquifer Thickness above Drain 50 300 15 r Pipe Radius 0.05 0.15 [0.05, 0.076, 0.102, 0.127, 0.15] Several regression models were tested to fit the relationship between x and F(x), including power law, exponential, rational, polynomial, and hybrid forms. The fitting process employed non-linear least squares regression , following the optimization techniques presented by Draper and Smith [14]. The performance of each model was evaluated using the coefficient of determination (R²), root mean square error (RMSE), and maximum relative error metrics. To ensure immediate usability by practitioners, the empirical formulation was translated into design nomographs and lookup tables for direct determination of F(x) and equivalent depth (d). This approach continues the ILRI tradition of visual design aids [3] and the practical framework of Smedema and Rycroft [8], while incorporating modern, data-driven contours based on the hybrid regression model. The integration of theoretical precision with empirical simplicity provides an effective, field-ready design methodology particularly suited for the preservation of heritage sites in dense urban environments. The models were fitted using non-linear least squares regression, and the best model was selected based on the coefficient of determination (R²), root mean squared error (RMSE), and maximum relative error. 3. Results The results of this two-part investigation are presented sequentially, beginning with the field-calibrated analysis that establishes the significance of the F(x) term, followed by the development and validation of the simplified empirical formulation. 3.1. F(x) Convergence and Impact Quantification The theoretical F(x) term is an infinite series, and the first step was to establish the number of terms required for mathematical convergence. As shown in Figure 1, the value of F(x) converges rapidly, with the value stabilizing after approximately 5-6 odd terms. Using a 10-term expansion ensures a fully converged value for all subsequent calculations. The practical impact of neglecting this converged term is substantial. A direct comparison of discharge calculations for the three mosque sites—one with the 10-term F(x) and one without—revealed that neglecting the term leads to an underestimation of the required discharge by over 40%. This finding from the field-calibrated models highlights a critical flaw in common engineering practice and serves as the primary motivation for developing a practical alternative. 3.2. Field-Calibrated Sensitivity and Discharge Analysis A sensitivity analysis was conducted to understand the influence of key design parameters on the required discharge, using the field-calibrated data from the three mosque sites. The results, shown in Figure 2, indicate that the discharge is highly sensitive to drain spacing (L) and hydraulic conductivity (K), reinforcing the need for accurate parameter determination. The analysis shows a strong, non-linear relationship between discharge and key parameters like drain spacing and hydraulic conductivity. Furthermore, the analysis of different permeability scenarios (Figure 3) demonstrates the dramatic impact of the underlying geology. A shift from a conservative (low permeability) to a moderate permeability assumption for the deeper soil layers can amplify the required discharge by a factor of over 100, underscoring the importance of thorough geotechnical investigation. The required discharge is significantly higher under moderate permeability assumptions for the deeper soil layers. 3.3. Uncertainty Quantification and Probabilistic Design Envelopes To account for the inherent uncertainty in the field-measured parameters, a 10,000-iteration Monte Carlo simulation was performed for each of the three mosque sites. The results (Figure 4) reveal a wide probability distribution for the required discharge, indicating that a deterministic, single-value design is insufficient. The analysis reveals a broad range of potential discharge values, highlighting the need for a probabilistic design approach. Based on this probabilistic analysis, we developed design discharge envelopes for each site (Figure 5). These envelopes provide P10, P50, and P90 percentile values, allowing for a risk-based design approach. For critical heritage sites, designing for the P90 (90th percentile) discharge provides a 90% confidence level that the drainage system will be adequate, representing a far more robust approach than a single deterministic calculation. The envelopes provide a range of design values (P10, P50, P90) for a risk-based approach to drainage system design. 3.4. Development and Validation of the Empirical F(x) Formulation While the field analysis confirmed the importance of F(x), it did not solve the problem of its computational complexity. To address this, we developed an empirical formulation from a 22,500-scenario dataset. Five different regression models were tested, and the results are summarized in Table 2. The Hybrid Formulation, designed to mimic the known asymptotic behavior of the F(x) function, demonstrated vastly superior performance. Table 2: Performance Comparison of Five Empirical F(x) Formulations Formulation R² RMSE MAE Max Relative Error (%) Power Law 0.2249 0.01097 0.00566 1.760 Exponential Decay 0.0640 0.01206 0.00527 1.956 Rational Function 0.5335 0.00851 0.00514 1.295 Polynomial (2nd Order) 0.1614 0.01141 0.00549 1.848 Hybrid (Recommended) 0.99998 0.00005 0.00002 0.008 Figure 6 visually confirms this result, showing the near-perfect fit of the Hybrid model compared to the other four approaches. The Hybrid model (bottom center) clearly outperforms all other approaches, with data points aligning almost perfectly along the 1:1 line. The recommended empirical equation for the F(x) correction term is: F(x) ≈ 5.432 + 2.832 * exp (-6.414 * x) ,where x = (Di - Dd) / L . This simple, three-parameter equation achieves a coefficient of determination (R²) greater than 0.99998 and a maximum relative error of less than 0.01%. The detailed performance analysis in Figure 7 confirms the exceptional accuracy of this formulation. The model accurately traces the theoretical curve, and the residual errors are symmetrically distributed around zero and are negligibly small. 3.5. Practical Design Tools To facilitate the immediate use of this empirical formulation in engineering practice, we have developed a set of design tools. Figure 8 presents two design nomographs that allow for the graphical determination of F(x) and the resulting equivalent depth ( d ) based on the site-specific drain spacing (L) and aquifer thickness (Di - Dd). Table 3: Lookup Table for F(x) vs. Dimensionless Parameter x x = (Di-Dd)/L 0.5 1.0 1.5 2.0 3.0 5.0 10.0 F(X) Value 5.546 5.437 5.432 5.432 5.432 5.432 5.432 Table 4: F(x) and Equivalent Depth (d) for Common Heritage Site Configurations (r = 0.1016m) L (m) Di-Dd (m) x F(X) d (m) 30 100 3.33 5.432 10.73 40 150 3.75 5.432 14.69 50 200 4.00 5.432 18.78 60 150 2.50 5.432 22.98 80 100 1.25 5.433 31.58 These tools collectively eliminate the need for any complex calculations, allowing practitioners to select appropriate design parameters with confidence in the accuracy of the underlying hydraulic model. 4. Discussion This study successfully bridges the gap between complex drainage theory and practical engineering design for heritage conservation. The integrated findings, from field-calibrated error quantification to the development of a simple and robust empirical solution, provide a comprehensive framework for improving the standard of practice. 4.1. The Significance of F(x) Precision in Heritage Contexts The field-calibrated analysis of the three Cairo mosques provides unequivocal evidence that the F(x) correction term is not a minor theoretical detail, but a critical component of accurate drainage design. The finding that neglecting F(x) can lead to discharge calculation errors exceeding 40% has profound implications for heritage conservation. An under-designed drainage system, based on an erroneously low discharge estimate, would fail to adequately lower the water table, leaving the historic foundations vulnerable to continued moisture-related decay. The cost of such a design error is not merely financial; it is the potential for irreversible loss of cultural heritage. This finding alone should compel a re-evaluation of current engineering practices that treat the F(x) term as optional or that rely on overly simplified approximations. 4.2. The Empirical Solution: A Practical Path to Accuracy Having established the problem, the second part of this study provides the solution. The development of the empirical equation, F(x) ≈ 5.432 + 2.832 * exp (-6.414 * x) , effectively eliminates the primary barrier to accurate F(x) computation: its complexity. With an R² value exceeding 0.9999 and a maximum relative error of less than 0.01%, this simple three-parameter equation provides the accuracy of the full theoretical series with the simplicity of a basic calculator function. This allows for the rapid iteration and optimization of drainage designs, a crucial capability in the complex, multi-stakeholder environment of heritage projects.The accompanying design nomographs and lookup tables further enhance the practical utility of this work, providing multiple entry points for practitioners of varying technical backgrounds to access this enhanced accuracy. 4.3. Comparison with and Contribution to Existing Approaches This work builds upon a long history of drainage research. The foundational nomographs published by the International Institute for Land Reclamation and Improvement (ILRI) [ 6 ] have been invaluable for decades, but our empirical approach offers a more flexible, precise, and digitally-native alternative. Our work is also distinct from, but complementary to, the research by Shokri & Bardsley (2015) [ 7 ], who focused on corrections related to soil-water retention properties. Our study provides the essential geometric correction, and a powerful future approach would be to integrate both corrections into a single, comprehensive design tool. The classical work of Wesseling (1964) [ 8 ] and more recent evaluations by Moriasi et al. (2013) [ 9 ] have compared various drainage spacing formulas, but none have provided the combination of field validation and simplified empirical formulation presented here. The primary contribution of our work is the creation of a complete narrative: we not only identify and quantify a critical error in current practice using real-world field data, but we also provide the direct, validated solution to that problem. 4.4. Application to Cairo Heritage Sites and Beyond The probabilistic design envelopes developed from the Monte Carlo analysis (Fig. 5 ) offer a new standard for drainage design at critical heritage sites. Moving away from a single deterministic value to a risk-based P90 design provides a quantifiable level of confidence that the system will perform as intended, even with the inherent uncertainties in the input parameters. This approach is directly applicable to the ongoing conservation efforts at the three mosques studied and can serve as a model for other heritage sites in Cairo.While the empirical formulation was specifically calibrated for Cairo-type hydrogeological conditions, the methodology itself is broadly applicable. The process of generating a large simulation dataset and fitting a robust empirical model can be replicated for other regions with different geological characteristics, creating a new generation of site-specific, highly accurate design tools. 4.5. Limitations It is important to acknowledge the limitations of this study. The empirical formulation is validated for the parameter space defined in Table 1 and should be used with caution outside of these ranges. Furthermore, this study, like the Hooghoudt equation itself, is based on the assumption of steady-state flow and homogeneous, isotropic soil. In highly complex, stratified geological settings, or where transient effects are significant, numerical modeling (e.g., with MODFLOW or HYDRUS) remains the most appropriate tool for final design verification. 5. Conclusion This study has successfully addressed a critical gap between subsurface drainage theory and practical engineering for heritage conservation. By integrating a field-calibrated analysis with the development of a simplified empirical model, we have created a comprehensive framework that both quantifies a significant problem and provides its direct solution. The key findings are twofold. First, the field-calibrated analysis of three historic mosques in Cairo demonstrates that neglecting the Hooghoudt F(x) correction term can lead to an underestimation of required drainage discharge by over 40%, posing a significant risk to the long-term preservation of these invaluable structures. Second, we have shown that the computational complexity of the F(x) term can be eliminated without sacrificing accuracy. The developed empirical equation, F(x) ≈ 5.432 + 2.832 * exp (-6.414 * x) , replicates the theoretical series with near-perfect accuracy (R² >0.9999) and is presented with a suite of practical design tools, including nomographs and lookup tables.This work provides a clear path for engineers and conservationists to move beyond overly simplified or inaccurate methods and adopt a design approach that is both theoretically sound and practically achievable. The ability to rapidly and reliably calculate the required drainage capacity is a critical component in the ongoing effort to protect our shared cultural heritage from the persistent threat of rising groundwater. The methodology and tools presented here offer a significant step forward in that endeavor. Declarations Funding This research received no external funding Data availability The article contains all of the data generated or analyzed throughout the research. Declarations Ethics approval and consent to participate Ethical approval was not sought for the present study because there are no human/animal subjects in this article. Consent to publication Author approved the manuscript for publication. Competing interests The author declares no competing interests Consent to Participate declaration Not applicable Author Contribution Writing original draft, Methodology, Formal analysis, Data curation, Conceptualization. Acknowledgement I would like to thank my colleagues at the Geotechnical Engineering Research Institute and the Greater Cairo Wastewater Company for their valuable support and collaboration. References Abdel-Shafy, H. I., & Aly, R. O. (2012). Water issues in Egypt: Resources, pollution, and protection endeavors. Environmental Science and Pollution Research , 19(6), 2163–2176. https://doi.org/10.1007/s11356-012-0877-8 UNESCO World Heritage Centre. (2021). Historic Cairo. Retrieved from https://whc.unesco.org/en/list/89/ Ritzema, H. P. (Ed.). (1994). Drainage Principles and Applications (Vol. 16). International Institute for Land Reclamation and Improvement (ILRI) , Wageningen, The Netherlands. FAO. (2002). Drainage Design Manual. FAO Irrigation and Drainage Paper 61. Rome: Food and Agriculture Organization of the United Nations. Hooghoudt, S. B. (1940). General consideration of the problem of field drainage by parallel drains, ditches, watercourses, and channels. Institute for Land and Water Management Research, Wageningen, The Netherlands. Van der Molen, W. H., & Wesseling, J. (1991). A solution in closed form and a series solution to replace the tables for the thickness of the equivalent layer in Hooghoudt’s drain spacing formula. Agricultural Water Management , 19(1), 1–16. https://doi.org/10.1016/0378-3774(91)90061-T Wesseling, J. (1964). A comparison of the steady state drain spacing formulas of Hooghoudt and Kirkham in connection with design practice. Journal of Hydrology , 2(1), 25–32. https://doi.org/10.1016/0022-1694(64)90004-6 Smedema, L. K., & Rycroft, D. W. (1983). Land Drainage: Planning and Design of Agricultural Drainage Systems. Ithaca, NY: Cornell University Press. Shokri, A., & Bardsley, W. E. (2015). Enhancement of the Hooghoudt drain-spacing equation. Journal of Irrigation and Drainage Engineering , 141(6), 04014082. https://doi.org/10.1061/(ASCE)IR.1943-4774.0000835 FAO & ILRI. (2002–1994). Standardized Field Methods for Subsurface Drainage Design. Joint Technical Framework, Rome–Wageningen. Moriasi, D. N., et al. (2013). Evaluation of the Hooghoudt and Kirkham tile drain equations in the Soil and Water Assessment Tool to simulate tile flow and nitrate-nitrogen. Journal of Environmental Quality , 42(6), 1699–1710. https://doi.org/10.2134/jeq2013.01.0018 Metropolis, N., & Ulam, S. (1949). The Monte Carlo Method. Journal of the American Statistical Association , 44(247), 335–341. https://doi.org/10.1080/01621459.1949.10483310 Saltelli, A., Chan, K., & Scott, E. M. (2010). Sensitivity Analysis in Practice: A Guide to Assessing Scientific Models. Chichester: John Wiley & Sons. Draper, N. R., & Smith, H. (1998). Applied Regression Analysis (3rd ed.). New York: Wiley-Interscience. U.S. Army Corps of Engineers (USACE). (2019). Engineering and Design: Risk-Informed Decision Making for Hydrologic Engineering (EM 1110-2-1619). Washington, D.C. Additional Declarations No competing interests reported. Cite Share Download PDF Status: Under Review Version 1 posted Editorial decision: Revision requested 14 Dec, 2025 Reviews received at journal 13 Dec, 2025 Reviewers agreed at journal 26 Nov, 2025 Reviewers agreed at journal 25 Nov, 2025 Reviews received at journal 22 Nov, 2025 Reviewers agreed at journal 11 Nov, 2025 Reviewers invited by journal 05 Nov, 2025 Editor invited by journal 04 Nov, 2025 Editor assigned by journal 27 Oct, 2025 Submission checks completed at journal 27 Oct, 2025 First submitted to journal 20 Oct, 2025 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-7906837","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":545452905,"identity":"dbcb0a4d-bf04-45f8-8237-b7077c26de97","order_by":0,"name":"Ahmed M. Noureldin","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABIElEQVRIiWNgGAWjYLACxgYgcQDMtGBsY2A+ABEhUosEUAtbAolaGhh4DPBq4Zduv8Dwc8c2Ob7bzcekKyokZPuke75J/NxhI8fAfvjoBixaJOecKWDsPXPbWPLOsTTJM2ckjNtkzm6T7D2TZszAk5Z2A4sWgxs5CQy8bbcTN9zIMZNsbJNIbJPI3SbB23Y4sUGCxwybFnugFsa/bbfrIVr+gbTkPJP8i0eLgUT6AWagLQkGYC0NYC1s0vhskbiRw8As23bbcOaNtGTLhmNAv0ikGVvLtqUZs+HwC/+M9AeMb9tuy/PdSD54s6HGRnb+jOSHN9+22cjxsx8+hk0LAwOP+Q90IRYJEMmGVTkIsD/AEGL+gFP1KBgFo2AUjEQAAF30au4Q4VVQAAAAAElFTkSuQmCC","orcid":"","institution":"Housing and Building National Research Center (HBRC)","correspondingAuthor":true,"prefix":"","firstName":"Ahmed","middleName":"M.","lastName":"Noureldin","suffix":""}],"badges":[],"createdAt":"2025-10-20 14:53:10","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-7906837/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-7906837/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":96247158,"identity":"1ffeb233-eddb-4794-aea9-1c62e14d84ad","added_by":"auto","created_at":"2025-11-19 07:27:13","extension":"docx","order_by":0,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":7523879,"visible":true,"origin":"","legend":"","description":"","filename":"FieldCalibratedSimplificationoftheFxCorrectioninSubsurfaceDrainageDesignforUrbanHeritageSites.261025.docx","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/aa576024c3b8a54c3e35b1a8.docx"},{"id":96246163,"identity":"4d3fd6e7-2f74-4ed5-a556-56541e3ca811","added_by":"auto","created_at":"2025-11-19 07:24:55","extension":"json","order_by":1,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":3253,"visible":true,"origin":"","legend":"","description":"","filename":"7b87d42866a64f3797a0924bd9fb2755.json","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/65bd260850f0ee519c49ac69.json"},{"id":96056997,"identity":"a9bed29f-4d81-43f8-9528-7f67211f21ef","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"xml","order_by":2,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":64209,"visible":true,"origin":"","legend":"","description":"","filename":"7b87d42866a64f3797a0924bd9fb27551enriched.xml","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/0ec3be319f56e17fd073cd6c.xml"},{"id":96057000,"identity":"cab7141d-7938-4ac1-b852-d3eff8e56a1e","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"png","order_by":3,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":459665,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage1.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/ea418364664afdc32b73b8db.png"},{"id":96057003,"identity":"b35bf8ff-b5ec-4fb7-9768-ab94a8e498c5","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"png","order_by":4,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":618877,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage2.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/7e51dff88f3019319dc2b78a.png"},{"id":96057009,"identity":"126bd78b-0f2d-4d8d-9f01-3c7295128b29","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"png","order_by":5,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":269568,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage3.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/9ed2e7e14ed2833e9ced02d9.png"},{"id":96057011,"identity":"cdd7e0a2-e206-4c83-a5fc-8fb110d46063","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"png","order_by":6,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":4592915,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage4.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/db01abe8ba566f7dcc928c45.png"},{"id":96246282,"identity":"34e820cc-2c77-4b31-bf90-11d7a5cecd0e","added_by":"auto","created_at":"2025-11-19 07:25:17","extension":"png","order_by":7,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":216576,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage5.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/516793d13374faf4b7420eec.png"},{"id":96246985,"identity":"cdedf2e2-9124-4624-894f-171988321456","added_by":"auto","created_at":"2025-11-19 07:26:59","extension":"png","order_by":8,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":429520,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage6.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/07dd7f8d95b647163d9b6d3f.png"},{"id":96248253,"identity":"d7a14d69-53e2-4532-a310-63789b3561c4","added_by":"auto","created_at":"2025-11-19 07:28:12","extension":"png","order_by":9,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":490335,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage7.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/7ff095244d46ac8b7b011e65.png"},{"id":96057006,"identity":"52c2530c-a6f0-4548-9eae-01c885dec949","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"png","order_by":10,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":409863,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage8.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/2a97f826d442e85bf1b60eed.png"},{"id":96247907,"identity":"aac13dc1-a989-45f9-a45d-c66a8027e899","added_by":"auto","created_at":"2025-11-19 07:27:50","extension":"png","order_by":11,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":191708,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage1.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/412a91a950614df904fd3c99.png"},{"id":96057020,"identity":"0e64bf99-1628-42f1-a143-9a4cf634f583","added_by":"auto","created_at":"2025-11-17 07:59:12","extension":"png","order_by":12,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":163005,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage2.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/d7275c4cd779f564bc4cc6c3.png"},{"id":96057021,"identity":"4e29cc1f-df0d-492e-93d1-2a05f000a130","added_by":"auto","created_at":"2025-11-17 07:59:12","extension":"png","order_by":13,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":110133,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage3.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/0ab0363138358795175bb5ff.png"},{"id":96057023,"identity":"f2f498a0-540f-4e6a-a85d-880ff938324c","added_by":"auto","created_at":"2025-11-17 07:59:12","extension":"png","order_by":14,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":940417,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage4.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/c235792048739dc49986bf15.png"},{"id":96057008,"identity":"a6d63428-47ba-4aeb-a247-296f0d2e9430","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"png","order_by":15,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":84927,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage5.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/98ef8c75a5f2d9d505835be1.png"},{"id":96057017,"identity":"16a198a5-5270-4f2e-a521-cbc46ffd89fa","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"png","order_by":16,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":231822,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage6.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/80c3f2811309e55452a1c010.png"},{"id":96057018,"identity":"5d2f2a51-dd1b-4954-9f2e-532834d9963b","added_by":"auto","created_at":"2025-11-17 07:59:12","extension":"png","order_by":17,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":160179,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage7.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/849ec0a8e0b501034ddb9f4d.png"},{"id":96246416,"identity":"e261d0e5-0d61-4be3-b209-456027f27745","added_by":"auto","created_at":"2025-11-19 07:25:56","extension":"png","order_by":18,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":112746,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage8.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/39709ef2013c6e9152f769f4.png"},{"id":96246372,"identity":"d0720548-58d7-44e1-a4d7-038e5fd17f5c","added_by":"auto","created_at":"2025-11-19 07:25:44","extension":"xml","order_by":19,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":63309,"visible":true,"origin":"","legend":"","description":"","filename":"7b87d42866a64f3797a0924bd9fb27551structuring.xml","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/ee45858df9429ad39beed3cd.xml"},{"id":96057004,"identity":"10591755-4d5b-46cb-bc3a-e61f8440a1eb","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"html","order_by":20,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":69688,"visible":true,"origin":"","legend":"","description":"","filename":"earlyproof.html","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/bc16f8da06f751188281dbc0.html"},{"id":96056995,"identity":"5c6a1895-41f6-443b-8f22-4c3c64b64543","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":459665,"visible":true,"origin":"","legend":"\u003cp\u003eF(X) Series Convergence and Impact on Equivalent Depth Accuracy.\u003c/p\u003e\n\u003cp\u003e(A) The F(x) value converges as the number of series terms increases. (B) The relative error in the equivalent depth \u003cu\u003ed\u003c/u\u003e decreases as more terms are included.\u003c/p\u003e","description":"","filename":"floatimage1.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/aa0bf95e8bf574672e5ac739.png"},{"id":96056994,"identity":"953016a3-e78f-4368-be49-3f4dd4ea02e2","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":618877,"visible":true,"origin":"","legend":"\u003cp\u003eDischarge Sensitivity to Design Parameters.\u003c/p\u003e","description":"","filename":"floatimage2.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/4e09823218a84b6fff349893.png"},{"id":96248362,"identity":"4ff9ea2b-2a8e-4f3f-a2ec-24d91ff6dfc3","added_by":"auto","created_at":"2025-11-19 07:28:22","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":269568,"visible":true,"origin":"","legend":"\u003cp\u003eComparative Discharge Analysis.\u003c/p\u003e","description":"","filename":"floatimage3.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/47fb97d34af8008ba8f7b2f5.png"},{"id":96056999,"identity":"609f53ed-98e1-44fd-b1f2-231b3275dcc6","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":4592915,"visible":true,"origin":"","legend":"\u003cp\u003eMonte Carlo Uncertainty Analysis for the Three Mosque Sites.\u003c/p\u003e","description":"","filename":"floatimage4.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/6b10404c594e8b3ea77c11d8.png"},{"id":96057002,"identity":"9355ed50-b7f4-428b-b2db-cb78abe5938c","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":216576,"visible":true,"origin":"","legend":"\u003cp\u003eProbabilistic Design Discharge Envelopes.\u003c/p\u003e","description":"","filename":"floatimage5.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/8eac56df9d90fdc3bb200c8e.png"},{"id":96057010,"identity":"3e19890b-6a52-4fe4-bb1b-76dff5944f58","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":429520,"visible":true,"origin":"","legend":"\u003cp\u003eComparison of Predicted vs. Theoretical F(x) for Five Empirical Formulations.\u003c/p\u003e","description":"","filename":"floatimage6.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/675de2ab1110bfffda71b561.png"},{"id":96057013,"identity":"fbfe11ad-9014-4b3c-bc68-8a1b5aa94247","added_by":"auto","created_at":"2025-11-17 07:59:11","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":490335,"visible":true,"origin":"","legend":"\u003cp\u003eDetailed Performance Analysis of the Recommended Hybrid Formulation.\u003c/p\u003e\n\u003cp\u003e(A) The empirical equation perfectly matches the theoretical values. (B) Residual errors are minimal and randomly distributed. (C) The distribution of relative error is centered at zero with a very small standard deviation. (D) The predicted vs. actual plot confirms the near-perfect fit.\u003c/p\u003e","description":"","filename":"floatimage7.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/8a9affe85f9b80def5b56791.png"},{"id":96057019,"identity":"b6560453-0643-4c77-974d-eb31126c1bcf","added_by":"auto","created_at":"2025-11-17 07:59:12","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":409863,"visible":true,"origin":"","legend":"\u003cp\u003eDesign Nomographs for F(x) and Equivalent Depth (d).\u003c/p\u003e\n\u003cp\u003e(A) F(x) can be determined from the aquifer thickness and drain spacing. (B) and (C) provide contour plots for F(x) and the resulting equivalent depth \u003cu\u003ed\u003c/u\u003e for a standard 8-inch pipe. Furthermore, we have compiled lookup tables for quick reference. Table 3 provides pre-calculated F(x) values for a range of \u003cu\u003ex\u003c/u\u003evalues, while Table 4 provides F(x) and \u003cu\u003ed\u003c/u\u003e values for common heritage site configurations in Cairo, allowing for rapid preliminary design.\u003c/p\u003e","description":"","filename":"floatimage8.png","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/0d1b0b7dadcf0a243d15a65d.png"},{"id":96256274,"identity":"2c804d80-8012-42c4-a170-f6f5e43d9d13","added_by":"auto","created_at":"2025-11-19 07:49:50","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":8381315,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-7906837/v1/6e4554da-c00d-4857-997f-3ceae3078f32.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":"Field-Calibrated Simplification of the F(x) Correction in Subsurface Drainage Design for Urban Heritage Sites","fulltext":[{"header":"1. Introduction","content":"\u003cdiv id=\"Sec2\" class=\"Section2\"\u003e\u003ch2\u003e1.1. The Heritage Conservation and Groundwater Challenge\u003c/h2\u003e\u003cp\u003eHistoric urban centers, such as Cairo, are repositories of invaluable cultural heritage. However, these sites are increasingly under threat from environmental pressures, most notably rising groundwater levels [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e, \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e]. The combination of urbanization, leaking infrastructure, and changes in regional hydrology has led to a persistent rise in the water table, which can cause catastrophic damage to the foundations and structures of ancient buildings. Effective management of the water table through engineered subsurface drainage is therefore not just a matter of geotechnical engineering, but a cornerstone of modern heritage conservation.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e\u003ch2\u003e1.2. The Hooghoudt Equation and the F(x) Problem\u003c/h2\u003e\u003cp\u003eThe Hooghoudt equation has long been a fundamental tool in drainage engineering, providing a steady-state approximation of the relationship between drain spacing, water table height, and soil hydraulic properties [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]. A critical, yet often overlooked, component of this equation is the F(x) correction term. This term accounts for the convergence of groundwater flow lines as they approach the drains, a phenomenon that is particularly pronounced in the constrained geometries typical of urban settings. The accurate determination of F(x) is essential for calculating the correct equivalent depth (\u003cspan type=\"Underline\" class=\"Underline\" name=\"Emphasis\"\u003ed\u003c/span\u003e) of the impermeable layer, which in turn governs the entire drainage calculation. The challenge lies in the fact that the F(x) term is defined by a complex, slowly converging infinite series. Its direct calculation is cumbersome and not conducive to the rapid iteration required in practical design workflows. This computational barrier has led to a significant gap between drainage theory and common engineering practice [\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e, \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e].\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec4\" class=\"Section2\"\u003e\u003ch2\u003e1.3. Current Practice and the Knowledge Gap\u003c/h2\u003e\u003cp\u003eIn practice, designers often resort to simplified methods to deal with the F(x) term. These include using pre-computed nomographs [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e], which may not be accurate for site-specific conditions, or, more alarmingly, neglecting the term altogether. While this simplifies the calculation, it introduces significant and often unquantified errors into the design, potentially leading to under-designed drainage systems that fail to protect the heritage structure.\u003c/p\u003e\u003cp\u003eThere is a clear knowledge gap in the literature regarding the magnitude of the error introduced by these simplifications, particularly in the context of urban heritage sites. Furthermore, there is a lack of practical, validated alternatives to the full series calculation that would allow engineers to maintain accuracy without sacrificing simplicity.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec5\" class=\"Section2\"\u003e\u003ch2\u003e1.4. Research Objectives\u003c/h2\u003e\u003cp\u003eThis study aims to bridge this gap through a comprehensive, two-part investigation. The primary objectives are: (1) \u003cb\u003eTo quantify the impact of the F(x) term\u003c/b\u003e on drainage design accuracy through a field-calibrated analysis of three historic mosques in Cairo, using a fully converged F(x) series calculation. (2) \u003cb\u003eTo develop and validate a simplified, robust, and highly accurate empirical formulation\u003c/b\u003e for the F(x) term that eliminates the need for complex series calculations. (3) \u003cb\u003eTo provide a set of practical design tools\u003c/b\u003e, including nomographs and lookup tables, to facilitate the immediate application of the empirical formulation in engineering practice. By pursuing these objectives, this paper presents a complete narrative, from problem identification and quantification in the field to the development and delivery of a practical, validated solution.\u003c/p\u003e\u003c/div\u003e"},{"header":"2. Methodology","content":"\u003cp\u003eThe methodological framework of this study integrates field-based calibration, analytical derivation, probabilistic modeling, and empirical regression to ensure both scientific rigor and practical applicability. The work was structured into four main components: (1) field data collection and site characterization, (2) theoretical computation of the converged F(x) correction term, (3) field-calibrated discharge and uncertainty analysis, and (4) development of a simplified empirical formulation.\u003c/p\u003e\n\u003cp\u003e2.1. Study Sites and Field Data Collection\u003c/p\u003e\n\u003cp\u003eThree heritage mosques in Cairo\u0026mdash;Amr Ibn al-As, Sayyida Ruqayya, and Hassan Pasha Taher were selected as representative case studies due to their historical and architectural significance, and their documented vulnerability to rising groundwater levels [1,2].\u003cbr\u003eField investigations included \u003cstrong\u003epiezometric monitoring\u003c/strong\u003e of the groundwater table and \u003cstrong\u003egeotechnical characterization\u003c/strong\u003e to determine soil stratigraphy and hydraulic conductivity. The procedures followed standard ILRI and FAO methodologies for subsurface drainage evaluation [3,4,10]. The resulting dataset was used to calibrate the hydraulic parameters (K, L, Di, Dd) for each site. This field-calibration process aligns with the approach adopted by Moriasi et al. [11] for validating drainage equations through observed field responses.\u003c/p\u003e\n\u003cp\u003e2.2. Theoretical Framework: Converged F(x) Computation\u003c/p\u003e\n\u003cp\u003eThe theoretical foundation is based on the \u003cstrong\u003eHooghoudt equation\u003c/strong\u003e [5], which provides a steady-state relationship between drain spacing (L), hydraulic conductivity (K), and water table height (m). A critical component of this equation is the \u003cstrong\u003eF(x) correction term\u003c/strong\u003e, which accounts for the convergence of flow lines near the drain. The term is defined by an infinite series that converges slowly, as documented in early analytical works by Hooghoudt [5] and later improved by Van der Molen and Wesseling [6].In this study, a \u003cstrong\u003e10-term series expansion\u003c/strong\u003e was adopted to ensure full numerical convergence of F(x), following the convergence criteria proposed by Van der Molen and Wesseling [6] and the design recommendations of ILRI [3]. This converged F(x) serves as the theoretical benchmark for calculating equivalent depth (d), in line with the comparative analyses by Wesseling [7] and the improved nomographic approaches of Ritzema [3].\u003c/p\u003e\n\u003cp\u003e2.3. Field-Calibrated Discharge and Uncertainty Analysis\u003c/p\u003e\n\u003cp\u003eUsing the field-calibrated hydraulic parameters obtained from the three mosque sites, the required drainage discharge was computed for two conditions:\u003cbr\u003e\u0026nbsp;(1) including the converged F(x) correction term, and (2) neglecting it entirely.\u003cbr\u003eThe resulting difference was used to quantify the design error introduced by simplified practices that ignore F(x).A \u003cstrong\u003esensitivity analysis\u003c/strong\u003e was performed to assess the effect of key design parameters (L, K, and d) on calculated discharge, consistent with the analytical procedures outlined by Smedema and Rycroft [8]. The approach was validated against empirical observations and prior numerical comparisons of Hooghoudt-type equations [11]. To incorporate the uncertainty associated with field-measured parameters, a \u003cstrong\u003eMonte Carlo simulation\u003c/strong\u003e with 10,000 iterations was conducted, applying the probabilistic framework introduced by Metropolis and Ulam [12] and refined through modern sensitivity analysis methods [13]. The resulting discharge distributions were summarized as \u003cstrong\u003eP10, P50, and P90 percentiles\u003c/strong\u003e, enabling risk-based design envelopes that define expected, median, and conservative design values. This risk-informed design concept aligns with the guidelines of the U.S. Army Corps of Engineers for hydraulic systems under uncertainty [15], thereby introducing a statistically robust framework for drainage design in heritage environments.\u003c/p\u003e\n\u003cp\u003e2.4. Empirical F(x) Formulation Development\u003c/p\u003e\n\u003cp\u003eTo avoid the gap between theoretical accuracy and practical usability, an extensive\u003cstrong\u003e\u0026nbsp;\u003cstrong\u003esimulation dataset of 22,500 unique scenarios\u003c/strong\u003e\u003c/strong\u003e was generated. Key geometric and hydraulic parameters drain depth (Dd), spacing (L), and aquifer thickness above the drain (Di\u0026minus;Dd) were varied systematically across ranges representative of Cairo-type conditions, as summarized in Table 1.\u003c/p\u003e\n\u003cp\u003eTable 1: Parameter Ranges for Simulation Dataset Generation\u0026nbsp;\u003c/p\u003e\n\u003ctable border=\"1\" cellspacing=\"0\" cellpadding=\"0\" width=\"100%\" class=\"fr-table-selection-hover\"\u003e\n \u003cthead\u003e\n \u003ctr\u003e\n \u003ctd\u003e\n \u003cp\u003eParameter\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003eDescription\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003eMin\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003eMax\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003ePoints/Values\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/thead\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd\u003e\n \u003cp\u003eDd\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003eDrain Depth\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003e0.5\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003e3.0\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003e15\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd\u003e\n \u003cp\u003eL\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003eDrain Spacing\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003e20\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003e100\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003e20\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd\u003e\n \u003cp\u003eDi - Dd\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003eAquifer Thickness above Drain\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003e50\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003e300\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003e15\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd\u003e\n \u003cp\u003er\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003ePipe Radius\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003e0.05\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003e0.15\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd\u003e\n \u003cp\u003e[0.05, 0.076, 0.102, 0.127, 0.15]\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003eSeveral regression models were tested to fit the relationship between x and F(x), including power law, exponential, rational, polynomial, and hybrid forms. The fitting process employed \u003cstrong\u003enon-linear least squares regression\u003c/strong\u003e\u003cstrong\u003e,\u003c/strong\u003e following the optimization techniques presented by Draper and Smith [14]. The performance of each model was evaluated using the coefficient of determination (R\u0026sup2;), root mean square error (RMSE), and maximum relative error metrics.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eTo ensure immediate usability by practitioners, the empirical formulation was translated into \u003cstrong\u003edesign nomographs and lookup tables\u003c/strong\u003e for direct determination of F(x) and equivalent depth (d). This approach continues the ILRI tradition of visual design aids [3] and the practical framework of Smedema and Rycroft [8], while incorporating modern, data-driven contours based on the hybrid regression model. The integration of theoretical precision with empirical simplicity provides an effective, field-ready design methodology particularly suited for the preservation of heritage sites in dense urban environments. The models were fitted using non-linear least squares regression, and the best model was selected based on the coefficient of determination (R\u0026sup2;), root mean squared error (RMSE), and maximum relative error.\u003c/p\u003e"},{"header":"3. Results","content":"\u003cp\u003eThe results of this two-part investigation are presented sequentially, beginning with the field-calibrated analysis that establishes the significance of the F(x) term, followed by the development and validation of the simplified empirical formulation.\u003c/p\u003e\n\u003cp\u003e3.1. F(x) Convergence and Impact Quantification\u003c/p\u003e\n\u003cp\u003eThe theoretical F(x) term is an infinite series, and the first step was to establish the number of terms required for mathematical convergence. As shown in Figure 1, the value of F(x) converges rapidly, with the value stabilizing after approximately 5-6 odd terms. Using a 10-term expansion ensures a fully converged value for all subsequent calculations. The practical impact of neglecting this converged term is substantial. A direct comparison of discharge calculations for the three mosque sites\u0026mdash;one with the 10-term F(x) and one without\u0026mdash;revealed that neglecting the term leads to an underestimation of the required discharge by over 40%. This finding from the field-calibrated models highlights a critical flaw in common engineering practice and serves as the primary motivation for developing a practical alternative.\u003c/p\u003e\n\u003cp\u003e3.2. Field-Calibrated Sensitivity and Discharge Analysis\u003c/p\u003e\n\u003cp\u003eA sensitivity analysis was conducted to understand the influence of key design parameters on the required discharge, using the field-calibrated data from the three mosque sites. The results, shown in Figure 2, indicate that the discharge is highly sensitive to drain spacing (L) and hydraulic conductivity (K), reinforcing the need for accurate parameter determination.\u003c/p\u003e\n\u003cp\u003eThe analysis shows a strong, non-linear relationship between discharge and key parameters like drain spacing and hydraulic conductivity. Furthermore, the analysis of different permeability scenarios (Figure 3) demonstrates the dramatic impact of the underlying geology. A shift from a conservative (low permeability) to a moderate permeability assumption for the deeper soil layers can amplify the required discharge by a factor of over 100, underscoring the importance of thorough geotechnical investigation.\u003c/p\u003e\n\u003cp\u003eThe required discharge is significantly higher under moderate permeability assumptions for the deeper soil layers.\u003c/p\u003e\n\u003cp\u003e3.3. Uncertainty Quantification and Probabilistic Design Envelopes\u003c/p\u003e\n\u003cp\u003eTo account for the inherent uncertainty in the field-measured parameters, a 10,000-iteration Monte Carlo simulation was performed for each of the three mosque sites. The results (Figure 4) reveal a wide probability distribution for the required discharge, indicating that a deterministic, single-value design is insufficient.\u003c/p\u003e\n\u003cp\u003eThe analysis reveals a broad range of potential discharge values, highlighting the need for a probabilistic design approach. Based on this probabilistic analysis, we developed design discharge envelopes for each site (Figure 5). These envelopes provide P10, P50, and P90 percentile values, allowing for a risk-based design approach. For critical heritage sites, designing for the P90 (90th percentile) discharge provides a 90% confidence level that the drainage system will be adequate, representing a far more robust approach than a single deterministic calculation.\u003c/p\u003e\n\u003cp\u003eThe envelopes provide a range of design values (P10, P50, P90) for a risk-based approach to drainage system design.\u003c/p\u003e\n\u003cp\u003e3.4. Development and Validation of the Empirical F(x) Formulation\u003c/p\u003e\n\u003cp\u003eWhile the field analysis confirmed the importance of F(x), it did not solve the problem of its computational complexity. To address this, we developed an empirical formulation from a 22,500-scenario dataset. Five different regression models were tested, and the results are summarized in Table 2. The Hybrid Formulation, designed to mimic the known asymptotic behavior of the F(x) function, demonstrated vastly superior performance.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eTable 2: Performance Comparison of Five Empirical F(x) Formulations\u0026nbsp;\u003c/p\u003e\n\u003ctable border=\"1\" cellspacing=\"0\" cellpadding=\"0\" width=\"100%\"\u003e\n \u003cthead\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 30px;\"\u003e\n \u003cp\u003eFormulation\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003eR\u0026sup2;\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003eRMSE\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003eMAE\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 33px;\"\u003e\n \u003cp\u003eMax Relative Error (%)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/thead\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 30px;\"\u003e\n \u003cp\u003ePower Law\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.2249\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.01097\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.00566\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 33px;\"\u003e\n \u003cp\u003e1.760\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 30px;\"\u003e\n \u003cp\u003eExponential Decay\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.0640\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.01206\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.00527\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 33px;\"\u003e\n \u003cp\u003e1.956\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 30px;\"\u003e\n \u003cp\u003eRational Function\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.5335\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.00851\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.00514\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 33px;\"\u003e\n \u003cp\u003e1.295\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 30px;\"\u003e\n \u003cp\u003ePolynomial (2nd Order)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.1614\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.01141\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.00549\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 33px;\"\u003e\n \u003cp\u003e1.848\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 30px;\"\u003e\n \u003cp\u003eHybrid (Recommended)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.99998\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.00005\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 11px;\"\u003e\n \u003cp\u003e0.00002\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 33px;\"\u003e\n \u003cp\u003e0.008\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003eFigure 6 visually confirms this result, showing the near-perfect fit of the Hybrid model compared to the other four approaches.\u003c/p\u003e\n\u003cp\u003eThe Hybrid model (bottom center) clearly outperforms all other approaches, with data points aligning almost perfectly along the 1:1 line. The recommended empirical equation for the F(x) correction term is: F(x) \u0026asymp; 5.432 + 2.832 * exp (-6.414 * x) ,where \u003cu\u003ex = (Di - Dd) / L\u003c/u\u003e. This simple, three-parameter equation achieves a coefficient of determination (R\u0026sup2;) greater than 0.99998 and a maximum relative error of less than 0.01%. The detailed performance analysis in Figure 7 confirms the exceptional accuracy of this formulation. The model accurately traces the theoretical curve, and the residual errors are symmetrically distributed around zero and are negligibly small.\u003c/p\u003e\n\u003cp\u003e3.5. Practical Design Tools\u003c/p\u003e\n\u003cp\u003eTo facilitate the immediate use of this empirical formulation in engineering practice, we have developed a set of design tools. Figure 8 presents two design nomographs that allow for the graphical determination of F(x) and the resulting equivalent depth (\u003cu\u003ed\u003c/u\u003e) based on the site-specific drain spacing (L) and aquifer thickness (Di - Dd).\u003c/p\u003e\n\u003cp\u003eTable 3: Lookup Table for F(x) vs. Dimensionless Parameter \u003cstrong\u003e\u003cu\u003ex\u003c/u\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003ctable border=\"1\" cellspacing=\"0\" cellpadding=\"0\" width=\"100%\"\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 22px;\"\u003e\n \u003cp\u003e\u003cstrong\u003ex = (Di-Dd)/L\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e0.5\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e1.0\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e1.5\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e2.0\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e3.0\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e5.0\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e10.0\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" style=\"width: 22px;\"\u003e\n \u003cp\u003e\u003cstrong\u003eF(X) Value\u003c/strong\u003e\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e5.546\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e5.437\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e5.432\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e5.432\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e5.432\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e5.432\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" style=\"width: 11px;\"\u003e\n \u003cp\u003e5.432\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eTable 4: F(x) and Equivalent Depth (d) for Common Heritage Site Configurations\u003c/p\u003e\n\u003cp\u003e(r = 0.1016m)\u003c/p\u003e\n\u003ctable border=\"1\" cellspacing=\"0\" cellpadding=\"0\" width=\"100%\"\u003e\n \u003cthead\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 19px;\"\u003e\n \u003cp\u003eL (m)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 30px;\"\u003e\n \u003cp\u003eDi-Dd (m)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 14px;\"\u003e\n \u003cp\u003ex\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 17px;\"\u003e\n \u003cp\u003eF(X)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 18px;\"\u003e\n \u003cp\u003ed (m)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/thead\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 19px;\"\u003e\n \u003cp\u003e30\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 30px;\"\u003e\n \u003cp\u003e100\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 14px;\"\u003e\n \u003cp\u003e3.33\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 17px;\"\u003e\n \u003cp\u003e5.432\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 18px;\"\u003e\n \u003cp\u003e10.73\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 19px;\"\u003e\n \u003cp\u003e40\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 30px;\"\u003e\n \u003cp\u003e150\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 14px;\"\u003e\n \u003cp\u003e3.75\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 17px;\"\u003e\n \u003cp\u003e5.432\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 18px;\"\u003e\n \u003cp\u003e14.69\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 19px;\"\u003e\n \u003cp\u003e50\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 30px;\"\u003e\n \u003cp\u003e200\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 14px;\"\u003e\n \u003cp\u003e4.00\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 17px;\"\u003e\n \u003cp\u003e5.432\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 18px;\"\u003e\n \u003cp\u003e18.78\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 19px;\"\u003e\n \u003cp\u003e60\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 30px;\"\u003e\n \u003cp\u003e150\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 14px;\"\u003e\n \u003cp\u003e2.50\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 17px;\"\u003e\n \u003cp\u003e5.432\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 18px;\"\u003e\n \u003cp\u003e22.98\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd style=\"width: 19px;\"\u003e\n \u003cp\u003e80\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 30px;\"\u003e\n \u003cp\u003e100\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 14px;\"\u003e\n \u003cp\u003e1.25\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 17px;\"\u003e\n \u003cp\u003e5.433\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd style=\"width: 18px;\"\u003e\n \u003cp\u003e31.58\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003eThese tools collectively eliminate the need for any complex calculations, allowing practitioners to select appropriate design parameters with confidence in the accuracy of the underlying hydraulic model.\u003c/p\u003e"},{"header":"4. Discussion","content":"\u003cp\u003eThis study successfully bridges the gap between complex drainage theory and practical engineering design for heritage conservation. The integrated findings, from field-calibrated error quantification to the development of a simple and robust empirical solution, provide a comprehensive framework for improving the standard of practice.\u003c/p\u003e\u003cdiv id=\"Sec18\" class=\"Section2\"\u003e\u003ch2\u003e4.1. The Significance of F(x) Precision in Heritage Contexts\u003c/h2\u003e\u003cp\u003eThe field-calibrated analysis of the three Cairo mosques provides unequivocal evidence that the F(x) correction term is not a minor theoretical detail, but a critical component of accurate drainage design. The finding that neglecting F(x) can lead to discharge calculation errors exceeding 40% has profound implications for heritage conservation. An under-designed drainage system, based on an erroneously low discharge estimate, would fail to adequately lower the water table, leaving the historic foundations vulnerable to continued moisture-related decay. The cost of such a design error is not merely financial; it is the potential for irreversible loss of cultural heritage. This finding alone should compel a re-evaluation of current engineering practices that treat the F(x) term as optional or that rely on overly simplified approximations.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec19\" class=\"Section2\"\u003e\u003ch2\u003e4.2. The Empirical Solution: A Practical Path to Accuracy\u003c/h2\u003e\u003cp\u003eHaving established the problem, the second part of this study provides the solution. The development of the empirical equation, \u003cb\u003eF(x)\u0026thinsp;\u0026asymp;\u0026thinsp;5.432\u0026thinsp;+\u0026thinsp;2.832 * exp (-6.414 * x)\u003c/b\u003e, effectively eliminates the primary barrier to accurate F(x) computation: its complexity. With an R\u0026sup2; value exceeding 0.9999 and a maximum relative error of less than 0.01%, this simple three-parameter equation provides the accuracy of the full theoretical series with the simplicity of a basic calculator function. This allows for the rapid iteration and optimization of drainage designs, a crucial capability in the complex, multi-stakeholder environment of heritage projects.The accompanying design nomographs and lookup tables further enhance the practical utility of this work, providing multiple entry points for practitioners of varying technical backgrounds to access this enhanced accuracy.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec20\" class=\"Section2\"\u003e\u003ch2\u003e4.3. Comparison with and Contribution to Existing Approaches\u003c/h2\u003e\u003cp\u003eThis work builds upon a long history of drainage research. The foundational nomographs published by the International Institute for Land Reclamation and Improvement (ILRI) [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e] have been invaluable for decades, but our empirical approach offers a more flexible, precise, and digitally-native alternative. Our work is also distinct from, but complementary to, the research by Shokri \u0026amp; Bardsley (2015) [\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e], who focused on corrections related to soil-water retention properties. Our study provides the essential geometric correction, and a powerful future approach would be to integrate both corrections into a single, comprehensive design tool. The classical work of Wesseling (1964) [\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e] and more recent evaluations by Moriasi et al. (2013) [\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e] have compared various drainage spacing formulas, but none have provided the combination of field validation and simplified empirical formulation presented here. The primary contribution of our work is the creation of a complete narrative: we not only identify and quantify a critical error in current practice using real-world field data, but we also provide the direct, validated solution to that problem.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec21\" class=\"Section2\"\u003e\u003ch2\u003e4.4. Application to Cairo Heritage Sites and Beyond\u003c/h2\u003e\u003cp\u003eThe probabilistic design envelopes developed from the Monte Carlo analysis (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e) offer a new standard for drainage design at critical heritage sites. Moving away from a single deterministic value to a risk-based P90 design provides a quantifiable level of confidence that the system will perform as intended, even with the inherent uncertainties in the input parameters. This approach is directly applicable to the ongoing conservation efforts at the three mosques studied and can serve as a model for other heritage sites in Cairo.While the empirical formulation was specifically calibrated for Cairo-type hydrogeological conditions, the methodology itself is broadly applicable. The process of generating a large simulation dataset and fitting a robust empirical model can be replicated for other regions with different geological characteristics, creating a new generation of site-specific, highly accurate design tools.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec22\" class=\"Section2\"\u003e\u003ch2\u003e4.5. Limitations\u003c/h2\u003e\u003cp\u003eIt is important to acknowledge the limitations of this study. The empirical formulation is validated for the parameter space defined in Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e and should be used with caution outside of these ranges. Furthermore, this study, like the Hooghoudt equation itself, is based on the assumption of steady-state flow and homogeneous, isotropic soil. In highly complex, stratified geological settings, or where transient effects are significant, numerical modeling (e.g., with MODFLOW or HYDRUS) remains the most appropriate tool for final design verification.\u003c/p\u003e\u003c/div\u003e"},{"header":"5. Conclusion","content":"\u003cp\u003eThis study has successfully addressed a critical gap between subsurface drainage theory and practical engineering for heritage conservation. By integrating a field-calibrated analysis with the development of a simplified empirical model, we have created a comprehensive framework that both quantifies a significant problem and provides its direct solution. The key findings are twofold. First, the field-calibrated analysis of three historic mosques in Cairo demonstrates that neglecting the Hooghoudt F(x) correction term can lead to an underestimation of required drainage discharge by over 40%, posing a significant risk to the long-term preservation of these invaluable structures. Second, we have shown that the computational complexity of the F(x) term can be eliminated without sacrificing accuracy. The developed empirical equation, \u003cb\u003eF(x)\u0026thinsp;\u0026asymp;\u0026thinsp;5.432\u0026thinsp;+\u0026thinsp;2.832 * exp (-6.414 * x)\u003c/b\u003e, replicates the theoretical series with near-perfect accuracy (R\u0026sup2; \u0026gt;0.9999) and is presented with a suite of practical design tools, including nomographs and lookup tables.This work provides a clear path for engineers and conservationists to move beyond overly simplified or inaccurate methods and adopt a design approach that is both theoretically sound and practically achievable. The ability to rapidly and reliably calculate the required drainage capacity is a critical component in the ongoing effort to protect our shared cultural heritage from the persistent threat of rising groundwater. The methodology and tools presented here offer a significant step forward in that endeavor.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003eFunding\u003c/p\u003e\n\u003cp\u003eThis research received no external funding\u003c/p\u003e\n\u003cp\u003eData availability\u003c/p\u003e\n\u003cp\u003eThe article contains all of the data generated or analyzed throughout the research.\u003c/p\u003e\n\u003cp\u003eDeclarations Ethics approval and consent to participate\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;Ethical approval was not sought for the present study because there are no human/animal subjects in this article.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eConsent to publication\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eAuthor approved the manuscript for publication.\u003c/p\u003e\n\u003cp\u003eCompeting interests\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eThe author declares no competing interests\u003c/p\u003e\n\u003cp\u003eConsent to Participate declaration\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;Not applicable\u003c/p\u003e\n\u003cp\u003eAuthor Contribution\u003c/p\u003e\n\u003cp\u003eWriting original draft, Methodology, Formal analysis, Data curation, Conceptualization.\u003c/p\u003e\n\u003cp\u003eAcknowledgement\u003c/p\u003e\n\u003cp\u003eI would like to thank my colleagues at the Geotechnical Engineering Research Institute and the Greater Cairo Wastewater Company for their valuable support and collaboration.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eAbdel-Shafy, H. I., \u0026amp; Aly, R. O. (2012). \u003cem\u003eWater issues in Egypt: Resources, pollution, and protection endeavors.\u003c/em\u003e \u003cstrong\u003eEnvironmental Science and Pollution Research\u003c/strong\u003e, 19(6), 2163\u0026ndash;2176. https://doi.org/10.1007/s11356-012-0877-8\u003c/li\u003e\n\u003cli\u003eUNESCO World Heritage Centre. (2021). \u003cem\u003eHistoric Cairo.\u003c/em\u003e Retrieved from https://whc.unesco.org/en/list/89/\u003c/li\u003e\n\u003cli\u003eRitzema, H. P. (Ed.). (1994). \u003cem\u003eDrainage Principles and Applications\u003c/em\u003e (Vol. 16). \u003cstrong\u003eInternational Institute for Land Reclamation and Improvement (ILRI)\u003c/strong\u003e, Wageningen, The Netherlands.\u003c/li\u003e\n\u003cli\u003eFAO. (2002). \u003cem\u003eDrainage Design Manual.\u003c/em\u003e \u003cstrong\u003eFAO Irrigation and Drainage Paper 61.\u003c/strong\u003e Rome: Food and Agriculture Organization of the United Nations.\u003c/li\u003e\n\u003cli\u003eHooghoudt, S. B. (1940). \u003cem\u003eGeneral consideration of the problem of field drainage by parallel drains, ditches, watercourses, and channels.\u003c/em\u003e Institute for Land and Water Management Research, Wageningen, The Netherlands.\u003c/li\u003e\n\u003cli\u003eVan der Molen, W. H., \u0026amp; Wesseling, J. (1991). \u003cem\u003eA solution in closed form and a series solution to replace the tables for the thickness of the equivalent layer in Hooghoudt\u0026rsquo;s drain spacing formula.\u003c/em\u003e \u003cstrong\u003eAgricultural Water Management\u003c/strong\u003e, 19(1), 1\u0026ndash;16. https://doi.org/10.1016/0378-3774(91)90061-T\u003c/li\u003e\n\u003cli\u003eWesseling, J. (1964). \u003cem\u003eA comparison of the steady state drain spacing formulas of Hooghoudt and Kirkham in connection with design practice.\u003c/em\u003e \u003cstrong\u003eJournal of Hydrology\u003c/strong\u003e, 2(1), 25\u0026ndash;32. https://doi.org/10.1016/0022-1694(64)90004-6\u003c/li\u003e\n\u003cli\u003eSmedema, L. K., \u0026amp; Rycroft, D. W. (1983). \u003cem\u003eLand Drainage: Planning and Design of Agricultural Drainage Systems.\u003c/em\u003e Ithaca, NY: Cornell University Press.\u003c/li\u003e\n\u003cli\u003eShokri, A., \u0026amp; Bardsley, W. E. (2015). \u003cem\u003eEnhancement of the Hooghoudt drain-spacing equation.\u003c/em\u003e \u003cstrong\u003eJournal of Irrigation and Drainage Engineering\u003c/strong\u003e, 141(6), 04014082. https://doi.org/10.1061/(ASCE)IR.1943-4774.0000835\u003c/li\u003e\n\u003cli\u003eFAO \u0026amp; ILRI. (2002\u0026ndash;1994). \u003cem\u003eStandardized Field Methods for Subsurface Drainage Design.\u003c/em\u003e Joint Technical Framework, Rome\u0026ndash;Wageningen.\u003c/li\u003e\n\u003cli\u003eMoriasi, D. N., et al. (2013). \u003cem\u003eEvaluation of the Hooghoudt and Kirkham tile drain equations in the Soil and Water Assessment Tool to simulate tile flow and nitrate-nitrogen.\u003c/em\u003e \u003cstrong\u003eJournal of Environmental Quality\u003c/strong\u003e, 42(6), 1699\u0026ndash;1710. https://doi.org/10.2134/jeq2013.01.0018\u003c/li\u003e\n\u003cli\u003eMetropolis, N., \u0026amp; Ulam, S. (1949). \u003cem\u003eThe Monte Carlo Method.\u003c/em\u003e \u003cstrong\u003eJournal of the American Statistical Association\u003c/strong\u003e, 44(247), 335\u0026ndash;341. https://doi.org/10.1080/01621459.1949.10483310\u003c/li\u003e\n\u003cli\u003eSaltelli, A., Chan, K., \u0026amp; Scott, E. M. (2010). \u003cem\u003eSensitivity Analysis in Practice: A Guide to Assessing Scientific Models.\u003c/em\u003e Chichester: John Wiley \u0026amp; Sons.\u003c/li\u003e\n\u003cli\u003eDraper, N. R., \u0026amp; Smith, H. (1998). \u003cem\u003eApplied Regression Analysis\u003c/em\u003e (3rd ed.). New York: Wiley-Interscience.\u003c/li\u003e\n\u003cli\u003eU.S. Army Corps of Engineers (USACE). (2019). \u003cem\u003eEngineering and Design: Risk-Informed Decision Making for Hydrologic Engineering\u003c/em\u003e (EM 1110-2-1619). Washington, D.C.\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"
[email protected]","identity":"discover-water","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"diwa","sideBox":"Learn more about [Discover Water](https://www.springer.com/43832)","snPcode":"","submissionUrl":"","title":"Discover Water","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Discover Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Subsurface Drainage, Hooghoudt Equation, F(x) Correction, Empirical Modeling, Heritage Site Conservation, Field Calibration, Monte Carlo Simulations, and Design Nomograph","lastPublishedDoi":"10.21203/rs.3.rs-7906837/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-7906837/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eRising groundwater levels pose a significant threat to urban heritage sites, highlighting the need for efficient subsurface drainage systems. The Hooghoudt equation is a foundational tool for such designs, yet its practical use is hindered by the complexity of computing the correction term F(x), which is critical in urban environments with irregular geometries. This study presents a two-phase investigation. The first phase involves a field-calibrated analysis of three historic mosques in Cairo, applying a refined 10-term series for F(x). Findings reveal that omitting this term can lead to discharge errors exceeding 40%. In the second phase, a simplified empirical formula for F(x) is developed using data from over 22,000 simulation scenarios. A hybrid exponential decay model delivers outstanding accuracy (R\u0026sup2; \u0026gt;0.9999; maximum error\u0026thinsp;\u0026lt;\u0026thinsp;0.01%). The resulting three-parameter equation, F(x)\u0026thinsp;\u0026asymp;\u0026thinsp;5.432\u0026thinsp;+\u0026thinsp;2.832\u0026sdot;exp(\u0026minus;\u0026thinsp;6.414x), offers a practical alternative to complex series computations. Accompanying tools, such as nomographs and lookup tables enable rapid and reliable drainage planning to support the preservation of heritage sites vulnerable to rising groundwater.\u003c/p\u003e","manuscriptTitle":"Field-Calibrated Simplification of the F(x) Correction in Subsurface Drainage Design for Urban Heritage Sites","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-11-17 07:59:06","doi":"10.21203/rs.3.rs-7906837/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2025-12-14T12:48:26+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-12-13T19:37:13+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"185540493854113938040442226421377825345","date":"2025-11-26T17:05:05+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"96446315880925753701701935031382635587","date":"2025-11-25T12:33:21+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-11-23T04:35:40+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"184992062315575005273772745562582946010","date":"2025-11-11T13:07:53+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2025-11-05T08:31:50+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2025-11-04T15:29:33+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2025-10-27T07:18:04+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2025-10-27T07:17:54+00:00","index":"","fulltext":""},{"type":"submitted","content":"Discover Water","date":"2025-10-20T14:37:10+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"
[email protected]","identity":"discover-water","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"diwa","sideBox":"Learn more about [Discover Water](https://www.springer.com/43832)","snPcode":"","submissionUrl":"","title":"Discover Water","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Discover Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"1d0d03f6-71a3-401c-bc16-b678e3cec804","owner":[],"postedDate":"November 17th, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"under-review","subjectAreas":[],"tags":[],"updatedAt":"2026-02-11T04:23:48+00:00","versionOfRecord":[],"versionCreatedAt":"2025-11-17 07:59:06","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-7906837","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-7906837","identity":"rs-7906837","version":["v1"]},"buildId":"XKTyCvWXoU3ODBz1xrDgd","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.