Estimating Soil Organic Carbon Stocks and Spatial Patterns with Statistical and GIS-Based Methods

Accurately quantifying soil organic carbon (SOC) is considered fundamental to studying soil quality, modeling the global carbon cycle, and assessing global climate change. This study evaluated the uncertainties caused by up-scaling of soil properties from the county scale to the provincial scale and from lower-level classification of Soil Species to Soil Group, using four methods: the mean, median, Soil Profile Statistics (SPS), and pedological professional knowledge based (PKB) methods. For the SPS method, SOC stock is calculated at the county scale by multiplying the mean SOC density value of each soil type in a county by its corresponding area. For the mean or median method, SOC density value of each soil type is calculated using provincial arithmetic mean or median. For the PKB method, SOC density value of each soil type is calculated at the county scale considering soil parent materials and spatial locations of all soil profiles. A newly constructed 1∶50,000 soil survey geographic database of Zhejiang Province, China, was used for evaluation. Results indicated that with soil classification levels up-scaling from Soil Species to Soil Group, the variation of estimated SOC stocks among different soil classification levels was obviously lower than that among different methods. The difference in the estimated SOC stocks among the four methods was lowest at the Soil Species level. The differences in SOC stocks among the mean, median, and PKB methods for different Soil Groups resulted from the differences in the procedure of aggregating soil profile properties to represent the attributes of one soil type. Compared with the other three estimation methods (i.e., the SPS, mean and median methods), the PKB method holds significant promise for characterizing spatial differences in SOC distribution because spatial locations of all soil profiles are considered during the aggregation procedure.     

Carbon Quality and Stocks in Organic Horizons in Boreal Forest Soils

We investigated the mechanisms that determine the quality and quantity of organic carbon (C) stocks in boreal forest soils by analyzing both qualitative and quantitative changes in the organic fractions in the soil organic matter (OM) in a vertical gradient in the decomposition continuum of the organic horizon [litter layer (L), fermentation layer (F), and humus layer (H)] in forest soils using a sequential fractionation method at two forest types along a climatic gradient in Finland. We predicted that the concentrations of water-soluble (WSE) and non-polar (NPE) extractives should decrease and those of the acid-soluble (AS) fraction and acid-insoluble residue (AIR) should increase from the L to the F, and from the F to the H layers, but the C/N ratio of soil OM should stay constant after reaching the critical quotient. We also predicted that the AIR concentrations should be higher in the south than north boreal, and in sub-xeric than mesic forests. Consistent with our hypothesis, the concentrations of WSE and NPE fractions decreased and concentrations of AIR increased in the vertical soil gradient. The highest concentrations of the AS fraction were found in the F layer. The C/N ratio was lowest in the F layer, and the highest in the H layer, indicating that soil OM is depleted in N in relation to C along the vertical soil gradient. Concentrations of WSE and NPE were lower, and concentrations of AIR were higher in the south than in north boreal forests, which is in agreement with our hypothesis that higher soil temperatures may enhance accumulation of slowly decomposable OM in the soil. The concentrations of AIR were higher in the sub-xeric than mesic forests. Contrary to our expectations, however, the differences in the chemical quality in soil OM between the site types were amplified from the L to the H layer. The size of the C storage was significantly larger in south than north boreal sites, and larger in the mesic than in the sub-xeric sites.     

Generation of Spatial Patterns in Boreal Forest Landscapes

Boreal forests are composed of a few plant species with contrasting traits with respect to ecosystem functioning and spatial patterning. Early successional deciduous species, such as birch and aspen, disperse seeds widely, do not tolerate low light and nitrogen availabilities, have rapidly decaying litter, and are highly preferred by herbivores. These later succeed to conifers, such as spruce and fir, which disperse seeds locally, tolerate low light levels and low nitrogen availability, have litter that decays slowly, and are unpalatable to most mammalian herbivores. Although there are also early successional conifers, such as jack pine and Scots pine, the aspen-birch-spruce-fir successional sequence is the most common over much of North America, and (without fir) in Fennoscandia and Siberia. The course of succession in these forests is controlled partly by seed dispersal and selective foraging by mammalian herbivores. Both of these processes are spatially dynamic, but little is known about how their spatial dynamics may affect ecosystem processes, such as nitrogen cycling or productivity. We present spatially explicit models that demonstrate the following: (a) Spatially explicit seed dispersal results in more clumped distribution of tree species and persistence of greater paper birch biomass than uniform seed rain across the landscape. Such results are consistent with current spatially explicit population models of dispersal and coexistence. (b) With localized seed dispersal, the concentrations of available soil nitrogen are distributed in larger patches with sharp transitions from low to high nitrogen availability near patch edges. In contrast, with a uniform seed rain, the distribution of soil nitrogen availability was more uniform and “hotspots” were more localized. Thus, the spatial pattern of an ecosystem process (nitrogen cycling) is determined by seed dispersal and competition for light among competing populations. (c) A dispersing herbivore, such as moose, that selectively forages on early successional deciduous species with high quality litter, such as aspen or birch, and discriminates against late successional conifers, such as spruce or fir, imposes higher-order repeated patterns of plant species and biomass distribution on the landscape. Thus, seed dispersal and herbivore foraging correlate properties in adjacent patches but in different ways, and different spatial patterns emerge. Other processes, such as insect outbreaks, fire, and water flow, also may correlate properties between adjacent patches and result in additional patterns. Received 8 February 1999; accepted 28 May 1999.     

Simulating Carbon Stocks and Fluxes of an African Tropical Montane Forest with an Individual-Based Forest Model

Tropical forests are carbon-dense and highly productive ecosystems. Consequently, they play an important role in the global carbon cycle. In the present study we used an individual-based forest model (FORMIND) to analyze the carbon balances of a tropical forest. The main processes of this model are tree growth, mortality, regeneration, and competition. Model parameters were calibrated using forest inventory data from a tropical forest at Mt. Kilimanjaro. The simulation results showed that the model successfully reproduces important characteristics of tropical forests (aboveground biomass, stem size distribution and leaf area index). The estimated aboveground biomass (385 t/ha) is comparable to biomass values in the Amazon and other tropical forests in Africa. The simulated forest reveals a gross primary production of 24 t_cha^-1yr^-1. Modeling above- and belowground carbon stocks, we analyzed the carbon balance of the investigated tropical forest. The simulated carbon balance of this old-growth forest is zero on average. This study provides an example of how forest models can be used in combination with forest inventory data to investigate forest structure and local carbon balances.     

Patterns of Land-use Abandonment Control Tree-recruitment and Forest Dynamics in Mediterranean Mountains

Abstract Mediterranean ecosystems have been impacted for millennia by human practices, particularly agricultural and pastoral activities. Since the middle of the nineteenth century, land-use abandonment has lead to scrubland and forest expansion, especially in mountain areas of the northern Mediterranean basin. This study aimed at analyzing how grazing history affects subsequent forest dynamics at a site located in the limestone foothills of the Southern Alps (France). The approach combines archival documents and dendroecology to investigate the origin, establishment and development of forest following land-use abandonment. Scots pine (Pinus sylvestris) started to colonize quickly in the 1870s, with the recruitment rate increasing during the first decade of the 1900s, associated with a decline of the local human population and regional livestock. Since the 1960s, European beech (Fagus sylvatica) and silver fir (Abies alba) have regenerated in the understorey of Scots pines. Regeneration is controlled by a threshold of grazing pressure. Noticeably, the rate of reforestation differs according to the former land-use, with pastures being colonized more quickly than ploughed areas. Different previous land-uses leading to different times of grazing cessation, combined with variable herbaceous competition explain the contrasting micro-scale regeneration patterns. Agricultural land-use and abandonment are both significant driving forces of vegetation dynamics. Knowledge of these factors is thus necessary to understand present patterns and to predict future forest pathways in the Mediterranean mountains.     

Temporal-Spatial Pattern of Carbon Stocks in Forest Ecosystems in Shaanxi,Northwest China

The precise and accurate quantitative evaluation of the temporal and spatial pattern of carbon (C) storage in forest ecosystems is critical for understanding the role of forests in the global terrestrial C cycle and is essential for formulating forest management policies to combat climate change. In this study, we examined the C dynamics of forest ecosystems in Shaanxi, northwest China, based on four forest inventories (1989–1993, 1994–1998, 1999–2003, and 2004–2008) and field-sampling measurements (2012). The results indicate that the total C storage of forest ecosystems in Shaanxi increased by approximately 29.3%, from 611.72 Tg in 1993 to 790.75 Tg in 2008, partially as a result of ecological restoration projects. The spatial pattern of C storage in forest ecosystems mainly exhibited a latitude-zonal distribution across the province, increasing from north (high latitude) to south (low latitude) generally, which signifies the effect of environmental conditions, chiefly water and heat related factors, on forest growth and C sequestration. In addition, different data sources and estimation methods had a significant effect on the results obtained, with the C stocks in 2008 being considerably overestimated (864.55 Tg) and slightly underestimated (778.07 Tg) when measured using the mean C density method and integrated method, respectively. Overall, our results demonstrated that the forest ecosystem in Shaanxi acted as a C sink over the last few decades. However, further studies should be carried out with a focus on adaption of plants to environmental factors along with forest management for vegetation restoration to maximize the C sequestration potential and to better cope with climate change.     

Seasonal and Spatial Patterns in Ground Foraging Ants in a Rain Forest in the Western Ghats,India1

Seasonal variation and spatial distribution in ground foraging rain forest ants were studied in South Kannada–Kodagu District in Karnataka (India) between 1990 and 1991 by pit-fall trap sampling. All ant species showed marked seasonality. A total of 31 species were recorded from the primary forest over a period of two years. More species were recorded from the closed canopy forest than from tree fall gaps in primary forest. All ant species showed marked seasonality with fewer species and individuals sampled/plot during the wetter seasons. The numerically dominant species, Pheidole sp., was markedly lower in abundance during the wet seasons. Spatial patterns were also studied during a dry season both in the primary forest and an adjacent logged forest. More species were recorded from the logged forest than the primary forest. Community composition in primary forest was different from that in logged forest. Common species were more ubiquitous than rarer species. Species were distributed bimodally across sampling plots. Probable underlying processes behind these seasonal and spatial patterns have been discussed.     

Carbon Stocks in an African Woodland Landscape: Spatial Distributions and Scales of Variation

Current knowledge of Africa’s carbon (C) pools is limited despite its importance in the global C budget. To increase the understanding of C stocks in African woodlands, we asked how C stocks in soil and vegetation vary across a miombo woodland landscape and to what degree and at what scales are these stocks linked? We sampled along a 5-km transect using a cyclic sampling scheme to allow geostatistical analyses. Soil C stocks in the top 5?cm (12.1?±?0.6?Mg?C?ha^?1 (±?SE)) and 30?cm depths (40.1?±?2.5?Mg?C?ha^?1) varied significantly at scales of a few meters (autocorrelation distance 14?m in 0–5-cm and 26?m in 0–30-cm interval), and aboveground (AG) woody C stocks (20.7?±?1.8?Mg?C?ha^?1) varied significantly at kilometer scales (1,426?m). Soil textural distributions were linked to topography (r ²?=?0.54) as were large-tree AG C stocks (r ²?=?0.70). AG C stocks were constrained to an upper boundary by soil texture with greater AG C being associated with coarser textured soils. Vegetation and soil C stocks were coupled in the landscape in the top 5?cm of soil (r ²?=?0.24) but not with deeper soil C stocks, which were coupled to soil clay content (r ²?=?0.38). This study is one of the most complete transect studies in an African miombo woodland, and suggests that C stock distributions are strongly linked to topography and soil texture. To optimize sampling strategies for C stock assessments in miombo, soil C should be sampled at more than 26?m apart, and AG C should be sampled at more than 1,426?m apart in plots larger than 0.5?ha.     

Spatial Variability and Stocks of Soil Organic Carbon in the Gobi Desert of Northwestern China

Soil organic carbon (SOC) plays an important role in improving soil properties and the C global cycle. Limited attention, though, has been given to assessing the spatial patterns and stocks of SOC in desert ecosystems. In this study, we quantitatively evaluated the spatial variability of SOC and its influencing factors and estimated SOC storage in a region (40 km²) of the Gobi desert. SOC exhibited a log-normal depth distribution with means of 1.6, 1.5, 1.4, and 1.4 g kg⁻¹ for the 0–10, 10–20, 20–30, and 30–40 cm layers, respectively, and was moderately variable according to the coefficients of variation (37–42%). Variability of SOC increased as the sampling area expanded and could be well parameterized as a power function of the sampling area. Significant correlations were detected between SOC and soil physical properties, i.e. stone, sand, silt, and clay contents and soil bulk density. The relatively coarse fractions, i.e. sand, silt, and stone contents, had the largest effects on SOC variability. Experimental semivariograms of SOC were best fitted by exponential models. Nugget-to-sill ratios indicated a strong spatial dependence for SOC concentrations at all depths in the study area. The surface layer (0–10 cm) had the largest spatial dependency compared with the other layers. The mapping revealed a decreasing trend of SOC concentrations from south to north across this region of the Gobi desert, with higher levels close to an oasis and lower levels surrounded by mountains and near the desert. SOC density to depths of 20 and 40 cm for this 40 km² area was estimated at 0.42 and 0.68 kg C m⁻², respectively. This study provides an important contribution to understanding the role of the Gobi desert in the global carbon cycle.     

Significance of Over-Mature and Decaying Trees for Carbon Stocks in a Central European Natural Spruce Forest

Old-growth forests are important stores for carbon as they may accumulate C for centuries. The alteration of biomass and soil carbon pools across the development stages of a forest dynamics cycle has rarely been quantified. We studied the above- and belowground C stocks in the five forest development stages (regeneration to decay stage) of a montane spruce (Picea abies) forest of the northern German Harz Mountains, one of Central Europe’s few forests where the natural forest dynamics have not been disturbed by man for several centuries. The over-mature and decay stages had the largest total (up to 480 Mg C ha^?1) and aboveground biomass carbon pools (200 Mg C ha^?1) with biomass C stored in dead wood in the decay stage. The soil C pool (220–275 Mg C ha^?1, 0–60 cm) was two to three times larger than in temperate lowland spruce forests and remained invariant across the forest dynamics cycle. On the landscape level, taking into account the frequency of the five forest development stages, the total carbon pool was approximately 420 Mg C ha^?1. The results evidence the high significance of over-mature and decaying stages of temperate mountain forests not only for conserving specialized forest organisms but also for their large carbon storage potential.     

Landscape-Scale Controls on Aboveground Forest Carbon Stocks on the Osa Peninsula,Costa Rica

Tropical forests store large amounts of carbon in tree biomass, although the environmental controls on forest carbon stocks remain poorly resolved. Emerging airborne remote sensing techniques offer a powerful approach to understand how aboveground carbon density (ACD) varies across tropical landscapes. In this study, we evaluate the accuracy of the Carnegie Airborne Observatory (CAO) Light Detection and Ranging (LiDAR) system to detect top-of-canopy tree height (TCH) and ACD across the Osa Peninsula, Costa Rica. LiDAR and field-estimated TCH and ACD were highly correlated across a wide range of forest ages and types. Top-of-canopy height (TCH) reached 67 m, and ACD surpassed 225 Mg C ha^-1, indicating both that airborne CAO LiDAR-based estimates of ACD are accurate in tall, high-biomass forests and that the Osa Peninsula harbors some of the most carbon-rich forests in the Neotropics. We also examined the relative influence of lithologic, topoedaphic and climatic factors on regional patterns in ACD, which are known to influence ACD by regulating forest productivity and turnover. Analyses revealed a spatially nested set of factors controlling ACD patterns, with geologic variation explaining up to 16% of the mapped ACD variation at the regional scale, while local variation in topographic slope explained an additional 18%. Lithologic and topoedaphic factors also explained more ACD variation at 30-m than at 100-m spatial resolution, suggesting that environmental filtering depends on the spatial scale of terrain variation. Our result indicate that patterns in ACD are partially controlled by spatial variation in geologic history and geomorphic processes underpinning topographic diversity across landscapes. ACD also exhibited spatial autocorrelation, which may reflect biological processes that influence ACD, such as the assembly of species or phenotypes across the landscape, but additional research is needed to resolve how abiotic and biotic factors contribute to ACD variation across high biomass, high diversity tropical landscapes.     

Changes in Biomass Carbon and Soil Organic Carbon Stocks following the Conversion from a Secondary Coniferous Forest to a Pine Plantation

The objectives of this study were to estimate changes of tree carbon (C) and soil organic carbon (SOC) stock following a conversion in land use, an issue that has been only insufficiently addressed. For this study, we examined a chronosequence of 2 to 54-year-old Pinus kesiya var. langbianensis plantations that replaced the original secondary coniferous forest (SCF) in Southwest China due to clearing. C stocks considered here consisted of tree, understory, litter, and SOC (0–1 m). The results showed that tree C stocks ranged from 0.02±0.001 Mg C ha^-1 to 141.43±5.29 Mg C ha^-1, and increased gradually with the stand age. Accumulation of tree C stocks occurred in 20 years after reforestaion and C stock level recoverd to SCF. The maximum of understory C stock was found in a 5-year-old stand (6.74±0.7 Mg C ha^-1) with 5.8 times that of SCF, thereafter, understory C stock decreased with the growth of plantation. Litter C stock had no difference excluding effects of prescribed burning. Tree C stock exhibited a significant decline in the 2, 5-year-old stand following the conversion to plantation, but later, increased until a steady state-level in the 20, 26-year-old stand. The SOC stocks ranged from 81.08±10.13 Mg C ha^-1 to 160.38±17.96 Mg C ha^-1. Reforestation significantly decreased SOC stocks of plantation in the 2-year-old stand which lost 42.29 Mg C ha^-1 in the 1 m soil depth compared with SCF by reason of soil disturbance from sites preparation, but then subsequently recovered to SCF level. SOC stocks of SCF had no significant difference with other plantation. The surface profile (0–0.1 m) contained s higher SOC stocks than deeper soil depth. C stock associated with tree biomass represented a higher proportion than SOC stocks as stand development proceeded.     

关帝山次生杨桦林种群结构与立木的空间点格局

选择关帝山典型天然次生杨桦林样地进行种群结构分析,并运用空间点格局分析方法O-ring统计系统分析了次生杨桦林群落内不同尺度下种群空间分布格局及种间关联性.结果表明,在个体组成上,演替先锋树种白桦(Betula platyphylla)、红桦(B.albo-sinensis)和山杨(Populus davidiana)是目前次生杨桦林样地的优势种.但从种群径级结构看,白桦、红桦和山杨天然更新不良,而云杉(Picea spp.)和华北落叶松(Larix principis-rupprechtii)幼苗呈现良好的更新状态.随着演替的进行,云杉、落叶松将逐渐进入林冠层并取代白桦、红桦和山杨成为该森林的优势种;从种群空间分布格局来看,红桦、白桦表现为小尺度上明显聚集,华北落叶松、皂柳(Salix wallichiana)表现为小尺度上强聚集性分布格局.种群空间分布格局受空间尺度、群落结构及种群径级结构综合作用的影响;空间关联性分析表明,红桦与华北落叶松、红桦与皂柳、皂柳与华北落叶松在小尺度上正相关,其它树种间没有表现出明显关联性.     

Holocene Carbon Stocks and Carbon Accumulation Rates Altered in Soils Undergoing Permafrost Thaw

Permafrost soils are a significant global store of carbon (C) with the potential to become a large C source to the atmosphere. Climate change is causing permafrost to thaw, which can affect primary production and decomposition, therefore affecting ecosystem C balance. To understand future responses of permafrost soils to climate change, we inventoried current soil C stocks, investigated ∆¹⁴C, C:N, δ¹³C, and δ¹⁵N depth profiles, modeled soil C accumulation rates, and calculated decadal net ecosystem production (NEP) in subarctic tundra soils undergoing minimal, moderate, and extensive permafrost thaw near Eight Mile Lake (EML) in Healy, Alaska. We modeled decadal and millennial soil C inputs, decomposition constants, and C accumulation rates by plotting cumulative C inventories against C ages based on radiocarbon dating of surface and deep soils, respectively. Soil C stocks at EML were substantial, over 50 kg C m⁻² in the top meter, and did not differ much among sites. Carbon to nitrogen ratio, δ¹³C, and δ¹⁵N depth profiles indicated most of the decomposition occurred within the organic soil horizon and practically ceased in deeper, frozen horizons. The average C accumulation rate for EML surface soils was 25.8 g C m⁻² y⁻¹ and the rate for the deep soil accumulation was 2.3 g C m⁻² y⁻¹, indicating these systems have been C sinks throughout the Holocene. Decadal net ecosystem production averaged 14.4 g C m⁻² y⁻¹. However, the shape of decadal C accumulation curves, combined with recent annual NEP measurements, indicates soil C accumulation has halted and the ecosystem may be becoming a C source. Thus, the net impact of climate warming on tundra ecosystem C balance includes not only becoming a C source but also the loss of C uptake capacity these systems have provided over the past ten thousand years.     

Effects of Land-Use Change on Carbon Stocks in Switzerland

We assessed how consequences of future land-use change may affect size and spatial shifts of C stocks under three potential trends in policy—(a) business-as-usual: continuation of land-use trends observed during the past 15 years; (b) extensification: full extensification of open-land; and (c) liberalization: full reforestation potential. The build-up times for the three scenarios are estimated at 30, 80 and 100 years, respectively. Potential C-stock change rates are derived from the literature. Whereas the business-as-usual scenario would cause marginal changes of 0.5%, liberalization would provoke a 13% increase in C stocks (+62 MtC). Gains of 24% would be expected for forests (+95 MtC), whereas open-land C stock would decrease 27% (−33 MtC). Extensification would lead to a C stock decrease of 3% (−12 MtC). Whereas forest C is expected to increase 12% (+36.5 MtC) at high elevations, stocks of open-land C would decline 38.5% (−48.5 MtC). Most affected are unfavorable grasslands, which increase in area (+59%) but contribute only 14.5% to the C stocks. C sinks would amount to 0.6 MtC y⁻¹ assuming a build-up time of 100 years for the liberalization scenario. C stocks on the current forest area are increasing by 1 MtC y⁻¹. The maximal total C sink of 1.6 MtC might thus suffice to compensate for agricultural greenhouse gases (2004: 1.4 Mt CO₂–C equivalents), but corresponds only to 11–13% of the anthropogenic greenhouse gas emission in Switzerland. Thus, even the largest of the expected terrestrial C stocks under liberalization will be small in comparison with current emissions of anthropogenic greenhouse gases.     

Understanding the Influences of Spatial Patterns on N Availability Within the Brazilian Amazon Forest

Nitrogen variations at different spatial scales and integrated across functional groups were addressed for lowland tropical forests in the Brazilian Amazon as follows: (1) how does N availability vary across the region over different spatial scales (regional × landscape scale); (2) how are these variations in N availability integrated across plant functional groups (legume × non-legume trees). Leaf N, P, and Ca concentrations as well the leaf N isotope ratios (δ¹⁵N) from a large set of legume and non-legume tree species were measured. Legumes had higher foliar N/Ca ratios than non-legumes, consistent with the high energetic costs in plant growth associated with higher foliar P/Ca ratios found in legumes than in non-legumes. At the regional level, foliar δ¹⁵N decreased with increasing rainfall. At the landscape level, N availability was higher in the forests on clayey soils on the plateau than in forests on sandier soils. The isotope as well as the non-isotope data relationships here documented, explain to a large extent the variation in δ¹⁵N signatures across gradients of rainfall and soil. Although at the regional level, the precipitation regime is a major determinant of differences in N availability, at the landscape level, under the same precipitation regime, soil type seems to be a major factor influencing the availability of N in the Brazilian Amazon forest.     

Personality and Situation Predictors of Consistent Eating Patterns

Introduction

A consistent eating style might be beneficial to avoid overeating in a food-rich environment. Eating consistency entails maintaining a similar dietary pattern across different eating situations. This construct is relatively under-studied, but the available evidence suggests that eating consistency supports successful weight maintenance and decreases risk for metabolic syndrome and cardiovascular disease. Yet, personality and situation predictors of consistency have not been studied.

Methods

A community-based sample of 164 women completed various personality tests, and 139 of them also reported their eating behaviour 6 times/day over 10 observational days. We focused on observations with meals (breakfast, lunch, or dinner). The participants indicated if their momentary eating patterns were consistent with their own baseline eating patterns in terms of healthiness or size of the meal. Further, participants described various characteristics of each eating situation.

Results

Eating consistency was positively predicted by trait self-control. Eating consistency was undermined by eating in the evening, eating with others, eating away from home, having consumed alcohol and having undertaken physical exercise. Interactions emerged between personality traits and situations, including punishment sensitivity, restraint, physical activity and alcohol consumption.

Conclusion

Trait self-control and several eating situation variables were related to eating consistency. These findings provide a starting point for targeting interventions to improve consistency, suggesting that a focus on self-control skills, together with addressing contextual factors such as social situations and time of day, may be most promising. This work is a first step to provide people with the tools they need to maintain a consistently healthy lifestyle in a food-rich environment.