Modeling coupled interactions of carbon,water, and ozone exchange between terrestrial ecosystems and the atmosphere. I: model description

A new biophysical model (FORFLUX) is presented to study the simultaneous exchange of ozone, carbon dioxide, and water vapor between terrestrial ecosystems and the atmosphere. The model mechanistically couples all major processes controlling ecosystem flows trace gases and water implementing recent concepts in plant eco-physiology, micrometeorology, and soil hydrology. FORFLUX consists of four interconnected modules-a leaf photosynthesis model, a canopy flux model, a soil heat-, water- and CO2- transport model, and a snow pack model. Photosynthesis, water-vapor flux and ozone uptake at the leaf level are computed by the LEAFC3 sub-model. The canopy module scales leaf responses to a stand level by numerical integration of the LEAFC3model over canopy leaf area index (LAI). The integration takes into account (1) radiative transfer inside the canopy, (2) variation of foliage photosynthetic capacity with canopy depth, (3) wind speed attenuation throughout the canopy, and (4) rainfall interception by foliage elements. The soil module uses principles of the diffusion theory to predict temperature and moisture dynamics within the soil column, evaporation, and CO2 efflux from soil. The effect of soil heterogeneity on field-scale fluxes is simulated employing the Bresler-Dagan stochastic concept. The accumulation and melt of snow on the ground is predicted using an explicit energy balance approach. Ozone deposition is modeled as a sum of three fluxes- ozone uptake via plant stomata, deposition to non-transpiring plant surfaces, and ozone flux into the ground. All biophysical interactions are computed hourly while model projections are made at either hourly or daily time step. FORFLUX represents a comprehensive approach to studying ozone deposition and its link to carbon and water cycles in terrestrial ecosystems.     

Carbon dioxide fluxes over a grazed prairie and seeded pasture in the Northern Great Plains

Temperate grasslands are vast terrestrial ecosystems that may be an important component of the global carbon (C) cycle; however, annual C flux data for these grasslands are limited. The Bowen ratio/energy balance (BREB) technique was used to measure CO2 fluxes over a grazed mixed-grass prairie and a seeded western wheatgrass [Pascopyrum smithii (Rybd) L?ve] site at Mandan, ND from 24 April to 26 October in 1996, 1997, and 1998. Above-ground biomass and leaf area index (LAI) were measured about every 21 days throughout the season. Root biomass and soil organic C and N content were determined to 110 cm depth in selected increments about mid-July each year. Peak above-ground biomass and LAI coincided with peak fluxes and occurred between mid-July to early August. Biomass averaged 1227 and 1726 kg ha(-1) and LAI 0.44 and 0.59, for prairie and western wheatgrass, respectively. Average CO2 flux for the growing season was 279 g CO2 m(-2) for prairie and 218 g CO2 m(-2) for western wheatgrass (positive flux is CO2 uptake and negative flux is CO2 loss to the atmosphere). Using prior measured dormant season CO2 fluxes from the prairie sites gave annual flux estimates that ranged from -131 to 128 g CO2 m(-2) for western wheatgrass and from -70 to 189 g CO2 m(-2) for the prairie. This wide range in calculated annual fluxes suggests that additional research is required concerning dormant season flux measurements to obtain accurate estimates of annual CO2 fluxes. These results suggest Northern Great Plains mixed-grass prairie grasslands can either be a sink or a source for atmospheric CO2 or near equilibrium, depending on the magnitude of the dormant season flux.     

Modeling above-canopy CO2 flux and evapotranspiration in wheat

Simulations of above-canopy water vapor and CO2 fluxes were calculated by the USGF linked model of canopy gas exchange and subsurface processes for the 1996-1997 winter wheat season at the AmeriFlux Wheat Site, Oklahoma. Soil surface CO2 flux plus canopy gas exchange and transpiration plus soil evaporation modeled the CO2 and water vapor fluxes, respectively. Parameter values for net photosynthesis, respiration and transpiration were obtained from published sources, generated from Wheat Site data, or estimated by minimizing standard deviation between model and data. The mean measured downward flux of CO2 during rapid growth and maturity of the crop was -0.45 mg m(-2) s(-1) compared to simulated flux of -0.47. Simulated downward CO2 flux exceeded measured values during rapid growth of the crop but underestimated the flux during maturity. For the entire 285-day period, the mean measured upward CO2 flux at night was 0.06 and simulated flux was 0.05.     

Seasonal air-water exchange fluxes of polychlorinated biphenyls in the Hudson River Estuary

Polychlorinated biphenyls (PCBs) were measured in the air and water over the Hudson River Estuary during six intensive field campaigns from December 1999 to April 2001. Over-water gas-phase SigmaPCB concentrations averaged 1100 pg/m3 and varied with temperature. Dissolved-phase SigmaPCB concentrations averaged 1100 pg/L and displayed no seasonal trend. Uncertainty analysis of the results suggests that PCBs with 5 or fewer chlorines exhibited net volatilization. The direction of net air/water exchange could not be determined for PCBs with 6 or more chlorines. Instantaneous net fluxes of SigmaPCBs ranged from +0.2 to +630 ng m(-2) d(-1). Annual fluxes of SigmaPCBs were predicted from modeled gas-phase concentrations, measured dissolved-phase concentrations, daily surface water temperatures and wind speeds. The net volatilization flux was +62 microg m(-2) yr(-1), corresponding to an annual loss of +28 kg/yr of SigmaPCBs from the Hudson River Estuary for the year of 2000.     

Atmospheric deposition of polychlorinated dibenzo-p-dioxins, polychlorinated dibenzofurans, and dioxin-like polychlorinated biphenyls in the Kanto Region, Japan

The atmospheric bulk (dry and wet) deposition of dioxins was investigated at four locations (Tokyo, Yokohama, Tsukuba, and Tanzawa) in the Kanto region (in Japan) over one year using a stainless-steel pot. Annual average polychlorinated dibenzo-p-dioxins/polychlorinated dibenzofurans (PCDD/PCDF) deposition fluxes were estimated to be from 450 to 1300 ng/m2/yr, and the annual average TEQ fluxes from 5.7 to 17 ng-TEQ/m2/yr at the four locations. The PCDD/PCDF deposition flux was higher in winter than in summer. The deposition flux could be related to ambient temperature, particularly for less chlorinated PCDDs/PCDFs, while the deposition flux is not necessarily related to the amount of precipitation. The PCDD/PCDF deposition flux increased as the particle deposition flux increased, for the winter samples. Based on the ratio of the PCDD/PCDF deposition fluxes to the particle deposition fluxes, the contribution of the reentrainment of soil particles to the TEQ of PCDD/PCDF deposition was considered to be negligible in this region. Based on the air concentrations monitored near our deposition sampling points by the municipalities, the ratio of the annual deposition flux to the annual average air concentration was roughly estimated to be 0.082 cm/s. The range of deposition flux in the Kanto region was estimated to be from 1.5 to 31 (median: 9.8) ng-TEQ/m2/yr based on the range of air concentration data measured by the municipalities. The total annual deposition flux in the entire Kanto region was estimated to range from 50 to 900 g-TEQ/yr (median 320 g-TEQ/yr). This estimated flux was of the same order as the sum of estimated emissions from municipal solid waste incinerators and industrial waste incinerators in the Kanto region. The contributions of dioxin-like PCBs in Yokohama, Tsukuba, and Tanzawa depositions were less than 10% of the total TEQ; however, in Tokyo it was almost equal to or more than 50%.     

Summer and autumn ozone fluxes to a forest in the Czech Republic Brdy Mountains

In an effort to examine ozone (O3) deposition over a forest site in the Czech Republic, a low cost eddy flux experiment using slow response ozone and temperature sensors was implemented in July 1993 within the Brdy Mountains. Half-hour 2-Hz ozone and sensible heat measurements made at the Brdy Mountains for 98 days during the period 7 July 1994-20 October 1994 are analyzed and reported. While the Czech Brdy Mountains AOT40 level for the overall 104 day period was 7.6 ppm h (15.1 ppm h for the full 24-h summation), indicating a slight potential for 03 injury, the 1994 summer to autumn'measured forest O3 uptake was 2.4 (+/- 0.9) g m(-2), not unusually high compared to other studies. Average summer midday 03 fluxes and depositidn velocities were -1.0 (+/- 0.6) microg m(-2) s(-1) and 1.1 (+/- 0.7) cm s(-1). and autumn values were -0.36 (+/- 0.4) microg m(-2) s(-1) and 0.7 (+/- 0.5) cm s(-1) respectively. A unique contribution of this study is the first time demonstrated use of slow responding sensors for eddy covariance flux measurements at heights of 20 m above a forest.     

Assessment of ozone effects on nitrate export from Hubbard Brook Watershed 6

The impact of the air pollution ozone on soil N dynamics and temporal and spatial patterns of streamflow nitrate flux at the Hubbard Brook Experimental Forest Watershed 6 during the 1964-1994 period was assessed using aggregated (one-cell) and spatially explicit (208-cell) versions of the SImple NItrogen Cycle (SINIC) model. Simulated ozone effects included reductions in stomatal conductance and plant N demand. Model uncertainty was evaluated using Monte Carlo simulations. Ambient ozone was estimated to cause an additional 0.042 gN/m2 per year of nitrate export, 12% of the mean annual streamflow nitrate flux. The 95% credible interval of this estimate was 0.002-0.083 gN/m2 per year, or 0.72-27.3% of the annual flux. The large uncertainty in this estimate suggests that it may be difficult to identify ozone effects on nitrate export utilizing long term data from a single site.     

Carbon budgets for catchments across a managed landscape mosaic in southeast Sweden: contributing to the safety assessment of a nuclear waste repository

Ecosystem budgets of matter contribute to the assessment of transport and accumulation of bioavailable contaminants in a landscape, since flows of matter and energy ultimately determine the rates at which contaminants will be partitioned in the environment. This study compares ecosystem properties, such as net primary production (NPP), sequestration of matter and fluxes to food sources for humans, which are of potential interest to describe fluxes and accumulation of bioavailable radionuclides in 14 catchments within a larger catchment area in southeast Sweden. The carbon budgets, used as a proxy for organic matter, are mainly based on local estimates of pools and fluxes, which have been distributed across a landscape mosaic of different vegetation types and management regimes using a geographical information system (GIS). NPP varied by a factor close to two (432 - 709 g x Cx m(-2)x y(-1)), while net ecosystem production ranged between -124 and 159 gx C x m(-2) x y (-1) for the different catchments. Carbon sequestration mainly occurred in the vegetation while the soil organic carbon pool was mainly a source of carbon. Large herbivores consumed on average 4.5 % of the above-ground green tissue production. When arable land was present in the catchment, the flux of carbon to humans was highest from crops and, in decreasing order, milk and beef, followed by the flux from hunting and berry/fungus picking. The results can be used to estimate the potential assimilation of radionuclides in vegetation and the potential exposure to humans of bioavailable radionuclides.     

Monitoring the multi-year carbon balance of a subarctic palsa mire with micrometeorological techniques

This article reports a dataset on 8 years of monitoring carbon fluxes in a subarctic palsa mire based on micrometeorological eddy covariance measurements. The mire is a complex with wet minerotrophic areas and elevated dry palsa as well as intermediate sub-ecosystems. The measurements document primarily the emission originating from the wet parts of the mire dominated by a rather homogenous cover of Eriophorum angustifolium. The CO(2)/CH(4) flux measurements performed during the years 2001-2008 showed that the areas represented in the measurements were a relatively stable sink of carbon with an average annual rate of uptake amounting to on average -46 g C m(-2) y(-1) including an equally stable loss through CH(4) emissions (18-22 g CH(4)-C m(-2) y(-1)). This consistent carbon sink combined with substantial CH(4) emissions is most likely what is to be expected as the permafrost under palsa mires degrades in response to climate warming.     

Sedimentation and seasonal variation of hexachlorocyclohexanes in sediments in a eutrophic lake, China

Hexachlorocyclohexane (HCH) concentrations in sediments and sediment trap fluxes of particulate organic carbon and HCHs were measured bi-weekly from March 31 to October 18, 2006 in an urban eutrophic lake in Tianjin, China, in order to investigate sedimentation and seasonal variation of HCHs in sediments. HCH concentrations (dry weight basis) ranged from 2.2 to 20.2 ng/g (mean 7.7 ng/g) in surface sediments and from 26.6 to 972.7 ng/g (mean 187.0 ng/g) in settling particles, respectively. A clear seasonal variation in HCH sedimentation and HCH concentrations in sediments was observed. The maximal HCH deposition occurred following a spring phytoplankton bloom. The average flux of HCHs to sediment was approximately 21-fold higher in April to mid-June as compared to late June to October. This was attributed to the high vertical fluxes at the end of the spring phytoplankton bloom. The maximum values of HCH concentrations in sediments were observed in mid-June to late July. Concentrations of HCHs in sediments from the eutrophic lake were well-correlated with organic carbon contents in sediments. The annual sediment trap flux of HCHs in the eutrophic lake, which was estimated using data obtained in the eutrophic lake, was 117 microg/m2 yr, about 72% of which was attributed to the sedimentation corresponding to spring bloom phytoplankton deposition in late May to mid-June. The high sediment trap flux of HCHs in the eutrophic lake was related to serious local contamination.     

Modelling annual pasture dynamics: Application to stomatal ozone deposition

Modelling ozone (O₃) deposition for impact risk assessment is still poorly developed for herbaceous vegetation, particularly for Mediterranean annual pastures. High inter-annual climatic variability in the Mediterranean area makes it difficult to develop models characterizing gas exchange behaviour and air pollutant absorption suitable for risk assessment. This paper presents a new model to estimate stomatal conductance (g_s) of Trifolium subterraneum, a characteristic species of dehesa pastures. The MEDPAS (MEDiterranean PAStures) model couples 3 modules estimating soil water content (SWC), vegetation growth and g_s. The g_s module is a reparameterized version of the stomatal component of the EMEP DO₃SE O₃ deposition model. The MEDPAS model was applied to two contrasting years representing typical dry and humid springs respectively and with different O₃ exposures. The MEDPAS model reproduced realistically the g_s seasonal and inter-annual variations observed in the field. SWC was identified as the major driver of differences across years. Despite the higher O₃ exposure in the dry year, meteorological conditions favoured 2.1 times higher g_s and 56 day longer growing season in the humid year compared to the dry year. This resulted in higher ozone fluxes absorbed by T. subterraneum in the humid year. High inter-family variability was found in gas exchange rates, therefore limiting the relevance of single species O₃ deposition flux modelling for dehesa pastures. Stomatal conductance dynamics at the canopy level need to be considered for more accurate O₃ flux modelling for present and future climate scenarios in the Mediterranean area.     

Transpiration of gaseous elemental mercury through vegetation in a subtropical wetland in Florida

Four seasonal sampling campaigns were carried out in the Florida Everglades to measure elemental Hg vapor (Hg°) fluxes over emergent macrophytes using a modified Bowen ratio gradient approach. The predominant flux of Hg° over both invasive cattail and native sawgrass stands was emission; mean day time fluxes over cattail ranged from ∼20 (winter) to ∼40 (summer) ng m⁻² h⁻¹. Sawgrass fluxes were about half those over cattail during comparable periods. Emission from vegetation significantly exceeded evasion of Hg° from the underlying water surface (∼1–2 ng m⁻² h⁻¹) measured simultaneously using floating chambers. Among several environmental factors (e.g. CO₂ flux, water vapor flux, wind speed, water, air and leaf temperature, and solar radiation), water vapor exhibited the strongest correlation with Hg° flux, and transpiration is suggested as an appropriate term to describe this phenomenon. The lack of significant Hg° emissions from a live, but uprooted (floating) cattail stand suggests that a likely source of the transpired Hg° is the underlying sediments. The pattern of Hg° fluxes typically measured indicated a diel cycle with two peaks, possibly related to different gas exchange dynamics: one in early morning related to lacunal gas release, and a second at midday related to transpiration; nighttime fluxes approached zero.     

Regional estimation of current and future forest biomass

The 90,674 wildland fires that burned 2.9 million ha at an estimated suppression cost of $1.6 billion in the United States during the 2000 fire season demonstrated that forest fuel loading has become a hazard to life, property, and ecosystem health as a result of past fire exclusion policies and practices. The fire regime at any given location in these regions is a result of complex interactions between forest biomass, topography, ignitions, and weather. Forest structure and biomass are important aspects in determining current and future fire regimes. Efforts to quantify live and dead forest biomass at the local to regional scale has been hindered by the uncertainty surrounding the measurement and modeling of forest ecosystem processes and fluxes. The interaction of elevated CO2 with climate, soil nutrients, and other forest management factors that affect forest growth and fuel loading will play a major role in determining future forest stand growth and the distribution of species across the southern United States. The use of satellite image analysis has been tested for timely and accurate measurement of spatially explicit land use change and is well suited for use in inventory and monitoring of forest carbon. The incorporation of Landsat Thematic Mapper data coupled with a physiologically based productivity model (PnET), soil water holding capacity, and historic and projected climatic data provides an opportunity to enhance field plot based forest inventory and monitoring methodologies. We use periodic forest inventory data from the USDA Forest Service's Forest Inventory and Analysis (FIA) project to obtain estimates of forest area and type to generate estimates of carbon storage for evergreen, deciduous, and mixed forest classes for use in an assessment of remotely sensed forest cover at the regional scale for the southern United States. The displays of net primary productivity (NPP) generated from the PnET model show areas of high and low forest carbon storage potential and their spatial relationship to other landscape features for the southern United States. At the regional scale, predicted annual NPP in 1992 ranged from 836 to 2181 g/m2/year for evergreen forests and 769-2634 g/m2/year for deciduous forests with a regional mean for all forest land of 1448 g/m2/year. Prediction of annual NPP in 2050 ranged from 913 to 2076 g/m2/year for evergreen forest types to 1214-2376 g/m2/year for deciduous forest types with a regional mean for all forest land of 1659 g/m2/year. The changes in forest productivity from 1992 to 2050 are shown to display potential areas of increased or decreased forest biomass. This methodology addresses the need for spatially quantifying forest carbon in the terrestrial biosphere to assess forest productivity and wildland fire fuels.     

The exchange of ozone between vegetation and atmosphere: micrometeorological measurement techniques and models

The European critical levels (CLs) to protect vegetation are expressed as an accumulative exposure over a threshold of 40 ppb (nl l(-1)). In view of the fact that these chamber-derived CLs are based on ozone (O(3)) concentrations at the top of the canopy the correct application to ambient conditions presupposes the application of Soil-Vegetation-Atmosphere-Transfer (SVAT) models for quantifying trace gas exchange between phytosphere and atmosphere. Especially in the context of establishing control strategies based on flux-oriented dose-response relationships, O(3) flux measurements and O(3) exchange simulations are needed for representative ecosystems. During the last decades several micrometeorological methods for quantifying energy and trace gas exchange were developed, as well as models for the simulation of the exchange of trace gases between phytosphere and atmosphere near the ground. This paper is a synthesis of observational and modeling techniques which discusses measurement methods, assumptions, and limitations and current modeling approaches. Because stomatal resistance for trace gas exchange is parameterized as a function of water vapor or carbon dioxide (CO(2)) exchange, the most important micrometeorological techniques especially for quantifying O(3), water vapor and CO(2) flux densities are discussed. A comparison of simulated and measured O(3) flux densities shows good agreement in the mean.     

Atmospheric deposition of polycyclic aromatic hydrocarbons to Izmit Bay, Turkey

Wet deposition and dry deposition samples were collected in an urban/industrialized area of Izmit Bay, North-eastern Marmara Sea, Turkey, from September 2002 to July 2003. The samples were analyzed for sixteen polycyclic aromatic hydrocarbon (PAH) compounds by using HPLC-UV technique. Wet and dry deposition concentrations and fluxes of PAHs were determined. The results showed that PAH concentrations were high because of industrial processes, heavy traffic and residential areas next to the sampling site. Total dry deposition flux of the fifteen 3-6 ring PAHs was 8.30 microg m(-2)day(-1), with a range of 0.034-1.77 microg m(-2)day(-1). The total wet deposition flux of the fifteen 3-6 ring PAHs was 1716 microg m(-2) 11 month(-1), with a range of 10-440 microg m(-2) 11 month(-1). Significant seasonal differences were observed in both types of deposition samples. The winter fluxes of total PAHs were 1.5 and 2.5 times greater than those of the warm period for wet and dry deposition samples, respectively. Factor analysis of dry deposition samples and back trajectory analysis of wet deposition samples were also used to characterize and identify the PAH emission sources in this study.     

Quantifying ozone uptake at the canopy level of spruce,pine and larch trees at the alpine timberline: an approach based on sap flow measurement

Micro-climatic and ambient ozone data were combined with measurements of sap flow through tree trunks in order to estimate whole-tree ozone uptake of adult Norway spruce (Picea abies), cembran pine (Pinus cembra), and European larch (Larix decidua) trees. Sap flow was monitored by means of the heat balance approach in two trees of each species during the growing season of 1998. In trees making up the stand canopy, the ozone uptake by evergreen foliages was significantly higher than by deciduous ones, when scaled to the ground area. However, if expressed per unit of whole-tree foliage area, ozone flux through the stomata into the needle mesophyll was 1.09, 1.18 and 1.40 nmol m(-2) s(-1) in Picea abies, Pinus cembra and Larix decidua, respectively. These fluxes are consistent with findings from measurements of needle gas exchange, published from the same species at the study site. It is concluded that the sap flow-based approach offers an inexpensive, spatially and temporally integrating way for estimating ozone uptake at the whole-tree and stand level, intrinsicly covering the effect of boundary layers on ozone flux.     

Dry deposition of nitrogen and sulfur to Ponderosa and Jeffrey pine in the San Bernardino national forest in Southern California

Little is known about the concentrations, deposition rates, and effects of nitrogenous and sulfurous compounds in photochemical smog in the San Bernardino National Forest (SBNF) in southern California. Dry deposition of NO(3)(-) and NH(4)(+) to foliage of ponderosa pine (Pinus ponderosa Laws.) and Jeffrey pine (Pinus jeffreyi Grev. & Balf.) was correlated (R = 0.83-0.88) with historical average hourly O(3) concentations at 10 sites across an O(3) gradient in the SBNF. Mean deposition fluxes of NO(3)(-) to ponderosa and Jeffrey pine branches were 0.82 nmol M(-2)s(-1) at Camp Paivika (CP), a high-pollution site, and 0.19 nmol m(-2) s(-1) at Camp Osceola (CAO), a low-pollution site. Deposition fluxes of NH(4)(+) were 0.32 nmol m(-2) s(-1) at CP and 0.17 nmol m(-2) s(-1) at CAO, while mean values for SO(4)(2-) were 0.03 at CP and 0.02 nmol m(-2) s(-1) at CAO. Deposition fluxes to paper and nylon filters were higher in most cases than fluxes to pine branches at the same site. The results of this study suggest that an atmospheric concentration and deposition gradient of N and S compounds occurs along with the west-east O(3) gradient in the SBNF. Annual stand-level dry deposition rates for S and N at CP and CAO were estimated. Further studies are needed to determine if high N deposition loads in the SBNF significantly affect plant/soil nutrient relations, tree health, and the response of ponderosa pine to ozone.     

Measurement of the dry deposition flux of NH3 on to coniferous forest

The dry deposition flux of NH3 to coniferous forest was determined by the micrometeorological gradient method using a 36 m high tower. Aerodynamic characteristics of the site were studied, using a second tower erected in the forest 100 m from the first. Fluxes and gradients of heat and momentum measured on both towers indicated a fairly homogeneous turbulent flow field over the studied area of the forest. Site specific flux profile functions for heat were derived from continuous measurements of turbulent fluxes and gradients. These functions were used to derive fluxes from the observed gradients of NH3. In total, eighty 90-min NH3 flux runs were performed. The results indicate a strong nonstomatal uptake of NH3 by the forest. A representative dry deposition velocity for NH3 of 3.6 cm(-1) s was derived. The annual average flux was roughly estimated to be equivalent to 50 kg N ha(-1) yr, significantly higher than the critical load for coniferous forest.     

Ozone dynamics and deposition processes at a deforested site in the Amazon basin

The goal of this study is to investigate the impact of deforestation on ozone dynamics and deposition in the Brazilian Amazon basin. This goal is accomplished through i) analyses of ozone levels and deposition rates at a deforested site during the rainy season; and ii) comparisons of these data with similar information derived at a forest. At the pasture site maximum ozone mixing ratios reach 20 parts per billion on a volume basis (ppbv) but about 6 ppbv prevail over the forest. Maximum ozone deposition velocities for pastures can reach 0.7 cm s-1, which is about threefold lower than values derived for forests. Combining ozone abundance and deposition velocities, pasture maximum ozone fluxes reach approximately 0.2 microgram (ozone) m-2 s-1. This flux represents approximately 70% of the deposition rates measured over the forest. Hence, this study suggests that conversion of rainforests to pastures could lead to a net reduction (30%) in the ozone sink in the Amazon.