تحلیل حساسیت پارامترهای هیدرولیکی انتقال فلزات سنگین کادمیم، نیکل و روی در ستون‌های خاک لومی دست‌خورده و دست‌نخورده

تحلیل حساسیت برای تحلیل رفتار مدل تحت شرایط مختلف کاربرد دارد. نتایج تحلیل حساسیت به‌منظور تعیین دقت قابل‌قبول در اندازه‌گیری داده‌های ورودی مهم است. در این پژوهش، انتقال فلزات سنگین کادمیم، نیکل و روی در ستون‌های خاک لومی دست‌خورده و دست‌نخورده با غلظت‌های اولیه‌ 50، 100 و 150 میلی‌گرم بر لیتر با استفاده از مدل hydrus-1d شبیه‌سازی و پارامترهای انتقال املاح نظیر ضریب پخش (d)، ضریب توزیع (kd) و ضریب پراکندگی (β) با استفاده از پارامترهای هیدرولیکی خاک و داده‌های مربوط به غلظت فلز کادمیم، نیکل و روی به روش مدل‌سازی معکوس برآورد شد. بررسی منحنی‌های رخنه‌ حاصل از شبیه‌سازی مدل hydrus-1d و داده‌های اندازه‌گیری ‌شده نشان‌دهنده‌ تطابق بالای منحنی‌های شبیه‌سازی با داده‌های اندازه‌گیری بود. برازش مدل در خاک دست‌خورده نسبت به خاک دست‌نخورده بهتر بود که این امر ممکن است به‌دلیل بهم‌خوردگی ساختمان خاک و افزایش سطح تماس ذرات خاک در خاک دست‌خورده و وجود ناهمگنی ذرات در ستون خاک دست‌نخورده باشد. بر اساس نتایج تحلیل حساسیت، پارامترهای پراکندگی (β) و ضریب پخش (d) به‌ترتیب با دامنه‌ حساسیت بیش‌تر از 5/1 و کم‌تر از 3/0 بیش‌ترین و کم‌ترین حساسیت را نسبت به تغییرات در هر دو نوع خاک دست‌خورده و دست‌نخورده نشان داد. در کل روند حساسیت پارامترهای مدل به‌صورت زیر بود: پارامتر پراکندگی> ضریب توزیع> ضریب پخش. به‌عبارتی با توجه به تاثیر قابل‌توجه تغییرات β بر مقادیر خروجی مدل، در صورت تعیین عملی و آزمایشگاهی این پارامتر بایستی با دقت بیش‌تری اندازه‌گیری شود و در مقابل می‌توان از خطاهای اندازه‌گیری پارامتر d چشم‌پوشی کرد. درجه‌ حساسیت پارامترهای مستقل از میزان غلظت اولیه‌ عناصر بود.

sensitivity analysis of hydraulic parameters of transport heavy metals cd, ni, and zn in disturbed and undisturbed loamy soil columns

introduction soil contamination due to heavy metals is a global environmental issue. one vital aspect for understanding the impact of a contaminant in porous media is to describe their transport behavior using appropriate models. the governing equations for solute transport in soil consist of the convection–dispersion equation (cde) and the mobile–immobile model (mim). mathematical models are usually used to evaluate solute transport in porous media. the first model used to express the transport of solutes and pollutants in porous media is cde it provides acceptable and satisfactory results in homogeneous soils in laboratory tests. hydrus-1d is a modeling environment for simulating water, heat, and solute movement in one-dimensional variably saturated media. sensitivity analyses and model identification are standard approaches in modeling applications to investigate the relative importance of model components that control the system’s behavior. the sensitivity analysis is applied to identify the parameters that influence the model performance most. the sensitivity analysis is defined as the rate of variation in the model outputs due to changes in the input parameters. this study is a fundamental practice for analyzing the behavior of a model under different conditions of an application. the sensitivity analysis could be a practical and powerful tool for investigating the role and importance of model components, such as parameters and forcing data on the model responses. materials and methods the loamy soil samples were collected in both disturbed and undisturbed forms from a farm in the qaramalek area with appropriate humidity located in western tabriz, iran, at 38º 5’ 59.89’’ north and 45º 12’ 38.57’’ east. to determine and present breakthrough curves, concentration values are required throughout the laboratory columns at different times. to simulate the cde model, hydrus software was used. solute transport parameters such as diffusion coefficient (d), distribution coefficient (kd), and dispersion coefficient (β) were estimated using soil hydraulic parameters and data related to the metal concentration of cadmium, nickel, and zinc by an inverse modeling method. a sensitivity analysis was carried out for the identification of the most influential factors on the model output. this method examines the impact of input data on a given model and its actual conditions. in line with this purpose, in each run, one input data was changed to a value equal to positive and negative five to 15%, and the other input data was kept constant. to identify the effect of the input parameters of a given model on its output, the sensitivity analysis for the hydrus model was utilized. the parameters of hydrodynamic dispersion coefficient (d), distribution coefficient (kd), and spreading parameter (β) were changed between five to 15 %. sensitivity analysis was carried out on cadmium, nickel, and zinc metals with densities equal to 50, 100, and 150 mg.l-1 in two disturbed and undisturbed soils. results and discussion examining the breakthrough curves of cadmium in disturbed and undisturbed soils shows that the fitted curves using the hydrus model and the measured curve almost coincide with each other, which is more obvious in disturbed soils. it should be noted that the model fits better in the disturbed soil than in the undisturbed soil. this may be due to the disruption of the structure the increase in the contact surface of the particles in the disturbed soil and the presence of heterogeneity in the undisturbed soil column. the simulation results show the transport of heavy metals (zn, ni, cd) and hydrus output have the highest and the lowest sensitivity to dispersion coefficient β and diffusion coefficient (d), respectively. in general, the impact of input parameters can be reported as follows: spreading parameter (β) > distribution coefficient (kd) > dispersion coefficient (d). therefore, it can be observed that d has a negligible effect on the model results; and consequently, measurement errors can be ignored. conclusion sensitivity analysis is used to analyze model behavior under different conditions. this analysis is used to investigate the relative importance of model components that control the system’s behavior. in this research, the transfer of hydraulic parameters of heavy metals cd, ni, and zn in disturbed and undisturbed loam soil columns with initial concentrations of 50, 100, and 150 mg.l-1 was performed under the simulation of the hydrus-1d model. the comparison of the simulated btcs of the hydrus-1d model and the measured data indicates a high agreement between the simulation curves and the measured data. solute transport parameters such as hydrodynamic dispersion coefficient (d), distribution coefficient (kd), and spreading parameter (β) were estimated using soil hydraulic parameters and data related to cd, ni, and zn metal concentration by inverse modeling method. based on the results of sensitivity analysis, the spreading parameter (β) and hydrodynamic dispersion coefficient (d) had the highest and lowest sensitivity, respectively. in other words, due to the significant effect of β changes on the output values of the model, this parameter should be measured more accurately and on the other hand, the measurement errors of parameter d can be ignored. the degree of sensitivity of the parameters was independent of the initial concentration of the elements.