A mathematical model is developed for investigating time dependent surface deformations of a hydrostatic water volume, when it is subjected to a sudden partial collapse or rise of the sea bottom. The model solves two-dimensional Navier-Stokes Equations on a vertical plane numerically by using Marker and Cell Method (MAC) for viscous and compressible fluid including all the nonlinear effects in the solution. For demonstration, a vertical motion was given to a section in a hypothetical reservoir bed within a short time period and the resulting velocity and pressure fields and the surface profile of the water body are obtained. Computational and physical aspects are discussed.

Volcanic subsidence, caused by partial emptying of magma in the subsurface reservoir has long been observed by spaceborne radar interferometry. Monitoring long-term crustal deformation at the most notable type of volcanic subsidence, caldera, gives us insights of the spatial and hazard-related information of subsurface reservoir. Several subsiding calderas, such as volcanoes on the Galapagos islands have shown a complex ground deformation pattern, which is often composed of a broad deflation signal affecting the entire edifice and a localized subsidence signal focused within the caldera floor. Although numerical or analytical models with multiple reservoirs are proposed as the interpretation, geologically and geophysically evidenced ring structures in the subsurface are often ignored. Therefore, it is still debatable how deep mechanisms relate to the observed deformation patterns near the surface. We aim to understand what kind of activities can lead to the complex deformation. Using two complementary approaches, we study the three-dimensional geometry and kinematics of deflation processes evolving from initial subsidence to later collapse of calderas. Firstly, the analog experiments analyzed by structure-from-motion photogrammetry (SfM) and particle image velocimetry (PIV) helps us to relate the surface deformation to the in-depth structures. Secondly, the numerical modeling using boundary element method (BEM) simulates the characteristic deformation patterns caused by a sill-like source and a ring-fault. Our results show that the volcano-wide broad deflation is primarily caused by the emptying of the deep magma reservoir, whereas the localized deformation on the caldera floor is related to ring-faulting at a shallower depth. The architecture of the ring-fault to a large extent determines the deformation localization on the surface. Since series evidence for ring-faulting at several volcanoes are provided, we highlight that it is vital to include ring-fault activity in numerical or analytical deformation source formulation. Ignoring the process of ring-faulting in models by using multiple point sources for various magma reservoirs will result in erroneous, thus meaningless estimates of depth and volume change of the magmatic reservoir(s).

Subsurface constructions, such as tunnels, create hydrogeological challenges in mitigating risk of subsidence due to groundwater drawdown. Presenting readily made precautionary mitigation plans, such as strategically planned artificial recharge applications, can help effectivise the mitigation process. The Bypass Stockholm project comprises of several subsurface constructions which may lower the surrounding groundwater level through tunnel leakage. Risk of land subsidence persists in the nearby urban area of Vinsta, Stockholm, where a groundwater drawdown may cause the clays in the area to experience land subsidence. A hydrogeological modelling approach was used in the area to create strategic artificial infiltration plans that could be employed as a mitigative response to the drop in groundwater head. In order to simulate the potential tunnel drainage, a steady state hydrogeological model was built using MODFLOW. A 220 l/s tunnel leakage was then simulated. Four different artificial groundwater infiltration scenarios were then conceptualized and simulated to observe effects on groundwater heads. The groundwater levels of the baseline model of the area fit the calibration targets with average absolute deviation of 0.18 m. The tunnel drainage scenario lowered the groundwater level in the till aquifer and bedrock by 0 - 1.5 m and 0.5 - 5 m respectively, with higher drawdowns observed closer to the tunnel. The infiltration scenarios mitigate the groundwater drawdown with different efficacies; proximity to the recharge point, and discharge into the till aquifer were observed to have the highest effect on groundwater recharge in the model. The model could have been improved by improving the data quality surrounding the hydraulic conductivity of the bedrock, as it had the highest effect according to the parameter sensitivity analysis. / Konstruktioner under mark kan skapa hydrogeologiska utmaningar, såsom sättningsrisk orsakade av grundvattenavsänkning. Ett sätt att effektivisera åtgärdsprocessen är att förbereda för eventuell artificiell grundvatteninfiltration. Vägprojektet Förbifart Stockholm innefattar konstruktioner under mark och riskerar, genom inläckage, att sänka grundvattennivån i omgivningen. Ett potentiellt problemområde är stadsdelen Vinsta, delar av vars är byggd på sättningskänslig lera som kan påverkas av en grundvattenavsänkning. För att kunna motverka en grundvattensänkning i Vinsta har hydrogeologisk modellering utförts för att strategiskt planera artificiell grundvatteninfiltration. Ett tunnelläckage på 220 l/s har simulerats genom en hydrogeologisk steady state-modell i MODFLOW. Fyra olika scenarier för grundvatteninfiltration har konceptualiserats och simulerats för att observera påverkan på grundvattennivån. Den spatialt variabla grundvattennivån i grundmodellen nådde kalibreringsmålen med en genomsnittlig absolutavvikelse på 0,18 m. Modellen för tunnelläckage resulterade i att grundvattennivån i moränakvifären och berget sjönk med 0 – 1,5 resp. 0,5 – 5 m, med större grundvattensänkning närmare tunneln. Scenarierna för infiltration motverkade grundvattensänkningen i olika grad. Närhet tilltunneln, eller platsen för inläckage, samt den hydrauliska konduktiviteten mellan infiltrationen och akvifären visade störst påverkan på resultatet för att motverka grundvattensänkningen. Känslighetsanalysen för parametrarna i modellen visade att berget och dess hydrauliska konduktivitet hade störst påverkan på resultatet. Tillgång till bättre data för berget möjliggör förbättrat modelleringsresultat.

No description available.

This thesis develops a two dimensional, axisymmetric finite element package for solving continuum problems including rubbleization subsidence and nonlinear fracture mechanics. This package includes both a user friendly preprocessor, PRESCRUBS.BYU, and a versatile analysis code SCRUBS.BYU. PRESCRUBS.BYU systematically creates the data file necessary to run SCRUBS.BYU. SCRUBS.BYU provides many options of nonlinear static analysis using either linear or quadratic isoparametric finite elements. Sample problems are presented that demonstrate the capabilities of this package.

Thermokarst is a typical process that indicates widespread permafrost degradation in yedoma landscapes. The Lena-Aldan interfluvial area in Central Yakutia in eastern Siberia is now facing extensive landscape changes with surface subsidence due to thermokarst development during the past few decades. To clarify the spatial extent and rate of subsidence, multiple spatial datasets, including GIS and remote sensing observations, were used to analyze the Churapcha rural locality, which has a typical yedoma landscape in Central Yakutia. Land cover classification maps for 1945 and 2009 provide basic information on anthropogenic disturbance to the natural landscape of boreal forest and dry grassland. Interferometric synthetic aperture radar (InSAR) with ALOS-2/PALSAR-2 data revealed activated surface subsidence of 2 cm/yr in the disturbed area, comprising mainly abandoned agricultural fields. Remote sensing with an unmanned aerial system also provided high-resolution information on polygonal relief formed by thermokarst development at a disused airfield where InSAR analysis exhibited extensive subsidence. It is worth noting that some historically deforested areas have likely recovered to the original landscape without further thermokarst development. Spatial information on historical land-use change is helpful because most areas with thermokarst development correspond to locations where land was used by humans in the past. Going forward, the integrated analysis of geospatial information will be essential for assessing permafrost degradation.

The overall goal of this dissertation is to determine and develop optimal strategies for inversely calibrating transmissivities (T), elastic and inelastic skeletal storage coefficients (Ske and Skv) of the developed-zone aquifer and conductance (CR) of the basin-fill faults for the entire Las Vegas basin, and to investigate future trends of land subsidence in Las Vegas Valley. This dissertation consists of three separate stand-alone chapters. Chapter 2 presents a discrete adjoint parameter estimation (APE) algorithm for automatically identifying suitable hydraulic parameter zonations from hydraulic head and subsidence measurements. Chapter 3 compares three different inversion strategies to determine the most accurate and computationally efficient method for estimating T and Ske and Skv at the basin scale: the zonation method (ZM), the adaptive multi-scale method and the Differential Evolution Adaptive Metropolis Markov chain Monte Carlo scheme (DREAM MCMC). Chapter 4 outlines a fine-scale numerical model capable of capturing far more hydrologic detail than any previously developed model of Las Vegas Valley The new model is calibrated using high-resolution InSAR data and hydraulic head data from 1912 to 2010. The calibrated model is used to investigate the influence of faults and their potential role on influencing clay thicknesses and land subsidence distributions, and to investigate future trends of land subsidence in Las Vegas Valley. / Ph. D.

台灣地區地層下陷問題肇始於六十年代迄今已逾三十餘載，持續下陷面積已達1,194平方公里，最大累積下陷量達到三公尺以上。而目前地層下陷地表監測所採用的傳統水準測量耗費人力、時間較多，且不易獲得連續和全面性之地層下陷資料，目前國內水利單位限於人力時間，無法針對所有監測區域每年皆施測一次。近年來由於空載光達測量技術興起，其具有短時間內獲取大區域高密度、高精度高程資料的特性，因此本研究之目的在探討如何利用空載光達測量技術快速獲取高精度之三維點雲資訊，進行大區域的地層下陷監測及其成效。 研究方法係首先將監測區內掃瞄的光達點雲資料進行網格化分群；接著，計算網格區域內所有光達點擬合平面的中心高程；然後，以人工或自動方法萃取出平坦、穩固的網格區域做為監測面；最後，進行不同時期網格監測面高程差異之統計測試分析，以求出地層下陷量。 實驗結果顯示改善點雲高程精度至5公分以內後，經由網格監測面的精度、坡度、坡向、反射強度、道路範圍等為門檻值，可萃取出80%以上正確率的穩固監測面，且其高差成果與長期監測成果的平均值差異在1.3公分至2.9公分之間，由此成果可以說明本研究成果對建立一套省時省力的監測模式，進而達到地層下陷監測自動化的目的有相當幫助。 / The issue of land subsidence in Taiwan has been concerned for over 30 years since 1970. Land subsidence area has been already over 1194 km2, the maximum amount of accumulative subsidence is more than 3 meters. The conventional leveling for the land subsidence monitoring is labor-intensive and time-consuming, so that the Water Resources Agency could not monitor all the subsidence area every year. Airborne LIDAR technology was developed in recent years, it has the characteristics of collecting 3-D point data at the high density and high elevation accuracy in short time. The purpose of this study, therefore, is to discuss how to utilize the airborne LIDAR technology to monitor the land subsidence. The proposed approach, therefore, is first to divide the collecting DSM points in the monitor area into regular grids. Secondly, all the points in the regular grids are fitted to one set planar parameters by least squares principle and the centric elevation of each grid is calculated. Third, the flatness and well-defined planar grids are selected as the monitoring surfaces with the manual or automatic method. Finally, the difference of centric elevation in each monitoring surfaces at different period is calculated and analyzed with statistical approach. This study shows that after refining the elevation accuracy of point clouds within 5 cm, our approach can extract stable monitoring surfaces by limiting planar fitting accuracy, flatness, slope, intensity, or by using road information. The extracted correct rate can be more than 80%. The discrepancy of elevation difference between this study and long-term monitoring result is between 1.3 cm and 2.9 cm. It proves the proposed approach is helpful on constructing the monitoring model in timesaving and efficient way, and our proposed approach has the potential for developing automatic land subsidence monitoring method.

Com a finalidade de estudar as perturbações médias sobre a Plataforma Continental Sudeste (PCSE) causadas pela passagem de uma frente fria durante o verão, obteve-se um estado básico médio de correntes. Esse campo médio foi obtido com o POM (Princeton Ocean Model) através de simulação diagnóstica, seguida por curta rodada prognostica. As correntes modeladas para o verão na Plataforma Continental Sudeste assemelharam-se em alguns aspectos a correntes obtidas em fundeios de correntógrafos realizados pela Diretoria de Hidrografia e Navegação (DHN), e pelo Instituto Oceanográfico da Universidade de São Paulo, durante os projetos Ciculação Oceânica na Região Oeste do Atlântico Sul (COROAS) e Dinâmica do Ecossistema de Plataforma da Região Oeste do Atlântico Sul (DEPROAS), validando os resultados do modelo Durante a passagem da frente fria as correntes invertem em cerca de seis a nove horas após a passagem do sistema frontal por cada região da PCSE. As etapas da subsidência costeira foram modeladas, ocorrendo transporte de deriva do vento em direção à costa, na camada superficial (0-10m), e transporte de compensação ao longo da coluna d\'água, havendo empilhamento de água na região costeira. A força de gradiente de pressão barotrópica desenvolvida implicou na formação e intensificação do jato costeiro em direção a menores latitudes. Durante esse estado transiente associado à passagem do sistema frontal, a Água Tropical aproxima-se da costa em nível sub-superficial, misturando-se com as águas da plataforma continental A porção da Plataforma Continental Sudeste localizada ao sul da Ilha de São Sebastião é a que responde com maior intensidade à passagem do sistema meteorológico frontal. Na direção transversal à costa, predomina o balanço geostrófico, o qual é intensificado durante a passagem da frente fria, enquanto que ao longo da costa as águas são desviadas do referido balanço pela tensão de cisalhamento do vento e pela aceleração local. Em Cananéia, ocorre a maior amplitude de variação do nível do mar, que pode estar associada a um mecanismo de ressonância entre os ventos do sistema frontal e Ondas de Plataforma Continental geradas pelo vento, uma vez que possuem velocidades de propagação próximas. A plataforma continental adquire vorticidade relativa ciclônica durante a passagem da frente fria. / In order to study the pertubations caused on the South Brazil Bigth (SBB) caused by the passage of a cold front during Summer, a mean basic state for currents was obtained using the POM (Princeton Ocean Model) through diagnostic simulation, followed by a short round prognostic. The mean currents modeled are in certain instancies similar to the currents measured by meter moorings deployed by the Directorate of Hydrography and Navigation (DHN), and by Oceanographic Institute of the University São Paulo during the Projects Oceanic Circulation in the Western South Atlantic (COROAS) and Dynamic of the Ecosystem of the Shelf of the Western South Atlantic Region (DEPROAS), validating results model. During the cold front passage, the currents invert six to nine hours after the passage of the frontal system through each region of the SBB. The stages of coastal subsidence were modeled, with coastward wind drift transport at superficial layer and offshore compensation transport throughout the water column, with piling of water in the coastal region. The strength of the barotropic pressure gradient developed caused the formation and intensification of a coastal jet. During this transient stage associated with the passage of the frontal system, the Tropical Water approaches the coast within sub-surface layer, mixing with continental shelf waters. The portion of SBB located to the South of São Sebastião Island is the portion which is most influenced by the passage of the frontal system. In a normal direction to the coast, the geostrophic balance predominates, being intensified during the passage of the cold front, while along the coast, the waters is diverted from this balance by the wind stress and by the local acceleration, which tend to have opposite signs. Near Cananéia, there is the highest sea level elevation, which may be associated with a mechanism of resonance between the winds of the frontal system and the Continental Shelf Waves generated by the wind, since they have the same speed of propagation. The continental shelf acquires cyclonic relative vorticity during the passage of the cold front.

Com a finalidade de estudar as perturbações médias sobre a Plataforma Continental Sudeste (PCSE) causadas pela passagem de uma frente fria durante o verão, obteve-se um estado básico médio de correntes. Esse campo médio foi obtido com o POM (Princeton Ocean Model) através de simulação diagnóstica, seguida por curta rodada prognostica. As correntes modeladas para o verão na Plataforma Continental Sudeste assemelharam-se em alguns aspectos a correntes obtidas em fundeios de correntógrafos realizados pela Diretoria de Hidrografia e Navegação (DHN), e pelo Instituto Oceanográfico da Universidade de São Paulo, durante os projetos Ciculação Oceânica na Região Oeste do Atlântico Sul (COROAS) e Dinâmica do Ecossistema de Plataforma da Região Oeste do Atlântico Sul (DEPROAS), validando os resultados do modelo Durante a passagem da frente fria as correntes invertem em cerca de seis a nove horas após a passagem do sistema frontal por cada região da PCSE. As etapas da subsidência costeira foram modeladas, ocorrendo transporte de deriva do vento em direção à costa, na camada superficial (0-10m), e transporte de compensação ao longo da coluna d\'água, havendo empilhamento de água na região costeira. A força de gradiente de pressão barotrópica desenvolvida implicou na formação e intensificação do jato costeiro em direção a menores latitudes. Durante esse estado transiente associado à passagem do sistema frontal, a Água Tropical aproxima-se da costa em nível sub-superficial, misturando-se com as águas da plataforma continental A porção da Plataforma Continental Sudeste localizada ao sul da Ilha de São Sebastião é a que responde com maior intensidade à passagem do sistema meteorológico frontal. Na direção transversal à costa, predomina o balanço geostrófico, o qual é intensificado durante a passagem da frente fria, enquanto que ao longo da costa as águas são desviadas do referido balanço pela tensão de cisalhamento do vento e pela aceleração local. Em Cananéia, ocorre a maior amplitude de variação do nível do mar, que pode estar associada a um mecanismo de ressonância entre os ventos do sistema frontal e Ondas de Plataforma Continental geradas pelo vento, uma vez que possuem velocidades de propagação próximas. A plataforma continental adquire vorticidade relativa ciclônica durante a passagem da frente fria. / In order to study the pertubations caused on the South Brazil Bigth (SBB) caused by the passage of a cold front during Summer, a mean basic state for currents was obtained using the POM (Princeton Ocean Model) through diagnostic simulation, followed by a short round prognostic. The mean currents modeled are in certain instancies similar to the currents measured by meter moorings deployed by the Directorate of Hydrography and Navigation (DHN), and by Oceanographic Institute of the University São Paulo during the Projects Oceanic Circulation in the Western South Atlantic (COROAS) and Dynamic of the Ecosystem of the Shelf of the Western South Atlantic Region (DEPROAS), validating results model. During the cold front passage, the currents invert six to nine hours after the passage of the frontal system through each region of the SBB. The stages of coastal subsidence were modeled, with coastward wind drift transport at superficial layer and offshore compensation transport throughout the water column, with piling of water in the coastal region. The strength of the barotropic pressure gradient developed caused the formation and intensification of a coastal jet. During this transient stage associated with the passage of the frontal system, the Tropical Water approaches the coast within sub-surface layer, mixing with continental shelf waters. The portion of SBB located to the South of São Sebastião Island is the portion which is most influenced by the passage of the frontal system. In a normal direction to the coast, the geostrophic balance predominates, being intensified during the passage of the cold front, while along the coast, the waters is diverted from this balance by the wind stress and by the local acceleration, which tend to have opposite signs. Near Cananéia, there is the highest sea level elevation, which may be associated with a mechanism of resonance between the winds of the frontal system and the Continental Shelf Waves generated by the wind, since they have the same speed of propagation. The continental shelf acquires cyclonic relative vorticity during the passage of the cold front.