Geological structures often exhibit smooth characteristics away from sharp discontinuities. One aim of geophysical inversion is to recover information about the smooth structures as well as about the sharp discontinui...Geological structures often exhibit smooth characteristics away from sharp discontinuities. One aim of geophysical inversion is to recover information about the smooth structures as well as about the sharp discontinuities. Because no specific operator can provide a perfect sparse representation of complicated geological models, hyper-parameter regularization inversion based on the iterative split Bregman method was used to recover the features of both smooth and sharp geological structures. A novel preconditioned matrix was proposed, which counteracted the natural decay of the sensitivity matrix and its inverse matrix was calculated easily. Application of the algorithm to synthetic data produces density models that are good representations of the designed models. The results show that the algorithm proposed is feasible and effective.展开更多
Distinguishing geochemical anomalies from background is a basic task in exploratory geochemistry. The derivation of geochemical anomalies from stream sediment geochemical data and the decomposition of these anomalies ...Distinguishing geochemical anomalies from background is a basic task in exploratory geochemistry. The derivation of geochemical anomalies from stream sediment geochemical data and the decomposition of these anomalies into their component patterns were described. A set of stream sediment geochemical data was obtained for 1 880 km 2 of the Pangxidong area, which is in the southern part of the recently recognized Qinzhou-Hangzhou joint tectonic belt. This belt crosses southern China and tends to the northwest (NE) direction. The total number of collected samples was 7 236, and the concentrations of Ag, Au, Cu, As, Pb and Zn were measured for each sample. The spatial combination distribution law of geochemical elements and principal component analysis (PCA) were used to construct combination models for the identification of combinations of geochemical anomalies. Spectrum-area (S-A) fractal modeling was used to strengthen weak anomalies and separate them from the background. Composite anomaly modeling was combined with fractal filtering techniques to process and analyze the geochemical data. The raster maps of Au, Ag, Cu, As, Pb and Zn were obtained by the multifractal inverse distance weighted (MIDW) method. PCA was used to combine the Au, Ag, Cu, As, Pb, and Zn concentration values. The S-A fractal method was used to decompose the first component pattern achieved by the PCA. The results show that combination anomalies from a combination of variables coincide with the known mineralization of the study area. Although the combination anomalies cannot reflect local anomalies closely enough, high-anomaly areas indicate good sites for further exploration for unknown deposits. On this basis, anomaly and background separation from combination anomalies using fractal filtering techniques can provide guidance for later work.展开更多
Total concentrations of arsenic, lead, cadmium, mercury, nickel, chromium, and copper in the soils from near a coal mine area in southwest Guizhou, China, were measured to evaluate the level of contamination, and the ...Total concentrations of arsenic, lead, cadmium, mercury, nickel, chromium, and copper in the soils from near a coal mine area in southwest Guizhou, China, were measured to evaluate the level of contamination, and the potential ecological risks posed by the heavy metals were quantitatively estimated. Results reveal that all heavy metals/metalloid exceeded the background values for soil environmental quality of heavy metals in Guizhou area. Geo-accumulation index(I_(geo)) showed that arsenic had the highest contamination level(I_(geo)=4) among the seven heavy metals/metalloid, and the contamination levels of mercury and lead were also relatively high(I_(geo)=3). Pearson correlation and cluster analysis identified that mercury, copper and arsenic had a relationship, and their presence might be mainly related to mining activity, coal and oil combustion, and vehicle emissions. Improved Nemerow index indicated that the overall level of heavy metal contamination in the studied area ranged from moderately–heavily contaminated to heavily contaminated level. Potential ecological risk index(R_I) analysis manifested that the whole ecological risk level ranged from high degree to very high degree(325.30≤R_I≤801.02) in the studied soil samples, and the potential ecological risk factors (E_r^i) of heavy metals/metalloid were as follows: Hg > As > Cd > Pb > Cu > Ni > Cr, and the E_r^i of Hg and As reached very high risk grade.展开更多
This paper focuses on understanding and evaluating the dynamic effect of the heavy-haul train system on the seismic performance of a long-span railway bridge. A systematic study on the effect of heavy-haul trains on b...This paper focuses on understanding and evaluating the dynamic effect of the heavy-haul train system on the seismic performance of a long-span railway bridge. A systematic study on the effect of heavy-haul trains on bridge seismic response has been conducted, considering the influence of vehicle modeling strategies and dynamic characteristics of the seismic waves. For this purpose, the performance of a long-span cable-stayed railway bridge is assessed with stationary trains atop it, where the heavy-haul vehicles are modeled in two different ways: the multi-rigid body model with suspension system and additional mass model. Comparison of the bridge response in the presence or absence of the train system has been conducted, and the vehicle loading situation, which includes full-load and no-load, is also discussed. The result shows that during the earthquake, the peak moment of the main girder and peak stress of stay cables increase by 80% and by 40% in the presence of fully loaded heavy-haul trains, respectively. At the same time, a considerable decrease appears in the peak acceleration of the main girder. This proves the existence of the damping effect of the heavy-haul train system, and this effect is more obvious for the fully loaded vehicles. Finally, this paper proposes an efficient vehicle modeling method with 2 degrees of freedom(DOF) for simplifying the treatment of the train system in bridge seismic checking.展开更多
The quality changes of shallow groundwater from karst hilly areas in Guiyang region of China impacted by the urbanization were investigated. The results show that the major ions in shallow groundwater from the karst h...The quality changes of shallow groundwater from karst hilly areas in Guiyang region of China impacted by the urbanization were investigated. The results show that the major ions in shallow groundwater from the karst hilly areas are mainly composed of HCO3^-, SO4^2-, Ca^2＋ and Mg^2＋, and the concentrations scopes of NO3^-, Cl^-, K^＋ and Na^＋ of the groundwater in agricultural, residential and industrial areas are 4.5-9.6, 2.8-7.1, 3.9-6.3 and 2.5-4.9 times higher than those in the forest areas, respectively. The concentrations of As, Pb and Cd of shallow groundwater in the industrial areas are also significantly enhanced, followed by those in the residential areas and the agricultural areas. The concentrations of NO3^-, SO4^2- , As, Pb and Cd of the groundwater in the industrial areas and those of NO3^- , SO4^2-, As and Cd of shallow groundwater in the residential areas reach grade Ⅲ of the Groundwater Quality Standard of China （GB/T 14848-93）, while the concentration of NO3^- in the groundwater from the industrial areas exceeds grade V. With the process of urbanization, NO3^- is the key factor to influence the groundwater quality in karst hilly areas, followed by SO4^2-, As, Pb and Cd.展开更多
The main principle and mathematical model of GOCE kinematic orbit adjustment for Earth gravity field model （EGM） validation and accelerometer calibration are presented. Based on 60 days GOCE kinematic orbits with 1-...The main principle and mathematical model of GOCE kinematic orbit adjustment for Earth gravity field model （EGM） validation and accelerometer calibration are presented. Based on 60 days GOCE kinematic orbits with 1-2 cm accuracy and accelerometer data from 2009-11-02 to 2009-12-31, the RMS-of-fit （ROF） of them using EGM2008, EIGEN-SC, ITG- GRACE2010S and GOCO01S up to 120, 150 and 180 degree and order （d/o） are evaluated and compared. The scale factors and biases of GOCE accelerometer data are calibrated and the energy balance method （EBM） is performed to test the accuracy of accelerometer calibration. The results show that GOCE orbits are also sensitive to EGM from 120 to 150 d/o. The ROFs of EGMs with 150 and 180 d/o are obviously better than those of EGMs with 120 d/o. The ROFs of GOCO01S and ITG-GRACE2010S are almost the same up to 120 and 150 d/o, which are about 3.3 cm and 1.8 cm, respectively. They are far better than those of EGM2008 and EIGEN-SC with the same d/o. The ROF of GOCO01S with 180 d/o is about 1.6 em, which is the best one among those EGMs. The accelerometer calibration accuracies （ACAs） of ITG-GRACE2010S and GOCO01S are obviously higher that those of EGM2008 and EIGEN-SC. The ACA of GOCO01S with 180 d/o is far higher than that of EGMs with 120 d/o, and a little higher than that of ITG-GRACE2010S with 150 d/o. I t is suggested that the newest released EGM such as GOCO01S or GOCO02S till at least 150 d/o should be chosen in GOCE precise orbit determination （POD） and accelerometer calibration.展开更多
In order to find the dynamic response laws of retaining walls affected by certain earthquake loads,the influence of the seismic wave characteristics and sub-grade fill parameters(including the foundation surface slope...In order to find the dynamic response laws of retaining walls affected by certain earthquake loads,the influence of the seismic wave characteristics and sub-grade fill parameters(including the foundation surface slope) were focused on,and a series of tests were performed.The results show that the maximum stress of the retaining wall decreases as internal friction angle,foundation slope,filled soil cohesion and the biggest dynamic elastic modulus increase,while it increases with the seismic frequency and seismic input peak dropping.The addition value of dynamics earth pressure increases when seismic frequency and seismic input peak are reduced,while it decreases when the filled soil cohesion and internal friction angle rise.Meanwhile,dynamic elastic modulus and foundation slope have no obvious influences on addition value of dynamics earth pressure.The slope will be instable if the seismic input peak exceeds 0.5g and be disruptive if seismic frequency is larger than 2.5 Hz.The mid-lower parts of retaining walls are in most heavy and obvious response to these factors,which reveals the mechanism of "belly burst" in retaining wall that appears commonly in practical projects.展开更多
Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique （SIRT） inversion ...Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique （SIRT） inversion is used to deduce the stress redistribution around the longwall mining panel. The mining-induced microseismic events were recorded by mounting an array of receivers on the surface, above the active panel. After processing and filtering the seismic data, the three-dimensional tomography images of the p-wave velocity variations by SIRT passive seismic velocity tomography were provided. To display the velocity changes on coal seam level and subsequently to infer the stress redistribution, these three-dimensional tomograms into the coal seam level were sliced. In addition, the boundary element method （BEM） was used to simulate the stress redistribution. The results show that the inferred stresses from the passive seismic tomograms are conformed to numerical models and theoretical concept of the stress redistribution around the longwall panel. In velocity tomograms, the main zones of the stress redistribution arotmd the panel, including front and side abutment pressures, and gob stress are obvious and also the movement of stress zones along the face advancement is evident. Moreover, the effect of the advance rate of the face on the stress redistribution is demonstrated in tomography images. The research result proves that the SIRT passive seismic velocity tomography has an ultimate potential for monitoring the changes of stress redistribution around the longwall mining panel continuously and subsequently to improve safety of mining operations.展开更多
Although extensive analyses of road segments and intersections located in urban road networks have examined the role of many factors that contribute to the frequency and severity of crashes, the explicit relationship ...Although extensive analyses of road segments and intersections located in urban road networks have examined the role of many factors that contribute to the frequency and severity of crashes, the explicit relationship between street pattern characteristics and traffic safety remains underexplored. Based on a zone-based Hong Kong database, the Space Syntax was used to quantify the topological characteristics of street patterns and investigate the role of street patterns and zone-related factors in zone-based traffic safety analysis. A joint probability model was adopted to analyze crash frequency and severity in an integrated modeling framework and the maximum likelihood estimation method was used to estimate the parameters. In addition to the characteristics of street patterns, speed, road geometry, land-use patterns, and temporal factors were considered. The vehicle hours was also included as an exposure proxy in the model to make crash frequency predictions. The results indicate that the joint probability model can reveal the relationship between zone-based traffic safety and various other factors, and that street pattern characteristics play an important role in crash frequency prediction.展开更多
The effects of soil texture, initial water content and bulk density on diesel oil infiltration in fine sand and silty clay loam materials were evaluated. Three physical and two empirical equations express diesel oil i...The effects of soil texture, initial water content and bulk density on diesel oil infiltration in fine sand and silty clay loam materials were evaluated. Three physical and two empirical equations express diesel oil infiltration through soils with time, with coefficients of determination greater than 0.99. Diesel oil infiltrates more quickly in the fine sand than in the silty clay loam material. Diesel oil infiltration rates are found to decrease with increasing initial water content and bulk density for the silty clay loam material. The infiltration rate of diesel oil in the fine sand material increases slightly with increasing initial water content. The diesel oil saturated conductivity(Kdiesel) decreases with increasing bulk density for the silty clay loam column. Diesel oil sorptivity(S) decreases linearly with increased initial water content and bulk density of the silty clay loam material. Changes in empirical parameters relative to initial water content and bulk density are similar to the parameter S.展开更多
基金Projects(41174061,41374120)supported by the National Natural Science Foundation of China
文摘Geological structures often exhibit smooth characteristics away from sharp discontinuities. One aim of geophysical inversion is to recover information about the smooth structures as well as about the sharp discontinuities. Because no specific operator can provide a perfect sparse representation of complicated geological models, hyper-parameter regularization inversion based on the iterative split Bregman method was used to recover the features of both smooth and sharp geological structures. A novel preconditioned matrix was proposed, which counteracted the natural decay of the sensitivity matrix and its inverse matrix was calculated easily. Application of the algorithm to synthetic data produces density models that are good representations of the designed models. The results show that the algorithm proposed is feasible and effective.
基金Project(1212010071012) supported by Guangdong Pangxidong Mineral Prospect Investigation, ChinaProject(41004051) supported by the National Natural Science Foundation of ChinaProject ([2007]038-01-18) supported by Nationwide Mineral Resource Potential Evaluation Projects of Ministry of Land and Resources, China
文摘Distinguishing geochemical anomalies from background is a basic task in exploratory geochemistry. The derivation of geochemical anomalies from stream sediment geochemical data and the decomposition of these anomalies into their component patterns were described. A set of stream sediment geochemical data was obtained for 1 880 km 2 of the Pangxidong area, which is in the southern part of the recently recognized Qinzhou-Hangzhou joint tectonic belt. This belt crosses southern China and tends to the northwest (NE) direction. The total number of collected samples was 7 236, and the concentrations of Ag, Au, Cu, As, Pb and Zn were measured for each sample. The spatial combination distribution law of geochemical elements and principal component analysis (PCA) were used to construct combination models for the identification of combinations of geochemical anomalies. Spectrum-area (S-A) fractal modeling was used to strengthen weak anomalies and separate them from the background. Composite anomaly modeling was combined with fractal filtering techniques to process and analyze the geochemical data. The raster maps of Au, Ag, Cu, As, Pb and Zn were obtained by the multifractal inverse distance weighted (MIDW) method. PCA was used to combine the Au, Ag, Cu, As, Pb, and Zn concentration values. The S-A fractal method was used to decompose the first component pattern achieved by the PCA. The results show that combination anomalies from a combination of variables coincide with the known mineralization of the study area. Although the combination anomalies cannot reflect local anomalies closely enough, high-anomaly areas indicate good sites for further exploration for unknown deposits. On this basis, anomaly and background separation from combination anomalies using fractal filtering techniques can provide guidance for later work.
基金Project(21467005)supported by the National Natural Science Foundation of China
文摘Total concentrations of arsenic, lead, cadmium, mercury, nickel, chromium, and copper in the soils from near a coal mine area in southwest Guizhou, China, were measured to evaluate the level of contamination, and the potential ecological risks posed by the heavy metals were quantitatively estimated. Results reveal that all heavy metals/metalloid exceeded the background values for soil environmental quality of heavy metals in Guizhou area. Geo-accumulation index(I_(geo)) showed that arsenic had the highest contamination level(I_(geo)=4) among the seven heavy metals/metalloid, and the contamination levels of mercury and lead were also relatively high(I_(geo)=3). Pearson correlation and cluster analysis identified that mercury, copper and arsenic had a relationship, and their presence might be mainly related to mining activity, coal and oil combustion, and vehicle emissions. Improved Nemerow index indicated that the overall level of heavy metal contamination in the studied area ranged from moderately–heavily contaminated to heavily contaminated level. Potential ecological risk index(R_I) analysis manifested that the whole ecological risk level ranged from high degree to very high degree(325.30≤R_I≤801.02) in the studied soil samples, and the potential ecological risk factors (E_r^i) of heavy metals/metalloid were as follows: Hg > As > Cd > Pb > Cu > Ni > Cr, and the E_r^i of Hg and As reached very high risk grade.
基金Project(51678576) supported by the National Natural Science Foundation of ChinaProject(2017YFB1201204) supported by the National Key R&D Program of China。
文摘This paper focuses on understanding and evaluating the dynamic effect of the heavy-haul train system on the seismic performance of a long-span railway bridge. A systematic study on the effect of heavy-haul trains on bridge seismic response has been conducted, considering the influence of vehicle modeling strategies and dynamic characteristics of the seismic waves. For this purpose, the performance of a long-span cable-stayed railway bridge is assessed with stationary trains atop it, where the heavy-haul vehicles are modeled in two different ways: the multi-rigid body model with suspension system and additional mass model. Comparison of the bridge response in the presence or absence of the train system has been conducted, and the vehicle loading situation, which includes full-load and no-load, is also discussed. The result shows that during the earthquake, the peak moment of the main girder and peak stress of stay cables increase by 80% and by 40% in the presence of fully loaded heavy-haul trains, respectively. At the same time, a considerable decrease appears in the peak acceleration of the main girder. This proves the existence of the damping effect of the heavy-haul train system, and this effect is more obvious for the fully loaded vehicles. Finally, this paper proposes an efficient vehicle modeling method with 2 degrees of freedom(DOF) for simplifying the treatment of the train system in bridge seismic checking.
基金Project(2006CB403205) supported by the National Key Basic Research Program of ChinaProject(40961030) supported by the National Natural Science Foundation of China
文摘The quality changes of shallow groundwater from karst hilly areas in Guiyang region of China impacted by the urbanization were investigated. The results show that the major ions in shallow groundwater from the karst hilly areas are mainly composed of HCO3^-, SO4^2-, Ca^2＋ and Mg^2＋, and the concentrations scopes of NO3^-, Cl^-, K^＋ and Na^＋ of the groundwater in agricultural, residential and industrial areas are 4.5-9.6, 2.8-7.1, 3.9-6.3 and 2.5-4.9 times higher than those in the forest areas, respectively. The concentrations of As, Pb and Cd of shallow groundwater in the industrial areas are also significantly enhanced, followed by those in the residential areas and the agricultural areas. The concentrations of NO3^-, SO4^2- , As, Pb and Cd of the groundwater in the industrial areas and those of NO3^- , SO4^2-, As and Cd of shallow groundwater in the residential areas reach grade Ⅲ of the Groundwater Quality Standard of China （GB/T 14848-93）, while the concentration of NO3^- in the groundwater from the industrial areas exceeds grade V. With the process of urbanization, NO3^- is the key factor to influence the groundwater quality in karst hilly areas, followed by SO4^2-, As, Pb and Cd.
基金Project(41174008)supported by the National Natural Science Foundation of ChinaProject(SKLGED2013-4-2-EZ)supported by the Open Foundation of State Key Laboratory of Geodesy and Earth’s Dynamics,ChinaProject(2007B51)supported by the Foundation for the Author of National Excellent Doctoral Dissertation of China
文摘The main principle and mathematical model of GOCE kinematic orbit adjustment for Earth gravity field model （EGM） validation and accelerometer calibration are presented. Based on 60 days GOCE kinematic orbits with 1-2 cm accuracy and accelerometer data from 2009-11-02 to 2009-12-31, the RMS-of-fit （ROF） of them using EGM2008, EIGEN-SC, ITG- GRACE2010S and GOCO01S up to 120, 150 and 180 degree and order （d/o） are evaluated and compared. The scale factors and biases of GOCE accelerometer data are calibrated and the energy balance method （EBM） is performed to test the accuracy of accelerometer calibration. The results show that GOCE orbits are also sensitive to EGM from 120 to 150 d/o. The ROFs of EGMs with 150 and 180 d/o are obviously better than those of EGMs with 120 d/o. The ROFs of GOCO01S and ITG-GRACE2010S are almost the same up to 120 and 150 d/o, which are about 3.3 cm and 1.8 cm, respectively. They are far better than those of EGM2008 and EIGEN-SC with the same d/o. The ROF of GOCO01S with 180 d/o is about 1.6 em, which is the best one among those EGMs. The accelerometer calibration accuracies （ACAs） of ITG-GRACE2010S and GOCO01S are obviously higher that those of EGM2008 and EIGEN-SC. The ACA of GOCO01S with 180 d/o is far higher than that of EGMs with 120 d/o, and a little higher than that of ITG-GRACE2010S with 150 d/o. I t is suggested that the newest released EGM such as GOCO01S or GOCO02S till at least 150 d/o should be chosen in GOCE precise orbit determination （POD） and accelerometer calibration.
基金Project(2006-318-740-20) supported by the West Project from the Department of Transportation of China
文摘In order to find the dynamic response laws of retaining walls affected by certain earthquake loads,the influence of the seismic wave characteristics and sub-grade fill parameters(including the foundation surface slope) were focused on,and a series of tests were performed.The results show that the maximum stress of the retaining wall decreases as internal friction angle,foundation slope,filled soil cohesion and the biggest dynamic elastic modulus increase,while it increases with the seismic frequency and seismic input peak dropping.The addition value of dynamics earth pressure increases when seismic frequency and seismic input peak are reduced,while it decreases when the filled soil cohesion and internal friction angle rise.Meanwhile,dynamic elastic modulus and foundation slope have no obvious influences on addition value of dynamics earth pressure.The slope will be instable if the seismic input peak exceeds 0.5g and be disruptive if seismic frequency is larger than 2.5 Hz.The mid-lower parts of retaining walls are in most heavy and obvious response to these factors,which reveals the mechanism of "belly burst" in retaining wall that appears commonly in practical projects.
文摘Mining operation, especially underground coal mining, always has the remarkable risks of ground control. Passive seismic velocity tomography based on simultaneous iterative reconstructive technique （SIRT） inversion is used to deduce the stress redistribution around the longwall mining panel. The mining-induced microseismic events were recorded by mounting an array of receivers on the surface, above the active panel. After processing and filtering the seismic data, the three-dimensional tomography images of the p-wave velocity variations by SIRT passive seismic velocity tomography were provided. To display the velocity changes on coal seam level and subsequently to infer the stress redistribution, these three-dimensional tomograms into the coal seam level were sliced. In addition, the boundary element method （BEM） was used to simulate the stress redistribution. The results show that the inferred stresses from the passive seismic tomograms are conformed to numerical models and theoretical concept of the stress redistribution around the longwall panel. In velocity tomograms, the main zones of the stress redistribution arotmd the panel, including front and side abutment pressures, and gob stress are obvious and also the movement of stress zones along the face advancement is evident. Moreover, the effect of the advance rate of the face on the stress redistribution is demonstrated in tomography images. The research result proves that the SIRT passive seismic velocity tomography has an ultimate potential for monitoring the changes of stress redistribution around the longwall mining panel continuously and subsequently to improve safety of mining operations.
基金Project(71301083)supported by the National Natural Science Foundation of ChinaProject(2012AA112305)supported by the National High-Tech Research and Development Program of China+1 种基金Project(2012CB725405)supported by the National Basic Research Program of ChinaProject(17208614)supported by the Research Grants Council of the Hong Kong Special Administrative Region,China
文摘Although extensive analyses of road segments and intersections located in urban road networks have examined the role of many factors that contribute to the frequency and severity of crashes, the explicit relationship between street pattern characteristics and traffic safety remains underexplored. Based on a zone-based Hong Kong database, the Space Syntax was used to quantify the topological characteristics of street patterns and investigate the role of street patterns and zone-related factors in zone-based traffic safety analysis. A joint probability model was adopted to analyze crash frequency and severity in an integrated modeling framework and the maximum likelihood estimation method was used to estimate the parameters. In addition to the characteristics of street patterns, speed, road geometry, land-use patterns, and temporal factors were considered. The vehicle hours was also included as an exposure proxy in the model to make crash frequency predictions. The results indicate that the joint probability model can reveal the relationship between zone-based traffic safety and various other factors, and that street pattern characteristics play an important role in crash frequency prediction.
基金Projects(40272108,41402208)supported by the National Natural Science Foundation of ChinaProjects(ZR2012DL05,ZR2015EL044)supported by Shandong Provincial Natural Science Foundation,China+1 种基金Project(4072-114017)supported by Young Teachers’ Development of Shandong University of Technology,ChinaProject(J12LC51)supported by Shandong Province Higher Educational Science and Technology Program,China
文摘The effects of soil texture, initial water content and bulk density on diesel oil infiltration in fine sand and silty clay loam materials were evaluated. Three physical and two empirical equations express diesel oil infiltration through soils with time, with coefficients of determination greater than 0.99. Diesel oil infiltrates more quickly in the fine sand than in the silty clay loam material. Diesel oil infiltration rates are found to decrease with increasing initial water content and bulk density for the silty clay loam material. The infiltration rate of diesel oil in the fine sand material increases slightly with increasing initial water content. The diesel oil saturated conductivity(Kdiesel) decreases with increasing bulk density for the silty clay loam column. Diesel oil sorptivity(S) decreases linearly with increased initial water content and bulk density of the silty clay loam material. Changes in empirical parameters relative to initial water content and bulk density are similar to the parameter S.
