Dense Deformation Field Estimation for Atlas Registration using the Active Contour Framework

In this paper, we propose a new paradigm to carry outthe registration task with a dense deformation fieldderived from the optical flow model and the activecontour method. The proposed framework merges differenttasks such as segmentation, regularization, incorporationof prior knowledge and registration into a singleframework. The active contour model is at the core of ourframework even if it is used in a different way than thestandard approaches. Indeed, active contours are awell-known technique for image segmentation. Thistechnique consists in finding the curve which minimizesan energy functional designed to be minimal when thecurve has reached the object contours. That way, we getaccurate and smooth segmentation results. So far, theactive contour model has been used to segment objectslying in images from boundary-based, region-based orshape-based information. Our registration technique willprofit of all these families of active contours todetermine a dense deformation field defined on the wholeimage. A well-suited application of our model is theatlas registration in medical imaging which consists inautomatically delineating anatomical structures. Wepresent results on 2D synthetic images to show theperformances of our non rigid deformation field based ona natural registration term. We also present registrationresults on real 3D medical data with a large spaceoccupying tumor substantially deforming surroundingstructures, which constitutes a high challenging problem.

Minimum price variations, time priority and quotes dynamics

We analyze the impact of a minimum price variation (tick) and timepriority on the dynamics of quotes and the trading costs when competitionfor the order flow is dynamic. We find that convergence to competitiveoutcomes can take time and that the speed of convergence is influencedby the tick size, the priority rule and the characteristics of the orderarrival process. We show also that a zero minimum price variation is neveroptimal when competition for the order flow is dynamic. We compare thetrading outcomes with and without time priority. Time priority is shownto guarantee that uncompetitive spreads cannot be sustained over time.However it can sometimes result in higher trading costs. Empiricalimplications are proposed. In particular, we relate the size of thetrading costs to the frequency of new offers and the dynamics of theinside spread to the state of the book.

Numerical modeling of the effect of variation of boundary conditions on vadose zone hydraulic properties

An accurate estimation of hydraulic fluxes in the vadose zone is essential for the prediction of water, nutrient and contaminant transport in natural systems. The objective of this study was to simulate the effect of variation of boundary conditions on the estimation of hydraulic properties (i.e. water content, effective unsaturated hydraulic conductivity and hydraulic flux) in a one-dimensional unsaturated flow model domain. Unsaturated one-dimensional vertical water flow was simulated in a pure phase clay loam profile and in clay loam interlayered with silt loam distributed according to the third iteration of the Cantor Bar fractal object Simulations were performed using the numerical model Hydrus 1D. The upper and lower pressure heads were varied around average values of -55 cm for the near-saturation range. This resulted in combinations for the upper and lower constant head boundary conditions, respectively, of -50 and -60 cm, -40 and -70 cm, -30 and -80 cm, -20 and -90 cm, and -10 and -100 cm. For the drier range the average head between the upper and lower boundary conditions was set to -550 cm, resulting in the combinations -500 and -600 cm, -400 and -700 cm, -300 and -800 cm, -200 and -900 cm, and -100 and -1,000 cm, for upper and lower boundary conditions, respectively. There was an increase in water contents, fluxes and hydraulic conductivities with the increase in head difference between boundary conditions. Variation in boundary conditions in the pure phase and interlayered one-dimensional profiles caused significant deviations in fluxes, water contents and hydraulic conductivities compared to the simplest case (a head difference between the upper and lower constant head boundaries of 10 cm in the wetter range and 100 cm in the drier range).

n=1/4 domain-growth universality class: Crossover to the n=1/2 class

The kinetic domain-growth exponent is studied by Monte Carlo simulation as a function of temperature for a nonconserved order-parameter model. In the limit of zero temperature, the model belongs to the n=(1/4 slow-growth unversality class. This is indicative of a temporal pinning in the domain-boundary network of mixed-, zero-, and finite-curvature boundaries. At finite temperature the growth kinetics is found to cross over to the Allen-Cahn exponent n=(1/2. We obtain that the pinning time of the zero-curvature boundary decreases rapidly with increasing temperature.

Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities.

Quantifying the spatial configuration of hydraulic conductivity (K) in heterogeneous geological environments is essential for accurate predictions of contaminant transport, but is difficult because of the inherent limitations in resolution and coverage associated with traditional hydrological measurements. To address this issue, we consider crosshole and surface-based electrical resistivity geophysical measurements, collected in time during a saline tracer experiment. We use a Bayesian Markov-chain-Monte-Carlo (McMC) methodology to jointly invert the dynamic resistivity data, together with borehole tracer concentration data, to generate multiple posterior realizations of K that are consistent with all available information. We do this within a coupled inversion framework, whereby the geophysical and hydrological forward models are linked through an uncertain relationship between electrical resistivity and concentration. To minimize computational expense, a facies-based subsurface parameterization is developed. The Bayesian-McMC methodology allows us to explore the potential benefits of including the geophysical data into the inverse problem by examining their effect on our ability to identify fast flowpaths in the subsurface, and their impact on hydrological prediction uncertainty. Using a complex, geostatistically generated, two-dimensional numerical example representative of a fluvial environment, we demonstrate that flow model calibration is improved and prediction error is decreased when the electrical resistivity data are included. The worth of the geophysical data is found to be greatest for long spatial correlation lengths of subsurface heterogeneity with respect to wellbore separation, where flow and transport are largely controlled by highly connected flowpaths.

Approximate polytope ensemble for one-class classification

In this work, a new one-class classification ensemble strategy called approximate polytope ensemble is presented. The main contribution of the paper is threefold. First, the geometrical concept of convex hull is used to define the boundary of the target class defining the problem. Expansions and contractions of this geometrical structure are introduced in order to avoid over-fitting. Second, the decision whether a point belongs to the convex hull model in high dimensional spaces is approximated by means of random projections and an ensemble decision process. Finally, a tiling strategy is proposed in order to model non-convex structures. Experimental results show that the proposed strategy is significantly better than state of the art one-class classification methods on over 200 datasets.

Structural, metamorphic and geochronological relations between the Zanskar Shear Zone and the Miyar Shear Zone (NW Indian Himalaya): Evidence for two distinct tectonic structures and implications for the evolution of the High Himalayan Crystalline of Zanskar

The geologic structures and metamorphic zonation of the northwestern Indian Himalaya contrast significantly with those in the central and eastern parts of the range, where the high-grade metamorphic rocks of the High Himalayan Crystalline (HHC) thrust southward over the weakly metamorphosed sediments of the Lesser Himalaya along the Main Central Thrust (MCT). Indeed, the hanging wall of the MCT in the NW Himalaya mainly consists of the greenschist facies metasediments of the Chamba zone, whereas HHC high-grade rocks are exposed more internally in the range as a large-scale dome called the Gianbul dome. This Gianbul dome is bounded by two oppositely directed shear zones, the NE-dipping Zanskar Shear Zone (ZSZ) on the northern flank and the SW-dipping Miyar Shear Zone (MSZ) on the southern limb. Current models for the emplacement of the HHC in NW India as a dome structure differ mainly in terms of the roles played by both the ZSZ and the MSZ during the tectonothermal evolution of the HHC. In both the channel flow model and wedge extrusion model, the ZSZ acts as a backstop normal fault along which the high-grade metamorphic rocks of the HHC of Zanskar are exhumed. In contrast, the recently proposed tectonic wedging model argues that the ZSZ and the MSZ correspond to one single detachment system that operates as a subhorizontal backthrust off of the MCT. Thus, the kinematic evolution of the two shear zones, the ZSZ and the MSZ, and their structural, metamorphic and chronological relations appear to be diagnostic features for discriminating the different models. In this paper, structural, metamorphic and geochronological data demonstrate that the MSZ and the ZSZ experienced two distinct kinematic evolutions. As such, the data presented in this paper rule out the hypothesis that the MSZ and the ZSZ constitute one single detachment system, as postulated by the tectonic wedging model. Structural, metamorphic and geochronological data are used to present an alternative tectonic model for the large-scale doming in the NW Indian Himalaya involving early NE-directed tectonics, weakness in the upper crust, reduced erosion at the orogenic front and rapid exhumation along both the ZSZ and the MSZ.

Timing of Capacity Change: Models for Capital Intensive Industry

The suitable timing of capacity investments is a remarkable issue especially in capital intensive industries. Despite its importance, fairly few studies have been published on the topic. In the present study models for the timing of capacity change in capital intensive industry are developed. The study considers mainly the optimal timing of single capacity changes. The review of earlier research describes connections between cost, capacity and timing literature, and empirical examples are used to describe the starting point of the study and to test the developed models. The study includes four models, which describe the timing question from different perspectives. The first model, which minimizes unit costs, has been built for capacity expansion and replacement situations. It is shown that the optimal timing of an investment can be presented with the capacity and cost advantage ratios. After the unit cost minimization model the view is extended to the direction of profit maximization. The second model states that early investments are preferable if the change of fixed costs is small compared to the change of the contribution margin. The third model is a numerical discounted cash flow model, which emphasizes the roles of start-up time, capacity utilization rate and value of waiting as drivers of the profitable timing of a project. The last model expands the view from project level to company level and connects the flexibility of assets and cost structures to the timing problem. The main results of the research are the solutions of the models and analysis or simulations done with the models. The relevance and applicability of the results are verified by evaluating the logic of the models and by numerical cases.

Stability Evaluation of a Paper Machine Wet End

Työssä analysoidaanprosessin vaikutusta paperikoneen stabiiliuteen. Kaksi modernia sanomalehtipaperikonetta analysoitiin ja sen perusteella molemmista prosesseista rakennettiin fysiikan lakeihin perustuvat simulointimallit APROS Paper simulointiohjelmistolla. Työn tavoitteena on selvittää, miten kyseisten koneiden prosessit eroavat toisistaan ja arvioida, miten havaitut erot vaikuttavat prosessien stabiiliuteen. Työssä tarkastellaan periodisten häiriöiden vaimenemista prosessissa. Simuloinnissa herätteenä käytettiin puhdasta valkoista kohinaa, jonka avulla eri taajuistenperiodisten häiriöiden vaimenemista analysoitiin. Prosessien häiriövasteet esitetään taajuuskoordinaatistossa. Suurimmat erot prosessien välillä löytyivät viirakaivosta ja sen sekoitusdynamiikasta. Perinteisen viirakaivon todettiin muistuttavan käyttäytymiseltään sarjaan kytkettyjä ideaalisekoittimia, kun taas pienempitilavuuksisen fluumin todettiin käyttäytyvän lähes kuin putkiviive. Vaikka erotprosessitilavuudessa sekä viirakaivon sekoitusdynamiikassa olivat hyvin selkeät, havaittiin vain marginaalinen ero prosessin välillä periodisten häiriöiden vaimenemisessa, koska erot viiraretentiotasoissa vaikuttivat eniten simulointituloksia. Matalammalla viiraretentiolla operoivan paperikoneen todettiin vaimentavan tehokkaammin prosessihäiriöitä. Samalla retentiotasolla pienempitilavuuksisen prosessin todettiin vaimentavan hitaita prosessihäiriöitä marginaalisesti paremmin. Tutkituista paperikoneista toisella simuloitiin viiraosan vedenpoistomuutoksenvaikutusta viiraretentioon ja paperin koostumukseen. Lisäksi arvioitiin viiraretention säädön toimivuutta. Viiraosan listakengän vedenpoiston todettiin aiheuttavan merkittäviä sakeus- ja retentiohäiriöitä, mikäli sen avulla poistettavan kiintoaineen virtaus tuplaantuisi. Viiraretention säädön todettiin estävän häiriöiden kierron prosessissa, mutta siirtävän ne suoraan rainaan. Retention säädön eikuitenkaan todettu olevan suoranainen häiriön lähde.

Rahoitusalan yrityksen arvonmäärittäminen

Tässä tutkimuksessa tutkittiin rahoitusalan yritysten arvonmäärittämistä. Tutkimuksen päätavoitteena oli perehtyä yrityksen arvonmäärityksen teoriaan ja testata erilaisten arvonmääritysmallien sopivuutta esimerkkiyrityksiin. Tätä kautta pyrittiin selvittämään perinteisten arvonmääritysmallien soveltuvuus rahoitusalan yritysten arvon määrittämiseksi.Tämän lisäksi oltiin kiinnostuneita menetelmien soveltuvuuseroista erityyppisten rahoitusalan yritysten arvonmäärityksessä. Teoriaosuudessa esiteltiin yrityksen arvonmääritysprosessi ja käytiin läpi perinteiset arvonmääritysmallit. Empiriaosuudessa näistä teoreettisista lähtökohdista pyrittiin havainnollistamaan erilaisten esimerkkiyritysten avulla rahoitusalan yritysten arvonmääritysprosessia. Tämän tutkimuksen mukaanperinteiset arvonmääritysmallit tuntuvat sopivan pankkeja ja vakuutusyhtiöitä lukuun ottamatta kohtuullisen hyvin rahoitusalan yritysten arvonmäärittämiseen. Pankkien arvon määrittäminen on erittäin vaikeata johtuen lukuisista niihin kohdistuvista erityissäädöksistä ja raportointimääräyksistä.

Value creation and capital market valuation in the forest product industry

Tämä työ tähtää löytämään mahdollisia poikkeamia metsäteollisuusyritysten markkina-arvoissa ja tunnistaa tekijöitä jotka ovat vaikuttaneet pääomamarkkinoiden odotuksiin yritysten tulevaisuuden suorituskyvystä sekä yrityksen markkina-arvoon. Tämän työn päätavoitteena on kehittää diskontattuihin kassavirtoihin perustuva arvonmääritysmalli jolla mitataan metsäteollisuusyritysten tosiasiallista arvoa yritysten suorituskyvyn ja arvoajureiden perusteella. Lisäksi tavoitteena on löytää selittäviä tekijöitä havaituille eroille yritysten tosiasiallisesten arvojen ja markkinaperusteisten arvojen välisillä. Teoreettisessa osassa esitellään rahoitusteorian pääpiirteet arvonmäärityksen kannalta, aikasempia tutkimuksia sekä metsäteollisuuden toimialakohtaisia tekijöitä. Empiirisessä osassa kehittetään diskontattuihin kassavirtoihin perustuvaa arvonmääritymalli. Otos koostuu 32 suurimmasta Pohjoismaisesta ja Pohjoisamerikkalaisesta metsäteollisuusyrityksestä vuonna 2000. Tutkimuksen aikavälien 1991 -2000. Tulokset tukevat aikaisempia tutkimuksia jonka mukaan kasvuinvestoinnit eivät luo positiivisia odotuksia yrityksen tulevaisuuden kassavirroista. Tarkemmat löydöt ovat, että arvon luominen tutkimusajanjakson aikana ei vaikuttanut yhtä merkittävästi pääomamarkkinoiden odotuksiin yrityksen tulevaisuuden suorituskyvystä kuin mitatut tosiasialliset arvot. Tulokset viittaavat siihen, että metsäteollisuusyritysten markkina- arvot olivat keskimäärin riippuvaisempia itse yrityksestä, kuin sen toiminnasta.

Yrityksen kasvu ja omistaja-arvo - Diskontatun kassavirran periaatteella lasketun omistaja-arvon aikasarjapohjainen tarkastelu

Tutkielman päätavoitteena oli selvittää yrityksen kasvun vaikutusta omistaja-arvoon. Omistaja-arvon laskentaa selvitetään diskontatun kassavirran periaatteella. Omistaja-arvot esitettiin ja niiden kehitystä tutkittiin vuosittain aikasarjana. Kasvun vaikutuksia tutkittiin kolmella kannattavuustasolla ja kahdella pääoman tuottovaatimusprosentilla. Omistaja-arvon kehitystä verrattiin sijoitetun pääoman tuottoprosenttiin ja pääoman kiertonopeuteen. Tutkimusta varten rakennettiin yksinkertainen malliyritys, jonka rakennetta muuttamalla kyettiin tutkimaan kasvun vaikutuksia ja omistaja-arvoon vaikuttavia tekijöitä. Yrityksen kasvaessa kannattavuudella on suuri merkitys siihen, kuinka omistaja-arvo kehittyy. Tutkimus osoittaa, että yrityksen kasvu vaikuttaa omistaja-arvoa lisäävästi tiettyyn rajaan saakka, jonka jälkeen kasvu alkaa tuhota omistaja-arvoa. Tutkimuksessa löydettiin se liikevaihdon kasvun taso, jolla omistaja-arvo on korkeimmillaan. Tämä arvo ei riipu pääoman tuottovaatimustasosta. Alhaisempi pääoman tuottovaatimus luo suuremman omistaja-arvon kuin korkeampi tuottovaatimus. Tutkimus osoittaa, että vaikka kannattavuus sijoitetun pääoman tuottovaatimusprosentilla mitattuna paranee edelleen yrityksen omistaja-arvo alkaa tuhoutua jollakin tietyllä liikevaihdon kasvuntasolla. Pääoman kiertonopeus käyttäytyy sijoitetun pääoman tuottoprosentin kanssa samalla tavalla. Tutkimus osoittaa, että omistaja-arvon tuhoutuminen johtuu yrityksen tuottaman vapaan kassavirran pienentymisestä. Vapaa kassavirta alkaa pienentyä, kun kasvun vaatimat investoinnit ja käyttöpääoman lisäys kasvavat nopeammin suhteessa kasvun tuottoon verrattuna.

Yrityksen arvon määrittäminen

Yrityksen arvoa joudutaan tänä päivänä määrittämään monissa eri tilanteissa ja arvonmäärittäjinä voivat toimia niin yksityiset osakkeenomistajat kuin yhtiöiden ylimmät johtajat. Arvonmäärityksen yleistymisen taustalla ei ole kuitenkaan sen yksinkertaisuus ja helppokäyttöisyys, vaan sen strateginen hyöty käyttäjälleen. Lisäksi arvonmäärittämisen vaikeus ja monipuolisuus tekevät aiheesta entistäkin kiehtovamman. Tutkimuksen lähtökohtana oli selvittää käytetyimmät ja suomalaiseen käytäntöön parhaiten sopivat arvonmääritysmenetelmät. Työssä käydään läpi yrityksen arvonmäärityksen motiivit, arvonmäärityksen edellytykset sekä maailmalla paljon käytetyt arvonmääritysmenetelmät. Lisäksi johtopäätöksissä tulkitaan jokaista arvonmääritysmenetelmää sekä niiden soveltumista suomalaisten yhtiöiden arvostamiseen. Yrityksen arvonmäärittäminen on siis hyvin monipuolinen yrityksen taloudellista tilaa tutkiva ja ennustava menetelmä. Aiheesta on kirjoitettu paljon kotimaista ja ulkomaista kirjallisuutta, jotka käsittelevät aihetta erilaisista näkökulmista, mutta harvoissa teoksissa on annettu konkreettisia esimerkkejä monista eri arvonmääritysmenetelmistä sekä tutkittu niiden yhteyksiä analyyttisesti. Tässä kandidaatintyössä on käytetty kahta hyvin erilaista pörssiyhtiötä, Fortumia ja M-realia, esimerkkeinä ja pyritty selvittämään niiden avulla eri menetelmien käyttökelpoisuutta ja problematiikkaa yrityksen arvonmäärityksessä.

Reaction Kinetics and Viscosity Modelling in the Fusion Syntheses of Ca- and Ca/Mg-resinates

Rosin is a natural product from pine forests and it is used as a raw material in resinate syntheses. Resinates are polyvalent metal salts of rosin acids and especially Ca- and Ca/Mg- resinates find wide application in the printing ink industry. In this thesis, analytical methods were applied to increase general knowledge of resinate chemistry and the reaction kinetics was studied in order to model the non linear solution viscosity increase during resinate syntheses by the fusion method. Solution viscosity in toluene is an important quality factor for resinates to be used in printing inks. The concept of critical resinate concentration, c crit, was introduced to define an abrupt change in viscosity dependence on resinate concentration in the solution. The concept was then used to explain the non-inear solution viscosity increase during resinate syntheses. A semi empirical model with two estimated parameters was derived for the viscosity increase on the basis of apparent reaction kinetics. The model was used to control the viscosity and to predict the total reaction time of the resinate process. The kinetic data from the complex reaction media was obtained by acid value titration and by FTIR spectroscopic analyses using a conventional calibration method to measure the resinate concentration and the concentration of free rosin acids. A multivariate calibration method was successfully applied to make partial least square (PLS) models for monitoring acid value and solution viscosity in both mid-infrared (MIR) and near infrared (NIR) regions during the syntheses. The calibration models can be used for on line resinate process monitoring. In kinetic studies, two main reaction steps were observed during the syntheses. First a fast irreversible resination reaction occurs at 235 °C and then a slow thermal decarboxylation of rosin acids starts to take place at 265 °C. Rosin oil is formed during the decarboxylation reaction step causing significant mass loss as the rosin oil evaporates from the system while the viscosity increases to the target level. The mass balance of the syntheses was determined based on the resinate concentration increase during the decarboxylation reaction step. A mechanistic study of the decarboxylation reaction was based on the observation that resinate molecules are partly solvated by rosin acids during the syntheses. Different decarboxylation mechanisms were proposed for the free and solvating rosin acids. The deduced kinetic model supported the analytical data of the syntheses in a wide resinate concentration region, over a wide range of viscosity values and at different reaction temperatures. In addition, the application of the kinetic model to the modified resinate syntheses gave a good fit. A novel synthesis method with the addition of decarboxylated rosin (i.e. rosin oil) to the reaction mixture was introduced. The conversion of rosin acid to resinate was increased to the level necessary to obtain the target viscosity for the product at 235 °C. Due to a lower reaction temperature than in traditional fusion synthesis at 265 °C, thermal decarboxylation is avoided. As a consequence, the mass yield of the resinate syntheses can be increased from ca. 70% to almost 100% by recycling the added rosin oil.

Methods for assessing the sustainability of integrated municipal waste management and

The general striving to bring down the number of municipal landfills and to increase the reuse and recycling of waste-derived materials across the EU supports the debates concerning the feasibility and rationality of waste management systems. Substantial decrease in the volume and mass of landfill-disposed waste flows can be achieved by directing suitable waste fractions to energy recovery. Global fossil energy supplies are becoming more and more valuable and expensive energy sources for the mankind, and efforts to save fossil fuels have been made. Waste-derived fuels offer one potential partial solution to two different problems. First, waste that cannot be feasibly re-used or recycled is utilized in the energy conversion process according to EU’s Waste Hierarchy. Second, fossil fuels can be saved for other purposes than energy, mainly as transport fuels. This thesis presents the principles of assessing the most sustainable system solution for an integrated municipal waste management and energy system. The assessment process includes: · formation of a SISMan (Simple Integrated System Management) model of an integrated system including mass, energy and financial flows, and · formation of a MEFLO (Mass, Energy, Financial, Legislational, Other decisionsupport data) decision matrix according to the selected decision criteria, including essential and optional decision criteria. The methods are described and theoretical examples of the utilization of the methods are presented in the thesis. The assessment process involves the selection of different system alternatives (process alternatives for treatment of different waste fractions) and comparison between the alternatives. The first of the two novelty values of the utilization of the presented methods is the perspective selected for the formation of the SISMan model. Normally waste management and energy systems are operated separately according to the targets and principles set for each system. In the thesis the waste management and energy supply systems are considered as one larger integrated system with one primary target of serving the customers, i.e. citizens, as efficiently as possible in the spirit of sustainable development, including the following requirements: · reasonable overall costs, including waste management costs and energy costs; · minimum environmental burdens caused by the integrated waste management and energy system, taking into account the requirement above; and · social acceptance of the selected waste treatment and energy production methods. The integrated waste management and energy system is described by forming a SISMan model including three different flows of the system: energy, mass and financial flows. By defining the three types of flows for an integrated system, the selected factor results needed in the decision-making process of the selection of waste management treatment processes for different waste fractions can be calculated. The model and its results form a transparent description of the integrated system under discussion. The MEFLO decision matrix has been formed from the results of the SISMan model, combined with additional data, including e.g. environmental restrictions and regional aspects. System alternatives which do not meet the requirements set by legislation can be deleted from the comparisons before any closer numerical considerations. The second novelty value of this thesis is the three-level ranking method for combining the factor results of the MEFLO decision matrix. As a result of the MEFLO decision matrix, a transparent ranking of different system alternatives, including selection of treatment processes for different waste fractions, is achieved. SISMan and MEFLO are methods meant to be utilized in municipal decision-making processes concerning waste management and energy supply as simple, transparent and easyto- understand tools. The methods can be utilized in the assessment of existing systems, and particularly in the planning processes of future regional integrated systems. The principles of SISMan and MEFLO can be utilized also in other environments, where synergies of integrating two (or more) systems can be obtained. The SISMan flow model and the MEFLO decision matrix can be formed with or without any applicable commercial or free-of-charge tool/software. SISMan and MEFLO are not bound to any libraries or data-bases including process information, such as different emission data libraries utilized in life cycle assessments.