L'hétérogénéité spatiale est aujourd'hui reconnue comme un facteur primordial pour la diversité et le fonctionnement des écosystèmes prairiaux, mais reste souvent négligée dans l'analyse de l'impact de différents modes d'utilisation pastorale. Les communautés prairiales subalpines sont des milieux hautement diversifiés, au cœur d'enjeux écologiques, socio-économiques et culturels. Depuis plusieurs décennies, des changements d'usage, notamment l'extensification conduisant parfois à l'abandon du pâturage, affectent ces prairies dans de vastes régions des Alpes. Connaitre la réponse spatiale de ces communautés à ces changements d'usage peut constituer un enjeu important pour mieux prédire les conséquences en termes de valeurs écologiques et agronomiques. L'objectif général de ce travail est d'identifier l'organisation spatiale de composantes spécifiques et fonctionnelles de communautés prairiales subalpines dans des situations contrastées d'usage : pâturage traditionnel et abandon sur le moyen terme (~20 ans), de comprendre les mécanismes liés au pâturage influençant cette organisation spatiale et d'explorer les conséquences potentielles de l'organisation spatiale des composantes spécifiques et fonctionnelles sur les valeurs d'usage de ces milieux. Ce travail révèle que l'arrêt du pâturage conduit à des modifications importantes de l'organisation spatiale à la fois des espèces et de traits fonctionnels aériens. La variabilité spatiale de ces composantes est systématiquement augmentée à différentes échelles spatiales en situation d'abandon révélant un grain plus grossier d'hétérogénéité spatiale, lié à une plus faible coexistence locale des espèces et des traits et à une convergence fonctionnelle plus forte que sous l'effet diversificateur à échelle fine du pâturage traditionnel. D'autres part, même si une réponse en termes de patrons spatiaux est également observée, celle-ci n'est pas systématique, indiquant l'absence d'une action structurante forte du pâturage qui agit donc surtout sur la variabilité. Une expérimentation in situ montre en effet que la formation de patrons spatiaux marqués en réponse à une action de défoliation sélective semble limitée dans ces prairies relativement peu productives. De plus, une faible structuration spatiale du couvert en termes de valeurs d'usage fourragère a été observée, pilotée dès des échelles très fines par la variabilité spatiale à un grain fin des assemblages d'espèces et de traits dans les zones pâturées. Ces résultats suggèrent ainsi l'expression d'un pâturage peu sélectif exprimé par les troupeaux ovins dans ces milieux peu productifs, et la prépondérance du filtre abiotique dans la détermination des assemblages spatiaux des espèces et des traits. La formation de patchs très contrastés et très structurés spatialement en réponse au pâturage est plus probable dans des milieux très productifs où la compétition est plus intense et la repousse de la végétation favorisée. La prise en compte de la variabilité autour de la réponse moyenne des communautés se révèle donc un élément essentiel pour analyser l'évolution des prairies subalpines en réponse au changement d'usages pastoraux. Une des perspectives fortes suite à ce travail est l'analyse de la synergie entre la réponse de l'hétérogénéité spatiale du couvert et celle du compartiment souterrain (traits racinaires, mycorhization, diversité microbienne…). / Spatial heterogeneity is now widely recognized as a key component for diversity and functioning of grazed plant communities but is still often neglected for the analysis of the impact of different pastoral uses. Subalpine grassland communities are species-rich habitats of high conservative, economic and cultural values. In large regions of the Alps, these grasslands have experienced a decrease in pastoral use and even abandonment, in the last decades. The investigation of the spatial response of these communities to these changes in pastoral use can be a key issue to better predict ecological and agronomical consequences of such changes. The main objective of this study is to identify the spatial organization of species and traits of subalpine grassland communities in contrasted situations of pastoral use: traditional grazing and abandonment (~20 years), to understand grazing-induced mechanisms influencing this spatial organisation and to explore the potential consequences of the spatial organization of species and traits on the pastoral value of these grasslands. This work reveals that grazing cessation leads to strong modifications of the spatial organization of both species and aerial functional traits. The spatial variability of these components are systematically increased at different spatial scales following grazing abandonment which revealed a coarser grain of spatial heterogeneity in relation to a lower species and trait local coexistence and to a stronger functional convergence than under the fine-scale diversifying action of traditional grazing. Besides, despite we also observed a response in terms of spatial patterns, this one is not systematic which indicates the absence of strongly structuring action of grazing that especially acts on variability. An in situ experiment indeed showed that spatial pattern creation in response to selective defoliation is limited in these relatively low productive grassland communities. In addition, we also evidenced that pastoral values (forage quality) exhibited poor spatial heterogeneity in these grasslands in relation to the fine grained spatial variability of species and plant traits in grazed areas. These results suggest the expression of a low selective grazing promoted by sheep herds grazing and low productivity of subalpine communities and the predominance of the abiotic filter in the determination of spatial assemblages of species and traits. Pronounced spatial patterns in response to grazing may be expected in more productive environments where competitive interactions between species are more intense and species regrowth promoted. Therefore, the consideration of the spatial variability in addition to the mean response of community components appeared to be essential in the analysis of subalpine grassland community dynamics in response to land use changes. One of the strong perspective following this work is to analyse the association between spatial responses of vegetation components and those of belowground components (root traits, microbial diversity…).

It is well known that methane (CH4) is emitted from soil, water and wetlands under anaerobic conditions through methanogenesis. CH4 is the final product of the anaerobic respiration of the microorganism methanogen. More recently, it has been shown that CH4 is also emitted by trees and if only the emissions from soil and water are measured the fluxes of CH4 in the ecosystem will be underestimated. Considering the emission from trees, the Amazon region greatly contributes to global emissions. To investigate if there is need for method development for measuring CH4 fluxes, the aim in this study was to statistically test the spatial and the short temporal variability of CH4 emissions from trees. This was done within and between two different seasons in three different plots in the Amazon basin during the year 2017. Samples of CH4 were collected using semi rigid chambers placed on tree stems. The samples were later analyzed in a laboratory environment using the Los Gatos Ultraportable Greenhouse Gas Analyzer (UGGA). For the statistical analysis non-parametric test were used, due to the non-parametric data. In this study, the result shows that the short temporal variability is not statistically significant in any of the three plots, but the short temporal variability is statistically significant between the two seasons. This tells us that it is of importance to collect samples during different seasons of the year when measuring CH4 emissions from trees. The spatial variability is statistically significant on all the three plots in both seasons. This tell us that it is important to collect samples from different heights of the tree stems when collecting CH4 samples regardless of the season. / The Global Methane Budget

Dans le cadre d’un projet d’aménagement, la reconnaissance du sous-sol est une étape prépondérante pour la construction d’un bâtiment. Le géotechnicien se doit de définir le comportement mécanique des sols à partir d’un nombre limité de sondages. La problématique soulevée par les prospections géotechniques est qu’il est difficile d’une part d’étendre une information ponctuelle dans l’environnement proche des sondages, et d’autre part d’estimer la variabilité latérale des matériaux. Dans cette optique, cette thèse développe une approche novatrice de combinaison de méthodes géophysiques aidant à la quantification de la variabilité spatiale des structures géoélectriques du sous-sol. Dans un second temps, au regard des limites rencontrées dans la modélisation 3D des structures géoélectriques du sous-sol, nous avons mis au point un processus innovant de modélisation 3D des propriétés physiques du sous-sol par combinaison de méthodes géophysiques. Les travaux menés ont permis de caractériser les paramètres régissant le positionnement et la quantité de données sources nécessaires.Enfin, un dispositif innovant de mesure de Tomographie de Résistivité Electrique en Cinématique (TREC) a été développé. La validation du dispositif sur trois sites différents a permis de mettre en évidence sa capacité à détecter les structures géoélectriques des sols sur un linéaire de plusieurs centaines de mètres. Ce dispositif constitue un outil précieux pour l’étude des structures géoélectriques du sous-sol dans le cadre de prospections géophysiques à grand rendement. / The geotechnical study of a building project is a decisive step for the construction of it. The geotechnical engineer must define the mechanical behavior of soils from a limited number of drills. The problem of the geotechnical surveys is, on the one hand that difficult to extend a punctual data in the near space to the drills,and on the other hand to estimate the lateral variability of materials. In this context, this thesis develops an innovative approach by combination of geophysical methods to quantify the spatial variability of geoelectrical structures. In asecond time, against the known limits of the 3D modeling of soil geoelectrical structures, we have developed an innovative process of 3D modeling of the physical properties of soil by geophysical methods combination. These works allow characterizing the parameters who governing the positioning and the number of necessary primary data.Finally, an innovative measuring device of Kinematic Electrical Resistivity Tomography (KERT) has been developed. The validation of the device on three different site has allowed to highlight its ability to detect the soil geoelectrical structures within a linear of hundreds meters. This device is a valuable tool for studying soil geoelectrical structures in the context of high yield geophysical prospection.

<p>The transport of elements by streams from headwater regions to the sea is influenced by landscape characteristics. This thesis focuses on the influence of landscape characteristics (e.g. proportion of wetland/forest coverage) on temporal and spatial variations of Fe, Mn, S and trace elements (As, Co, Pb) in streams located in northern Sweden, a boreal region characterized by coniferous forests and peat wetlands.</p><p>Water samples from a network of 15 streams revealed a different hydrogeochemistry in forested catchments compared to wetland catchments. The temporal variation was dominated by spring flood, when concentrations of Fe, Mn and trace elements increased in forested headwaters. However, in streams of wetland catchments concentrations decreased, but Pb concentrations were higher in comparison to other streams. Both Fe and Pb showed positive correlations with wetland area, while Co correlated with forest coverage. The anthropogenic contribution of As and Pb appear to be larger than the supply from natural sources.</p><p>During spring flood SO₄^2- decreased in most streams, although concentrations increased in streams of wetland catchments. Concentrations of SO₄^2-were higher in streams of forested catchments than in wetland dominated streams, the former being net exporters of S and the latter net accumulators. Isotope values of stream water SO₄^2- (δ³⁴S_SO4) were close to that of precipitation during spring flood, indicating that the major source of S is from deposition. The results show that, although emissions of anthropogenic S have been reduced, there is still a strong influence of past and current S deposition on runoff in this region.</p><p>In conclusion, wetlands are key areas for the hydrogeochemistry in this boreal landscape. The findings emphasize the importance of understanding stream water chemistry and element cycling from a landscape perspective. This may be important for predicting how boreal regions respond to environmental disturbances such as climate change.</p>

The focus of this study is on potential damages to buildings resulting from expansive soils in Tanzania, particularly clay soils in Kibaha. For the fact that most of the affected structures are founded on expansive soils, a clear understanding of the behaviour of soils and their interaction with structures has been of interest to the study in order to evaluate properly the source of the problem.The geotechnical behaviour of expansive clay soils is investigated by looking into the geomorphologic, geological and climatic conditions and mineralogical composition of the soils in the study area.Two sites, representative of known problem-areas in Kibaha were selected for geotechnical tests. Geotechnical site investigation consisted of open trial pits, profile description and the collection of both disturbed and undisturbed samples. To extend and amplify the findings, supplementary samples were collected from the environs of the two sites.The collected samples were submitted to soil laboratories at KTH, ARU, SEAMIC and DIT for mineralogical composition tests, natural water content, density, Atterberg limits and swell tests. The results of this investigation indicate that soils in Kibaha contains clay (31%), have high liquid limit (59%) and plastic limit (37%) which indicate high potential swell.Since swell pressure, free swell and swell percent are key properties of expansive soils, the swell properties were measured by free swell tests and one-dimensional oedometer swell tests. The free swell ranged from 100% to 150% and the swell pressure was in the region of 45 kPa. The coefficient of linear extensibility (COLE) was determined for characterizing expansive clays. For all tested samples, COLE ranged from 0.09 to 0.14 indicating that soils fall in the region of high to very high expansion potential rating. The properties of expansive soils were confirmed by the x-ray diffraction test which showed the presence of smectite in the soil. Furthermore, total suction measurement technique using filter paper method indicated that the soils have high suction values, signifying that they have a tendency to swell upon wetting depending on plasticity of particular soil.The depth of the active zone was measured as a function of moisture variations in the profiles during two extreme weather conditions. The active zone depth was found to be between 1.0 and 2.0 m deep. Procedures to assess models to predict swell in the case study were outlined together with their validity.Vertical and horizontal spatial variability in selected soil properties was defined using geostatistical techniques through the fitting of variogram. The indicator semivariograms of both clay contents and free swell gave a range of 20 m horizontally and 1.0 m vertically, with the horizontal variograms exhibiting greater ranges than the dipping variograms.Physical conditions of the surveyed properties in the area confirmed that building damages are associated with poor building materials triggered by expansive soils. In support of the obtained data, the actual behaviour of the foundations was supplemented with prototypes of strip foundations whose performances were monitored over a period of four months. Finally, suggested are the ways forward to solve the problem of foundation on expansive soil / QC 20100824

Die Berücksichtigung der flächenhaften bodenhydraulischen Variabilität gilt bei der Modellierung von Wasser- und Stofftransportprozessen als problematisch. Dies liegt vorrangig an ihrer Erfassung, die kosten- und zeitintensiv ist. Die vorliegende Arbeit untersucht verschiedene Messstrategien, die zur Abbildung der flächenhaften Bodenhydraulik mit wenigen, einfach zu bestimmenden und physikalisch begründeten Bodenparametern führen. Die Vorgehensweise erfolgt mit der Anwendung eines Ähnlichkeitskonzeptes, das die Böden in bodenhydraulisch ähnliche Klassen unterteilt. Innerhalb einer Klasse kann die Variabilität der Retentions- und hydraulischen Leitfähigkeitcharakteristik auf einen freien Parameter (Skalierungsparameter) reduziert werden. Die Analyse der Zusammenhänge zwischen Boden- und Skalierungsparametern führt letztendlich zu den geeigneten Parametern die eine flächenhafte Abbildung möglich machen. Diese Untersuchungen bilden die Grundlage für die weitere Entwicklung eines stochastischen Modellansatzes, der die Variabilität der Bodenhydraulik bei der Modellierung des Bodenwassertransportes im Feldmaßstab berücksichtigen kann. An Hand von drei Datensätzen unterschiedlicher Skalenausbreitung konnte dieses Konzept angewendet werden. Die Ergebnisse zeigen, dass die Beschreibung der hydraulischen Variabilität nur für die vertikale (Profil) nicht aber für die flächenhafte Ausbreitung mit einfachen Bodenparametern möglich ist. Mit einer ersten Modellanwendung konnte gezeigt werden, dass über die Variabilität der Bodenparameter Trockenrohdichte und Tongehalt auch die Variabilität der Bodenhydraulik und damit die Berechnung des Bodenfeuchteverlaufs am Standort darstellbar ist. / The consideration of the spatial variability of the unsaturated soil hydraulic characteristics still remains an unsolved problem in the modelling of the water and matter transport in the vadose zone. This can be mainly explained by the rather cumbersome measurement of this variability, which is both, time-consuming and cost-intensive. The presented thesis analyses various measurement strategies which aim at the description of the soil-hydraulic heterogeneity by a small number of proxy-parameters, which should be easily measurable and still have a soil-physical meaning. The developed approach uses a similarity concept, which groups soils into similar soil hydraulic classes. Within a class, the variability of the retention and hydraulic conductivity curves can be explained by a single parameter (scaling parameter). The analysis of the correlation between the soil parameters and the scaling parameters can eventually indicate which soil parameters can be used for describing the soil hydraulic variability in a given area. This investigation forms the basis for the further development of a stochastic model, which can integrate the soil-hydraulic variability in the modelling of the soil water transport. Three data sets, all covering different scales, were subsequently used in the application of the developed concept. The results show that depth development of the soil-hydraulic variability in a soil profile can be explained by a single soil parameter. Contrarily, the explanation of the horizontal variability of the soil-hydraulic properties was not possible with the given data sets. First model applications for a soil profile showed that including the variability of the soil parameters bulk density and clay fraction in the water transport simulations could describe the variability of the soil-hydraulic variability and thus, the dynamics of the soil water content at the investigated profile.

This thesis aims to increase our understanding of mechanisms that influence larval dispersal in marine benthic invertebrates, particularly in the absence of strong oceanographic features (e.g. estuarine plumes, upwelling events, or markedly different water masses). Laboratory experiments identified behavioural mechanisms that regulate the vertical distribution of larvae in response to thermal stratification, and field studies in St. George’s Bay, Nova Scotia (NS), Canada, examined the relationship between larval abundance and physical variables (temperature, salinity, fluorescence, etc) and identified mechanisms that regulate larval distributions in situ. In the laboratory, I demonstrated that thermal stratification affects the vertical distribution of larvae by acting as a barrier to migration, or through temperature-dependent vertical swimming velocities. I also developed a random walk based model which highlighted that the key to successfully simulating larval response to temperature was 1) determining the temperature-dependent distribution of vertical swimming velocities and 2) the temporal autocorrelation in these velocities. In the field, the most striking pattern was that the larval distributions for species with similar swimming abilities were significantly correlated to one another at all scales (0.5 to 40 km). This suggests a common mechanism, related to larval swimming ability, which greatly influences the horizontal larval distribution. I found that the spatial scale of variability in larval distributions (~ 3 km) matches that in both the environmental variables and of coherent structures in current velocities (i.e. the tidal excursion). Results from an aggregation-diffusion model suggest that horizontal larval swimming could not be responsible for the observed level of aggregation in the larval horizontal distributions. I suggest that these horizontal patterns are the result of 1) an aggregative process (i.e. larvae swimming against a vertical current and maintaining their vertical position) and 2) a diffusive process which scales the aggregations to the scale of the coherent structures in current velocity (i.e. tidal excursion). In conclusion, this thesis increases our understanding of larval behaviour and its effects on larval dispersal. The results will be particularly useful to those who are interested in mechanisms regulate population connectivity, particularly those using bio-physical models to model dispersal trajectories.

The consequences of land use choices on soil water and gas transport properties are significant for gas and water flux in agricultural environments. Spatial and temporal patterns and associations of soil water and soil gas characteristics and processes in different land uses are not well understood. The objectives of this study were to 1) characterize soil structure under crop and grass systems, 2) quantify spatial patterns and associations of soil physical characteristics in crop and grass systems, and 3) quantify spatial and temporal patterns and associations of CO2 and N2O fluxes. The research was conducted in a 60 by 80 m field divided into grass and crop systems. Sixty sampling points were distributed in four transects with 5- and 1-m spatial intervals between measurement points. Gas fluxes were measured, at two-week time intervals, 22 times during a year. Pore size distribution was more homogeneous and more continuous pores were found in the grass than in the crop system. The spatial variability of most selected soil physical characteristics was more structured in the crop than in the grass system, which reflected the impact of land use and soil structure on their spatial patterns. CO2 flux was dependent for a longer distance in the grass than in the crop system, however, the two land-use systems exhibited similar spatial ranges of N2O flux. Gas fluxes were temporally dependent for a longer period in the grass than in the crop system. The spatial associations between CO2 and N2O fluxes and selected biochemical and physical factors depended on the flux sampling season and land use. Soil temperature was the dominant controlling factor on the temporal variability of CO2 and N2O fluxes but not on the spatial behavior. Considering the spatial and temporal ranges and dependency strength of soil variables helps identify efficient sampling designs that can result in better time and resource management. Spatial and temporal relationships between the selected soil variables also improve understanding soil management and sampling soil variables. This study provides the baseline and recommendations for future investigations specifically for sampling designs, soil management, and predictions of different soil processes related to gas fluxes.

Frequent use of phenylurea herbicide isoproturon (IPU) in agricultural fields has resulted not only in the contamination of the natural resources including soil and water but also in the adaptation of the soil microflora to its rapid degradation. However, up to now, the mechanisms underlying this microbial adaptation are not well elucidated. The aim of this study was to explore the processes and factors implicated in IPU degradation from the agricultural field to the genes coding for catabolic genes. The study carried out at the experimental field of Epoisses cropped with a winter wheat / barley / rape seed crop rotation indicated that as a result of its periodically repeated use, the soil microflora adapted to IPU mineralization activity. Further analysis using exploratory and geostatistical tools demonstrated the existence of spatial variability in IPU mineralization activity at the field scale which was correlated not only with several soil physico-chemical parameters like organic matter content, CEC and C/N ratio but also with the pesticide application plan over a three year crop rotation. In order to get further insight into underlying mechanisms, an IPU mineralizing bacterial culture and strain Sphingomonas sp. SH were isolated through enrichment cultures performed from two different adapted soils. Both had the catabolic activities highly specific for the mineralization of IPU and its metabolites but none of other structurally related phenylurea herbicides. IPU metabolic activity of both the mixed culture and the strain SH was found to be affected by pH with optimal activity taking place at pH 7.5. Based on the accumulation of different known metabolites during mineralization kinetics, IPU metabolic pathway was proposed to be initiated by two successive demethylations, followed by cleavage of the urea side chain resulting in the accumulation of 4-isopropylaniline, and ultimately the mineralization of the phenyl ring. In order to identify the genes involved in IPU degradation, BAC clone library was established from the genomic DNA of the bacterial culture. Although, the functional screening did not yield in identifying any BAC clone able to degrade IPU or its known metabolites, the PCR based screening led us to identify a cat gene cluster involved in ortho-cleavage of the phenyl ring of catechol through beta-ketoadipate pathway. Based on this finding, it was hypothesized that phenyl ring of 4-isopropylaniline formed during IPU transformation might be mineralized through ortho-cleavage of catechol. This finding allowed us to propose the lower IPU metabolic pathway which was not yet described.

Masters Research - Master of Philospohy (MPhil) / It has been shown that masonry material properties, in particular, unit flexural bond strength (ft), vary significantly throughout masonry structures, despite the fact that often only one type of brick and mortar are used. Unit flexural bond strength was previously identified as one of the most important material parameters contributing to the strength of clay brick unreinforced masonry (URM) walls in flexure. It was the objectives of this research, in the context of clay brick URM walls subject to vertical bending, to examine how unit flexural bond strength varied spatially in a clay brick URM wall, determine a best fit probability distribution function which can describe expected variability in unit flexural bond strength and determine how this variability and other factors affect wall behaviour and failure load using 3D non-linear finite element analysis (FEA). It was hoped that modelling a full sized clay brick URM wall subject to vertical bending using a 3D non-linear FEA model would more accurately predict wall failure load (compared to current analytical methods) and allow the examination of crack pattern development as the wall progresses to failure upon being laterally loaded. The first part of the research project was to conduct an experimental program to examine unit-to-unit spatial strength correlation within six full sized clay brick URM walls and to characterise a unit flexural bond strength probability distribution. It was observed that although weak correlation in unit flexural bond strength exists in some courses and between courses, these locations were difficult to predict and didn����t follow any particular pattern relating to for example, mortar batch. Therefore, although somewhat counter-intuitive, the results indicate that statistically significant correlation between adjacent unit flexural bond strengths is not likely to be observed. It was also observed that clay brick wall unit flexural bond strengths obtained for all of the walls tested best fit a truncated Normal probability distribution. Strength of the brick/mortar interface appeared to be governed by factors relating to workmanship (and therefore mortar quality and moisture content), weather (which can affect material characteristics like brick suction rate) and inherent material variability. It would appear that brick suction rate can significantly affect the overall strength of a URM wall. v Stochastic analysis was conducted for walls with and without uncorrelated spatial variability in unit flexural bond strength and associated tensile fracture energy (GfI ). It was found that the TNO DIANA 9.2 FEA package could be used to implement spatial variability of various material parameters and reasonably accurately model failure of clay brick URM walls in vertical bending. From the non-linear FEA model development stage, it was observed that because the brick/mortar bond has significantly more strength capacity in compression, it appears that the lateral load resistance of the wall comes from a combination of the ability of the brick/mortar bond to tensile soften while providing significant compressive resistance at the compressive edge. It was found for a spatial stochastic analysis with spatial variability in bond strength (referred to from now on as a spatial stochastic analysis), with COVs of 0.1, 0.3 and 0.5, that COV of wall failure loads were relatively small, being 0.02, 0.04 and 0.06 respectively. For the non-spatially varying stochastic analysis with fully correlated bond strength (now referred to as non-spatial stochastic analysis), with COVs of 0.1, 0.3 and 0.5, COV of wall failure loads were 0.07, 0.20 and 0.32 respectively. For the spatial stochastic analysis, it was found that with a bond strength COV increase from 0.1 to 0.5 the mean wall failure load dropped from 2.25 kPa to 2.0 kPa (an 11% reduction). Despite the relatively small drop in magnitude of the mean wall failure load with increase in bond strength COV, the mean wall failure loads were statistically different to one another. For the non-spatial stochastic analysis, mean failure load stayed relatively constant at 2.24-2.25 kPa. These results could be explained by examining the 3D wall progression to failure. For walls with spatial variability in bond strength, it is expected that wall failure load COVs would be smaller because those walls would consistently be composed of smaller valued bond strengths which would consistently contribute to weakness in the wall. For the non-spatial wall simulations, this effect would not occur as failure load is determined by one uniform weak or strong bond strength. It was proposed that failure of a clay brick URM wall is not governed by one course only cracking, but rather, instability in the wall is governed by several courses in the vicinity of locations of large bending moment. It was shown that various current stochastic approximations which employ a unit failure hypotheses in combination with a linear/elastic approximation for first cracking load all underestimated wall capacity significantly. The reason for this is suggested as being vi because all hypotheses only assume failure is governed by one course and linear/elastic theory only considers the tensile capacity of a joint and neglects strength capacity available as a result of joint tension softening and the resistance to failure provided by compressive strength on the compression side of the wall. The hypotheses also don’t take into consideration factors which affect overall wall bond strength mean which result from influences such as workmanship, weather and material variability factors, such as (for example), variation in brick suction rate due to weather conditions which can make the overall strength of the wall stronger or weaker. Based upon a comparison in wall failure load COV for the spatial and non-spatial stochastic wall analysis results, a more realistic approach for future modelling attempts of spatial variability in masonry material properties is suggested. This would address the issue of external factors such as workmanship and weather on the overall strength of the wall, as well as the inherent bond strength variability due to material variability. For walls with spatial variability in bond strength, upon examination of numerous wall simulation results, several crack patterns were witnessed and are discussed.