Modelling spatial decisión support systems focused on the development of a data warehouse: applying software engineering

Nowadays, organizations have plenty of data stored in DB databases, which contain invaluable information. Decision Support Systems DSS provide the support needed to manage this information and planning médium and long-term ?the modus operandi? of these organizations. Despite the growing importance of these systems, most proposals do not include its total evelopment, mostly limiting itself on the development of isolated parts, which often have serious integration problems. Hence, methodologies that include models and processes that consider every factor are necessary. This paper will try to fill this void as it proposes an approach for developing spatial DSS driven by the development of their associated Data Warehouse DW, without forgetting its other components. To the end of framing the proposal different Engineering Software focus (The Software Engineering Process and Model Driven Architecture) are used, and coupling with the DB development methodology, (and both of them adapted to DW peculiarities). Finally, an example illustrates the proposal.

Spatial analysis of geometric design consistency and road sight distance

Because of the high number of crashes occurring on highways, it is necessary to intensify the search for new tools that help in understanding their causes. This research explores the use of a geographic information system (GIS) for an integrated analysis, taking into account two accident-related factors: design consistency (DC) (based on vehicle speed) and available sight distance (ASD) (based on visibility). Both factors require specific GIS software add-ins, which are explained. Digital terrain models (DTMs), vehicle paths, road centerlines, a speed prediction model, and crash data are integrated in the GIS. The usefulness of this approach has been assessed through a study of more than 500 crashes. From a regularly spaced grid, the terrain (bare ground) has been modeled through a triangulated irregular network (TIN). The length of the roads analyzed is greater than 100 km. Results have shown that DC and ASD could be related to crashes in approximately 4% of cases. In order to illustrate the potential of GIS, two crashes are fully analyzed: a car rollover after running off road on the right side and a rear-end collision of two moving vehicles. Although this procedure uses two software add-ins that are available only for ArcGIS, the study gives a practical demonstration of the suitability of GIS for conducting integrated studies of road safety.

Spatial aggregation of time-variant stream water ages in urbanizing catchments

Date of Acceptance: 25/03/2015

Large vortex-like structure of dipole field in computer models of liquid water and dipole-bridge between biomolecules

We propose a framework to describe the cooperative orientational motions of water molecules in liquid water and around solute molecules in water solutions. From molecular dynamics (MD) simulation a new quantity “site-dipole field” is defined as the averaged orientation of water molecules that pass through each spatial position. In the site-dipole field of bulk water we found large vortex-like structures of more than 10 Å in size. Such coherent patterns persist more than 300 ps although the orientational memory of individual molecules is quickly lost. A 1-ns MD simulation of systems consisting of two amino acids shows that the fluctuations of site-dipole field of solvent are pinned around the amino acids, resulting in a stable dipole-bridge between side-chains of amino acids. The dipole-bridge is significantly formed even for the side-chain separation of 14 Å, which corresponds to five layers of water. The way that dipole-bridge forms sensitively depends on the side-chain orientations and thereby explains the specificity in the solvent-mediated interactions between biomolecules.

Impact of vaccination on the spatial correlation and persistence of measles dynamics.

The onset of measles vaccination in England and Wales in 1968 coincided with a marked drop in the temporal correlation of epidemic patterns between major cities. We analyze a variety of hypotheses for the mechanisms driving this change. Straightforward stochastic models suggest that the interaction between a lowered susceptible population (and hence increased demographic noise) and nonlinear dynamics is sufficient to cause the observed drop in correlation. The decorrelation of epidemics could potentially lessen the chance of global extinction and so inhibit attempts at measles eradication.

Applying reactive models to column experiments to assess the hydrogeochemistry of seawater intrusion: optimising ACUAINTRUSION and selecting cation exchange coefficients with PHREEQC

Three sets of laboratory column experimental results concerning the hydrogeochemistry of seawater intrusion have been modelled using two codes: ACUAINTRUSION (Chemical Engineering Department, University of Alicante) and PHREEQC (U.S.G.S.). These reactive models utilise the hydrodynamic parameters determined using the ACUAINTRUSION TRANSPORT software and fit the chloride breakthrough curves perfectly. The ACUAINTRUSION code was improved, and the instabilities were studied relative to the discretisation. The relative square errors were obtained using different combinations of the spatial and temporal steps: the global error for the total experimental data and the partial error for each element. Good simulations for the three experiments were obtained using the ACUAINTRUSION software with slight variations in the selectivity coefficients for both sediments determined in batch experiments with fresh water. The cation exchange parameters included in ACUAINTRUSION are those reported by the Gapon convention with modified exponents for the Ca/Mg exchange. PHREEQC simulations performed using the Gains-Thomas convention were unsatisfactory, with the exchange coefficients from the database of PHREEQC (or range), but those determined with fresh water – natural sediment allowed only an approximation to be obtained. For the treated sediment, the adjusted exchange coefficients were determined to improve the simulation and are vastly different from those from the database of PHREEQC or batch experiment values; however, these values fall in an order similar to the others determined under dynamic conditions. Different cation concentrations were simulated using two different software packages; this disparity could be attributed to the defined selectivity coefficients that affect the gypsum equilibrium. Consequently, different calculated sulphate concentrations are obtained using each type of software; a smaller mismatch was predicted using ACUAINTRUSION. In general, the presented simulations by ACUAINTRUSION and PHREEQC produced similar results, making predictions consistent with the experimental data. However, the simulated results are not identical to the experimental data; sulphate (total S) is overpredicted by both models, most likely due to such factors as the kinetics of gypsum, the possible variations in the exchange coefficients due to salinity and the neglect of other processes.

Transferability of Regional Permafrost Disturbance Susceptibility Modelling Using Generalized Linear and Generalized Additive Models

To effectively assess and mitigate risk of permafrost disturbance, disturbance-p rone areas can be predicted through the application of susceptibility models. In this study we developed regional susceptibility models for permafrost disturbances using a field disturbance inventory to test the transferability of the model to a broader region in the Canadian High Arctic. Resulting maps of susceptibility were then used to explore the effect of terrain variables on the occurrence of disturbances within this region. To account for a large range of landscape charac- teristics, the model was calibrated using two locations: Sabine Peninsula, Melville Island, NU, and Fosheim Pen- insula, Ellesmere Island, NU. Spatial patterns of disturbance were predicted with a generalized linear model (GLM) and generalized additive model (GAM), each calibrated using disturbed and randomized undisturbed lo- cations from both locations and GIS-derived terrain predictor variables including slope, potential incoming solar radiation, wetness index, topographic position index, elevation, and distance to water. Each model was validated for the Sabine and Fosheim Peninsulas using independent data sets while the transferability of the model to an independent site was assessed at Cape Bounty, Melville Island, NU. The regional GLM and GAM validated well for both calibration sites (Sabine and Fosheim) with the area under the receiver operating curves (AUROC) N 0.79. Both models were applied directly to Cape Bounty without calibration and validated equally with AUROC's of 0.76; however, each model predicted disturbed and undisturbed samples differently. Addition- ally, the sensitivity of the transferred model was assessed using data sets with different sample sizes. Results in- dicated that models based on larger sample sizes transferred more consistently and captured the variability within the terrain attributes in the respective study areas. Terrain attributes associated with the initiation of dis- turbances were similar regardless of the location. Disturbances commonly occurred on slopes between 4 and 15°, below Holocene marine limit, and in areas with low potential incoming solar radiation

3D models created from ROV videos and their corresponding subtransect dimensions

Underwater video transects have become a common tool for quantitative analysis of the seafloor. However a major difficulty remains in the accurate determination of the area surveyed as underwater navigation can be unreliable and image scaling does not always compensate for distortions due to perspective and topography. Depending on the camera set-up and available instruments, different methods of surface measurement are applied, which make it difficult to compare data obtained by different vehicles. 3-D modelling of the seafloor based on 2-D video data and a reference scale can be used to compute subtransect dimensions. Focussing on the length of the subtransect, the data obtained from 3-D models created with the software PhotoModeler Scanner are compared with those determined from underwater acoustic positioning (ultra short baseline, USBL) and bottom tracking (Doppler velocity log, DVL). 3-D model building and scaling was successfully conducted on all three tested set-ups and the distortion of the reference scales due to substrate roughness was identified as the main source of imprecision. Acoustic positioning was generally inaccurate and bottom tracking unreliable on rough terrain. Subtransect lengths assessed with PhotoModeler were on average 20% longer than those derived from acoustic positioning due to the higher spatial resolution and the inclusion of slope. On a high relief wall bottom tracking and 3-D modelling yielded similar results. At present, 3-D modelling is the most powerful, albeit the most time-consuming, method for accurate determination of video subtransect dimensions.

Age models, iron intensity, magnetic susceptibility records and dry bulk density of sediment cores from around the Canary Islands

A set of 43 sediment cores from around the Canary Islands is used to characterise this region, which intersects meridional climatic regimes and zonal productivity gradients in a high spatial resolution. Using rapid and nondestructive core logging techniques we carried out Fe intensity and magnetic susceptibility (MS) measurements and created a stack on the basis of five stratigraphic reference cores, for which a stratigraphic age model was available from d18O and 14C analyses on planktonic foraminifera. By correlation of the stack with the Fe and MS records of the other cores, we were able to develop age depth models at all investigated sites of the region. We present the bulk sediment accumulation rates (AR) of the Canary Islands region as an indicator of shifts in the upwelling-influenced areas for the Holocene (0-12 ky), the deglaciation (12-18 ky) and the last glacial (18-40 ky). General observations are an enhanced productivity during glacial times with highest values during the deglaciation. The main differences between the analysed time intervals we interpret as result of the sea-level effects, changes in the extent of high productivity areas, and current intensity.

Modelled spatial and seasonal distribution of Blue Whiting (Micromesistius poutassou) larvae in the North-East Atlantic (1951 to 2005)

Blue whiting (Micromesistius poutassou, http://www.marinespecies.org/aphia.php?p=taxdetails&id=126439) is a small mesopelagic planktivorous gadoid found throughout the North-East Atlantic. This data contains the results of a model-based analysis of larvae captured by the Continuous Plankton Recorder (CPR) during the period 1951-2005. The observations are analysed using Generalised Additive Models (GAMs) of the the spatial, seasonal and interannual variation in the occurrence of larvae. The best fitting model is chosen using the Aikaike Information Criteria (AIC). The probability of occurrence in the continous plankton recorder is then normalised and converted to a probability distribution function in space (UTM projection Zone 28) and season (day of year). The best fitting model splits the distribution into two separate spawning grounds north and south of a dividing line at 53 N. The probability distribution is therefore normalised in these two regions (ie the space-time integral over each of the two regions is 1). The modelled outputs are on a UTM Zone 28 grid: however, for convenience, the latitude ("lat") and longitude ("lon") of each of these grid points are also included as a variable in the NetCDF file. The assignment of each grid point to either the Northern or Southern component (defined here as north/south of 53 N), is also included as a further variable ("component"). Finally, the day of year ("doy") is stored as the number of days elapsed from and included January 1 (ie doy=1 on January 1) - the year is thereafter divided into 180 grid points.

Determining Spatial Distribution and Physical Properties of the Vashon Advance Outwash near Mountlake Terrace, Washington

This study is aimed at determining the spatial distribution, physical properties, and groundwater conditions of the Vashon advance outwash (Qva) in the Mountlake Terrace, WA area. The Qva is correlative with the Esperance Sand, as defined at its type section; however, local variations in the Qva are not well-characterized (Mullineaux, 1965). While the Qva is a dense glacial unit with low compressibility and high frictional shear strength (Gurtowski and Boirum, 1989), the strength of this unit can be reduced when it becomes saturated (Tubbs, 1974). This can lead to caving or flowing in excavations, and on a larger scale, can lead to slope failures and mass-wasting when intersected by steep slopes. By studying the Qva, we can better predict how it will behave under certain conditions, which will be beneficial to geologists, hydrogeologists, engineers, and environmental scientists during site assessments and early phases of project planning. In this study, I use data from 27 geotechnical borings from previous field investigations and C-Tech Corporation’s EnterVol software to create three-dimensional models of the subsurface geology in the study area. These models made it possible to visualize the spatial distribution of the Qva in relation to other geologic units. I also conducted a comparative study between data from the borings and generalized published data on the spatial distribution, relative density, soil classification, grain-size distribution, moisture content, groundwater conditions, and aquifer properties of the Qva. I found that the elevation of the top of the Qva ranges from 247 to 477 ft. I found that the Qva is thickest where the modern topography is high, and is thinnest where the topography is low. The thickness of the Qva ranges from absent to 242 ft. Along the northern, east-west trending transect, the Qva thins to the east as it rises above a ridge composed of Pre- Vashon glacial deposits. Along the southern, east-west trending transect, the Qva pinches out against a ridge composed of pre-Vashon interglacial deposits. Two plausible explanations for this ridge are paleotopography and active faulting associated with the Southern Whidbey Fault Zone. Further investigations should be done using geophysical methods and the modeling methods described in this study to determine the nature of this ridge. The relative density of the Qva in the study area ranges from loose to very dense, with the loose end of the spectrum probably relating to heave in saturated sands. I found subtle correlations between density and depth. Volumetric analysis of the soil groups listed in the boring logs indicate that the Qva in the study area is composed of approximately 9.5% gravel, 89.3% sand, and 1.2% silt and clay. The natural moisture content ranges from 3.0 to 35.4% in select samples from the Qva. The moisture content appears to increase with depth and fines content. The water table in the study area ranges in elevation from 231.9 to 458 ft, based on observations and measurements recorded in the boring logs. The results from rising-head and falling-head slug tests done at a single well in the study area indicate that the geometric mean of hydraulic conductivity is 15.93 ft/d (5.62 x 10-03 cm/s), the storativity is 3.28x10-03, and the estimated transmissivity is 738.58 ft2/d in the vicinity of this observation well. At this location, there was 1.73 ft of seasonal variation in groundwater elevation between August 2014 and March 2015.

Reducing errors in simulated satellite views of clouds from large-scale models

Thesis (Ph.D.)--University of Washington, 2016-06

Space-Time Smoothing Models for Surveillance and Complex Survey Data

Thesis (Ph.D.)--University of Washington, 2016-06

Remotely sensed data for ecosystem analyses: Combining hierarchy and scene models

Remotely sensed data have been used extensively for environmental monitoring and modeling at a number of spatial scales; however, a limited range of satellite imaging systems often. constrained the scales of these analyses. A wider variety of data sets is now available, allowing image data to be selected to match the scale of environmental structure(s) or process(es) being examined. A framework is presented for use by environmental scientists and managers, enabling their spatial data collection needs to be linked to a suitable form of remotely sensed data. A six-step approach is used, combining image spatial analysis and scaling tools, within the context of hierarchy theory. The main steps involved are: (1) identification of information requirements for the monitoring or management problem; (2) development of ideal image dimensions (scene model), (3) exploratory analysis of existing remotely sensed data using scaling techniques, (4) selection and evaluation of suitable remotely sensed data based on the scene model, (5) selection of suitable spatial analytic techniques to meet information requirements, and (6) cost-benefit analysis. Results from a case study show that the framework provided an objective mechanism to identify relevant aspects of the monitoring problem and environmental characteristics for selecting remotely sensed data and analysis techniques.

Spatial and temporal variability in somatic growth of green sea turtles (Chelonia mydas) resident in the Hawaiian Archipelago

The somatic growth dynamics of green turtles ( Chelonia mydas) resident in five separate foraging grounds within the Hawaiian Archipelago were assessed using a robust non-parametric regression modelling approach. The foraging grounds range from coral reef habitats at the north-western end of the archipelago, to coastal habitats around the main islands at the southeastern end of the archipelago. Pelagic juveniles recruit to these neritic foraging grounds from ca. 35 cm SCL or 5 kg ( similar to 6 years of age), but grow at foraging-ground-specific rates, which results in quite different size- and age-specific growth rate functions. Growth rates were estimated for the five populations as change in straight carapace length ( cm SCL year) 1) and, for two of the populations, also as change in body mass ( kg year) 1). Expected growth rates varied from ca. 0 - 2.5 cm SCL year) 1, depending on the foraging-ground population, which is indicative of slow growth and decades to sexual maturity, since expected size of first-time nesters is greater than or equal to 80 cm SCL. The expected size- specific growth rate functions for four populations sampled in the southeastern archipelago displayed a non-monotonic function, with an immature growth spurt at ca. 50 - 53 cm SCL ( similar to 18 - 23 kg) or ca. 13 - 19 years of age. The growth spurt for the Midway atoll population in the northwestern archipelago occurs at a much larger size ( ca. 65 cm SCL or 36 kg), because of slower immature growth rates that might be due to a limited food stock and cooler sea surface temperature. Expected age-at-maturity was estimated to be ca. 35 - 40 years for the four populations sampled at the south-eastern end of the archipelago, but it might well be > 50 years for the Midway population. The Hawaiian stock comprises mainly the same mtDNA haplotype, with no differences in mtDNA stock composition between foraging-ground populations, so that the geographic variability in somatic growth rates within the archipelago is more likely due to local environmental factors rather than genetic factors. Significant temporal variability was also evident, with expected growth rates declining over the last 10 - 20 years, while green turtle abundance within the archipelago has increased significantly since the mid-1970s. This inverse relationship between somatic growth rates and population abundance suggests a density-dependent effect on somatic growth dynamics that has also been reported recently for a Caribbean green turtle stock. The Hawaiian green turtle stock is characterised by slow growth rates displaying significant spatial and temporal variation and an immature growth spurt. This is consistent with similar findings for a Great Barrier Reef green turtle stock that also comprises many foraging-ground populations spanning a wide geographic range.