Interlinking ecosystem services and Ostrom’s framework through orientation in sustainability research

Structuring integrated social-ecological systems (SES) research remains a core challenge for achieving sustainability. Numerous concepts and frameworks exist, but there is a lack of mutual learning and orientation of knowledge between them. We focus on two approaches in particular: the ecosystem services concept and Elinor Ostrom’s diagnostic SES framework. We analyze the strengths and weaknesses of each and discuss their potential for mutual learning. We use knowledge types in sustainability research as a boundary object to compare the contributions of each approach. Sustainability research is conceptualized as a multi-step knowledge generation process that includes system, target, and transformative knowledge. A case study of the Southern California spiny lobster fishery is used to comparatively demonstrate how each approach contributes a different lens and knowledge when applied to the same case. We draw on this case example in our discussion to highlight potential interlinkages and areas for mutual learning. We intend for this analysis to facilitate a broader discussion that can further integrate SES research across its diverse communities.

Modelling of water removal during a paper vacuum dewatering process using a Level-Set method

Water removal in paper manufacturing is an energy-intensive process. The dewatering process generally consists of four stages of which the first three stages include mechanical water removal through gravity filtration, vacuum dewatering and wet pressing. In the fourth stage, water is removed thermally, which is the most expensive stage in terms of energy use. In order to analyse water removal during a vacuum dewatering process, a numerical model was created by using a Level-Set method. Various different 2D structures of the paper model were created in MATLAB code with randomly positioned circular fibres with identical orientation. The model considers the influence of the forming fabric which supports the paper sheet during the dewatering process, by using volume forces to represent flow resistance in the momentum equation. The models were used to estimate the dry content of the porous structure for various dwell times. The relation between dry content and dwell time was compared to laboratory data for paper sheets with basis weights of 20 and 50 g/m2 exposed to vacuum levels between 20 kPa and 60 kPa. The comparison showed reasonable results for dewatering and air flow rates. The random positioning of the fibres influences the dewatering rate slightly. In order to achieve more accurate comparisons, the random orientation of the fibres needs to be considered, as well as the deformation and displacement of the fibres during the dewatering process.

Shapes and Dynamics of Blood Cells in Poiseuille and Shear Flows

The dynamics, shape, deformation, and orientation of red blood cells in microcirculation affect the rheology, flow resistance and transport properties of whole blood. This leads to important correlations of cellular and continuum scales. Furthermore, the dynamics of RBCs subject to different flow conditions and vessel geometries is relevant for both fundamental research and biomedical applications (e.g drug delivery). In this thesis, the behaviour of RBCs is investigated for different flow conditions via computer simulations. We use a combination of two mesoscopic particle-based simulation techniques, dissipative particle dynamics and smoothed dissipative particle dynamics. We focus on the microcapillary scale of several μm. At this scale, blood cannot be considered at the continuum but has to be studied at the cellular level. The connection between cellular motion and overall blood rheology will be investigated. Red blood cells are modelled as viscoelastic objects interacting hydrodynamically with a viscous fluid environment. The properties of the membrane, such as resistance against bending or shearing, are set to correspond to experimental values. Furthermore, thermal fluctuations are considered via random forces. Analyses corresponding to light scattering measurements are performed in order to compare to experiments and suggest for which situations this method is suitable. Static light scattering by red blood cells characterises their shape and allows comparison to objects such as spheres or cylinders, whose scattering signals have analytical solutions, in contrast to those of red blood cells. Dynamic light scattering by red blood cells is studied concerning its suitability to detect and analyse motion, deformation and membrane fluctuations. Dynamic light scattering analysis is performed for both diffusing and flowing cells. We find that scattering signals depend on various cell properties, thus allowing to distinguish different cells. The scattering of diffusing cells allows to draw conclusions on their bending rigidity via the effective diffusion coefficient. The scattering of flowing cells allows to draw conclusions on the shear rate via the scattering amplitude correlation. In flow, a RBC shows different shapes and dynamic states, depending on conditions such as confinement, physiological/pathological state and cell age. Here, two essential flow conditions are studied: simple shear flow and tube flow. Simple shear flow as a basic flow condition is part of any more complex flow. The velocity profile is linear and shear stress is homogeneous. In simple shear flow, we find a sequence of different cell shapes by increasing the shear rate. With increasing shear rate, we find rolling cells with cup shapes, trilobe shapes and quadrulobe shapes. This agrees with recent experiments. Furthermore, the impact of the initial orientation on the dynamics is studied. To study crowding and collective effects, systems with higher haematocrit are set up. Tube flow is an idealised model for the flow through cylindric microvessels. Without cell, a parabolic flow profile prevails. A single red blood cell is placed into the tube and subject to a Poiseuille profile. In tube flow, we find different cell shapes and dynamics depending on confinement, shear rate and cell properties. For strong confinements and high shear rates, we find parachute-like shapes. Although not perfectly symmetric, they are adjusted to the flow profile and maintain a stationary shape and orientation. For weak confinements and low shear rates, we find tumbling slippers that rotate and moderately change their shape. For weak confinements and high shear rates, we find tank-treading slippers that oscillate in a limited range of inclination angles and strongly change their shape. For the lowest shear rates, we find cells performing a snaking motion. Due to cell properties and resultant deformations, all shapes differ from hitherto descriptions, such as steady tank-treading or symmetric parachutes. We introduce phase diagrams to identify flow regimes for the different shapes and dynamics. Changing cell properties, the regime borders in the phase diagrams change. In both flow types, both the viscosity contrast and the choice of stress-free shape are important. For in vitro experiments, the solvent viscosity has often been higher than the cytosol viscosity, leading to a different pattern of dynamics, such as steady tank-treading. The stress-free state of a RBC, which is the state at zero shear stress, is still controversial, and computer simulations enable direct comparisons of possible candidates in equivalent flow conditions.

Implications of TiO2 surface functionalization on polycrystalline mixed halide perovskite films and photovoltaic devices

We exploit TiO2 surface functionalization as a tool to induce the crystallization process of CH3NH3PbI3xClx perovskite thin films resulting in a reduction of the degree of orientation of the (110) crystallographic planes. Notably, the variation of the film crystalline orientational order does not affect the photovoltaic performances of the perovskite-based devices, whose efficiency remains mostly unchanged. Our findings suggest that other factors are more significant in determining the device efficiency, such as the non-homogenous coverage of the TiO2 surface causing charge recombination at the organic/TiO2 interface, defect distribution on the perovskite bulk and at the interfaces, and transport in the organic or TiO2 layer. This observation represents a step towards the comprehension of the perovskite film peculiarities influencing the photovoltaic efficiency for high performance devices.

Superlattice structures in twisted bilayers of folded graphene

The electronic properties of bilayer graphene strongly depend on relative orientation of the two atomic lattices. Whereas Bernal-stacked graphene is most commonly studied, a rotational mismatch between layers opens up a whole new field of rich physics, especially at small interlayer twist. Here we report on magnetotransport measurements on twisted graphene bilayers, prepared by folding of single layers. These reveal a strong dependence on the twist angle, which can be estimated by means of sample geometry. At small rotation, superlattices with a wavelength in the order of 10 nm arise and are observed by friction atomic force microscopy. Magnetotransport measurements in this small-angle regime show the formation of satellite Landau fans. These are attributed to additional Dirac singularities in the band structure and discussed with respect to the wide range of interlayer coupling models.

Report of the Working Group on Introductions and Transfers of Marine Organisms (WGITMO). 16-18 March 2016 Olbia, Italy

The 42th meeting of the ICES Working Group on Introductions and Transfers of Marine Organisms (WGITMO) was held in Olbia, Italy, 16–18 March 2016, with Anna Occhipinti-Ambrogi as host and Henn Ojaveer as chairperson. Representatives from 19 countries participated in the meeting. Attendants were from Belgium, Canada, Dennark, Estonia, Finland, France, Germany, Ireland, Israel, Italy, Lithuania, Norway, Poland, Portugal, Russia, Spain, Sweden, United Kingdom and United States. Sweden contributed by cor-respondence. The objectives of the meeting were to update information and discuss several aspects related to the introductions and transfers of non-indigenous aquatic species. Data and information management were two of the discussion topics of the meeting, with special focus on the better exploitation of the ‘Information system on aquatic non-indigenous and cryptogenic species’ (AquaNIS). The WGITMO also dedicated time for addressing the MSFD D2 issues: indicator on new non-indigenous species introduced by human activities, and opportunities and problems related to cross-regional comparison of non-indigenous species indicators. Preparation of the manuscript of the alert report on sea squirt Didemnum vexillum, which is to be published in ICES CRR series, was discussed and the steps to be taken to finalise the report were agreed. As usual, adequate time was devoted to discuss national reports, to exchange of information on the management of NIS and to review ongoing and planned research activities. The approach taken during the meeting facilitated presentations and discussions on the issues of relevance related to the Terms of References as well as on a few generic and strategically-important issues of general relevance to bioinvasions. The meeting began with a full-day joint meeting with the Working Group on Ballast and Other Ship Vectors (WGBOSV), which provided an opportunity to discuss and address issues of common interest, such as shipping and biofouling as introduction vectors. The proposed ICES demonstration advice on ‘Risk management of non-indigenous species associated with shipping in the Arctic’ was discussed, and edits were suggested for both the orientation of the demonstration advice as well as for the exact questions to be asked. Both working groups agreed that the practice of conducting back-to-back meetings with one joint day is useful and will continue in 2017. All Terms of References to be addressed for 2016 were discussed. For some Terms of Ref-erence, more detailed presentations were given, and a short overview of the information and subsequent discussion is provided herein at the end of each section. This report is structured so that each Term of Reference is dealt with in sequential order. The main body of the report contains summaries of the presentations and discussions with the more detailed documents being contained in the Annexes. WGITMO progressed each of the Terms of Reference by either completing the task or clearly identifying and agreeing on the inter-sessional activities required to still finalise the work in 2016. From 2017, WGITMO will be shifted to multi-annual management.

The impact of mechanical compression on cortical microtubules in Arabidopsis: a quantitative pipeline

Exogenous mechanical perturbations on living tissues are commonly used to investigate whether cell effectors can respond to mechanical cues. However, in most of these experiments, the applied mechanical stress and/or the biological response are described only qualitatively. We developed a quantitative pipeline based on microindentation and image analysis to investigate the impact of a controlled and prolonged compression on microtubule behaviour in the Arabidopsis shoot apical meristem, using microtubule fluorescent marker lines. We found that a compressive stress, in the order of magnitude of turgor pressure, induced apparent microtubule bundling. Importantly, that response could be reversed several hours after the release of compression. Next, we tested the contribution of microtubule severing to compression-induced bundling: microtubule bundling seemed less pronounced in the katanin mutant, in which microtubule severing is dramatically reduced. Conversely, some microtubule bundles could still be observed 16 hours after the release of compression in the spiral2 mutant, in which severing rate is instead increased. To quantify the impact of mechanical stress on anisotropy and orientation of microtubule arrays, we used the nematic tensor based FibrilTool ImageJ/Fiji plugin. To assess the degree of apparent bundling of the network, we developed several methods, some of which were borrowed from geostatistics. The final microtubule bundling response could notably be related to tissue growth velocity that was recorded by the indenter during compression. Because both input and output are quantified, this pipeline is an initial step towards correlating more precisely the cytoskeleton response to mechanical stress in living tissues.

TAKING THE PERSPECTIVE OF A SELLER AND A BUYER: IMPLICATIONS FOR PRICE ELICITATION AND PRICE FRAMING

This dissertation consists of two essays which investigate how assuming the role of a seller or a buyer affects valuations in a price elicitation task (essay I) and how different presentations of an equivalent price affect evaluations when a consumer plays the dual roles of a buyer and a seller in transactions involving trade-ins (essay II). Sellers’ willingness to accept (WTA) to give up a good is typically higher than buyers' willingness to pay (WTP) to obtain the good. Essay I proposes that valuation processes of sellers and buyers are guided by a motivational orientation of “getting the best.” For a seller (buyer) indicating WTA (WTP), getting the best implies receiving as much as possible to give up a specific good (giving up as little as possible to get the specific good). Results of six studies suggest that the WTA-WTP elicitation task activates different directional goals, leading to the WTA-WTP disparity. The different directional goals lead sellers and buyers to focus on different aspects and bias their cognitive reasoning and interpretation of information. By connecting the valuation process to the general motivation of getting the best, this research provides a unifying framework to explain the disparate interpretations of the WTA-WTP disparity. Many new purchases and replacement decisions involve consumers’ trading in their old products. In such transactions, the overall exchange may be priced either as separate transactions (partitioned) with price tags for the payment and the receipt or as a single net price (consolidated) which takes into account the value of the trade-in. Essay II examines whether consumers prefer a partitioned price versus a consolidated price presentation. The findings suggest that when consumers are trading in a product which has a low value relative to the price of a new product, they prefer a consolidated price. In contrast, when trading in a product which has high value, they prefer a partitioned price. The results suggest that consumers use the price of the new product as an anchor to evaluate the trade-in value, and the perception of the trade-in value influences the overall evaluation especially when the transaction is partitioned.

MULTISCALE MODELING OF LIQUID CRYSTALLINE/NANOTUBE COMPOSITES

The objective of this research is to synthesize structural composites designed with particular areas defined with custom modulus, strength and toughness values in order to improve the overall mechanical behavior of the composite. Such composites are defined and referred to as 3D-designer composites. These composites will be formed from liquid crystalline polymers and carbon nanotubes. The fabrication process is a variation of rapid prototyping process, which is a layered, additive-manufacturing approach. Composites formed using this process can be custom designed by apt modeling methods for superior performance in advanced applications. The focus of this research is on enhancement of Young's modulus in order to make the final composite stiffer. Strength and toughness of the final composite with respect to various applications is also discussed. We have taken into consideration the mechanical properties of final composite at different fiber volume content as well as at different orientations and lengths of the fibers. The orientation of the LC monomers is supposed to be carried out using electric or magnetic fields. A computer program is modeled incorporating the Mori-Tanaka modeling scheme to generate the stiffness matrix of the final composite. The final properties are then deduced from the stiffness matrix using composite micromechanics. Eshelby's tensor, required to calculate the stiffness tensor using Mori-Tanaka method, is calculated using a numerical scheme that determines the components of the Eshelby's tensor (Gavazzi and Lagoudas 1990). The numerical integration is solved using Gaussian Quadrature scheme and is worked out using MATLAB as well. . MATLAB provides a good deal of commands and algorithms that can be used efficiently to elaborate the continuum of the formula to its extents. Graphs are plotted using different combinations of results and parameters involved in finding these results

Nuclear Magnetic Resonance Spectroscopic and Computational Investigations of Chloroperoxidase Catalyzed Regio- and Enantio-Selective Transformations

Chloroperoxidase (CPO) is the most versatile heme-containing enzyme that catalyzes a broad spectrum of reactions. The remarkable feature of this enzyme is the high regio- and enantio-selectivity exhibited in CPO-catalyzed oxidation reactions. The aim of this dissertation is to elucidate the structural basis for regio- and enantio-selective transformations and investigate the application of CPO in biodegradation of synthetic dyes. To unravel the mechanism of CPO-catalyzed regioselective oxidation of indole, the dissertation explored the structure of CPO-indole complex using paramagnetic relaxation and molecular modeling. The distances between the protons of indole and the heme iron revealed that the pyrrole ring of indole is oriented toward the heme with its 2-H pointing directly at the heme iron. This provides the first experimental and theoretical explanation for the "unexpected" regioselectivity of CPO-catalyzed indole oxidation. Furthermore, the residues including Leu 70, Phe 103, Ile 179, Val 182, Glu 183, and Phe 186 were found essential to the substrate binding to CPO. These results will serve as a lighthouse in guiding the design of CPO mutants with tailor-made activities for biotechnological applications. To understand the origin of the enantioselectivity of CPO-catalyzed oxidation reactions, the interactions of CPO with substrates such as 2-(methylthio)thiophene were investigated by nuclear magnetic resonance spectroscopy (NMR) and computational techniques. In particular, the enantioselectivity is partly explained by the binding orientation of substrates. In third facet of this dissertation, a green and efficient system for degradation of synthetic dyes was developed. Several commercial dyes such as orange G were tested in the CPO-H2O2-Cl- system, where degradation of these dyes was found very efficient. The presence of halide ions and acidic pH were found necessary to the decomposition of dyes. Significantly, the results revealed that this degradation of azo dyes involves a ferric hypochlorite intermediate of CPO (Fe-OCl), compound X.

Analysis of the failure of cracked biscuits

Cracks or checks in biscuits weaken the material and cause the product to break at low load levels that are perceived as injurious to product quality. In this work, the structural response of circular digestive biscuits, with diameter 72 mm and thickness 7.2 mm, simply supported around the circumference and loaded by a central concentrated force was investigated by experiment and theory. Tests were conducted to quantify the distribution in breakage strength for structurally sound biscuits, biscuits with natural checks and biscuits with a single known part-through crack. For sound biscuits the breakage force is Normally distributed with a mean of 12.5 N and standard deviation of 1.2 N. For biscuits with checks, the corresponding statistics are 9.6 N ± 2.62 N respectively. The presence of a crack weakens the biscuit and strength, as measured by breakage force falls almost linearly with crack length and crack depth. The orientation of the crack, whether radial or tangential, and its location (i.e. position of the crack mid-point on the biscuit surface) are also important. Deep, radial, cracks located close to the biscuit centre can reduce the strength by up to 50%. Two separate failure criteria were examined for sound and cracked biscuits respectively. The results from these tests were in good accord with theory. For a biscuit without defects, breakage occurred when maximum biscuit stress reached or exceeded the failure stress of 420 kPa. For a biscuit with cracks, breakage occurred as above or alternatively when its critical stress intensity factor of 18 kPam0.5 was reached.

Influence of kaolin addition on the dynamics of oxygen mass transport in polyvinyl alcohol dispersion coatings

The permeability of dispersion barriers produced from polyvinyl alcohol (PVOH) and kaolin clay blends coated onto polymeric supports has been studied by employing two different measurement methods: the oxygen transmission rate (OTR) and the ambient oxygen ingress rate (AOIR). Coatings with different thicknesses and kaolin contents were studied. Structural information of the dispersion-barrier coatings was obtained by Fourier transform infrared spectroscopy (FTIR) spectroscopy and scanning electron microscopy (SEM). These results showed that the kaolin content influences both the orientation of the kaolin and the degree of crystallinity of the PVOH coating. Increased kaolin content increased the alignment of the kaolin platelets to the basal plane of the coating. Higher kaolin content was accompanied by higher degree of crystallinity of the PVOH. The barrier thickness proved to be less important in the early stages of the mass transport process, whereas it had a significant influence on the steady-state permeability. The results from this study demonstrate the need for better understanding of how permeability is influenced by (chemical and physical) structure.

Near-infrared polarimetry in the galactic center. On properties of Sagittarius A* and a dusty S-cluster object

Near-infrared polarimetry observation is a powerful tool to study the central sources at the center of the Milky Way. My aim of this thesis is to analyze the polarized emission present in the central few light years of the Galactic Center region, in particular the non-thermal polarized emission of Sagittarius~A* (Sgr~A*), the electromagnetic manifestation of the super-massive black hole, and the polarized emission of an infrared-excess source in the literature referred to as DSO/G2. This source is in orbit about Sgr~A*. In this thesis I focus onto the Galactic Center observations at $\lambda=2.2~\mu m$ ($K_\mathrm{s}$-band) in polarimetry mode during several epochs from 2004 to 2012. The near-infrared polarized observations have been carried out using the adaptive optics instrument NAOS/CONICA and Wollaston prism at the Very Large Telescope of ESO (European Southern Observatory). Linear polarization at 2.2 $\mu m$, its flux statistics and time variation, can be used to constrain the physical conditions of the accretion process onto the central super-massive black hole. I present a statistical analysis of polarized $K_\mathrm{s}$-band emission from Sgr~A* and investigate the most comprehensive sample of near-infrared polarimetric light curves of this source up to now. I find several polarized flux excursions during the years and obtain an exponent of about 4 for the power-law fitted to polarized flux density distribution of fluxes above 5~mJy. Therefore, this distribution is closely linked to the single state power-law distribution of the total $K_\mathrm{s}$-band flux densities reported earlier by us. I find polarization degrees of the order of 20\%$\pm$10\% and a preferred polarization angle of $13^o\pm15^o$. Based on simulations of polarimetric measurements given the observed flux density and its uncertainty in orthogonal polarimetry channels, I find that the uncertainties of polarization parameters under a total flux density of $\sim 2\,{\mathrm{mJy}}$ are probably dominated by observational uncertainties. At higher flux densities there are intrinsic variations of polarization degree and angle within rather well constrained ranges. Since the emission is most likely due to optically thin synchrotron radiation, the obtained preferred polarization angle is very likely reflecting the intrinsic orientation of the Sgr~A* system i.e. an accretion disk or jet/wind scenario coupled to the super-massive black hole. Our polarization statistics show that Sgr~A* must be a stable system, both in terms of geometry, and the accretion process. I also investigate an infrared-excess source called G2 or Dusty S-cluster Object (DSO) moving on a highly eccentric orbit around the Galaxy's central black hole, Sgr~A*. I use for the first time the near-infrared polarimetric imaging data to determine the nature and the properties of DSO and obtain an improved $K_\mathrm{s}$-band identification of this source in median polarimetry images of different observing years. The source starts to deviate from the stellar confusion in 2008 data and it does not show a flux density variability based on our data set. Furthermore, I measure the polarization degree and angle of this source and conclude based on the simulations on polarization parameters that it is an intrinsically polarized source with a varying polarization angle as it approaches Sgr~A* position. I use the interpretation of the DSO polarimetry measurements to assess its possible properties.

Chiroptical spectroscopy of biomolecules using chiral plasmonic nanostructures

This thesis explores the potential of chiral plasmonic nanostructures for the ultrasensitive detection of protein structure. These nanostructures support the generation of fields with enhanced chirality relative to circularly polarised light and are an extremely incisive probe of protein structure. In chapter 4 we introduce a nanopatterned Au film (Templated Plasmonic Substrate, TPS) fabricated using a high through-put injection moulding technique which is a viable alternative to expensive lithographically fabricated nanostructures. The optical and chiroptical properties of TPS nanostructures are found to be highly dependent on the coupling between the electric and magnetic modes of the constituent solid and inverse structures. Significantly, refractive index based measurements of strongly coupled TPSs display a similar sensitivity to protein structure as previous lithographic nanostructures. We subsequently endeavour to improve the sensing properties of TPS nanostructures by developing a high through-put nanoscale chemical functionalisation technique. This process involves a chemical protection/deprotection strategy. The protection step generates a self-assembled monolayer (SAM) of a thermally responsive polymer on the TPS surface which inhibits protein binding. The deprotection step exploits the presence of nanolocalised thermal gradients in the water surrounding the TPS upon irradiation with an 8ns pulsed laser to modify the SAM conformation on surfaces with high net chirality. This allows binding of biomaterial in these regions and subsequently enhances the TPS sensitivity levels. In chapter 6 an alternative method for the detection of protein structure using TPS nanostructures is introduced. This technique relies on mediation of the electric/magnetic coupling in the TPS by the adsorbed protein. This phenomenon is probed through both linear reflectance and nonlinear second harmonic generation (SHG) measurements. Detection of protein structure using this method does not require the presence of fields of enhanced chirality whilst it is also sensitive to a larger array of secondary structure motifs than the measurements in chapters 4 and 5. Finally, a preliminary investigation into the detection of mesoscale biological structure is presented. Sensitivity to the mesoscale helical pitch of insulin amyloid fibrils is displayed through the asymmetry in the circular dichroism (CD) of lithographic gammadions of varying thickness upon adsorption of insulin amyloid fibril spherulites and fragmented fibrils. The proposed model for this sensitivity to the helical pitch relies on the vertical height of the nanostructures relative to this structural property as well as the binding orientation of the fibrils.

Detection of Glass pattern configurations by the primary visual cortex

The neurons in the primary visual cortex that respond to the orientation of visual stimuli were discovered in the late 1950s (Hubel, D.H. & Wiesel, T.N. 1959. J. Physiol. 148:574-591) but how they achieve this response is poorly understood. Recently, experiments have demonstrated that the visual cortex may use the image processing techniques of cross or auto-correlation to detect the streaks in random dot patterns (Barlow, H. & Berry, D.L. 2010. Proc. R. Soc. B. 278: 2069-2075). These experiments made use of sinusoidally modulated random dot patterns and of the so-called Glass patterns - where randomly positioned dot pairs are oriented in a parallel configuration (Glass, L. 1969. Nature. 223: 578-580). The image processing used by the visual cortex could be inferred from how the threshold of detection of these patterns in the presence of random noise varied as a function of the dot density in the patterns. In the present study, the detection thresholds have been measured for other types of patterns including circular, hyperbolic, spiral and radial Glass patterns and an indication of the type of image processing (cross or auto-correlation) by the visual cortex is presented. As a result, it is hoped that this study will contribute to an understanding of what David Marr called the ‘computational goal’ of the primary visual cortex (Marr, D. 1982. Vision: A Computational Investigation into the Human Representation and Processing of Visual Information. New York: Freeman.)