Calculation Of Rate Research Articles

A high-resolution map of diffuse groundwater recharge rates for Australia

Abstract. Estimating groundwater recharge rates is important to understand and manage groundwater. Numerous studies have used collated recharge datasets to understand and project regional- or global-scale groundwater recharge rates. However, recharge estimation methods all have distinct assumptions, quantify different recharge components and operate over different temporal scales. We use over 200 000 groundwater chloride measurements to estimate groundwater recharge rates using an improved chloride mass balance (CMB) method across Australia. Groundwater recharge rates were produced stochastically using gridded chloride deposition, runoff and precipitation datasets. After filtering out groundwater recharge rates where the assumptions of the method may have been compromised, 98 568 estimates of recharge were produced. The resulting groundwater recharge rates and 17 spatial datasets were integrated into a random forest regression algorithm, generating a high-resolution (0.05°) model of groundwater recharge rates across Australia. The regression reveals that climate-related variables, including precipitation, rainfall seasonality and potential evapotranspiration, exert the most significant influence on groundwater recharge rates, with vegetation (the normalised difference vegetation index or NDVI) also contributing significantly. Importantly, the mean values of both the recharge point dataset (43.5 mm yr−1) and the spatial recharge model (22.7 mm yr−1) are notably lower than those reported in previous studies, underscoring the prolonged timescale of the CMB method, the potential disparities arising from distinct recharge estimation methodologies and limited averaging across climate zones. This study presents a robust and automated approach to estimate recharge using the CMB method, offering a unified model based on a single estimation method. The resulting datasets, the Python script for recharge rate calculation and the spatial recharge models collectively provide valuable insights for water resource management across the Australian continent, and similar approaches can be applied globally.

Innovative recovery of matrix layered double hydroxide from simulated acid mine wastewater for the removal of copper and cadmium from wastewater.

This study innovatively added biochar to optimize regulation in the neutralization process of simulated acid mine drainage (AMD) and recovered a new type of matrix layered double hydroxides (MLDH), which can be used to remove copper (Cu(II)) and cadmium (Cd(II)) from wastewater. A series of batch experiments show that MLDH with strong selective removal ability of Cu(II) and Cd(II) can be successfully obtained by adding biochar (BC) at pH = 5 end in the neutralization process. Kinetic and isotherm modeling studies indicated that the removal of Cu(II) and Cd(II) by the MLDH was a chemical multilayer adsorption process. The removal mechanism of Cu(II) and Cd(II) was further analyzed through related characterization analysis with contribution rate calculation: the removal rates of Cu(II) and Cd(II) by ion exchange were 42.7% and 26%, while that by precipitation were 34.5% and 49.9%, respectively. This study can provide a theoretical reference and experimental basis for the recovery and utilization of valuable by-products in AMD and the treatment of heavy metal wastewater.

Occupational prevalence rates in metallurgists in the Sverdlovsk Region

Introduction. The previous analysis conducted for 2003–2022 showed ferrous industry employing over twenty five thousand people to make the most contribution to occupational disease incidence rates in the Sverdlovsk Region since 68.6 % of its workers are affected by industrial exposures above occupational safety standards. The purpose of the study was the analysis of occupational disease incidence in workers of a ferrous metal processing plant of the Sverdlovsk Region and assessment of the impact of adverse factors of the work environment and industrial processes on its rates.
 Materials and methods. Long-term occupational disease incidence in workers of the ferrous metal processing plant has been analyzed by the causative agent, nosology, age, gender, and job. Data of the reporting form No. 1-T (working conditions) for the study period have been scrutinized; comparative data evaluation by working conditions and occupational disease incidence rates has been made.
 Results. High incidence rates of work-related diseases in the metallurgists indicate to a significant proportion (70 %) of them to be exposed to occupational hazards. The established 20-year rate at the enterprise was 5.6 and 1.6 times higher than the regional and Russian averages in this industry, respectively. Increased gas and dust levels in the workplace air remain the leading health risk factors in ferrous metallurgy.
 Limitations. The limitation of this study was the lack of annual data on workers by job preventing calculation of job-specific disease incidence rates per 10 thousand workers.
 Conclusion. Despite the substantial modernization of the metallurgical complex, metal processing technology still does not comply with safety standards, thus requiring the development and implementation of science-based preventive measures aimed at maintaining workers’ health.

The role of the Holocene transgression in the environmental changes of lagoons and marshes of the Mediterranean coast

The aim of this paper is to evaluate the position of the Holocene transgression on the coast of the Valencia Gulf. To achieve this goal, a sedimentological and micropaleontological analysis of samples of fossil benthic foraminifera in cores from six boreholes drilled in three different wetlands (Almenara marsh, Moro marsh and Valencia lagoon) was carried out. In order to assess the extent of sea level change, at least three factors must be taken into account: 1) eustatic rise (global scale) 2) neotectonics (regional scale) and 3) the sedimentary factor (local and regional scale). Regarding the first factor, some authors have noted that the Holocene transgression reached +1 m above the current sea level in this area. As for the second factor, while the Valencia lagoon is located on a subsiding coast the northernmost marshes are located on a tectonically stable coast. Finally, sedimentation rates are highly variable depending on the area; whereas in the Valencia lagoon the shoreline generally shifts seawards (progradation) in the northernmost areas the coast is currently undergoing retrogradation, although the anthropic factor clearly interferes with this pattern. One of the main conclusions is that in the Valencian coast the Holocene maximum relative sea level did not exceed the present-day one, except in the case of the subsiding area of the Valencia lagoon. The calculation of sedimentation rates, the relative sea-level curve and the evolution of the water column in the different wetlands seem to corroborate it. Radiocarbon and AAR analyses allowed us to date the maximum extent of the Holocene transgression on the Valencian coast at around 5500 cal yr BP in all cases. Despite the different subsidence the three studied wetlands reflected the predominant marine influence until 5500–5200 BP, which was later followed by a growing fluvial or continental incidence, when glacioeustatic adjustments defined the patterns of coastal progradation.

Comprehensive evaluation of natural ventilation potential of buildings in urban areas under the influence of multiple environment-related factors

Realistic natural ventilation potential estimation is critical for energy-efficient building design. For buildings in urban areas, evaluations using design expectations tend to overestimate potential. This study presents a comprehensive framework and corresponding method to predict the natural ventilation potential of buildings, integrating multiple key environmental factors of climate features, local wind status, and air pollution level. It addresses their complex interplay. Climate impacts are evaluated through simulating using a pre-set natural ventilation rate. Computational fluid dynamics models typical neighbourhood wind fields, giving empirical wind pressure coefficient distributions by layout and wind direction to enable real-time ventilation rate calculation. Substituting real-time for pre-set rates, the natural ventilation potential considering the community layout can be obtained. Based on the data-driven rapid pollutant prediction model established in the previous study, the predictions of the temporal and spatial distributions of particulate matter can be obtained. Filtering excessive pollution periods yields air quality-integrated potentials. Natural ventilation potentials are compared under different analysis scenarios for a case-study building. Hotter climates and more frequent heatwaves lower the potential. Morphology-based potentials are significantly lower than design expectations using constant ventilation rates. Deteriorating urban air quality also directly limits ventilation use, substantially decreasing potentials when filtering for poor conditions. By capturing complex environmental interactions, this research provides a robust platform to reveal localization potentials and limitations.

Calculation of Reaction Rate and Uncertainty for d(p,γ)3He and 3He(d,p)4He Using Bayesian Statistical Approach in Astrophysics

Calculation of Reaction Rate and Uncertainty for d(p,γ)3He and 3He(d,p)4He Using Bayesian Statistical Approach in Astrophysics

An improved analytical solution to outcrossing rate for scalar nonstationary and non-gaussian processes

For time-dependent reliability methods, there is a lack of analytical solutions to outcrossing rates of nonstationary and non-Gaussian processes. This paper aims to propose an analytical method to determine this type of outcrossing rate with improved accuracy and efficiency. The novelty of this proposed method is that it can consider the probabilistic properties of non-Gaussian processes by combining a two-component parallel system model with higher-order moments-based reliability indexes. The main contributions of this method are: (1) compared with PHI2+ method, it is more accurate for the calculation of the outcrossing rate of nonstationary and non-Gaussian processes with nonlinear limit state functions; (2) compared with MPHI2 method, it is analytical and not only applicable for the nonstationary stochastic processes but also insensitive to time increments; and (3) it is not only more accurate than the existing methods for nonstationary and non-Gaussian processes but also more computation efficient than the Monte Carlo simulation method. The proposed method has shown its advantages for practical structures with neither Gaussian processes nor linear limit state functions, which are beneficial for both researchers and asset managers to evaluate the time-dependent reliability of structures accurately and to develop risk-informed maintenance schemes with a view to prolonging their service life.

Trend of using cementless total knee arthroplasty: a nationwide analysis from 2015 to 2021

BackgroundModern cementless total knee arthroplasty (TKA) fixation has shown comparable long-term outcomes to cemented TKA, but the trend of using cementless TKA remains unclear. This study aimed to investigate the trend of using cementless TKA based on a national database.MethodsThe patients undergoing cementless TKA between 2015 and 2021 were retrospectively extracted from the PearlDiver (Mariner dataset) Database. The annual percentage of cementless TKA was calculated using the following formula: annual number of cementless TKA/annual number of TKA. The trend of the number of patients undergoing cementless TKA was created according to a compounded annual growth rate (CAGR) calculation of annual percentages. Patient age, comorbidity, region, insurance type, etc., were also investigated. Differences were considered statistically significant at P < 0.05.ResultsOf the 574,848 patients who received TKA, 546,731 (95%) underwent cemented fixation and 28,117 (5%) underwent cementless fixation. From 2015 to 2021, the use of cementless TKA significantly increased by 242% from 3 to 9% (compounded annual growth rate (CAGR): + 20%; P < 0.05). From 2015 to 2021, we observed a CAGR greater than 15% for all age groups (< 50, 50–59, 60–69, 70–74, 75 +), insurance types (cash, commercial, government, Medicare, Medicaid), regions (Midwest, Northeast, South, West), sex (male and female), and certain comorbidities (osteoporosis, diabetes mellitus, tobacco use, underweight (BMI < 18.5), rheumatoid arthritis) (P < 0.05 for all). Patients undergoing TKA with chronic kidney disease, prior fragility fractures, and dementia demonstrated a CAGR of + 9%–13% from 2015 to 2021 (P < 0.05).ConclusionFrom 2015 to 2021, the use of cementless TKA saw a dramatic increase in all patient populations. However, there is still no consensus on when to cement and in whom. Clinical practice guidelines are needed to ensure safe and effective use of cementless fixation.

Prioritization of control measures in leakage scenario using Hendershot theory and FBWM-TOPSIS.

Currently, there is increasing concern about the safety and leakage of process industries. Therefore, the present study aims to prioritize control measures before and after the leakage scenario by using the Hendershot theory and MCDM techniques. In this study, two proactive and reactive layers were selected before and after leakage of tanks, respectively. Then, criteria and alternatives were selected to perform fuzzy TOPSIS (FTOPSIS) and find the best alternative based on the literature review and Hendershot approach. The linear model of the fuzzy Best-Worst method (FBWM) was constructed and resolved using Lingo 17 software. Subsequently, criteria were assigned weights based on thorough calculations of the inconsistency rate. The weight of study experts was equal to 0.25. The results of FBWM showed that the reliability index with a weight of 0.3727 was ranked first and the inconsistency rate ([Formula: see text]) was calculated to be equal to 0.040. Inherent Safety Design (ISD) (0.899) and passive safety (0.767) also ranked first before and after tank leaks, respectively. Using the FBWM method leads to fewer pairwise comparisons and at the same time more stability. Although ISD and passive strategies are more valid and strict, elements of all strategies are necessary for a comprehensive process safety management program.

NLTE modelling of water-rich exoplanet atmospheres. Cooling and heating rates

The hydrogen and water molecules respond very differently to the collisional-radiative processes taking place in planetary atmospheres. Naturally, the question arises whether H2O-rich atmospheres are more (or less) resilient to long-term mass loss than H2-dominated ones if they radiate away the incident stellar energy more (or less) efficiently. If confirmed, the finding would have implications on our understanding of the evolution of sub-Neptune exoplanets. As a key step towards answering this question, we present a non-local thermodynamic equilibrium (NLTE) model of H2O for the atmospheric region where the gas accelerates to escape the planet and conditions relevant to close-in sub-Neptunes. Our exploratory calculations for isothermal gas composed of H2, H2O and e− reveal that: (1) In the pressure region ∼10−2–10−4 dyn cm−2 where the stellar extreme-ultraviolet (XUV) photons are typically deposited in the atmosphere, H2O is in rotational LTE but vibrational NLTE. Vibrational LTE is facilitated by high H2O abundances and fractional ionizations, and we report critical densities for the LTE-NLTE transition; (2) Vibrational cooling may locally dominate over rotational cooling, partly because of the comparatively small opacities of ro-vibrational lines; (3) Even low H2O abundances notably enhance the cooling, foreseeably offsetting some of the stellar heating; (4) Heating due to the deposition of stellar infrared (IR) photons is significant at pressures ≳0.1 dyn cm−2. We estimate the contribution of H2O excitation to the internal energy of the gas and speculate on the photodissociation from the excited vibrational states. Ultimately, our findings motivate the consideration of NLTE in the mass loss rate calculations of H2O-rich atmospheres.

From video to vital signs: using personal device cameras to measure pulse rate and predict blood pressure using explainable AI

In this study, we address the need for low-cost, rapid, and reliable methods for blood pressure measurement to facilitate timely diagnosis and treatment of hypertension. We aim to enhance understanding and trust in the novel remote photoplethysmography technology, specifically as deployed by Lifelight®\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^\\circledR$$\\end{document}, for vital sign monitoring. Lifelight®\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^\\circledR$$\\end{document} utilizes remote photoplethysmography to extract physiologically meaningful features from video recordings on personal devices. This study describes the accurate calculation of pulse rates and the extraction of morphological and time series features from the remote photoplethysmography signal for blood pressure prediction. Unlike common approaches, distinct model types are employed for predicting systolic and diastolic blood pressure. Evaluation on an independent clinical laboratory’s dataset demonstrates the efficacy of the proposed method. The pulse rate accuracy, measured as root mean square error, is 1.9 beats per minute. The model for systolic blood pressure prediction shows a mean error of 1.2 mmHg and a standard deviation of 15.9 mmHg. The model for diastolic prediction shows a mean error of -\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$-$$\\end{document}3.2 mmHg and a standard deviation of 13.1 mmHg. This research establishes the validity of remote photoplethysmography technology by showcasing accurate pulse rate calculations and effective feature extraction for blood pressure prediction. The results indicate the potential for widespread adoption of remote photoplethysmography in measuring blood pressure and other vital signs.

High Accuracy Photodissociation Rates in VULCAN with MCRT

Photochemical processing is a key influence on the chemical compositions of exoplanet atmospheres. We develop a 1D plane-parallel Monte Carlo radiative-transfer (MCRT) path length method to calculate actinic fluxes for photochemical modeling, which we couple to the VULCAN photochemical model. Due to the flexible nature of MCRT methods, this enables accurate calculation of photodissociation rates for complex multiple scattering scenarios encountered in exoplanet atmospheres such as high altitude haze layers.

Excitation and recombination studies with silicon and sulphur ions at an EBIT

Measurements of electron-impact excitation and recombination rate coefficients of highly charged Si and S ions at the Stockholm electron beam ion trap are reported. The experimental method was a combination of photon detection from the trapped ions during probing and subsequently extraction and time-of-flight (TOF) charge analysis of these ions. The TOF technique allows to measure recombination rate coefficients separately for every charge state, and together with the photon spectra of these ions also the excitation rate coefficients. In this paper, we present more details of the experimental procedure and summarize the experimental results in comparison with two different state-of-the-art calculations of recombination and excitation rates for Si10+–Si13+ and S12+–S15+ ions. One of these uses a relativistic configuration interaction approach (flexible atomic code) and the other is a relativistic many-body perturbation theory. A good to excellent agreement with both of them is found in energy and resonance strength for the investigated ions.

Towards the computation of inclusive decay rates using lattice QCD

We present a non-perturbative computation of inclusive rates of semileptonic decays of heavy mesons from lattice QCD simulations. The calculation is based on the extraction of smeared spectral functions obtained from four-point Euclidean correlation functions computed on configuration ensembles of the JLQCD and ETM collaborations. We compare our results for the inclusive decay rates with analytical predictions from the operator-product expansion, finding a good agreement for the calculation of the inclusive decay rate. This opens the path to the theoretical determination of the magnitude of the Cabibbo-Kobayashi-Maskawa matrix element V_{cb}Vcb to a level of precision competitive with the present experimental uncertainty.

Study on backfire characteristics of port fuel injection single-cylinder hydrogen internal combustion engine

Hydrogen energy is one of the most promising sources of carbon neutrality from well to wheel. Among various applications of hydrogen energy, hydrogen internal combustion engines serve as the critical technical pathway. Nevertheless, several challenges persist, primarily related to abnormal combustion resulting from excessive flame propagation velocity. Of these challenges, backfire emerges as the most recurrent and pressing issue to be resolved in the port fuel-injected hydrogen internal combustion engine (PFI-H2ICE). In order to delve into the characteristics and mechanisms of backfire, this study concentrates on a critical internal parameter, namely the start of injection (SOI), and comprehensively assesses the backfire occurrences under varying lambda and engine speed conditions. Experimental findings reveal that SOI encompasses one of the pivotal factors influencing backfire. Backfire predominantly arises within two ranges of SOI conditions: the region before −430 degCA ATDC (Zone I) and the stretch between −360 degCA ATDC and − 280 degCA ATDC (Zone II). In the process of employing EGR to depress the effects of backfire on subsequent tests, it was incidentally found that the Universal Exhaust Gas Oxygen Sensor (UEGO) installed on the intake manifold for executing EGR rate calculations acted as a hotspot (with temperatures exceeding 100 °C). This not only intensified the backfire strength and the scope in Zone I but also caused the most severe backfire events to occur after −270 degCA ATDC (in Zone III). Despite this, the EGR strategy still demonstrated a significant backfire suppression effect, only exhibiting almost no influence on backfire caused by hotspots within the intake manifold. With an escalation in engine speed, the intensity of backfire diminishes, and the range of backfire occurrence decreases accordingly. By exercising control over SOI within the operating range from -430degCA to -360degCA ATDC (Zone A), the risk of backfire can be dramatically reduced, thereby achieving nearly backfire-free operation. Detailed analysis of cylinder pressure indicates that backfire may also give rise to other abnormal combustions, such as misfire, pre-ignition, and knock, in subsequent cycles. Moreover, this effect amplifies in tandem with an increase in mixture concentration. A certain degree of mutual promotion among diverse abnormal combustions is also observed. This study furnishes the most direct and practical foundation for internal and external strategies in the realm of backfire control.

Optical penetration models for practical prediction of femtosecond laser ablation of dental hard tissue.

To develop and practically test high-precision femtosecond laser ablation models for dental hard tissue that are useful for detailed planning of automated laser dental restorative treatment. Analytical models are proposed, derived, and demonstrated for practical calculation of ablation rates, ablation efficiency and ablated morphology of human dental enamel and dentin using femtosecond lasers. The models assume an effective optical attenuation coefficient for the irradiated material. To achieve ablation, it is necessary for the local energy density of the attenuated pulse in the hard tissue to surpass a predefined threshold that signifies the minimum energy density required for material ionization. A 1029 nm, 40 W carbide 275 fs laser was used to ablate sliced adult human teeth and generate the data necessary for testing the models. The volume of material removed, and the shape of the ablated channel were measured using optical profilometry. The models fit with the measured ablation efficiency curve against laser fluence for both enamel and dentin, correctly capturing the fluence for optimum ablation and the volume of ablated material per pulse. The detailed shapes of a 400-micrometer wide channel and a single-pulse width channel are accurately predicted using the superposition of the analytical result for a single pulse. The findings have value for planning automated dental restorative treatment using femtosecond lasers. The measurements and analysis give estimates of the optical properties of enamel and dentin irradiated with an infrared femtosecond laser at above-threshold fluence and the proposed models give insight into the physics of femtosecond laser processing of dental hard tissue.

How many motor units is enough? An assessment of the influence of the number of motor units on firing rate calculations

The number of motor units included in calculations of mean firing rates varies widely in the literature. It is unknown how the number of decomposed motor units included in the calculation of firing rate per participant compares to the total number of active motor units in the muscle, and if this is different for males and females. Bootstrapped distributions and confidence intervals (CI) of mean motor unit firing rates decomposed from the tibialis anterior were used to represent the total number of active motor units for individual participants in trials from 20 to 100 % of maximal voluntary contraction. Bootstrapped distributions of mean firing rates were constructed using different numbers of motor units, from one to the maximum number for each participant, and compared to the CIs. A probability measure for each number of motor units involved in firing rate was calculated and then averaged across all individuals. Motor unit numbers required for similar levels of probability increased as contraction intensity increased (p < 0.001). Increased levels of probability also required higher numbers of motor units (p < 0.001). There was no effect of sex (p ≥ 0.97) for any comparison. This methodology should be repeated in other muscles, and aged populations.

Uplifted Pleistocene Marine Terraces at Active Margins: Modeling Reveals the Effects of Sea Reoccupation and Coseismic Uplift on Uplift Rate Calculation

AbstractUplifted Pleistocene marine terrace sequences are used to quantify uplift rates along active margins by knowing terrace age and elevation, and sea level (SL) position at the time of terrace formation. When terraces are undated, ages are assigned by correlating terraces at progressively higher elevations with progressively older highstands. Uplift at convergent margins can be constant over time or occur coseismically during upper plate earthquakes. We explore the formation of terrace sequences under conditions of constant and earthquake‐driven uplift by using a forward numerical model. The modeling reveals that terraces are generally abandoned at SL highstands but they are carved during all stands, depending on the time spent within the sea erosional‐depth‐range. Therefore sea reoccupation of a same platform after formation is a common occurrence that decreases with increasing uplift rates, suggesting that most platforms in nature may be in fact polygenetic. Furthermore, the model run time influences the terrace sequences: terraces formed at the beginning of longer runs constitute an ‘inherited morphology’ affecting subsequent sequences. When coseismic uplift is applied, the formation and preservation of terraces for a given average uplift rate depend stochastically on the coseismic displacement ‐ recurrence interval combination in relation to the SL position at the time of the earthquake. These factors significantly contribute to a higher likelihood of non‐preserved terraces along a terrace sequence, which may affect age correlation and, consequently, the resulting uplift rates. Further research is needed to explore the effect of the full seismic cycle in shaping a terrace sequence.

Transformation, leaching and plant uptake simulations of 6:2 and 8:2 polyfluoroalkyl phosphate diesters (diPAPs) and related transformation products under near-natural conditions

In response to the growing concern over PFAS contamination, employing models to simulate PFAS behavior in the environment becomes necessary. This facilitates evaluating risks tied to leaching into groundwater, adsorption in soil, plant uptake, entry into the food chain, and the conversion of precursors into persistent PFAS. We utilized the MACRO model to simulate the behavior of the precursors 6:2 diPAP and 8:2 diPAP using data from a 2-year lysimeter experiment, key compound parameters were optimized via the caRamel evolutionary algorithm. We assumed that the transformation of both diPAP precursors into stable PFAAs is influenced by temperature and soil moisture, similar to pesticide degradation by microorganisms. Results reveal that the model accurately represents transformation, leaching, soil retention, and plant uptake of diPAP and transformation products. A comparison with a lab-based soil column study supports the slower natural degradation of precursors, affirming our modeling approach. Temperature and soil moisture could indicate that a worst-case scenario for transformation product leaching into groundwater could occur during a mild summer with moderate evapotranspiration and heavy rainfall. Plant uptake involves multiple elements: PFAS availability in the root zone depends on prior degradation or presence. Increased moisture in the root zone favors PFAS uptake combined with temperatures high enough for prior biotransformation. The calculation of temperature and moisture-based conversion rates was adopted directly from MACRO. It is recommended to further investigate these effects to validate and possibly modify them.

The uncertainty of the calculative value of the volumetric flow rate in open channel

The well-known Manning formula is usually used for the calculation of the calculative volumetric flow rate in a river or open canal. The discharge depends on the geometry of the channel, i.e. the water area, the wetted perimeter and the slope, as well as on the roughness coefficients. All these quantities are determined with some uncertainty. The article proposes a methodology for calculating the uncertainty of the roughness coefficients of the riverbed and the floodplain as well as the uncertainty of the geometric dimensions of the riverbed. Then, the method of calculating the uncertainty of the calculative discharge is then given. If these uncertainties are taken into consideration in the process of discharge calculation, then, as has been demonstrated for a hypothetical river channel, the ratio of the uncertainty to the calculated value of the discharge will change from several dozen percent in case of small flows to about ten percent in case of big, flood flows. It has also been shown that the uncertainty of the roughness coefficients has the biggest influence on the uncertainty of the flow rate. The presented calculations show that in order to take into account the influence of uncertainty of linear dimensions and roughness coefficients, the engineer designing the riverbed should assume for the calculations the flow rate increased by 10% then design flow. The obtained results can be used for homogeneous flows only, which is usually assumed in practical engineering calculations.