LBA-ECO CD-10 CO2 and H2O Eddy Flux Data at km 67 Tower Site, Tapajos National Forest
简介
该数据集包含一个文本文件,报告了帕拉西部(圣塔伦)地区 67 公里处原始森林塔点的二氧化碳和水交换涡流通量测量值以及相关气象测量值。该地点位于巴西中北部的塔帕若斯国家森林内。测量时间跨度为 2002 年 1 月至 2006 年 1 月。
使用塔式闭路 Licor 6262 气体分析仪和 Campbell CSAT3 声波风速计,在两个高度(58 米和 47 米 )测量了二氧化碳和水的涡旋通量(图 1)。涡旋通量测量的采样率为 8 赫兹,并以 1 小时为间隔进行平均。此外,还测量了一系列气象参数(气温、气压、光合有效辐射、净辐射、降水等)。
同地测量包括第三台 Licor 气体分析仪,用于测量 (a) 冠层内及上方 8 个高度的二氧化碳和水浓度分布(每 2 分钟测量 1 个高度);以及 (b) 冠层二氧化碳和水的瞬时总储量, 测量时采用的设计是同时从所有 8 个进气口吸入空气(每 20 分钟一次)。参见相关数据集。
摘要
The flux and meteorological data are reported in one comma separated ASCII text file, km67_eddyflux_2002_2006.txt .
Data File Documentation:
column variable description
1 "hours" Hour of measurement (GMT, continuous from 1/1/00)
2 "JDstart.GMT" Decimal day (GMT, continuous from 1/1/00)
NOTE 1: Tapajos Forest Local time (LT) = GMT - 4 hours
NOTE 2: these times are time at the BEGINNING of the hour-long
data aggregation interval, i.e., data at 12:00 are from
aggregating measurements between 12:00 and 13:00
3 "ws" Sonic wind speed, rotated u (m/sec, @ 57.8 m)
4 "wdir" Wind direction (degrees (0 - 360 degrees))
5 "Tamb" Ambient temp at 58m (level 1 eddy) from chilled mirror sensor (deg C)
6 "Tdew" Ambient dew point temperature @ 58m (deg C)
7 "Tson" Sonic temp (deg C, unadjusted for moisture @ 57.8m)
8 "Tasp" Distance weighted aspirated ambient temperature from thermistors (deg C)
9 "fheat" Heat flux (deg C m/sec)
10 "fmom" Momentum flux (m2/sec2)
11 "co2" CO2 concentration (mmol/mol)
12 "fco2" Eddy flux of CO2, (umol/m2/sec) [relative to dry air]
13 "h2o.mmol.m" H2O concentration (mmol/mol)
14 "fh2o" Eddy flux of H2O (mmol/m2/sec) [relative to dry air]
15 "h2o.mix" Mixing ratio (g H2O/kg (dry) air)
16 "Pamb.Pa" Ambient pressure (Pascals)
17 "H" Sensible heat flux (W/m2)
18 "LHdry" Latent heat flux (W/m2) [relative to dry air]
19 "ustar" Friction velocity, sqrt(-) (m/sec),
where w,u are the rotated wind components
20 "ppm2umol" Density conversion factor for fluxes (mol/m3) [relative to dry air]
21 "compiler" Which fortran compiler was used in processing [IDs ed1 vs. ed2]
22 "NetRad" Net Radiation @ 64.1 m (W/m2, corrected)
23 "sNetRad" Std deviation on hourly mean of NetRad
24 "par" PAR (umol/m2/sec) level 1 (63.6 m)
25 "sPAR1up" Std deviation on hourly mean of par
26 "PAR2up" PAR (umol/m2/sec) level 2 (15.09 m)
27 "sPAR2up" Std deviation on hourly mean of PAR2up
28 "PAR1dn" Downward PAR (umol/m2/sec) level 1 (63.6 m)
29 "sPAR1dn" Std deviation on hourly mean of PAR1dn
30 "Tair1" Temp (via thermistor) at profile level 1 (61.94 m) (deg C)
31 "Tair2" Temp (via thermistor) at profile level 2 (49.75 m) (deg C)
32 "WS1" Wind speed, Cup anemometer #1 (m/sec, 64.1 m)
33 "sWS1" Std deviation on hourly mean of WS1
34 "rain" Precipitation (mm in each hour) (tipping bucket @ 42.6m)
35 "co2col.wt" Mean column CO2 concentration (ppm)
36 "storage.wt" Storage flux (umol/m2/sec) below this level = d/dt (CO2 column avg)
based on a discrete integral approach using the profile concentration
measurements
37 "nee.wt" NEE (umol/m2/sec) = fco2 + storage.wt
38 "T.filled" Filled temperature for level 1 (filled Tasp). Filled based on adjusted
temperature measurements from the Jamaraqua met. station
(-2.80639, -55.03639). Data courtesy of D. Fitzjarrald (SUNY Albany)
39 "GMT" GMT time of day
40 "yr" Year
41 "NEE" Integrated NEE (umol/m2/sec) [filtered and integrated across gaps],
see status for filling information [u*>0.22]
42 "R" Respiration based on nighttime NEE with u*>0.22 m/s (umol/m2/sec)
see status for filling information
43 "GEE" Gross Ecosystem Exchange (umol/m2/sec)
see status for filling information
44 "R.light" Respiration based on light-curve interpolation to 0 PAR (umol/m2/sec)
45 "stg.filled" Filled storage (umol/m2/sec). This storage estimate has also been balanced
to 0 on 5 day intervals.
46 "par.filled" Filled PAR timeseries (umol/m2/sec) see par.status for filling codes
47 "par.status" Status 0 = unfilled data, measured with Wofsy sensor
Status 1 = filled based on relation between Wofsy PAR sensor & Wofsy NetRad
sensor over 10 day intervals
Status 2 = filled based on relation between D. Fitzjarrald PAR sensor &
Wofsy NetRad sensor over 10 day intervals
Status 3 = filled based on lookup table for a composite year, based on
entire time series.
48 "NetRad.filled" Filled NetRad timeseries (W/m2), see NetRad.status for filling codes
49 "NetRad.status" Status 0 = unfilled data, measured with Wofsy sensor
Status 1 = filled based on relation between Wofsy PAR sensor & Wofsy NetRad
sensor over 10 day intervals
Status 2 = filled based on relation between D. Fitzjarrald PAR sensor &
Wofsy NetRad sensor over 10 day intervals
Status 3 = filled based on lookup table for a composite year, based on
entire time series.
50 "hr.2" Hour of the day (GMT)
51 "Status" 6-bit status word to describe what, if any, filling technique was used
on the NEE, R, and GEE data
bit value description
0 1 u < 0.22 or u = NA; NEE = R + GEE
1 2 no PAR; PAR filled from lookup table for GEE calculation
2 4 no Tair available; filled in from Jamaraqua for R calculation
3 8 No measured NEE available; filled in with R+GEE
4 16 R calculated, not measured (daytime & calm nights)
5 32 Storage gap of >= 5 days; filled w/mean storage value
Missing value code is "NA"
Values are comma separated
Sample Data Records:
All of the flux and meteorological data are reported in km67_eddyflux_2002_2006.txt.
hours,JDstart.GMT.,ws,wdir,Tamb,Tdew,Tson,Tasp,fheat,fmom,co2,fco2,h2o.mmol.m.,fh2o,
h2o.mix,Pamb.Pa,H, LHdry,ustar,ppm2umol,compiler,NetRad,sNetRad,par,sPAR1up,PAR2up,sPAR2up,PAR1dn,
sPAR1dn,Tair1,Tair2, WS1,sWS1,rain,co2col.wt,storage.wt,nee.wt,T.filled,GMT,yr,NEE,R,
GEE, R.light,stg.filled,par.filled,par.status,NetRad.filled,NetRad.status,hr.2,Status
17610,733.75,1.95,179.804,29.183,23.13,32.475,29.079,0.043,-0.227,394.86,NA,30.834,4.402,
19.188,97867.5,47.953, 192.815,0.43,37.743,11,NA,0.128,541.576,36.264,138.926,13.245,24.118,
1.604,29.161,28.92, 2.799,0.94,0,373.085,0.539,NA,29.079,2.75,2002,-7.251,9.231,
-16.482, NA,0.437,541.576,0,239.15,3,18,24
17611,733.792,2.226,274.378,27.888,24.574,31.531,27.906,0.052,-0.197,377.74,-6.896,32.106,1.562,
19.98,97929.8,57.2, 68.492,0.41,37.805,11,NA,0.244,570.412,68.772,102.103,14.096,24.946,
2.716,27.978,27.768, 2.138,0.772,0,373.88,0.587,-6.309,27.906,2.792,2002,-6.411,9.231,
-15.642, NA,0.485,570.412,0,190.283,3,19,16
17612,733.833,1.856,275.332,28.21,24.312,31.909,28.337,0.022,-0.114,376.175,-3.952,31.252,1.011,
19.449,97965.3,24.009, 44.332,0.337,37.804,11,NA,0.13,337.28,38.094,67.313,6.905,15.59,
1.868,28.402,28.212, 1.746,0.684,0,374.584,0.518,-3.434,28.337,2.833,2002,-3.536,9.231,
-12.767, NA,0.416,337.28,0,59.45,3,20,16
...
**Line Breaks Added to Improve Readability
Site boundaries: (All latitude and longitude given in degrees and fractions)
Site (Region) Westernmost Longitude Easternmost Longitude Northernmost Latitude Southernmost Latitude Geodetic Datum
Para Western (Santarem) - km 67 Primary Forest Tower Site (Para Western (Santarem)) -54.959 -54.959 -2.857 -2.857 World Geodetic System, 1984 (WGS-84)
代码
!pip install leafmap
!pip install pandas
!pip install folium
!pip install matplotlib
!pip install mapclassify
import pandas as pd
import leafmap
url = "https://github.com/opengeos/NASA-Earth-Data"
df = pd.read_csv(url, sep="\t")
df
leafmap.nasa_data_login()
results, gdf = leafmap.nasa_data_search(
short_name="CD10_EddyFlux_Tapajos_860",
cloud_hosted=True,
bounding_box=(-55.21, -4.05, -54.91, -2.84),
temporal=("2002-01-01", "2006-01-18"),
count=-1, # use -1 to return all datasets
return_gdf=True,
)
gdf.explore()
leafmap.nasa_data_download(results[:5], out_dir="data")