Hydrological dynamics between a coastal aquifer and the adjacent estuarine system, Biscayne Bay, South Florida

Geochemical and geophysical approaches have been used to investigate the freshwater and saltwater dynamics in the coastal Biscayne Aquifer and Biscayne Bay. Stable isotopes of oxygen and hydrogen, and concentrations of Sr2+ and Ca2+ were combined in two geochemical mixing models to provide estimates of the various freshwater inputs (precipitation, canal water, and groundwater) to Biscayne Bay and the coastal canal system in South Florida. Shallow geophysical electromagnetic and direct current resistivity surveys were used to image the geometry and stratification of the saltwater mixing zone in the near coastal (less than 1km inland) Biscayne Aquifer. The combined stable isotope and trace metal models suggest a ratio of canal input-precipitation-groundwater of 38%–52%–10% in the wet season and 37%–58%–5% in the dry season with an error of 25%, where most (20%) of the error was attributed to the isotope regression model, while the remaining 5% error was attributed to the Sr2+/Ca2+ mixing model. These models suggest rainfall is the dominate source of freshwater to Biscayne Bay. For a bay-wide water budget that includes saltwater and freshwater mixing, fresh groundwater accounts for less than 2% of the total input. A similar Sr 2+/Ca2+ tracer model indicates precipitation is the dominate source in 9 out of 10 canals that discharge into Biscayne Bay. The two-component mixing model converged for 100% of the freshwater canal samples in this study with 63% of the water contributed to the canals coming from precipitation and 37% from groundwater inputs ±4%. There was a seasonal shift from 63% precipitation input in the dry season to 55% precipitation input in the wet season. The three end-member mixing model converged for only 60% of the saline canal samples possibly due to non-conservative behavior of Sr2+ and Ca2+ in saline groundwater discharging into the canal system. Electromagnetic and Direct Current resistivity surveys were successful at locating and estimating the geometry and depth of the freshwater/saltwater interface in the Biscayne Aquifer at two near coastal sites. A saltwater interface that deepened as the survey moved inland was detected with a maximum interpreted depth to the interface of 15 meters, approximately 0.33 km inland from the shoreline. ^

Statistical Models for Predicting College Success

Colleges base their admission decisions on a number of factors to determine which applicants have the potential to succeed. This study utilized data for students that graduated from Florida International University between 2006 and 2012. Two models were developed (one using SAT as the principal explanatory variable and the other using ACT as the principal explanatory variable) to predict college success, measured using the student’s college grade point average at graduation. Some of the other factors that were used to make these predictions were high school performance, socioeconomic status, major, gender, and ethnicity. The model using ACT had a higher R^2 but the model using SAT had a lower mean square error. African Americans had a significantly lower college grade point average than graduates of other ethnicities. Females had a significantly higher college grade point average than males.

Relative roles of dispersal dynamics and competition in determining the isotopic niche breadth of a wetland fish

1. The niche variation hypothesis predicts that among-individual variation in niche use will increase in the presence of intraspecific competition and decrease in the presence of interspecific competition. We sought to determine whether the local isotopic niche breadth of fish inhabiting a wetland was best explained by competition for resources and the niche variation hypothesis, by dispersal of individuals from locations with different prey resources or by a combination of the two. We analysed stable isotopes of carbon and nitrogen as indices of feeding niche and compared metrics of within-site spread to characterise site-level isotopic niche breadth. We then evaluated the explanatory power of competing models of the direct and indirect effects of several environmental variables spanning gradients of disturbance, competition strength and food availability on among-individual variation of the eastern mosquitofish (Gambusia holbrooki). 2. The Dispersal model posits that only the direct effect of disturbance (i.e. changes in water level known to induce fish movement) influences among-individual variation in isotopic niche. The Partitioning model allows for only direct effects of local food availability on among-individual variation. The Combined model allows for both hypotheses by including the direct effects of disturbance and food availability. 3. A linear regression of the Combined model described more variance than models limited to the variables of either the Dispersal or Partitioning models. Of the independent variables considered, the food availability variable (per cent edible periphyton) explained the most variation in isotopic niche breadth, followed closely by the disturbance variable (days since last drying event). 4. Structural equation modelling provided further evidence that the Combined model was best supported by the data, with the Partitioning and the Dispersal models only modestly less informative. Again, the per cent edible periphyton was the variable with the largest direct effect on niche variability, with other food availability variables and the disturbance variable only slightly less important. Indirect effects of heterospecific and conspecific competitor densities were also important, through their effects on prey density. 5. Our results support the Combined hypotheses, although partitioning mechanisms appear to explain the most diet variation among individuals in the eastern mosquitofish. The results also support some predictions of the niche variation hypothesis, although both conspecific and interspecific competition appeared to increase isotopic niche breadth in contrast to predictions that interspecific competition would decrease it. We think this resulted from high diet overlap of co-occurring species, most of which consume similar macroinvertebrates.

Evaluating the effectiveness of signal timing optimization based on microscopic simulation evaluation

The optimization of the timing parameters of traffic signals provides for efficient operation of traffic along a signalized transportation system. Optimization tools with macroscopic simulation models have been used to determine optimal timing plans. These plans have been, in some cases, evaluated and fine tuned using microscopic simulation tools. A number of studies show inconsistencies between optimization tool results based on macroscopic simulation and the results obtained from microscopic simulation. No attempts have been made to determine the reason behind these inconsistencies. This research investigates whether adjusting the parameters of macroscopic simulation models to correspond to the calibrated microscopic simulation model parameters can reduce said inconsistencies. The adjusted parameters include platoon dispersion model parameters, saturation flow rates, and cruise speeds. The results from this work show that adjusting cruise speeds and saturation flow rates can have significant impacts on improving the optimization/macroscopic simulation results as assessed by microscopic simulation models.

Sedimentology of the northwestern Weddel Sea continental shelf and slope

Geologie cores on two profiles oriented normaly to the continental shelf and slope, have been investigated to reconstruct the Quaternary sedimentary history of the southeast continental border of South Orkney (NW Weddell Sea). The sediments were described macroscopically and their fabric investigated by use of X-radiographs. Laboratory work comprised detailed grain-size analysis, determination of the watercontent, carbonate, organic carbon and sand fraction.composition. Stable oxygen and carbon isotopes have been measured On planktonic foraminifera. Palaeomagnetism, analysis of 230Th-content and detailed comparison of the lithlogic Parameters with the oxygen isotope stages (Martinson curve) were used for stratigraphic classification of the sediments. The sediment cores from the continental slope comprise a maximum age of 300,000 years B. P.. Bottom currents, ice rafting and biogenic input are the main sources of sediment. Based on lithologic parameters a distinction between glacial and interglacial facies is possible. Silty clays without microfossils and few bioturbation characterise the sediments of the glacial facies. Only small amounts of icerafted debris can be recognized. This type of sediment was accumulated during times of lower sea-level and drastically reduced rate of bottom water production. Based on grain-size distribution, bottom current velocities of 0.01 cmls were calculated. Thick sea-ice coverage reduced biogenic production in the surface water, and as consequence benthic communities were depleted. Because of the reduced benthic life, sediments are only slithly bioturbated. At the beginning of the interglacial Stage, the sea-level rised rapidly, and calving rate of icebergs, combined with input of ice-rafted material, increased considerably. Sediments of this transition facies are silty cliiys with a high proportion of coarse ice-rafted debris, but without microfossils. With the onset of bottom water production in connection with shelf ice water, sediments of interglacial facies were formed. They consist of silty clays to clayey silts with considerable content of sand and gravel. Sediments are strongly bioturbated. Based On the sediment caracteristics, current velocities of the bottom water were calculated to be of 0.96 cmls for interglacials. At the southern slope of a NW/SE-striking ridge, bottom water current is channelized, resulting in a drastic increase of current velocities. Current velocities up to 7.5 cm/s lead to formation of residual sediments. While the continental slope has predominantly fine sediments, the South Orkney shelf are mainly sandy silts and silty sands with a high proportion of gravel. These sediments were formed dominantly by ice-rafting during Brunhes- and Matuyama-Epoch. Currents removed the fine fraction of the sediments. Based on microfossil contents it was not possible to differentiate sediments from glacial to interglacial. In the upper Parts of the cores graded sequences truncated by erosion were observed. These sequences were formed during Brunhes-Epoch by strong currents with velocities decreasing periodically from about 7.5 cm/s to about 1 cm/s. Sediments with a high proportion of siliceous microfossils but barren of foraminifera compose the lower part of the shelf cores. These sediments have formed during the warmer Matuyama-Epoch.

Investigations on sediment coring site PS1823 in the Riiser-Larsen Sea, Southern Ocean

Physical and sedimentological investigations were carried out on a 14 m long gravity core and a 0.5 m long box core from 4440 m water depth off Queen Maud Land, East-Antarctica. Strongly bioturbated hemipelagic muds of predominantly terrigenous origin and a very small biogenic part build up the 'Normal-Facies'. Several sandy to silty layers are inserted in the 'Normal-Facies'. These layers are seperated by lithology, structure and the investigated parameters of this study and are interpreted as turbidites. The source area for the turbidity currents is supposed to be at the uppermost continental margin, close to the shelf break and there is evidenee for this gravity transport within the erosive Ritscher-Canyon, which extends close to the core position. The distribution of biogenic components indicates an age of 1.3 million years or more, with an average sedimentation rate of about 1 cm/1000 years. Early diagenetic proeesses caused water loss by compaction, errosion and dissolution of biogenic components and precipitation and recrystallization of manganese micronodules. Cyclic fluctuations of the sediment-parameters within the 'Normal-Facies' enable the distinction of a 'Glazial'- and an 'Interglazial'-Facies. The 'Glazial'-Facies reflects glacial sedimentary conditions and shows a dark olive gray colour, high susceptibility, low silt/clay-ratios, only a few biogenic components and the regular occurence of interrelated turbidite layers. In contrast, the 'Interglazial'-Facies is dominated by a light olive or olive-brown colour, low susceptibility, high silt/clay-ratios and an increased number of biogenic components. This facies corresponds to interglacial conditions. Three main processes are supposed to have been responsible for the observed facies changes: (1) the bottom water mass circulation, (2) the gravity transport by turbidity currents and (3) the biogenic surface production. These processes are related to the quaternary climatic changes. The extension of the ice shelves directed the gravity transport to the deep sea and the formation of Antarctic Bottom Water, which in turn influenced the silt/clay-ratios in the sediment record. Fluctuations in sea ice coverage controlled the biogenic surface production.

Sources and mixing state of size-resolved elemental carbon particles in a European megacity: Paris

An Aerosol Time-Of-Flight Mass Spectrometer (ATOFMS) was deployed to investigate the size-resolved chemical composition of single particles at an urban background site in Paris, France, as part of the MEGAPOLI winter campaign in January/February 2010. ATOFMS particle counts were scaled to match coincident Twin Differential Mobility Particle Sizer (TDMPS) data in order to generate hourly size-resolved mass concentrations for the single particle classes observed. The total scaled ATOFMS particle mass concentration in the size range 150–1067 nm was found to agree very well with the sum of concurrent High-Resolution Time-of-Flight Aerosol Mass Spectrometer (HR-ToF-AMS) and Multi-Angle Absorption Photometer (MAAP) mass concentration measurements of organic carbon (OC), inorganic ions and black carbon (BC) (R2 = 0.91). Clustering analysis of the ATOFMS single particle mass spectra allowed the separation of elemental carbon (EC) particles into four classes: (i) EC attributed to biomass burning (ECbiomass), (ii) EC attributed to traffic (ECtraffic), (iii) EC internally mixed with OC and ammonium sulfate (ECOCSOx), and (iv) EC internally mixed with OC and ammonium nitrate (ECOCNOx). Average hourly mass concentrations for EC-containing particles detected by the ATOFMS were found to agree reasonably well with semi-continuous quantitative thermal/optical EC and optical BC measurements (r2 = 0.61 and 0.65–0.68 respectively, n = 552). The EC particle mass assigned to fossil fuel and biomass burning sources also agreed reasonably well with BC mass fractions assigned to the same sources using seven-wavelength aethalometer data (r2 = 0.60 and 0.48, respectively, n = 568). Agreement between the ATOFMS and other instrumentation improved noticeably when a period influenced by significantly aged, internally mixed EC particles was removed from the intercomparison. 88% and 12% of EC particle mass was apportioned to fossil fuel and biomass burning respectively using the ATOFMS data compared with 85% and 15% respectively for BC estimated from the aethalometer model. On average, the mass size distribution for EC particles is bimodal; the smaller mode is attributed to locally emitted, mostly externally mixed EC particles, while the larger mode is dominated by aged, internally mixed ECOCNOx particles associated with continental transport events. Periods of continental influence were identified using the Lagrangian Particle Dispersion Model (LPDM) "FLEXPART". A consistent minimum between the two EC mass size modes was observed at approximately 400 nm for the measurement period. EC particles below this size are attributed to local emissions using chemical mixing state information and contribute 79% of the scaled ATOFMS EC particle mass, while particles above this size are attributed to continental transport events and contribute 21% of the EC particle mass. These results clearly demonstrate the potential benefit of monitoring size-resolved mass concentrations for the separation of local and continental EC emissions. Knowledge of the relative input of these emissions is essential for assessing the effectiveness of local abatement strategies.

Essays on Firm Heterogeneity with Empirical Applications in Economic History and Agricultural Economics

This thesis uses models of firm-heterogeneity to complete empirical analyses in economic history and agricultural economics. In Chapter 2, a theoretical model of firm heterogeneity is used to derive a statistic that summarizes the welfare gains from the introduction of a new technology. The empirical application considers the use of mechanical steam power in the Canadian manufacturing sector during the late nineteenth century. I exploit exogenous variation in geography to estimate several parameters of the model. My results indicate that the use of steam power resulted in a 15.1 percent increase in firm-level productivity and a 3.0-5.2 percent increase in aggregate welfare. Chapter 3 considers various policy alternatives to price ceiling legislation in the market for production quotas in the dairy farming sector in Quebec. I develop a dynamic model of the demand for quotas with farmers that are heterogeneous in their marginal cost of milk production. The econometric analysis uses farm-level data and estimates a parameter of the theoretical model that is required for the counterfactual experiments. The results indicate that the price of quotas could be reduced to the ceiling price through a 4.16 percent expansion of the aggregate supply of quotas, or through moderate trade liberalization of Canadian dairy products. In Chapter 4, I study the relationship between farm-level productivity and participation in the Commercial Export Milk (CEM) program. I use a difference-in-difference research design with inverse propensity weights to test for causality between participation in the CEM program and total factor productivity (TFP). I find a positive correlation between participation in the CEM program and TFP, however I find no statistically significant evidence that the CEM program affected TFP.

Evaluación del estado poblaciónal del Cangrejo violáceo Platyxanthus orbignyi (Milne Edwards y Lucas, 1843) del área de Lambayeque por medio de un modelo dinámico de Biomasa

Se hace un análisis del estado poblacional del cangrejo violáceo Platyxanthus orbignyi (Milne Edwards y Lucas, 1843) del litoral de Lambayeque – Perú para el periodo 2001-2010 por medio de: 1) El modelo dinámico de biomasa de Schaefer en su versión de error de observación; a este modelo se le introdujo la variable ambiental anomalía de la temperatura superficial del mar (ATSM) del área de San José (Lambayeque) y se obtuvo 2) el modelo dinámico con variable ambiental, ambos basados en datos de captura, esfuerzo y CPUE. Se utilizó el método de máxima verosimilitud en el proceso de ajuste y el bootstrap para determinar los intervalos de confianza de los parámetros. Los parámetros poblacionales y pesqueros estimados por el modelo dinámico de biomasa de Schaefer (MDB) fueron: K: 750 000 kg, r : 0,21 y q: 8,36 x 10-6 y por el modelo dinámico con variable ambiental (MDVA) los parámetros fueron K: 765 000 kg, r: 0,23 y q: 8,02 x 10-6. Con los valores de los parámetros estimados mediante el MDB y el MDVA se calcularon los principales puntos biológicos de referencia (PBR) los cuales fueron: MRS: 39 822 kg, BMRS: 375 000 kg, fMRS: 12 561 nasas, FMRS: 0,11, F0.1: 0,10 para el MDB; y MRS: 44 069 kg, BMRS: 382 500 kg, fMRS: 13 782 nasas, FMRS: 0,12, F0.1: 0,10 para el MDVA. Los resultados indican que el estado actual de la pesquería del cangrejo violáceo del Litoral de Lambayeque se encuentra muy cerca al nivel óptimo. En vista de que no se dispone de información de evaluaciones directas de este recurso que confirme o no los resultados del MDB y MDVA y en virtud de la calidad de datos, se sugiere que el manejo de la pesquería sea del tipo adaptativo alrededor del punto de referencia F0.1 y teniendo en cuenta las condiciones ambientales.

Intercity Bus Route Evaluation for Statewide Planning, 1981

The goal of this research was to evaluate the needs of the intercity common carrier bus service in Iowa. Within the framework of the overall goal, the objectives were to: (1) Examine the detailed operating cost and revenue data of the intercity carriers in Iowa; (2) Develop a model or models to estimate demand in cities and corridors served by the bus industry; (3) Develop a cost function model for estimating a carrier's operating costs; (4) Establish the criteria to be used in assessing the need for changes in bus service; (5) Outline the procedures for estimating route operating costs and revenues and develop a matrix of community and social factors to be considered in evaluation; and (6) Present a case study to demonstrate the methodology. The results of the research are presented in the following chapters: (1) Introduction; (2) Intercity Bus Research and Development; (3) Operating Characteristics of Intercity Carriers in Iowa; (4) Commuter Carriers; (5) Passenger and Revenue Forecasting Models; (6) Operating Cost Relationships; (7) Social and General Welfare Aspects of Intercity Bus Service; (8) Case Study Analysis; and (9) Additional Service Considerations and Recommendations.

Error Recovery Performance: The Impact of Leader Behavioral Integrity and Job Satisfaction

This article introduces the concept of error recovery performance, followed by the development and validation of an instrument to measure it. The first objective of this article is to broaden the current concept of service recovery to be relevant to the back-of-house operations. The second objective is to examine the influence of leader behavioral integrity (BI) on error recovery performance. Moreover, the study examines the mediating effect of job satisfaction between BI and error recovery performance. Finally, the study links error management performance with work-unit effectiveness. Data for Study 1 were collected from 369 hotel employees in Turkey. The same relationships were tested again in Study 2 to validate the findings of Study 1 with a different sample. Data for Study 2 were collected from 33 departmental managers from the same hotels. Linear regression analysis was used to test the direct effects. The mediating effects were tested using the mediation test suggested by Preacher and Hayes. In addition, in Study 2, general managers of the hotels were asked to rate the effectiveness of each manager and their respective department. Results from Study 1 indicate that BI drives error recovery performance, and this impact is mediated by employee job satisfaction. Results of Study 2 confirm this model and finds further that managers’ self-rated error recovery performance was associated with their general managers’ assessment of their deliverables and of their department’s overall performance.

Hoshin Kanri – the Japanese way of piloting : An exploratory study of a Japanese strategic management system

Strategy is a highly topical subject among managers and since the world is constantlychanging it is also an important subject for companies’ competitive advantage and survival.At the same time experts in the field of strategic management describe western techniques ascomplex and ineffective while the Japanese techniques have been seen as unambiguous andcharacterized by focus on quality, productivity and teamwork. This calls for greaterknowledge in the Japanese management systems. Hoshin Kanri is a collection of Japanesebest strategic management practices and therefore an interesting target for our study. Thus, onthe one hand this study investigates the theory of Hoshin Kanri in order to give structure to itand provide a way for practitioner into the management system. On the other hand this studyinvestigates Hoshin Kanri in order to reveal how Japanese subsidiaries based in Sweden haveimplemented this strategic management system. This is firstly done by reviewing the existingliterature on the subject and secondly by a collective case study with in-depth interviewsconducted with managers at Japanese owned subsidiaries based in Sweden. There are somelimitations in this study. One is that the results of the study do not include all Japanesesubsidiaries in Sweden as not all companies participated in the study. Moreover, the study islimited by one individuals’ knowledge and perception of Hoshin Kanri in each of the threecompanies. The study contributes to the existing literature on the topic of Hoshin Kanri by;(1) structuring the literature and the existing models under one of two categories, namelycyclical or sequential; (2) providing a model that aims at making it more understandable andattractive for practitioner to apply; (3) initiating the mapping of the spread of Hoshin Kanriamong Japanese subsidiaries in Sweden and (4) providing a Swedish model for theapplication of HK in Japanese subsidiaries.

Lifetime and distribution of ozone and related pollutants in the eastern United States

Most major cities in the eastern United States have air quality deemed unhealthy by the EPA under a set of regulations known as the National Ambient Air Quality Standards (NAAQS). The worst air quality in Maryland is measured in Edgewood, MD, a small community located along the Chesapeake Bay and generally downwind of Baltimore during hot, summertime days. Direct measurements and numerical simulations were used to investigate how meteorology and chemistry conspire to create adverse levels of photochemical smog especially at this coastal location. Ozone (O3) and oxidized reactive nitrogen (NOy), a family of ozone precursors, were measured over the Chesapeake Bay during a ten day experiment in July 2011 to better understand the formation of ozone over the Bay and its impact on coastal communities such as Edgewood. Ozone over the Bay during the afternoon was 10% to 20% higher than the closest upwind ground sites. A combination of complex boundary layer dynamics, deposition rates, and unaccounted marine emissions play an integral role in the regional maximum of ozone over the Bay. The CAMx regional air quality model was assessed and enhanced through comparison with data from NASA’s 2011 DISCOVER-AQ field campaign. Comparisons show a model overestimate of NOy by +86.2% and a model underestimate of formaldehyde (HCHO) by –28.3%. I present a revised model framework that better captures these observations and the response of ozone to reductions of precursor emissions. Incremental controls on electricity generating stations will produce greater benefits for surface ozone while additional controls on mobile sources may yield less benefit because cars emit less pollution than expected. Model results also indicate that as ozone concentrations improve with decreasing anthropogenic emissions, the photochemical lifetime of tropospheric ozone increases. The lifetime of ozone lengthens because the two primary gas-phase sinks for odd oxygen (Ox ≈ NO2 + O3) – attack by hydroperoxyl radicals (HO2) on ozone and formation of nitrate – weaken with decreasing pollutant emissions. This unintended consequence of air quality regulation causes pollutants to persist longer in the atmosphere, and indicates that pollutant transport between states and countries will likely play a greater role in the future.

Quantification of the human postural control system to perturbations

Human standing posture is inherently unstable. The postural control system (PCS), which maintains standing posture, is composed of the sensory, musculoskeletal, and central nervous systems. Together these systems integrate sensory afferents and generate appropriate motor efferents to adjust posture. The PCS maintains the body center of mass (COM) with respect to the base of support while constantly resisting destabilizing forces from internal and external perturbations. To assess the human PCS, postural sway during quiet standing or in response to external perturbation have frequently been examined descriptively. Minimal work has been done to understand and quantify the robustness of the PCS to perturbations. Further, there have been some previous attempts to assess the dynamical systems aspects of the PCS or time evolutionary properties of postural sway. However those techniques can only provide summary information about the PCS characteristics; they cannot provide specific information about or recreate the actual sway behavior. This dissertation consists of two parts: part I, the development of two novel methods to assess the human PCS and, part II, the application of these methods. In study 1, a systematic method for analyzing the human PCS during perturbed stance was developed. A mild impulsive perturbation that subjects can easily experience in their daily lives was used. A measure of robustness of the PCS, 1/MaxSens that was based on the inverse of the sensitivity of the system, was introduced. 1/MaxSens successfully quantified the reduced robustness to external perturbations due to age-related degradation of the PCS. In study 2, a stochastic model was used to better understand the human PCS in terms of dynamical systems aspect. This methodology also has the advantage over previous methods in that the sway behavior is captured in a model that can be used to recreate the random oscillatory properties of the PCS. The invariant density which describes the long-term stationary behavior of the center of pressure (COP) was computed from a Markov chain model that was applied to postural sway data during quiet stance. In order to validate the Invariant Density Analysis (IDA), we applied the technique to COP data from different age groups. We found that older adults swayed farther from the centroid and in more stochastic and random manner than young adults. In part II, the tools developed in part I were applied to both occupational and clinical situations. In study 3, 1/MaxSens and IDA were applied to a population of firefighters to investigate the effects of air bottle configuration (weight and size) and vision on the postural stability of firefighters. We found that both air bottle weight and loss of vision, but not size of air bottle, significantly decreased balance performance and increased fall risk. In study 4, IDA was applied to data collected on 444 community-dwelling elderly adults from the MOBILIZE Boston Study. Four out of five IDA parameters were able to successfully differentiate recurrent fallers from non-fallers, while only five out of 30 more common descriptive and stochastic COP measures could distinguish the two groups. Fall history and the IDA parameter of entropy were found to be significant risk factors for falls. This research proposed a new measure for the PCS robustness (1/MaxSens) and a new technique for quantifying the dynamical systems aspect of the PCS (IDA). These new PCS analysis techniques provide easy and effective ways to assess the PCS in occupational and clinical environments.

Three-dimensional kinetic simulations of active suspensions: effect of chemotaxis and steric interaction

In this work, the effects of chemotaxis and steric interactions in active suspensions are analyzed by extending the kinetic model proposed by Saintillan and Shelley [1, 2]. In this model, a conservation equation for the active particle configuration is coupled to the Stokes equation for the flow arising from the force dipole exerted by the particles on the fluid. The fluid flow equations are solved spectrally and the conservation equation is solved by second-order finite differencing in space and second-order Adams-Bashforth time marching. First, the dynamics in suspensions of oxytactic run-and-tumble bacteria confined in thin liquid films surrounded by air is investigated. These bacteria modify their tumbling behavior by making temporal comparisons of the oxygen concentration, and, on average, swim towards high concentrations of oxygen. The kinetic model proposed by Saintillan and Shelley [1, 2] is modified to include run-and-tumble effects and oxygentaxis. The spatio-temporal dynamics of the oxygen and bacterial concentration are analyzed. For small film thicknesses, there is a weak migration of bacteria to the boundaries, and the oxygen concentration is high inside the film as a result of diffusion; both bacterial and oxygen concentrations quickly reach steady states. Above a critical film thickness (approximately 200 micron), a transition to chaotic dynamics is observed and is characterized by turbulent-like 3D motion, the formation of bacterial plumes, enhanced oxygen mixing and transport into the film, and hydrodynamic velocities of magnitudes up to 7 times the single bacterial swimming speed. The simulations demonstrate that the combined effects of hydrodynamic interactions and oxygentaxis create collective three-dimensional instabilities which enhances oxygen availability for the bacteria. Our simulation results are consistent with the experimental findings of Sokolov et al. [3], who also observed a similar transition with increasing film thickness. Next, the dynamics in concentrated suspensions of active self-propelled particles in a 3D periodic domain are analyzed. We modify the kinetic model of Saintillan and Shelley [1, 2] by including an additional nematic alignment torque proportional to the local concentration in the equation for the rotational velocity of the particles, causing them to align locally with their neighbors (Doi and Edwards [4]). Large-scale three- dimensional simulations show that, in the presence of such a torque both pusher and puller suspensions are unstable to random fluctuations and are characterized by highly nematic structures. Detailed measures are defined to quantify the degree and direction of alignment, and the effects of steric interactions on pattern formation will be presented. Our analysis shows that steric interactions have a destabilizing effect in active suspensions.