3D magnetotelluric modelling including surface topography

An edge finite‐element method has been applied to compute magnetotelluric (MT) responses to three‐dimensional (3D) earth topography. The finite‐element algorithm uses a single edge shape function at each edge of hexahedral elements, guaranteeing the continuity of the tangential electric field while conserving the continuity of magnetic flux at boundaries. We solve the resulting system of equations using the biconjugate gradient method with a Jacobian preconditioner. The solution gives electric fields parallel to the slope of a surface relief that is often encountered in MT surveys. The algorithm is successfully verified by comparison with other numerical solutions for a 3D‐2 model for comparison of modelling methods for EM induction and a ridge model. We use a 3D trapezoidal‐hill model to investigate 3D topographic effects, which are caused mainly by galvanic effects, not only in the Zxy mode but also in the Zyx mode. If a 3D topography were approximated by a two‐dimensional topography therefore errors occurring in the transverse electric mode would be more serious than those in the transverse magnetic mode.     

A hybrid near-field/far-field thermal discharge model for coastal areas

A hybrid technique has been used to simulate the dispersion of heat from surface discharges in coastal areas. Characteristics of the near field thermal dispersion are described by the CORMIX3 model. A two-dimensional harmonic finite element hydrodynamic model (TEA) and a Eulerian–Lagrangian transport model (ELA) are applied for the far-field computation. A Gaussian puff algorithm in ELA, which represents the near field plume as a series of patches, is used to link the two regimes. The computed results are compared to available field measurements. Very reasonable agreement is observed.     

Development of updateable model output statistics (UMOS) system for air temperature over South Korea

In this study, a Updateable Model Output Statistics (UMOS) system has been developed for the forecast of 3-h temperature over South Korea using two significantly different models’ (Regional Data Assimilation and Prediction System (RDAPS) and Korea Meteorological Administration (KMA) Weather Research and Forecasting (WRF) model (KWRF)) outputs based on the Canadian UMOS system (Wilson and Vallee, 2002; 2003). The UMOS system is designed to consider the local climatology and the model’s forecasting skills. The 20 most frequently selected potential predictors for each season, station, and forecast projection time from the 67 potential predictors of the Model Output Statistics (MOS) system, were used as potential predictors of the UMOS system. The UMOS equations are developed by a weighted blending of the new and old model data, with weights chosen to emphasize the new model data while including enough old model data in the development to ensure stable equations and a smooth transition to dependency on the new model. The UMOS equations were updated regularly at a predefined time interval to consider the changes of covariance structure between the new model output and observations as the new model data increase. The validation results showed that seasonal mean bias, Root Mean Square Error (RMSE), and correlation coefficients for the total forecast projection times are ?0.379～0.055°C, 1.951～2.078°C, and 0.741～0.965, respectively. Although, the forecasting skills of UMOS system are very consistent without regard to the season and geographic location, the performance is slightly better in autumn and winter than in spring and summer, and better in coastal regions than in inland region. When we take into account the significant differences of the RDAPS and KWRF, the UMOS system can be used as a supplementary forecasting tool of the MOS system for 3-h temperature over South Korea. However, the UMOS system is very sensitive to the selected number and/or types of predictors. Therefore, more work is needed to enable the use of the UMOS system in operation, including tuning of the number and types of potential predictors and automation of the updating processes of the UMOS equations.     

National-scale assessment of landslide susceptibility to rank the vulnerability to failure of rock-cut slopes along expressways in Korea

The objective of this study is to perform a preliminary national-scale assessment of the landslide susceptibility of rock-cut slopes along expressways in Korea. A geographic information system (GIS) database was compiled based on data from topographical and geological maps, and rock-cut slope data, including the locations of past landslides. Seven factors (i.e., slope height, slope length, slope gradient, upper slope gradient, lithology, distance from nearest fault, and dip direction of slope) were extracted from the GIS database to assess the relationship between each factor and landslide events. Weight of evidence (WOE), analytic hierarchy process (AHP), and fuzzy logic methods, as well as hybrid methods, were used to establish the rating of classes for each factor, weightings for the factors, and to combine multiple factor layers into landslide-susceptibility maps. A comparison of the results obtained using several different methods, based on the area under curve technique, revealed that the WOE method showed the highest accuracy of 74%. The annual cost of traffic congestion resulting from slope failures was evaluated to identify those rock-cut slopes where detailed investigations and landslide warning systems are required.     

Three-dimensional inversion of magnetotelluric data including sea effects obtained in Pohang, Korea

Magnetotelluric (MT) surveys were conducted in Pohang, Korea, for low-temperature geothermal exploration in 2002 and 2003. Pohang is located in the southeastern part of the Korean Peninsula and close to the East Sea. In the interpretation of MT data from a coastal environment, sea effects must be correctly included because seawater is a strong conductor. We first constructed a five-layered earth model with a realistic coastline and bathymetry to investigate sea effects on MT data measured in Pohang. This model clearly shows that the Pohang data are significantly influenced by sea water at frequencies blow 1 Hz at the whole measurement sites. Next, we utilized a three-dimensional inversion algorithm based on the Gauss–Newton approach to produce a reliable resistivity model. Seawater is excluded from the inversion domain to fix the resistivity, while included in the modeling domain to simulate sea effects on MT responses. Blocks for the sub-seafloor are included in unknown parameters since they are sufficiently close to the survey area to affect MT responses in Pohang. Static shifts are also considered in inversion for more accurate interpretation. The rms data misfit is smoothly reduced from 11.2 to 1.87 after 7 iterations. The resulting resistivity model shows a pattern of low–high–low resistivity with depth. The model is compatible with resistivity logs obtained from four boreholes in the survey area, and can explain major geological features in Pohang.     

Separation of incident and reflected waves in wave–current flumes

A technique is developed to separate the incident and reflected waves propagating on a known current in a laboratory wave–current flume by analyzing wave records measured at two or more locations using a least squares method. It can be applied to both regular and irregular waves. To examine its performance, numerical tests are made for waves propagating on quiescent or flowing water. In some cases, to represent the signal noise and measurement error, white noise is superimposed on the numerically generated wave signal. For all the cases, good agreement is observed between target and estimation.     

An analytic solution to the mild slope equation for waves propagating over an axi-symmetric pit

An analytic solution to the mild slope equation is derived for waves propagating over an axi-symmetric pit located in an otherwise constant depth region. The water depth inside the pit decreases in proportion to an integer power of radial distance from the pit center. The mild slope equation in cylindrical coordinates is transformed into ordinary differential equations by using the method of separation of variables, and the coefficients of the equation in radial direction are transformed into explicit forms by using the direct solution for the wave dispersion equation by Hunt (Hunt, J.N., 1979. Direct solution of wave dispersion equation. J. Waterw., Port, Coast., Ocean Div., Proc. ASCE, 105, 457–459). Finally, the Frobenius series is used to obtain the analytic solution. Due to the feature of the Hunt's solution, the present analytic solution is accurate in shallow and deep waters, while it is less accurate in intermediate depth waters. The validity of the analytic solution is demonstrated by comparison with numerical solutions of the hyperbolic mild slope equations. The analytic solution is also used to examine the effects of the pit geometry and relative depth on wave transformation. Finally, wave attenuation in the region over the pit is discussed.     

Changes in spatial variations of sap flow in Korean pine trees due to environmental factors and their effects on estimates of stand transpiration

It is difficult to scale up measurements of the sap flux density (Js) for the characterization of tree or stand transpiration (E) due to spatial variations in JS and their temporal changes. To assess spatial variations in the sap flux density of Korean pine (Pinus koraiensis) and their effects on E estimates, we measured the Js using Granier-type sensors. Within trees, the Js decreased exponentially with the radial depth, and the Js of the east aspects were higher than those of the west aspects. Among trees, there was a positive relationship between Js and the tree diameter at breast height, and this positive relationship became stronger as the transpiration demand increased. The spatial variations that caused large errors in E estimates (i.e., up to 110.8 % when radial variation was ignored) had varied systematically with environmental factors systematic characteristics in relation to environmental factors. However, changes in these variations did not generate substantial errors in the E estimates. For our study periods, the differences in the daily E (E _D) calculated by ignoring radial, azimuthal and tree-to-tree variations and the measured E _D were fairly constant, especially when the daily vapor pressure deficit (D__D) was higher than 0.6 kPa. These results imply that the effect of spatial variations changes on sap flow can be a minor source of error compared with spatial variations (radial, azimuthal and tree-to-tree variations) when considering E estimates.