A peak-clustering method for MEG group analysis to minimise artefacts due to smoothness

Magnetoencephalography (MEG), a non-invasive technique for characterizing brain electrical activity, is gaining popularity as a tool for assessing group-level differences between experimental conditions. One method for assessing task-condition effects involves beamforming, where a weighted sum of field measurements is used to tune activity on a voxel-by-voxel basis. However, this method has been shown to produce inhomogeneous smoothness differences as a function of signal-to-noise across a volumetric image, which can then produce false positives at the group level. Here we describe a novel method for group-level analysis with MEG beamformer images that utilizes the peak locations within each participant's volumetric image to assess group-level effects. We compared our peak-clustering algorithm with SnPM using simulated data. We found that our method was immune to artefactual group effects that can arise as a result of inhomogeneous smoothness differences across a volumetric image. We also used our peak-clustering algorithm on experimental data and found that regions were identified that corresponded with task-related regions identified in the literature. These findings suggest that our technique is a robust method for group-level analysis with MEG beamformer images.

Mode-locked semiconductor lasers with optical injection

We perform characterization of the pulse shape and noise properties of quantum dot passively mode-locked lasers (PMLLs). We propose a novel method to determine the RF linewidth and timing jitter, applicable to high repetition rate PMLLs, through the dependence of modal linewidth on the mode number. Complex electric field measurements show asymmetric pulses with parabolic phase close to threshold, with the appearance of waveform instabilities at higher currents. We demonstrate that the waveform instabilities can be overcome through optical injection-locking to the continues wave (CW) master laser, leading to time-bandwidth product (TBP) improvement, spectral narrowing, and spectral tunability. We discuss the benefits of single- and dual-tone master sources and demonstrate that dual-tone optical injection can additionally improve the noise properties of the slave laser with RF linewidth reduction below instrument limits (1 kHz) and integrated timing jitter values below 300 fs. Dual-tone injection allowed slave laser repetition rate control over a 25 MHz range with reduction of all modal optical linewidths to the master source linewidth, demonstrating phase-locking of all slave modes and coherence improvement.

Mode-locked semiconductor lasers with optical injection

We perform characterization of the pulse shape and noise properties of quantum dot passively mode-locked lasers (PMLLs). We propose a novel method to determine the RF linewidth and timing jitter, applicable to high repetition rate PMLLs, through the dependence of modal linewidth on the mode number. Complex electric field measurements show asymmetric pulses with parabolic phase close to threshold, with the appearance of waveform instabilities at higher currents. We demonstrate that the waveform instabilities can be overcome through optical injection-locking to the continues wave (CW) master laser, leading to time-bandwidth product (TBP) improvement, spectral narrowing, and spectral tunability. We discuss the benefits of single- and dual-tone master sources and demonstrate that dual-tone optical injection can additionally improve the noise properties of the slave laser with RF linewidth reduction below instrument limits (1 kHz) and integrated timing jitter values below 300 fs. Dual-tone injection allowed slave laser repetition rate control over a 25 MHz range with reduction of all modal optical linewidths to the master source linewidth, demonstrating phase-locking of all slave modes and coherence improvement.

Geographical variation in carbon dioxide fluxes from soils in agro-ecosystems and its implications for life-cycle assessment

1. Exchange of carbon dioxide (CO2) from soils can contribute significantly to the global warming potential (GWP) of agro-ecosystems. Due to variations in soil type, climatic onditions and land management practices, exchange of CO2 can differ markedly in different geographical locations. The food industry is developing carbon footprints for their products necessitating integration of CO2 exchange from soils with other CO2 emissions along the food chain. It may be advantageous to grow certain crops in different geographical locations to minimize CO2 emissions from the soil, and this may provide potential to offset other emissions in the food chain, such as transport. 2. Values are derived for the C balance of soils growing horticultural crops in the UK, Spain and Uganda. Net ecosystem production (NEP) is firstly calculated from the difference in net primary production (NPP) and heterotrophic soil respiration (Rh). Both NPP and Rh were estimated from intensive direct field measurements. Secondly, net biome production (NBP) is calculated by subtracting the crop biomass from NEP to give an indication of C balance. The importance of soil exchange is discussed in the light of recent discussions on carbon footprints and within the context of food life-cycle assessment (LCA). 3. The amount of crop relative to the biomass and the Rh prevailing in the different countries were the dominant factors influencing the magnitude of NEP and NBP. The majority of the biomass for lettuce Lactuca sativa and vining peas Pisum sativum, was removed from the field as crop; therefore, NEP and NBP were mainly negative. This was amplified for lettuces grown in Uganda (-16·5 and -17 t C ha-1 year-1 compared to UK and Spain -4·8 to 7·4 and -5·1 to 6·3 t C ha-1 year-1 for NEP and NBP, respectively) where the climate elevated Rh. 4. Synthesis and applications. This study demonstrates the importance of soil emissions in the overall life cycle of vegetables. Variability in such emissions suggests that assigning a single value to food carbon footprints may not be adequate, even within a country. Locations with high heterotrophic soil respiration, such as Spain and Uganda (21·9 and 21·6 t C ha-1 year-1, respectively), could mitigate the negative effects of climate on the C costs of crop production by growth of crops with greater returns of residue to the soil. This would minimize net CO2 emissions from these agricultural ecosystems.

A tidal creek water budget: Estimation of groundwater discharge and overland flow using hydrologic modeling in the Southern Everglades

Taylor Slough is one of the natural freshwater contributors to Florida Bay through a network of microtidal creeks crossing the Everglades Mangrove Ecotone Region (EMER). The EMER ecological function is critical since it mediates freshwater and nutrient inputs and controls the water quality in Eastern Florida Bay. Furthermore, this region is vulnerable to changing hydrodynamics and nutrient loadings as a result of upstream freshwater management practices proposed by the Comprehensive Everglades Restoration Program (CERP), currently the largest wetland restoration project in the USA. Despite the hydrological importance of Taylor Slough in the water budget of Florida Bay, there are no fine scale (∼1 km2) hydrodynamic models of this system that can be utilized as a tool to evaluate potential changes in water flow, salinity, and water quality. Taylor River is one of the major creeks draining Taylor Slough freshwater into Florida Bay. We performed a water budget analysis for the Taylor River area, based on long-term hydrologic data (1999–2007) and supplemented by hydrodynamic modeling using a MIKE FLOOD (DHI,http://dhigroup.com/) model to evaluate groundwater and overland water discharges. The seasonal hydrologic characteristics are very distinctive (average Taylor River wet vs. dry season outflow was 6 to 1 during 1999–2006) with a pronounced interannual variability of flow. The water budget shows a net dominance of through flow in the tidal mixing zone, while local precipitation and evapotranspiration play only a secondary role, at least in the wet season. During the dry season, the tidal flood reaches the upstream boundary of the study area during approximately 80 days per year on average. The groundwater field measurements indicate a mostly upwards-oriented leakage, which possibly equals the evapotranspiration term. The model results suggest a high importance of groundwater contribution to the water salinity in the EMER. The model performance is satisfactory during the dry season where surface flow in the area is confined to the Taylor River channel. The model also provided guidance on the importance of capturing the overland flow component, which enters the area as sheet flow during the rainy season. Overall, the modeling approach is suitable to reach better understanding of the water budget in the mangrove region. However, more detailed field data is needed to ascertain model predictions by further calibrating overland flow parameters.

Computational evaluation of wind loads on low- and highrise buildings

Buildings and other infrastructures located in the coastal regions of the US have a higher level of wind vulnerability. Reducing the increasing property losses and causalities associated with severe windstorms has been the central research focus of the wind engineering community. The present wind engineering toolbox consists of building codes and standards, laboratory experiments, and field measurements. The American Society of Civil Engineers (ASCE) 7 standard provides wind loads only for buildings with common shapes. For complex cases it refers to physical modeling. Although this option can be economically viable for large projects, it is not cost-effective for low-rise residential houses. To circumvent these limitations, a numerical approach based on the techniques of Computational Fluid Dynamics (CFD) has been developed. The recent advance in computing technology and significant developments in turbulence modeling is making numerical evaluation of wind effects a more affordable approach. The present study targeted those cases that are not addressed by the standards. These include wind loads on complex roofs for low-rise buildings, aerodynamics of tall buildings, and effects of complex surrounding buildings. Among all the turbulence models investigated, the large eddy simulation (LES) model performed the best in predicting wind loads. The application of a spatially evolving time-dependent wind velocity field with the relevant turbulence structures at the inlet boundaries was found to be essential. All the results were compared and validated with experimental data. The study also revealed CFD's unique flow visualization and aerodynamic data generation capabilities along with a better understanding of the complex three-dimensional aerodynamics of wind-structure interactions. With the proper modeling that realistically represents the actual turbulent atmospheric boundary layer flow, CFD can offer an economical alternative to the existing wind engineering tools. CFD's easy accessibility is expected to transform the practice of structural design for wind, resulting in more wind-resilient and sustainable systems by encouraging optimal aerodynamic and sustainable structural/building design. Thus, this method will help ensure public safety and reduce economic losses due to wind perils.

Development of effective approaches to the large-scale aerodynamic testing of low-rise buildings

Low-rise buildings are often subjected to high wind loads during hurricanes that lead to severe damage and cause water intrusion. It is therefore important to estimate accurate wind pressures for design purposes to reduce losses. Wind loads on low-rise buildings can differ significantly depending upon the laboratory in which they were measured. The differences are due in large part to inadequate simulations of the low-frequency content of atmospheric velocity fluctuations in the laboratory and to the small scale of the models used for the measurements. A new partial turbulence simulation methodology was developed for simulating the effect of low-frequency flow fluctuations on low-rise buildings more effectively from the point of view of testing accuracy and repeatability than is currently the case. The methodology was validated by comparing aerodynamic pressure data for building models obtained in the open-jet 12-Fan Wall of Wind (WOW) facility against their counterparts in a boundary-layer wind tunnel. Field measurements of pressures on Texas Tech University building and Silsoe building were also used for validation purposes. The tests in partial simulation are freed of integral length scale constraints, meaning that model length scales in such testing are only limited by blockage considerations. Thus the partial simulation methodology can be used to produce aerodynamic data for low-rise buildings by using large-scale models in wind tunnels and WOW-like facilities. This is a major advantage, because large-scale models allow for accurate modeling of architectural details, testing at higher Reynolds number, using greater spatial resolution of the pressure taps in high pressure zones, and assessing the performance of aerodynamic devices to reduce wind effects. The technique eliminates a major cause of discrepancies among measurements conducted in different laboratories and can help to standardize flow simulations for testing residential homes as well as significantly improving testing accuracy and repeatability. Partial turbulence simulation was used in the WOW to determine the performance of discontinuous perforated parapets in mitigating roof pressures. The comparisons of pressures with and without parapets showed significant reductions in pressure coefficients in the zones with high suctions. This demonstrated the potential of such aerodynamic add-on devices to reduce uplift forces.

Submarine slope failures in gassy soils

Submarine slope stability has become an important concern and a subject of research with increasing demand for offshore developments and technological advancement for harsh and challenging environments. The consequences of submarine slope failure adjacent to oil and gas facilities would have a large financial, safety and regulatory impact. This current research work investigates potential failure of submarine gassy slopes triggered by tidal variations. Due to tidal variations, failure of an unsaturated slope may occur under specific combinations of increasing degree of saturation and soil permeability, and decreasing tidal period. Novel physical model tests in a geotechnical centrifuge were undertaken to examine submarine slope failure mechanisms containing gassy sediments. The model preparation techniques, measurement systems and results are presented. The response observed in the model test is discussed and further developments proposed. The buried PPT’s response of the submarine slope are comparable in terms of attenuation and phase lag with Nagaswaran (1983) and with field measurements of Atigh and Byrne (2004) in terms of phase lag.

Crystallographic and magnetic identification of secondary phase in orientated Bi 5Fe 0.5Co 0.5Ti 3O 15 ceramics

The fabrication of highly-oriented polycrystalline ceramics of Bi 5Fe 0.5Co 0.5Ti 3O 15, prepared via molten salt synthesis and uniaxial pressing of high aspect ratio platelets is reported. Electron backscatter images show a secondary phase within the ceramic which is rich in cobalt and iron. The concentration of the secondary phase obtained from scanning electron microscopy is estimated at less than 2% by volume, below the detection limit of x-ray diffraction (XRD). The samples were characterized by x-ray diffraction, polarization-electric field measurements, superconducting quantum interference device as a function of sample orientation and vibrating sample magnetometry as a function of temperature. It is inferred from the data that the observed ferromagnetic response is dominated by the secondary phase. This work highlights the importance of rigorous materials characterisation in the study of multiferroics as small amounts of secondary phase, below the limit of XRD, can lead to false conclusions.

Land cover classification of Samoylov Island and Landsat subpixel water cover of Lena River Delta, Siberia, with links to ESRI grid files

Ignoring small-scale heterogeneities in Arctic land cover may bias estimates of water, heat and carbon fluxes in large-scale climate and ecosystem models. We investigated subpixel-scale heterogeneity in CHRIS/PROBA and Landsat-7 ETM+ satellite imagery over ice-wedge polygonal tundra in the Lena Delta of Siberia, and the associated implications for evapotranspiration (ET) estimation. Field measurements were combined with aerial and satellite data to link fine-scale (0.3 m resolution) with coarse-scale (upto 30 m resolution) land cover data. A large portion of the total wet tundra (80%) and water body area (30%) appeared in the form of patches less than 0.1 ha in size, which could not be resolved with satellite data. Wet tundra and small water bodies represented about half of the total ET in summer. Their contribution was reduced to 20% in fall, during which ET rates from dry tundra were highest instead. Inclusion of subpixel-scale water bodies increased the total water surface area of the Lena Delta from 13% to 20%. The actual land/water proportions within each composite satellite pixel was best captured with Landsat data using a statistical downscaling approach, which is recommended for reliable large-scale modelling of water, heat and carbon exchange from permafrost landscapes.

(Table 1) Mean annual ground temperatures and description of permafrost boreholes in Russia

The results of the International Permafrost Association's International Polar Year Thermal State of Permafrost (TSP) project are presented based on field measurements from Russia during the IPY years (2007-09) and collected historical data. Most ground temperatures measured in existing and new boreholes show a substantial warming during the last 20 to 30 years. The magnitude of the warming varied with location, but was typically from 0.5°C to 2°C at the depth of zero annual amplitude. Thawing of Little Ice Age permafrost is ongoing at many locations. There are some indications that the late Holocene permafrost has begun to thaw at some undisturbed locations in northeastern Europe and northwest Siberia. Thawing of permafrost is most noticeable within the discontinuous permafrost domain. However, permafrost in Russia is also starting to thaw at some limited locations in the continuous permafrost zone. As a result, a northward displacement of the boundary between continuous and discontinuous permafrost zones was observed. This data set will serve as a baseline against which to measure changes of near-surface permafrost temperatures and permafrost boundaries, to validate climate model scenarios, and for temperature reanalysis.

Sediment budget and geomorphic process units of the Swedish Kärkevagge, link to ArcInfo E00 files

Geomorphic process units have been derived in order to allow quantification via GIS techniques at a catchment scale. Mass movement rates based on existing field measurements are employed in the budget calculations. In the Kärkevagge catchment, Northern Sweden, 80% of the area can be identified either as a source area for sediments or as a zone where sediments are deposited. The overall budget for the slopes beneath the rockwalls in the Kärkevagge is approximately 680 t/a whilst about 150 t a-1 are transported into the fluvial system.

Monthly air temperatures, snow density and snow and ice thickness at Malcolm Ramsay Lake, North America

In northern regions where observational data is sparse, lake ice models are ideal tools as they can provide valuable information on ice cover regimes. The Canadian Lake Ice Model was used to simulate ice cover for a lake near Churchill, Manitoba, Canada throughout the 2008/2009 and 2009/2010 ice covered seasons. To validate and improve the model results, in situ measurements of the ice cover through both seasons were obtained using an upward-looking sonar device Shallow Water Ice Profiler (SWIP) installed on the bottom of the lake. The SWIP identified the ice-on/off dates as well as collected ice thickness measurements. In addition, a digital camera was installed on shore to capture images of the ice cover through the seasons and field measurements were obtained of snow depth on the ice, and both the thickness of snow ice (if present) and total ice cover. Altering the amounts of snow cover on the ice surface to represent potential snow redistribution affected simulated freeze-up dates by a maximum of 22 days and break-up dates by a maximum of 12 days, highlighting the importance of accurately representing the snowpack for lake ice modelling. The late season ice thickness tended to be under estimated by the simulations with break-up occurring too early, however, the evolution of the ice cover was simulated to fall between the range of the full snow and no snow scenario, with the thickness being dependant on the amount of snow cover on the ice surface.

Global distributions of diazotrophs abundance, biomass and nitrogen fixation rates - Gridded data product (NetCDF) - Contribution to the MAREDAT World Ocean Atlas of Plankton Functional Types

The MAREDAT atlas covers 11 types of plankton, ranging in size from bacteria to jellyfish. Together, these plankton groups determine the health and productivity of the global ocean and play a vital role in the global carbon cycle. Working within a uniform and consistent spatial and depth grid (map) of the global ocean, the researchers compiled thousands and tens of thousands of data points to identify regions of plankton abundance and scarcity as well as areas of data abundance and scarcity. At many of the grid points, the MAREDAT team accomplished the difficult conversion from abundance (numbers of organisms) to biomass (carbon mass of organisms). The MAREDAT atlas provides an unprecedented global data set for ecological and biochemical analysis and modeling as well as a clear mandate for compiling additional existing data and for focusing future data gathering efforts on key groups in key areas of the ocean. This is a gridded data product about diazotrophic organisms . There are 6 variables. Each variable is gridded on a dimension of 360 (longitude) * 180 (latitude) * 33 (depth) * 12 (month). The first group of 3 variables are: (1) number of biomass observations, (2) biomass, and (3) special nifH-gene-based biomass. The second group of 3 variables is same as the first group except that it only grids non-zero data. We have constructed a database on diazotrophic organisms in the global pelagic upper ocean by compiling more than 11,000 direct field measurements including 3 sub-databases: (1) nitrogen fixation rates, (2) cyanobacterial diazotroph abundances from cell counts and (3) cyanobacterial diazotroph abundances from qPCR assays targeting nifH genes. Biomass conversion factors are estimated based on cell sizes to convert abundance data to diazotrophic biomass. Data are assigned to 3 groups including Trichodesmium, unicellular diazotrophic cyanobacteria (group A, B and C when applicable) and heterocystous cyanobacteria (Richelia and Calothrix). Total nitrogen fixation rates and diazotrophic biomass are calculated by summing the values from all the groups. Some of nitrogen fixation rates are whole seawater measurements and are used as total nitrogen fixation rates. Both volumetric and depth-integrated values were reported. Depth-integrated values are also calculated for those vertical profiles with values at 3 or more depths.

Investigating the potential magnetic origin of wind variability in OB stars

In this thesis, the origin of large-scale structures in hot star winds, believed to be responsible for the presence of discrete absorption components (DACs) in the absorption troughs of ultraviolet resonance lines, is constrained using both observations and numerical simulations. These structures are understood as arising from bright regions on the stellar surface, although their physical cause remains unknown. First, we use high quality circular spectropolarimetric observations of 13 well-studied OB stars to evaluate the potential role of dipolar magnetic fields in producing DACs. We perform longitudinal field measurements and place limits on the field strength using Bayesian inference, assuming that it is dipolar. No magnetic field was detected within this sample. The derived constraints statistically refute any significant dynamical influence from a magnetic dipole on the wind for all of these stars, ruling out such fields as a cause for DACs. Second, we perform numerical simulations using bright spots constrained by broadband optical photometric observations. We calculate hydrodynamical wind models using three sets of spot sizes and strengths. Co-rotating interaction regions are yielded in each model, and radiative transfer shows that the properties of the variations in the UV resonance lines synthesized from these models are consistent with those found in observed UV spectra, establishing the first consistent link between UV spectroscopic line profile variability and photometric variations and thus supporting the bright spot paradigm (BSP). Finally, we develop and apply a phenomenological model to quantify the measurable effects co-rotating bright spots would have on broadband optical photometry and on the profiles of photopheric lines in optical spectra. This model can be used to evaluate the existence of these spots, and, in the event of their detection, characterize them. Furthermore, a tentative spot evolution model is presented. A preliminary analysis of its output, compared to the observed photometric variations of xi Persei, suggests the possible existence of “active longitudes” on the surface of this star. Future work will expand the range of observational diagnostics that can be interpreted within the BSP, and link phenomenology (bright spots) to physical processes (magnetic spots or non-radial pulsations).