A three dimensional numerical investigation of hillslope flow processes

Physically based distributed models of catchment hydrology are likely to be made available as engineering tools in the near future. Although these models are based on theoretically acceptable equations of continuity, there are still limitations in the present modelling strategy. Of interest to this thesis are the current modelling assumptions made concerning the effects of soil spatial variability, including formations producing distinct zones of preferential flow. The thesis contains a review of current physically based modelling strategies and a field based assessment of soil spatial variability. In order to investigate the effects of soil nonuniformity a fully three dimensional model of variability saturated flow in porous media is developed. The model is based on a Galerkin finite element approximation to Richards equation. Accessibility to a vector processor permits numerical solutions on grids containing several thousand node points. The model is applied to a single hillslope segment under various degrees of soil spatial variability. Such variability is introduced by generating random fields of saturated hydraulic conductivity using the turning bands method. Similar experiments are performed under conditions of preferred soil moisture movement. The results show that the influence of soil variability on subsurface flow may be less significant than suggested in the literature, due to the integrating effects of three dimensional flow. Under conditions of widespread infiltration excess runoff, the results indicate a greater significance of soil nonuniformity. The recognition of zones of preferential flow is also shown to be an important factor in accurate rainfall-runoff modelling. Using the results of various fields of soil variability, experiments are carried out to assess the validity of the commonly used concept of `effective parameters'. The results of these experiments suggest that such a concept may be valid in modelling subsurface flow. However, the effective parameter is observed to be event dependent when the dominating mechanism is infiltration excess runoff.

Quantification of visual field loss in age-related macular degeneration

Presentation Purpose:To determine methods of quantifying the sensitivity loss in the central 10o visual field in a cross section of patients at various stages of age-related macular degeneration (AMD). Methods:Standard and short-wavelength automated perimetry (SAP and SWAP) visual fields were collected using program 10-2 of the Humphrey Field Analyzer, in 44 eyes of 27 patients with AMD and 41 eyes of 22 normal subjects. Stereoscopic fundus photographs were graded by two independent observers and the stage of disease determined. Global indices were compared for their ability to delineate the normal visual field from early stages of AMD and to differentiate between stages. Results:Mean Deviation (MD) and Pattern Standard Deviation (PSD) varied significantly with stage of disease in SAP (both p<0.001) and SWAP (both p<0.001), but post-hoc analysis revealed overlap of functional values between stages. Global indices of focal loss, PSD and local spatial variability (LSV) were the most sensitive to detecting differences between normal subjects and early stage AMD patients, in SAP and SWAP, respectively. Overall, defects were confined to the central 5°. SWAP defects were consistently greater in depth and area than those in SAP. The most vulnerable region of the 10° field to sensitivity loss with increasing stage of AMD was the central 1°, in which the sensitivity decline was -4.8dB per stage in SAP and -4.9dB per stage in SWAP. Based on the pattern deviation defect maps, a severity index of AMD visual field loss was derived. Threshold variability was considerably increased in late stage AMD eyes. Conclusions:Global indices of focal loss were more sensitive to the detection of early stage AMD from normal. The sensitivity decline with advancing stage of AMD was greater in SWAP compared to SAP, however the trend was not strong across all stages of disease. The less commonly used index LSV represents relatively statistically unmanipulated summary measure of focal loss. A new severity index is described which is sensitive to visual field change in AMD, measures visual field defects on a continuous scale and may serve as a useful measure of functional change in AMD in longitudinal studies. Keywords: visual fields • age-related macular degeneration • perimetry

Airborne laser quantification of Florida shoreline and beach volume change caused by hurricanes

This dissertation combines three separate studies that measure coastal change using airborne laser data. The initial study develops a method for measuring subaerial and subaqueous volume change incrementally alongshore, and compares those measurements to shoreline change in order to quantify their relationship in Palm Beach County, Florida. A poor correlation (R2 = 0.39) was found between shoreline and volume change before the hurricane season in the northern section of Palm Beach County because of beach nourishment and inlet dynamics. However, a relatively high R2 value of 0.78 in the southern section of Palm Beach County was found due to little disturbance from tidal inlets and coastal engineering projects. The shoreline and volume change caused by the 2004 hurricane season was poorly correlated with R 2 values of 0.02 and 0.42 for the north and south sections, respectively. The second study uses airborne laser data to investigate if there is a significant relationship between shoreline migration before and after Hurricane Ivan near Panama City, Florida. In addition, the relationship between shoreline change and subaerial volume was quantified and a new method for quantifying subaqueous sediment change was developed. No significant spatial relationship was found between shoreline migration before and after the hurricane. Utilization of a single coefficient to represent all relationships between shoreline and subaerial volume change was found to be problematic due to the spatial variability in the linear relationship. Differences in bathymetric data show only a small portion of sediment was transported beyond the active zone and most sediment remained within the active zone despite the occurrence of a hurricane. The third study uses airborne laser bathymetry to measure the offshore limit of change, and compares that location with calculated depth of closures and subaqueous geomorphology. There appears to be strong geologic control of the depth of closure in Broward and Miami-Dade Counties. North of Hillsboro Inlet, hydrodynamics control the geomorphology which in turn indicates the location of the depth of closure.

Phosphorus Availability and Salinity Control Productivity and Demography of the Seagrass Thalassia testudinum in Florida Bay

Biomass, net primary productivity (NPP), foliar elemental content, and demography of Thalassia testudinum were monitored in populations from five sites across Florida Bay beginning in January 2001. Sites were selected to take advantage of the spatial variability in phosphorus (P) availability and salinity climates across the bay. Aboveground biomass and NPP of T. testudinum were determined five to six times annually. Short-shoot demography, belowground biomass, and belowground NPP were assessed from a single destructive harvest at each site and short-shoot cohorts were estimated from leaf scar counts multiplied by site-specific leaf production rates. Biomass, relative growth rate (RGR), and overall NPP were positively correlated with P availability. Additionally, a positive correlation between P availability and the ratio of photosynthetic to non-photosynthetic biomass suggests that T. testudinum increases allocation to aboveground biomass as P availability increases. Population turnover increased with P availability, evident in positive correlations of recruitment and mortality rates with P availability. Departures from seasonally modeled estimates of RGR were found to be influenced by salinity, which depressed RGR when below 20 psu or above 40 psu. Freshwater management in the headwaters of Florida Bay will alter salinity and nutrient climates. It is becoming clear that such changes will affect T. testudinum, with likely feedbacks on ecosystem structure, function, and habitat quality.

Across-scale patterning of plant-soil-water interactions surrounding tree islands in Southern Everglades landscapes

The freshwater Everglades is a complex system containing thousands of tree islands embedded within a marsh-grassland matrix. The tree island-marsh mosaic is shaped and maintained by hydrologic, edaphic and biological mechanisms that interact across multiple scales. Preserving tree islands requires a more integrated understanding of how scale-dependent phenomena interact in the larger freshwater system. The hierarchical patch dynamics paradigm provides a conceptual framework for exploring multi-scale interactions within complex systems. We used a three-tiered approach to examine the spatial variability and patterning of nutrients in relation to site parameters within and between two hydrologically defined Everglades landscapes: the freshwater Marl Prairie and the Ridge and Slough. Results were scale-dependent and complexly interrelated. Total carbon and nitrogen patterning were correlated with organic matter accumulation, driven by hydrologic conditions at the system scale. Total and bioavailable phosphorus were most strongly related to woody plant patterning within landscapes, and were found to be 3 to 11 times more concentrated in tree island soils compared to surrounding marshes. Below canopy resource islands in the slough were elongated in a downstream direction, indicating soil resource directional drift. Combined multi-scale results suggest that hydrology plays a significant role in landscape patterning and also the development and maintenance of tree islands. Once developed, tree islands appear to exert influence over the spatial distribution of nutrients, which can reciprocally affect other ecological processes.

Does landscape context affect habitat value? The importance of seascape ecology in back-reef systems

Seascape ecology provides a useful framework from which to understand the processes governing spatial variability in ecological patterns. Seascape context, or the composition and pattern of habitat surrounding a focal patch, has the potential to impact resource availability, predator-prey interactions, and connectivity with other habitats. For my dissertation research, I combined a variety of approaches to examine how habitat quality for fishes is influenced by a diverse range of seascape factors in sub-tropical, back-reef ecosystems. In the first part of my dissertation, I examined how seascape context can affect reef fish communities on an experimental array of artificial reefs created in various seascape contexts in Abaco, Bahamas. I found that the amount of seagrass at large spatial scales was an important predictor of community assembly on these reefs. Additionally, seascape context had differing effects on various aspects of habitat quality for the most common reef species, White grunt Haemulon plumierii. The amount of seagrass at large spatial scales had positive effects on fish abundance and secondary production, but not on metrics of condition and growth. The second part of my dissertation focused on how foraging conditions for fish varied across a linear seascape gradient in the Loxahatchee River estuary in Florida, USA. Gray snapper, Lutjanus griseus, traded food quality for quantity along this estuarine gradient, maintaining similar growth rates and condition among sites. Additional work focused on identifying major energy flow pathways to two consumers in oyster-reef food webs in the Loxahatchee. Algal and microphytobenthos resource pools supported most of the production to these consumers, and body size for one of the consumers mediated food web linkages with surrounding mangrove habitats. All of these studies examined a different facet of the importance of seascape context in governing ecological processes occurring in focal habitats and underscore the role of connectivity among habitats in back-reef systems. The results suggest that management approaches consider the surrounding seascape when prioritizing areas for conservation or attempting to understand the impacts of seascape change on focal habitat patches. For this reason, spatially-based management approaches are recommended to most effectively manage back-reef systems.

Airborne Laser Quantification of Florida Shoreline and Beach Volume Change Caused by Hurricanes

This dissertation combines three separate studies that measure coastal change using airborne laser data. The initial study develops a method for measuring subaerial and subaqueous volume change incrementally alongshore, and compares those measurements to shoreline change in order to quantify their relationship in Palm Beach County, Florida. A poor correlation (R2 = 0.39) was found between shoreline and volume change before the hurricane season in the northern section of Palm Beach County because of beach nourishment and inlet dynamics. However, a relatively high R2 value of 0.78 in the southern section of Palm Beach County was found due to little disturbance from tidal inlets and coastal engineering projects. The shoreline and volume change caused by the 2004 hurricane season was poorly correlated with R2 values of 0.02 and 0.42 for the north and south sections, respectively. The second study uses airborne laser data to investigate if there is a significant relationship between shoreline migration before and after Hurricane Ivan near Panama City, Florida. In addition, the relationship between shoreline change and subaerial volume was quantified and a new method for quantifying subaqueous sediment change was developed. No significant spatial relationship was found between shoreline migration before and after the hurricane. Utilization of a single coefficient to represent all relationships between shoreline and subaerial volume change was found to be problematic due to the spatial variability in the linear relationship. Differences in bathymetric data show only a small portion of sediment was transported beyond the active zone and most sediment remained within the active zone despite the occurrence of a hurricane. The third study uses airborne laser bathymetry to measure the offshore limit of change, and compares that location with calculated depth of closures and subaqueous geomorphology. There appears to be strong geologic control of the depth of closure in Broward and Miami-Dade Counties. North of Hillsboro Inlet, hydrodynamics control the geomorphology which in turn indicates the location of the depth of closure.

Impacts of the Madden-Julian oscillation on intraseasonal precipitation over northeast Brazil

Este estudo tem como objetivo investigar os impactos da oscilação de Madden-Julian (OMJ) na precipitação da região Nordeste do Brasil (NEB). Para tanto foram utilizados dados diários de precipitação baseados em 492 pluviômetros distribuídos na região e cobrindo um período de 30 anos (1981 − 2010). As análises através de composições de anomalias de precipitação, radiação de onda longa e fluxo de umidade, foram obtidas com base no índice da OMJ desenvolvido por Jones-Carvalho. Para distinguir o sinal da OMJ de outros padrões de variabilidade climática, todos os dados diários foram filtrados na escala de 20 − 90 dias; portanto somente dias classificados como eventos da OMJ foram considerados nas composições. Uma análise preliminar baseada apenas nos dados de precipitação foi feita para uma pequena área localizada no interior semiárido do NEB, conhecida como Seridó. Essa microrregião é uma das áreas mais secas do NEB e foi reconhecida pela Convenção das Nações Unidas para o Combate à Desertificação e Mitigação dos Efeitos das Secas como particularmente vulnerável à desertificação. Composições de anomalias de precipitação foram feitas para cada uma das oito fases da OMJ durante Fevereiro-Maio (principal período chuvoso da microrregião). Os resultados mostraram a existência de variações significativas nos padrões de precipitação (de precipitação excessiva à deficiente) associados à propagação da OMJ. A combinação dos sinais de precipitação obtidos durantes as fases úmidas e secas da OMJ mostrou que a diferença corresponde cerca de 50 − 150% de modulação das chuvas na microrregião. Em seguida, uma investigação abrangente sobre o papel da OMJ sobre toda a região Nordeste foi feita considerando-se as quatro estações do ano. Os resultados mostraram que os impactos da OMJ na precipitação intrassazonal do NEB apresentam forte sazonalidade. A maior coerência espacial dos sinais de precipitação ocorreram durante o verão austral, quando cerca de 80% das estações pluviométricas apresentaram anomalias positivas de precipitação durante as fases 1 − 2 da OMJ e anomalias negativas de precipitação nas fases 5 − 6 da oscilação. Embora impactos da OMJ na precipitação intrassazonal tenham sido encontrados na maioria das localidades e em todas as estações do ano, eles apresentaram variações na magnitude dos sinais e dependem da fase da oscilação. As anomalias de precipitação do NEB observadas são explicadas através da interação existente entre as ondas de Kelvin-Rossby acopladas convectivamente e as características climáticas predominantes sobre a região em cada estação do ano. O aumento de precipitação observado sobre a maior parte do NEB durante o verão e primavera austrais encontra-se associado com o fluxo de umidade de oeste (regime de oeste), o qual favorece a atividade convectiva em amplas áreas da América do Sul tropical. Por outro lado, as anomalias de precipitação durante o inverno e outono austrais apresentaram uma variabilidade espacial mais complexa. Durante estas estações, as anomalias de precipitação observadas nas estações localizadas na costa leste do NEB dependem da intensidade do anticiclone do Atlântico Sul, o qual é modulado em grande parte por ondas de Rossby. As características topográficas do NEB parecem desempenhar um papel importante na variabilidade observada na precipitação, principalmente nestas áreas costeiras. A intensificação do anticiclone aumenta a convergência dos ventos alísios na costa contribuindo para a ocorrência de precipitação observada à barlavento do planalto da Borborema. Por outro lado, o aumento da subsidência parece ser responsável pelos déficits de precipitação observados à sotavento. Tais condições mostraram-se típicas durante o predomínio do regime de leste sobre a região tropical da América do Sul e o NEB, durante o qual ocorre uma diminuição no fluxo de umidade proveniente da Amazônia.

Age determination and sedimentology on profile Lz1020 from Lake Hoare, Antarctica

Up to 2.3 m long sediment sequences were recovered from the deepest part of Lake Hoare in Taylor Valley, southern Victoria Land, Antarctica. Sedimentological, biogeochemical, and mineralogical analyses revealed a high spatial variability of these parameters in Lake Hoare. Five distinct lithological units were recognized. Radiocarbon dating of bulk organic carbon samples from the sediment sequences yielded apparently too old ages and significant age reversals, which prevented the establishment of reliable age-depth models. However, cross correlation of the sedimentary characteristics with those of sediment records from neighbouring Lake Fryxell indicates that the lowermost two units of the Lake Hoare sediment sequences were probably deposited during the final phase of proglacial Lake Washburn, which occupied Taylor Valley during the late Pleistocene and early Holocene. High amounts of angular gravel and the absence of fine-grained material imply a complete desiccation with subaerial conditions in the Lake Hoare basin in the middle of the Holocene. The late Holocene (< c. 3300 calendar yr BP) is characterized by the establishment of environmental conditions similar to those existing today. A late Holocene desiccation event, such as proposed in former studies, is not indicated in the sediment sequences recovered.

Stable-isotope ratios and accumulation rates of four firn cores from the Dronning Maud Land, Antarctica

Four firn cores were retrieved in 2007 at two ridges in the area of the Ekström Ice Shelf, Dronning Maud Land, coastal East Antarctica, in order to investigate the recent regional climate variability and the potential for future extraction of an intermediate-depth core. Stable water-isotope analysis, tritium content and electrical conductivity were used to date the cores. For the period 1981-2006 a strong and significant correlation between the stable-isotope composition of firn cores in the hinterland and mean monthly air temperatures at Neumayer station was (r=0.54-0.71). No atmospheric warming or cooling trend is inferred from our stable-isotope data for the period 1962-2006. The stable-isotope record of the ice/firn cores could expand well beyond the meteorological record of the region. No significant temporal variation of accumulation rates was detected. However, decreasing accumulation rates were found from coast to hinterland, as well as from east (Halvfarryggen) to west (Søråsen). The deuterium excess (d) exhibits similar differences (higher d at Søråsen, lower d at Halvfarryggen), with a weak negative temporal trend on Halvfarryggen (0.04 per mil/a), probably implying increasing oceanic input. We conclude that Halvfarryggen acts as a natural barrier for moisture-carrying air masses circulating in the region from east to west.

Estimates of annual mean surface temperatures in the late Quaternary of the tropical Atlantic

At least two modes of glacial-interglacial climate change have existed within the tropical Atlantic Ocean during the last 20,000 years. The first mode (defined by cold glacial and warm interglacial conditions) occurred symmetrically north and south of the equator and dominated the eastern boundary currents and tropical upwelling areas. This pattern suggests that mode 1 is driven by a glacial modification of surface winds in both hemispheres. The second mode of oceanic climate change, defined by temperature extremes centered on the deglaciation, was hemispherically asymmetrical, with the northern tropical Atlantic relatively cold and the southern tropical Atlantic relatively warm during deglaciation. A likely cause for this pattern of variation is a reduction of the presently northward cross-equatorial heat flux during deglaciation. No single mechanism accounts for all the data. Potential contributors to oceanic climate changes are linkage to high-latitude climates, modification of monsoonal winds by ice sheet and/or insolation changes, atmospheric CO2 and greenhouse effects, indirect effects of glacial meltwater, and variations in thermohaline overturn of the oceans.

Sedimentology on surface sediments from Lake Ohrid

Lake Ohrid is likely of Pliocene age and thus commonly referred to as the oldest existing lake in Europe. In this study spatial variability of recent sediment composition is assessed using >50 basin wide distributed surface sediment samples. Analysis of biogeochemical bulk parameters, selected metals, pigment concentrations as well as grain size distributions revealed a significant spatial heterogeneity in surface sediment composition. It implies that sedimentation in Lake Ohrid is controlled by an interaction of multiple natural and anthropogenic factors and processes. Major factors controlling surface sediment composition are related to differences in geological catchment characteristics, anthropogenic land use, and a counterclockwise rotating surface water current. In some instances processes controlling sediment composition also seem to impact distribution patterns of biodiversity, which suggests a common interaction of processes responsible for both patterns.

Evaluation of the ECOSSE model for simulating soil organic carbon under Miscanthus and short rotation coppice-willow crops in Britain

Acknowledgements This work contributes to the ELUM (Ecosystem Land Use Modelling & Soil Carbon GHG Flux Trial) project, which was commissioned and funded by the Energy Technologies Institute (ETI). We acknowledge the E-OBS data set from the EU-FP6 project ENSEMBLES (http://ensembles-eu.metoffice.com) and the data providers in the ECA&D project (http://www.ecad.eu).

Biological and Physical Factors Affecting the Natural History and Evolution of Encapsulated Development

The evolution of reproductive strategies involves a complex calculus of costs and benefits to both parents and offspring. Many marine animals produce embryos packaged in tough egg capsules or gelatinous egg masses attached to benthic surfaces. While these egg structures can protect against environmental stresses, the packaging is energetically costly for parents to produce. In this series of studies, I examined a variety of ecological factors affecting the evolution of benthic development as a life history strategy. I used marine gastropods as my model system because they are incredibly diverse and abundant worldwide, and they exhibit a variety of reproductive and developmental strategies.

The first study examines predation on benthic egg masses. I investigated: 1) behavioral mechanisms of predation when embryos are targeted (rather than the whole egg mass); 2) the specific role of gelatinous matrix in predation. I hypothesized that gelatinous matrix does not facilitate predation. One study system was the sea slug Olea hansineensis, an obligate egg mass predator, feeding on the sea slug Haminoea vesicula. Olea fed intensely and efficiently on individual Haminoea embryos inside egg masses but showed no response to live embryos removed from gel, suggesting that gelatinous matrix enables predation. This may be due to mechanical support of the feeding predator by the matrix. However, Haminoea egg masses outnumber Olea by two orders of magnitude in the field, and each egg mass can contain many tens of thousands of embryos, so predation pressure on individuals is likely not strong. The second system involved the snail Nassarius vibex, a non-obligate egg mass predator, feeding on the polychaete worm Clymenella mucosa. Gel neither inhibits nor promotes embryo predation for Nassarius, but because it cannot target individual embryos inside an egg mass, its feeding is slow and inefficient, and feeding rates in the field are quite low. However, snails that compete with Nassarius for scavenged food have not been seen to eat egg masses in the field, leaving Nassarius free to exploit the resource. Overall, egg mass predation in these two systems likely benefits the predators much more than it negatively affects the prey. Thus, selection for environmentally protective aspects of egg mass production may be much stronger than selection for defense against predation.

In the second study, I examined desiccation resistance in intertidal egg masses made by Haminoea vesicula, which preferentially attaches its flat, ribbon-shaped egg masses to submerged substrata. Egg masses occasionally detach and become stranded on exposed sand at low tide. Unlike adults, the encased embryos cannot avoid desiccation by selectively moving about the habitat, and the egg mass shape has high surface-area-to-volume ratio that should make it prone to drying out. Thus, I hypothesized that the embryos would not survive stranding. I tested this by deploying individual egg masses of two age classes on exposed sand bars for the duration of low tide. After rehydration, embryos midway through development showed higher rates of survival than newly-laid embryos, though for both stages survival rates over 25% were frequently observed. Laboratory desiccation trials showed that >75% survival is possible in an egg mass that has lost 65% of its water weight, and some survival (<25%) was observed even after 83% water weight lost. Although many surviving embryos in both experiments showed damage, these data demonstrate that egg mass stranding is not necessarily fatal to embryos. They may be able to survive a far greater range of conditions than they normally encounter, compensating for their lack of ability to move. Also, desiccation tolerance of embryos may reduce pressure on parents to find optimal laying substrata.

The third study takes a big-picture approach to investigating the evolution of different developmental strategies in cone snails, the largest genus of marine invertebrates. Cone snail species hatch out of their capsules as either swimming larvae or non-dispersing forms, and their developmental mode has direct consequences for biogeographic patterns. Variability in life history strategies among taxa may be influenced by biological, environmental, or phylogenetic factors, or a combination of these. While most prior research has examined these factors singularly, my aim was to investigate the effects of a host of intrinsic, extrinsic, and historical factors on two fundamental aspects of life history: egg size and egg number. I used phylogenetic generalized least-squares regression models to examine relationships between these two egg traits and a variety of hypothesized intrinsic and extrinsic variables. Adult shell morphology and spatial variability in productivity and salinity across a species geographic range had the strongest effects on egg diameter and number of eggs per capsule. Phylogeny had no significant influence. Developmental mode in Conus appears to be influenced mostly by species-level adaptations and niche specificity rather than phylogenetic conservatism. Patterns of egg size and egg number appear to reflect energetic tradeoffs with body size and specific morphologies as well as adaptations to variable environments. Overall, this series of studies highlights the importance of organism-scale biotic and abiotic interactions in evolutionary patterns.

Inter-site and inter-scanner diffusion MRI data harmonization.

We propose a novel method to harmonize diffusion MRI data acquired from multiple sites and scanners, which is imperative for joint analysis of the data to significantly increase sample size and statistical power of neuroimaging studies. Our method incorporates the following main novelties: i) we take into account the scanner-dependent spatial variability of the diffusion signal in different parts of the brain; ii) our method is independent of compartmental modeling of diffusion (e.g., tensor, and intra/extra cellular compartments) and the acquired signal itself is corrected for scanner related differences; and iii) inter-subject variability as measured by the coefficient of variation is maintained at each site. We represent the signal in a basis of spherical harmonics and compute several rotation invariant spherical harmonic features to estimate a region and tissue specific linear mapping between the signal from different sites (and scanners). We validate our method on diffusion data acquired from seven different sites (including two GE, three Philips, and two Siemens scanners) on a group of age-matched healthy subjects. Since the extracted rotation invariant spherical harmonic features depend on the accuracy of the brain parcellation provided by Freesurfer, we propose a feature based refinement of the original parcellation such that it better characterizes the anatomy and provides robust linear mappings to harmonize the dMRI data. We demonstrate the efficacy of our method by statistically comparing diffusion measures such as fractional anisotropy, mean diffusivity and generalized fractional anisotropy across multiple sites before and after data harmonization. We also show results using tract-based spatial statistics before and after harmonization for independent validation of the proposed methodology. Our experimental results demonstrate that, for nearly identical acquisition protocol across sites, scanner-specific differences can be accurately removed using the proposed method.