Establishing a sampling protocol to estimate tobacco specific nitrosamines in growers bales: Results from the first year of a two year study

Similar documents
RFT SG Report. (Agrochemical Residue Field Trials SG) Berlin, Germany October 13 th 2016

PRESENTED BY HOP GROWERS OF AMERICA

Calculation of Predicted Environmental Concentrations in Soil -

2017 CORESTA JOINT STUDY GROUPS MEETINGS

Hydronic Systems Balance

Growth performance evaluation of genetically improved silver barb (Barbonymus gonionotus Bleeker) in different agro-ecological zones in Bangladesh

Competitive Performance of Elite Olympic-Distance Triathletes: Reliability and Smallest Worthwhile Enhancement

m v = 1.04 x 10-4 m 2 /kn, C v = 1.29 x 10-2 cm 2 /min

INTERSECTION SAFETY RESEARCH IN CRACOW UNIVERSITY OF TECHNOLOGY

Regulatory Mechanisms of Turkish Tobacco Sector

E. Agu, M. Kasperski Ruhr-University Bochum Department of Civil and Environmental Engineering Sciences

Standardized CPUE of Indian Albacore caught by Taiwanese longliners from 1980 to 2014 with simultaneous nominal CPUE portion from observer data

Agilent Dimension Software for ELSD User Manual

Chapter 2: Modeling Distributions of Data

CORESTA GUIDE N 12. May 2013 CONTROLLED ATMOSPHERE PARAMETERS FOR THE CONTROL OF CIGARETTE BEETLE AND TOBACCO MOTH

All AQA Unit 1 Questions Higher

A Cost Effective and Efficient Way to Assess Trail Conditions: A New Sampling Approach

Modeling Approaches to Increase the Efficiency of Clear-Point- Based Solubility Characterization

WindProspector TM Lockheed Martin Corporation

Case Study: How Misinterpreting Probabilities Can Cost You the Game. Kurt W. Rotthoff Seton Hall University Stillman School of Business

The Summer of 2007: A Look at Niagara

Modeling Fantasy Football Quarterbacks

Highway & Transportation (I) ECIV 4333 Chapter (4): Traffic Engineering Studies. Spot Speed

UK Trade Statistics 2017

Draft. Hiroki Yokoi, Yasuko Semba, Keisuke Satoh, Tom Nishida

The Aftermath of Global Financial Crises

TO INFINITY AND BEYOND? CAN RETAIL NON WOVENS KEEP GROWING? DANIEL GRIMSEY, RESEARCH ANALYST

More than half the world lives on less than $2 a day

Finite Element Analysis of an Aluminium Bike Frame

Confidence Intervals with proportions

QT2 - Coach s Reference

USING A PROBABILITY APPROACH TO RANK IN-LINE INSPECTION ANOMALIES FOR EXCAVATION AND FOR SETTING REINSPECTION INTERVALS

Saithe (Pollachius virens) in subareas 4 and 6, and in Division 3.a (North Sea, Rockall and West of Scotland, Skagerrak and Kattegat)

UK Trade Statistics 2016

Lab 3 Introduction to Quantitative Analysis: Pumps and Measurements of Flow

Is Ground Rolling For You? Jodi DeJong-Hughes Regional Extension Educator, Marshall cell

Q5. If x = 19, find (a) 1 (b) 2 (c) 5 (d) 4

Elements for the Biennial Review And Assessment of the World Timber Situation SELECTED DATA AND ANALYSIS.

An early warning system to predict house price bubbles

CORESTA RECOMMENDED METHOD N 6

WIND DATA REPORT. Bourne Water District

Smart Water Application Technologies (SWAT)

WATER OIL RELATIVE PERMEABILITY COMPARATIVE STUDY: STEADY VERSUS UNSTEADY STATE

The effect of the Gravity-Loading Countermeasure Skinsuit (GLCS) upon aerobic exercise performance

The Evolution of Transport Planning

Analysis of Variance. Copyright 2014 Pearson Education, Inc.

Requirements and Methods for Calibration of Equipment in Clinical Laboratories. NK Chan (HKRCBTS) 24 August 2007

Chapter 12 Practice Test

Update of trade weights data underlying the EERs and HCIs

API 10A Cooperative Testing Report Austin, TX January 2017

This Time is Different: Eight Centuries of Financial Folly

Cumulative Frequency Diagrams

Robust Imputation of Missing Values in Compositional Data Using the R-Package robcompositions

Cod (Gadus morhua) in subdivisions 24 32, eastern Baltic stock (eastern Baltic Sea) *

Geotechnical Engineering Laboratory CIVIL ENGINEERING VIRTUAL LABORATORY Experiment no 8 Standard Penetration Test

The Herzliya Indices. National Security Balance The Civilian Quantitative Dimension. Herzliya Conference Prof. Rafi Melnick, IDC Herzliya

WIND DATA REPORT. Ragged Mt Maine

AP 11.1 Notes WEB.notebook March 25, 2014

Forecasting sports tournaments by ratings of (prob)abilities

Prediction Market and Parimutuel Mechanism

Stat 139 Homework 3 Solutions, Spring 2015

Monthly Digest February 2016 No. 2016/04. Copyrights Statistics Botswana 2016

Fundamentals of Machine Learning for Predictive Data Analytics

COMPLETING THE RESULTS OF THE 2013 BOSTON MARATHON

Feeding fish to fish is this a responsible practice?

One-factor ANOVA by example

Commissioning an IMRT System for MLC Delivery. Gary A. Ezzell., Ph.D. Mayo Clinic Scottsdale

Bicycle Produced by IAR Team Focus Technology Co., Ltd

4-3 Rate of Change and Slope. Warm Up. 1. Find the x- and y-intercepts of 2x 5y = 20. Describe the correlation shown by the scatter plot. 2.

PGA Tour Scores as a Gaussian Random Variable

INSTITUTE AND FACULTY OF ACTUARIES. Curriculum 2019 AUDIT TRAIL

1. The data below gives the eye colors of 20 students in a Statistics class. Make a frequency table for the data.

Introduction to Scientific Notation & Significant Figures. Packet #6

CEE 345, Part 2, Winter 2012, Final Exam Solutions (Open Channel Flow)

Global Economic Briefing: Industrial Production

WIND DATA REPORT. Swan s Island, ME

Announcements. Unit 7: Multiple Linear Regression Lecture 3: Case Study. From last lab. Predicting income


ISyE 6414 Regression Analysis

Climatic Parameters Variability & Climate Change in Lebanon

WMO LABORATORY INTERCOMPARISON OF RAIN INTENSITY GAUGES

Measuring Relative Achievements: Percentile rank and Percentile point

Straw CO 2 Leak Testing

Optimizing Compressed Air Storage for Energy Efficiency

Individual Behavior and Beliefs in Parimutuel Betting Markets

Full file at

To Illuminate or Not to Illuminate: Roadway Lighting as It Affects Traffic Safety at Intersections

Lab 4: Transpiration

ECO 745: Theory of International Economics. Jack Rossbach Fall Lecture 6

NOBENDEM. Advanced Decompression Tables

Process Simulator Evaluates Blower and Valve Control Strategies for WWTP Aeration

Alternate Suction to Reduce Prolonged Air Leak

Global growth forecasts Key countries/regions,

A preliminary analysis of in-depth accident data for powered two-wheelers and bicycles in Europe

GUIDE TO THE ARCTIC SHIPPING CALCULATION TOOL

6. EXPERIMENTAL METHOD. A primary result of the current research effort is the design of an experimental

Market Correlations: Copper

Chart Collection for Morning Briefing

Vitamin C Market Research Report- Forecast to 2023

Unconventional Monetary Policy: Thoughts on the U.S. and Japanese Experiences

Transcription:

Establishing a sampling protocol to estimate tobacco specific nitrosamines in growers bales: Results from the first year of a two year study Colin Fisher 1 & Kristen McQuerry 2 1 Plant & Soil Science, University of Kentucky 2 Applied Statistics Laboratory, University of Kentucky 2016 CORESTA

Introduction Every attempt to reduce TSNA at all stages of production and manufacture Growing and curing conditions impact TSNA in farmer packages at sale. Point of sale is first opportunity for buyer to sample cured leaf Sampling may be required anytime before redrying.

Grower marketing packages Country Max wt. cu. M kg Indonesia 0.14 35 US 0.16 41 Poland 0.19 40 Brazil 0.19 55 France 0.20 25 Philippines 0.23 23 Thailand 0.28 60 Guatemala 0.32 62 Zimbabwe 0.32 90 Indonesia 0.32 100 Thailand 0.43 70 US (big) 1.16 230/300

CORESTA Sub-Group TSNA in Air-Cured and Fire-Cured Tobacco Draft Protocol for sample handling of post-cure tobacco (a.k.a. farmer bales) included: Section 1aii: 1. Collect 2 samples per 25 to 40 kg bale 2. Open bale and collect 2 samples of at least 0.7 kg per sample from center of bale. 3. For larger bales (250 kg), core sampler can be used. Collect at least 1.4 kg of leaf per bale (at least 14 cores if core is 100 g).

Quebec 2014 Aspects of protocol queried: 1. Total amount sample from each bale (1.4 kg) 2. Number samples from each size of bale Small bales - 2 grab samples of whole leaf Big bales 14 core samples 3. Effect of inclusion of midrib, especially in cored samples

Objective: Clarify sampling protocol to optimize sample size number of samples to best represent TSNA content of whole bale Year 1: Characterize variability of TSNA within individual bales

Assumptions and Limitations: Industry most interested in lamina chemistry, rather than whole leaf or midrib Number of bales tested limited by total number of samples for chemistry analysis Only burley in this study Only NNN addressed in this presentation Lamina+Midrib data = weighted concentration

Methods Small bales - 25-100 kg 6 bales Big bales - 250 kg 6 bales

Sampling Methods: Big Bales Three 45 cm cores diagonally across top and side of each bale (2 cores from same position if one core too little sample)

Sampling Methods: Small Bales Three methods: 1. Core Tip ¾ leaf length Mid-leaf Midrib ¼ leaf length x 4

Sampling Methods: Small Bales 2. Grab Sample (protocol) 30 40 leaves from each of 4 depths of each bale

Sampling Methods: Small Bales 3. Perfect Sample + 160 individual leaves randomly selected from throughout bale. then randomly divided into 4 samples of 40 leaves each

Sample preparation: 6 big bales X top & side X 3 cores each X lamina & midrib = 72 samples 6 small bales X 20 samples/bale X lamina & midrib = 240 samples

Big Bales - Variability of NNN in Lamina + Midrib between samples from same bale NNN, µg g- 1 30 25 20 15 10 5 0 Whole Leaf (weighted Lamina + Midrib) 1 2 3 4 5 6 Bale # Max Mean Min SD

NNN, ug/g NNN, ug/g Variability of NNN in Lamina and Midrib between samples from same bale 25 20 15 10 5 0 Lamina Max Mean Min SD 1 2 3 4 5 6 Bale # 80 60 40 20 0 Midrib Max Mean Min SD 1 2 3 4 5 6 Bale # Therefore: log transformed to normalize data for statistical analysis

Big Bales: Sampling from Side vs. Top

Ln single sample Lamina+Midrib Ln single sample Lamina+Midrib Big Bales: Sampling from Side vs. Top 3.5 Ln Lamina+Midrib Single SIDE Sample vs bale mean side samples (n=3) 3.5 Ln Lamina+Midrib Single TOP Sample vs. bale mean top samples (n=3) 3.0 3.0 2.5 2.0 1.5 1.0 0.5 Bottom Side R² = 0.984 R² = 0.955 Top R² = 0.991 0.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 Ln Mean Lamina+Midrib 2.5 2.0 1.5 1.0 0.5 Centre Left Right R² = 0.860 R² = 0.858 R² = 0.890 0.0 0.0 1.0 2.0 3.0 4.0 Ln Mean Lamina+Midrib

Lamina:Lam+Mid Ratio Big Bales: Ratio of Lamina NNN: Midrib NNN Bale # Lamina + Midrib Lamina Lamina/ Lam + Mid 1 21.9 18.4 0.84 2 0.8 0.6 0.78 3 8.6 6.6 0.77 4 12.9 8.7 0.67 5 5.4 4.4 0.81 6 5.4 4.4 0.82 1.0 0.8 0.6 0.4 0.2 0.0 Ratio of Lamina NNN:Midrib NNN 1 2 3 4 5 6 Bale #

NNN Lamina+Midrib Small Bales: Variability of NNN between samples in perfect samples 20 15 10 5 0 Mean Max Min SD 1 2 3 Bale # 4 5 6

Small Bale Sampling Method Burton, Dye, Bush. 1992. Distribution of tobacco constituents in tobacco leaf tissue. 1. Tobacco-specific nitrosamines, nitrate, nitrite, and alkaloids. J. Agric. Food Chem., 40 (6):1050 1055 0.8 NNN µg g -1 0.6 0.4 0.2 0.0 0 1 2 3 4 5 6 7 8 9 10 Tip Segment # Butt

Proportion of contribution to total NNN Fisher, Pearce. 2013. Assessing yield loss from hail in burley tobacco. CORESTA Agro-Phyto Meeting, Italy. NNN concentrations along length of leaf 30 25 20 15 10 5 0 0 2 4 6 8 10 Burton Segment # (1 = tip, 10 = butt) Lamina Midrib Sum L+M Straight laid small bale core samples: NNN of lamina + midrib at half way along length of leaf Lamina only NNN

Small Bales: NNN by sample method & position across all bales Sample Method & Position Indiv Lvs Perfect Sample Lamina NNN Midrib Lamina + Midrib 3.2 13.1 5.5 Hand Pull 3.5 10.9 5.3 Core Butt 5.2 13.8 8.7 Core Mid 3.0 8.0 3.9 Core Tip 2.5 7.3 2.9

Results of Single Sample NNN Position of Lamina+Midrib predicting Lamina of Individual Leaves The Rep the single sample originated Order of Significance per Single Sample Position p-value 1 1 Mid <.0001 2 Butt 0.0018 3 Tip 0.0033 4 Hand Pull 0.0505

Results of Single Sample NNN Position of Lamina+Midrib predicting Lamina of Individual Leaves The Rep the single sample originated Order of Significance per Single Sample Position p-value 1 1 Mid <.0001 2 Butt 0.0018 3 Tip 0.0033 4 Hand Pull 0.0505 2 1 Mid 0.0030 2 Tip 0.0030 3 Butt 0.0115 4 Hand Pull 0.0492

Results of Single Sample NNN Position of Lamina+Midrib predicting Lamina of Individual Leaves The Rep the single sample originated Order of Significance per Single Sample Position p-value 1 1 Mid <.0001 2 Butt 0.0018 3 Tip 0.0033 4 Hand Pull 0.0505 2 1 Mid 0.0030 2 Tip 0.0030 3 Butt 0.0115 4 Hand Pull 0.0492 3 1 Mid 0.0002 2 Tip 0.0004 3 Butt 0.0005 4 Hand Pull 0.0005

Results of Single Sample Position of Lamina+Midrib predicting Lamina of Individual Leaves The Rep the single sample originated Order of Significance per Single Sample Position p-value 1 1 Mid <.0001 2 Butt 0.0018 3 Tip 0.0033 4 Hand Pull 0.0505 2 1 Mid 0.0030 2 Tip 0.0030 3 Butt 0.0115 4 Hand Pull 0.0492 3 1 Mid 0.0002 2 Tip 0.0004 3 Butt 0.0005 4 Hand Pull 0.0005 4 1 Hand Pull 0.0007 2 Tip 0.0014 3 Mid 0.0024 4 Butt 0.0067

Predicting Whole Leaf NNN from single samples of Mid-leaf core Lamina+Midrib NNN Position Frequency Percent Hand Pull 4 16.67 Butt 4 16.67 Mid 12 50.00 Tip 4 16.67 Total Possible Iterations 24

Both bales sizes: Small data set n=6 Conclusions Large variability of NNN within single bale Gross variability diminishes as mean NNN concentration decreases Data normalized with log-transformation At approx. 3 µg g -1, range +1 2 µg g -1 Probably acceptable?

Big bales: Conclusions No possibility avoiding butts in sample NNN of lamina estimated from 0.8 x (lamina + midrib)? If not, lamina should be separated from midrib before analysis if NNN of lamina only required Single core could be sufficient, but more to reduce error

Small bales: Conclusions NNN concentration of mid-leaf core samples NNN concentration of whole leaf lamina One or two samples sufficient Sample through most of bale depth 2016 Study: 3 cores from mid-leaf and either side of mid-leaf

Acknowledgments: Sample Bales Financial Support Alliance One International British American Tobacco Philip Morris International Altria Client Services Kentucky Tobacco Research and Development Center RJ Reynolds