(4) What is the effect of the final surface temperature on the carbon dioxide flux data? 1. Introduction

Size: px
Start display at page:

Download "(4) What is the effect of the final surface temperature on the carbon dioxide flux data? 1. Introduction"

Transcription

1 Correcting aarent off-season CO 2 utake ue to surface heating of an oen ath gas analyzer: rogress reort of an ongoing stuy G.G. urba, D. J. nerson, Liukang Xu an D.K. McDermitt LI-COR iosciences, Lincoln, NE, US bstract number of Ey Coariance flux sites with oen ath analyzers hae obsere small aarent CO 2 utake outsie the growing season. In this stuy, we roose an aitional correction ue to instrument surface heating as a art of the Webb-Pearman- Leuning correction. Effects of the roose correction were examine on hourly an seasonal time scales. Signifit reuction in aarent CO 2 utake uring off-season erios was obsere, while fluxes uring the growing season were not noticeably affecte. The roose correction may be useful for future CO 2 flux research, an also be alie to re-existing ata. 1. Introuction Oen-ath infrare gas analyzers are low maintenance high-quality sensors wiely use in fiel research (meriflux, 26). s more yearroun carbon flux stuies become aailable, there is growing concern about aarent off-season CO 2 utake being measure when it is hysiologically unreasonable to exect it (Hirata et al., 25). Soling this roblem is an imortant ste in imroing CO 2 measurements, annual NEE estimates, an subsequent global carbon exchange an climate moeling efforts. One laboratory stuy an two fiel exeriments are currently unerway to examine the influence of the oen-ath infrare gas analyzer surface temerature on the measure CO 2 flux using the Ey Coariance technique. The original concet an some initial results were resente at the meriflux nnual Meeting in October, 25 (urba et al., 25a,b), an has been eeloe further in subsequent months ue to the interest an enthusiasm generate in the flux research community towars resoling this issue. Here we resent the latest uates from these exeriments: a theoretical correction ue to the surface heating is introuce, the relationshi of instrument surface temerature an ambient air temerature is inestigate, an the effect of an aitional term on the CO 2 fluxes is examine. The following secific questions are aresse: (1) Do the electronics insie the sensor hea contribute substantially to instrument surface heating? (2) Is the air in the otical ath of the oen-ath analyzer signifitly warmer then the ambient air ue to the heate instrument surface? (3) How o the electronics an ambient factors combine to affect the final surface temerature? (4) What is the effect of the final surface temerature on the carbon ioxie flux ata? 2. Concet of the correction ue to instrument surface heating The surface of the oen-ath infrare gas analyzers (e.g., LI-75, LI-COR Inc., Lincoln, NE) coul be warmer than the ambient air ue to heat generate by the electronics in an instrument s sensor hea (choer motor, heating of the choer housing, infrare source, thermoelectric coolers, etc.), an by solar raiation. The instrument may also be coole by emission of long-wae raiation an by win-inuce conectie exchange (Fig. 1). The heat flux from an instrument surface is in aition to the ambient heat flux: coler air continuously enters the oen otical ath, an is warme an exane by the warm instrument surface in the same way that the atmoshere is warme by the soil surface (urba et al., 25a). nother way to look at this rocess is on the scale of two measurement ss (Fig. 1) for a situation in which there is no actual flux an no change in actual CO 2 mole fraction. With a ecrease in horizontal win see (s1), air in the ath gets warmer than ambient because there is more time to heat it, an exansion takes lace; so, CO 2 number ensity is lower than ambient. t the same time, w is ositie to assure momentum flux is always towar the surface. With an increase in horizontal win see (s 2), the air is heate less, exansion is less, an CO 2 number ensity is closer to ambient. t the same time w is negatie (to kee momentum transfer towar the surface). Thus, from the coariance calculation ( w 'C'), a false utake woul be obsere when no actual flux is occurring.

2 2 S 1 U weaker S 2 U stronger w w ball bounary layer sars lumes urising warm lume bounary layers U win - u negatie: -w is ositie - more heating - more exansion - reas less CO 2 - u ositie - w is negatie - less heating - less exansion - reas more CO 2 In blue flow atterns of warme air (bounary layers & wakes) In otte re thermal heat issiation (longwae) w u - less CO 2 w own - more CO 2 result - reaing CO 2 utake Figure 1: () Schematics of the stanar LI-75 setu & relate heat issiation; () Visualization of the concet: two measurement ss when there is no actual flux an no change in CO2 mole fraction occurs. S 1: with a ecrease in a horizontal win see, air in the ath gets warmer (more time to heat it), exansion takes lace, CO2 number ensity is lower, an w is ositie. For s 2 the situation is reerse. Therefore, an aarent utake be obsere when there is no actual flux. 3. Deriation of the correction ue to instrument surface heating The effects of such heating coul be inclue in the Webb-Pearman-Leuning correction (WPL: Webb et al., 198) as an aitional heat flux term (urba et al., 25a). The WPL correction, in general, comensates for the fluctuations of temerature an water aor affecting measure fluctuations of CO 2 an H 2 O. In general for an oen-ath CO 2 measurement, utting all constants together, the flux equation coul be written in the following form: Fc = Fc o E qc H C qc Ta (1) where Fc is final correcte CO 2 flux; ry air ensity; Fc o initial uncorrecte CO 2 flux; water aor ensity; E eaotransiration; total air ensity; q c mean ambient CO 2 ensity; C secific heat of air; H sonic sensible heat flux; T a air temerature. In a similar manner for water aor, the flux equation coul be written as follows: E = ( ) E o V H + ( ) C Ta (2) Now, let us imagine that the instrument is warm, an we not istinguish between the multile rocesses shown in Fig 1. Yet, we coul assume some resultant exansion of the air by a sensible heat flux create by instrument surface heating. Such a flux be estimate from the temerature graient using resistance aroach (Montieth, 1963) as follows:

3 3 C ( Ts Ta) HI = (3) ra where H I sensible heat flux from the warm instrument to the ambient air, T s instrument surface temerature, an r a aeroynamic resistance. For an instrument oriente near ertically, only a fraction of H I woul affect the measurement because the win will remoe much of the warme air from the otical ath, carrying away most of H without affecting the measurements. Howeer, a ortion of the warme an exane air will remain in the bounary layer, reucing the measure air ensity an influencing the measurement (Fig. 1). The fraction fr of H I retaine in the otical ath is calle H P an be exresse as follows: H P = frhi (4) So, combining Eqs. 1,3, an 4, the WPL woul be transforme as follows: Fc new E qc = Fco ( H + fr( C ( Ts Ta) / ra )) qc + C Ta (5) The term E in Eq. 5 also nees an aitional heat correction as well (Eq.2), such that: E = ( )Eo + ( H + fr( C ( Ts Ta ) / ra )) V + ( ) C Ta (6) If oen-ath flux ata hae alreay been correcte in the traitional way (frequency resonse an WPL), then the aitional correction for E woul be (from Eqs. 2 an 6): E new = E + fr( ) r a ( T ( T s a Ta) ) (7) an the correction to Fc for ata alreay correcte in the traitional way, woul be (from Eqs. 1,5 an 7): Fc new ( Ts Ta) qc = Fc + fr ( ) (8) r ( Ta ) a where ( Ts Ta) qc fr reresents the instrument ra ( Ta ) heating correction for Fc. The aitional correction to E in Eq. 6 is reresente by ( ), which is tyically ery small, on the orer of.5-2% for a majority of ambient temerature an humiity ranges. For the articular case of the LI-75, each of the two major structural elements of the instrument ( an ball, Fig. 1) may contribute ifferently to the roose aitional term in WPL correction: Fc new ( Ts Ta) qc = Fc + fr ( ) + ra ( Ta ) ball ( Ts Ta) qc + frball ( ) (9) ball r ( Ta ) a Resistances to heat transfer be estimate assuming that force conection reails near the surfaces of the instrument (Cambell an Norman, 1998): r 7. (1) U a 4 r ball 7. (11) U a ball 4 where is the characteristic imension of each element (.42 m for the ball, an.133 m for the ); =.133 m was comute from the raius of the cure making the to surface of the ). oth structural elements are assume here to be cyliners an sheres, so no aitional multilier for is require. ssume that the fraction of the sensible heat flux affecting measurement fr (=H P /H I ) coul be aroximate by the ratio of athlengths containing air influence by H P iie by the total mechanical athlength of the instrument (l ath =.128 m). The

4 4 athlengths affecte by H P will be the sum of the thermal bounary layer thicknesses of the an the ball. Thus: l fr = (12) ath l ball l ball fr = (13) ath l an fr=fr +fr ball. The layers thickness coul be aroximate, assuming laminar bounary layer an ertical instrument osition, as follows: l * Re 2 win + tg( α) (14) ball ball 5.5 l Re (15) x U Re = µ (16) where Re is Reynols number (Gates, 198); =.33 m (raius of the ); α=2 o (angle between the bottom winow an metal shoulers aroun it); win =.12 m (raius of the flat ortion of the winow) ); ball =.22 m (raius of the ball); x is or ball ; µ is ynamic iscosity of the air. The air heate insie the sensor between winows an transucers is CO 2 free, an its temerature fluctuations o not contribute to the flux correction. 4. Results 4.1.Electronics contribution enironmental chamber exeriment To ealuate contribution of the electronics in the sensor hea to the heating of the instrument surface, an exeriment was conucte in an enironmental chamber (Moel F466PC-1, lue M Electronic Comany, lue Islan, Illinois) with a ertical-own recirculating conitioning stream to ensure uniform mixing of the air. Thermocoules were attache to arious arts of the instrument (sie of the cyliner, lower winow area, uer winow area, an sar) an coere with small ieces of foam to minimize irect heating an cooling of the thermocoule by chamber air. Figure 2 shows the instrument surface temeratures lotte ersus ambient temerature for the range from -25 to +5 C. The bottom winow area was consistently warmer than ambient chamber air with an aerage relationshi as follows: T s =.92 T a +4.83, r 2 =1., n>9 (17) Equation 17 hols for the original ersion of the instrument electronic boar, house in the white control box (ersion 1). Eq. 17 emonstrates that ifferences between ambient temerature an the surface temeratures of the lower winow area increase with a ecrease in ambient temerature: the instrument was warmer than ambient by +.8 C at +5 C ambient, an by +6.8 C at -25 C ambient. Other arts of the instrument i not islay signifit ifference from the ambient temerature, though all relationshis were ery tight, with little or no hysteresis. These obserations are corroborate by the location of the majority of the internal instrument electronics, an articularly the heater for the instrument source, near the bottom winow of the instrument. somewhat smaller influence of electronic heating was obsere for the newer ersion of the instrument electronic boar (ersion 2): T s =.93 T a +3.53, r 2 =1., n>9 (18) The instrument was heate by the electronics to +.3 C aboe ambient at +5 C, an by +5.3 C at - 25 C ue to higher efficiency an lower heat issiation of the new boar configuration. Relationshis were linear in all cases, with total number of samlings n>9, groue by ambient temeratures in 9 grous, an yieling r-square of 1.. oth of these equations (Eqs. 17 an 18) emonstrate a signifit influence of the sensor hea electronics on the instrument surface heating, esecially at low ambient temeratures. lso, both are likely to be secific only to the ersion of the electronic boar, an shoul hol for all instrument configurations an for arious exerimental setus. In rincile, these relationshis coul be use, with aroriate assumtions an caution, to ealuate the effect of the surface heating on the measure flux (urba et al., 25a,b).

5 5 5 Enironmental Chamber 1:1 Ts, C sloe offset r 2 Lower winow area er er Uer winow area Ta, C 25 5 Figure 2. Instrument surface temerature (Ts) is lotte ersus air temerature (Ta) in the climate controlle chamber for two ersions of electronic boars. Each oint on the lot is an aerage of 2-3 one-minute ata oints after temerature equilibrium between chamber air an instrument surface was achiee. Data gathere before equilibrium was achiee hae been remoe from the analysis. 4.2 mbient air warming near the instrument surface LI-COR facility exeriment Fiel conitions are ifferent from those in the enironmentally controlle chamber, rimarily ue to factors influencing the instrument surface irectly (not ia electronics heating): sunlight, win an raiatie cooling. fiel exeriment is being conucte at the LI-COR exerimental facility (Lincoln, NE) starting in October, 25 to ealuate the resulting surface-to-air temerature graient in fiel conitions with an without the influence of the electronics. During the exeriment instrument ower was turne off for seeral weeks an on for seeral months. Fine-wire thermocoules were attache to LI-75 analyzer at the center of the lower winow area aroximately.5 mm an 1. mm aboe the lower winow, an in the center of the uer winow area. mbient fine-wire thermocoule was installe 2 cm away from the oen-ath sensor. Figure 3 emonstrates that air coming into the instrument ath is warme u signifitly by the heate instrument surface in the ambient enironment. Temerature graients are ery stee in the first 1 mm aboe the lower winow area: for col air in calm atmosheric conitions, graients may excee 2 C er 1 mm, while uring winy conitions graients seem to ecrease (Fig. 3). During warm ays graients are consierably smaller than in col temerature, while some win effect is still obsere (Fig. 3). Figure 3 C emonstrates that een without electronics heating (with ower off) the instrument surface be seeral egrees warmer than the ambient air in the mile of the ay (likely ue to sunlight), an seeral egrees cooler at night (ue to raiatie cooling). This henomenon has imortant imlications for any oen-ath instrument measuring flux: unless the instrument surface temerature is controlle to be similar to the ambient, instrument surface heating/cooling may signifitly affect the measure flux. 4.3 Combination of all factors affecting instrument surface temerature - exeriment at Mea, NE fiel exeriment to irectly ealuate the influence of instrument heating on the measure CO 2 flux uring the off-season is being conucte (starting Noember, 25) at the Uniersity of Nebraska Carbon Sequestration Program research site near Mea, Nebraska (in collaboration with S.. Verma an. Suyker). This site is equie

6 6 Col ay: 12/9/5 winy -6 T, C calm Lower winow area.5 mm eg C 1. mm eg C Uer winow area Tair Time Warm ay: 1/26/6 12 winy T, C 6 calm T, C 26 27C Time Power off: 12/9/ : 338 6: 12: 18: 24: 339 Central Stanar Time Figure 3. Daily atterns of surface temerature of the lower an uer winow areas; air temerature was measure at aroximately. 5 mm an 1 mm aboe the lower winow, an in the oen air away from the instrument on selecte ays: () col ay, 12/9/5; () warm ay, 1/26/6; an (C) ay with instrument ower turne off, 12/9/5. LI-COR facility exeriment. with a comrehensie Ey Coariance station, which has an oen-ath CO 2 /H 2 O gas analyzer (LI- 75), an a comlete set of instruments for measuring weather an raiation (Suyker et al., 25). Fine-wire thermocoules were attache to an LI-75 oen-ath analyzer near the lower winow, on the sie of the cyliner, on a sar, an near the uer winow. ll thermocoules were ainte with the same titanium oxie aint as the boy of the instrument to assure a similar albeo. Relationshis between the air an instrument surface temeratures near the lower an uer winows are shown in Fig. 4. Similar to the enironmental chamber exeriment, the temerature of the uer winow was not signifitly ifferent from the ambient, at least for the majority of the ata (linear regression yiele sloe of.99 an offset of.23 C). The lower winow area, howeer, was warmer than ambient, but not as much as insie of the enironmental chamber. Linear fits from all exeriments to ate are shown in Table I. mbient exeriments (LI- COR facility an Mea) rouce consistently lower Ts for the same Ta, as comare to the enironmental chamber. We attribute this to better thermal exchange in ambient conitions ue to

7 7 4 Mea exeriment: Noember, 25-March, :1 2 Ts, C equation r 2 Lower winow area linear fit Ts=.9Ta olynomial Ts=.25Ta 2 +.9Ta Uer winow area x Ts=.99Ta Figure 4. Hourly alues of instrument surface temerature (Ts) at the lower an uer winow areas lotte ersus air temerature (Ta) uring the erio from Noember, 25 to March, 26. Ta, C win conection, an ue to raiatie cooling at night (as seen in Fig. 3C). This is somewhat corroborate by frequently obsere conensation an ew formation near the to winow in col an humi fiel enironments. nother imortant feature of the ambient relationshi between Ts an Ta is that the best fit may not actually be linear, because a linear fit woul imly the occurrence of the cooling of the surface at ambient temeratures aboe +2 to +3 C. It is ifficult to justify such a situation, in which the electronics woul cool the instrument below the ambient temerature. Exeriment Ver. Influencing factors lower winow area uer winow area sloe offset r 2 n sloe offset r 2 n Climate control chamber 1 electronics > Climate control chamber 2 electronics > >9 Li-Cor Facility* 2 electronics & ambient > >2 Mea exeriment* 2 electronics & ambient > >15 *non-linear may be better: regression fit will become more comlete after warm-season ata become aailable Table I. Linear fit arameters of the relationshi between ambient air temerature (Ta) an instrument surface temerature (Ts) for the lower an uer winow areas: Ts = sloe Ta+ offset. Influence of the heating by electronics (chamber exeriment) is noticeably offset by atmosheric exchange for the outsie exeriments. Uer winow Ts i not seem to be consistently affecte by Ta. (*) Non-linear may be better: regression fit will become more comlete after warm-season ata become aailable. With only coler temeratures aailable in the Mea exeriment to ate, a olynomial fit not be confiently extraolate beyon +2C. Howeer, a olynomial relationshi for the ambient temerature range from -25 to +2 C oes seem to fit the ata as well or better, than a linear fit (Fig. 4):

8 8 Ts=.25Ta 2 +.9Ta+2.7, r 2 =.99, n>15 (19) Temeratures for the relationshi resente in Eq. 19 were measure irectly, an all factors influencing Ts (e.g., electronics, sunlight, raiatie cooling, an win exchange) hae been incororate into the actual measurement of Ts, howeer, this relationshi may also be less general than those resente in Eqs. 17 an 18, an may ary from site to site an from one setu to another (because of ifferences in shaowing of the instrument from sun an win, temerature ranges, instrument orientation etc.), so it shoul be use with caution..5.5 Maize: October 6, 22 Maize: January 29, 23 Fc, mg m -2 s ol new : 12: : C : 12: : D.5.5 Soybean: October 9, 22 Soybean: February 8, 23. Fc, mg m -2 s : 12: : -.5 : 12: : Central Stanar Time Figure 5. Daily atterns of hourly carbon ioxie flux (Fc) as obtaine using close-ath LI-COR analyzer (6262), traitionally correcte LI-75 (75 ol), an flux correcte with a roose aitional member of WPL correction (75 new): () maize, warm ay; () maize, col ay; (C) soybean, warm ay; (D) soybean, col ay. Utake of CO2 is ositie an release is negatie lying corrections to re-existing ata To take aantage of exactly the same instrument configuration an setu, we hae use Ts relationshi obtaine from Mea exeriment in 25-6 (Eq. 19), an alie it ia Eq. 9 to the traitionally correcte Fc flux ata from 22-3 from the same stuy site uner maize, an from a similar neighboring site uner soybean (ata courtesy of S.. Verma an. Suyker: Suyker et al., 25). Fc from the LI-75 was comare to that from a close-ath instrument, the LI-6262, which was not affecte by instrument surface heating, uring off-season with no green oy, an uring the entire year. Fc from the LI-75 with the traitional form of WPL was comare to that from the LI-75 with an aitional correction ue to instrument surface heating. Figs.5-D comare aily atterns of hourly Fc from the LI-6262 an LI-75 (correcte traitionally an with roose aitional correction) uring warm an col erios without a green oy. Utake of CO 2 is ositie an release is negatie. n aitional term in WPL (Eq. 9) hele to reuce the aarent off-season utake an mae the LI-75 an LI-6262 match reasonably well for such a small oerall flux (fluctuations are in the orer of.5 mg CO 2 m -2 s -1 ). The magnitue

9 9 of the aitional correction itself was een smaller: -.7 mg CO 2 m -2 s -1. The relatie size of the correction increase uring col erios with Fc fluctuations near zero. oth the ost-maize fiel (with substantial resiue) an the soybean fiel (ery little resiue) hae shown similar results. To assess the yearly atterns of the roose correction, hourly alues are shown for the entire year, incluing growing season, for maize an soybean in Figs.6(-). s reicte from Eqs. 9 an 19, the largest corrections were obsere uring colest erios. The majority of col-season corrections range from.25 to.7 mg m -2 s -1. Corrections were signifitly smaller uring warm season, ranging rimarily between.5 an.2 mg m -2 s -1 in both sites..1 Maize: Fc, mg m -2 s /1/2 12/1/2 1/31/3 4/2/3 6/2/3 8/2/3 1/2/3.1 Soybean: Fc, mg m -2 s /1/2 12/1/2 1/31/3 4/2/3 6/2/3 8/2/3 1/2/3 Figure 6. Seasonal atterns of the hourly alues of the roose aitional member of WPL ue to instrument surface heating ( Fc): () irrigate maize; () rainfe soybean. Note the minimal magnitue of this correction - scale on the lot is from to.1 mg CO2 m-2 s-1. Utake of CO2 is ositie an release is negatie. The imact of the roose correction on the resultant hourly Fc is emonstrate in Figs.7(-D). The ratio of Fc from the LI-6262 to that from the traitionally correcte LI-75 is lotte ersus Ta for the erio from October 22 through October 23. In both maize (Fig. 7) an soybean (Fig. 7C), the ratio signifitly eiate from 1when T a was below 2C, an change sign below 1 C, inicating hysiologically unreasonable CO 2 utake. Using the roose correction (Figs. 7,D) eliminate most of the sign fliing, an returne the ratio to nearly 1. It is imortant to note that correction i not just increase the ratio, but greatly reuce ariability on both sies of zero (Figs.

10 1 1 Maize: Fc 75 ol 1 Maize: Fc 75 new Fc 6262/Fc C Soybean: Fc 75 ol Soybean: Fc 75 new D Fc 6262/Fc Ta, C Figure 7. Hourly ratio of carbon ioxie flux measure with LI-6262 (Fc 6262) to that measure with LI-75 (Fc 75) correcte in a traitional way (ol) an correcte with a roose aitional member of WPL correction (new): () maize, traitional WPL; () maize, WPL with aitional term; (C) soybean, traitional WPL; (D) soybean, WPL with aitional term. The erio is from October 22 to October, 23. itional correction imroe both means an istribution of the ata, an hae not just shifte all ata uwar. 7,D). It is also worth mentioning that ery small quantities are iie in these ratios, an therefore some ariability shoul be execte, esecially at low temeratures. The contrast between seemingly small hourly alues of the roose correction an its imact on the integrate fluxes be emonstrate by comaring Figs. 8 (, ) an Table II (-). Using the roose correction reuce the sloe from 3% to 2% an reuce the offset from.3 to.1 mg m -2 s -1 in maze (Fig. 8), an from.3 to zero mg m -2 s -1 in soybean (sloe was not affecte; Fig. 8). This shows that off-season fluxes are somewhat affecte, while growing season fluxes i not signifitly change. In contrast, the imact of the correction on off-season integrations of Fc was extremely imortant (Table II). The roose correction ramatically reuce the number of utake hours (Table II), from 716 to 252 hours in maize, an from 571 to 189 hours in soybean. s a result, the unerestimation of off-season losses of CO 2 (an resulting oerestimation of the yearly utake of CO 2 ) was reuce seeral fol (from -8 g CO 2 to -34 g in maize an from -32 g CO 2 to g in soybean (Table II: utake is ositie an release is negatie). fter alying the roose correction the LI-6262 an LI-75 ata matche to within few ercent, which is well within the accuracy of the Ey Coariance metho. 5. Discussion 5.1 Note of caution The metho escribe aboe for correcting Fc ue to instrument surface heating was relying, on the assumtion that the instrument was oriente nearly ertically. There are two ariables in Eq. 9 that woul be signifitly affecte by ifferent instrument configurations: fr an T s. The

11 11 Fc 75, mg m -2 s Maize: black - Fc 75 ol re - Fc 75 new y =.97x +.3 R 2 =.99 y =.98x +.1 R 2 = Fc 6262, mg m -2 s -1 Fc 75, mg m -2 s Soybean black - Fc 75 ol re - Fc 75 new y = 1.2x +.3 R 2 =.99 y = 1.2x +. R 2 = Fc 6262, mg m -2 s -1 Figure 8. Hourly carbon ioxie flux measure with LI-75 (Fc 75) lotte ersus that measure with LI-6262 (Fc 6262) in: () maize, an () soybean uring the erio is from October 22 to October, 23. 'Fc 75 ol' was correcte in a traitional way, an 'Fc 75 new' was correcte with a roose aitional member of WPL correction. itional correction was so small as comare to the hourly fluxes uring the growing season, that the latter were not noticeably affecte. Utake of CO2 is ositie an release is negatie. Off-season utake hours LI-6262 LI-75 ol LI-75 new Irrigate maize fallow Rainfe soybean fallow Off-season CO2 release, when both 6262 an 75 aailable, g C LI-6262 Irrigate maize fallow Rainfe soybean fallow LI-75 ol LI-75 new Table II. Comarison of the () off-season utake hours an () integrate carbon ioxie release between a close ath analyzer (LI-6262), traitionally correcte oen-ath ata (LI-75 ol) an ata correcte using the roose aitional member of WPL correction (LI-75 new). ll tables an lots inclue only hours with comlete ata: when actual reaings from the LI-62, LI-75, an other instruments use in the roose corrections were aailable. Non-stationary conitions, rain, snow, instrument malfunctions an other fille-in erios were exclue to assure roer comarison between Fc75 an Fc6262. Off-season erios are erios without green foliage area from October 1 to ril 3.

12 12 orientation of the instrument may imact the fr term, since the instrument woul not be ertical, an the ortion of the ath affecte by warm air may be higher than assume for a nearly-ertical instrument. lso, in such a situation, Ts shoul inclue the influence of the surface temerature that woul now also affect the ath, in aition to the winow area temerature. Win irection in a non-ertical configuration woul be an imortant factor as well, since a configuration with a signifitly incline sensor woul not be symmetrical for all win irections. In aition to instrument orientation, the lacement of other instruments an the tower might affect both fr an r a in Eq. 9, thus affecting the correction, esecially when win blows through the structures into the LI- 75 ath. Therefore, the finings resente here shoul be treate with caution, as only an uate on the ongoing inestigation of the imact of instrument surface heating on measurements mae by a nearly-ertical oen-ath CO 2 /H 2 O analyzer. 5.2 lternatie way to ealuate effect of arious sensor configurations In some cases there may be a way to ealuate the effects of a articular comlex sensor configuration on fr an r a in Eq.9, an thus aoi some assumtions (e.g., near-ertical sensor osition, realence of laminar bounary layer aroun sensor arts, absence of flow obstruction from other nearby sensors). If the correct alue of flux is known from close ath instrument or coul be confiently assume for a range of win sees an win irections, Eq.9 may be sole for fr/r a ratio. This woul establish a site-secific fr/r a ratio for ifferent win sees an irections, assuming these relationshis woul not change signifitly with temerature, humiity or ressure. Such solution may inclue the short-term use of the close-ath sensor (e.g., LI-6262 or LI-7) as stanar, or establishing erios when flux alue shoul be hysiologically close to zero (e.g., extremely low temeratures below -25 or -3 C). The rest of the ata coul be correcte for arious sensor configurations using these ratios in Eq.9 in lace of the Eqs This aroach shoul be treate with secial care ue to ossible statistical errors inole in estimating the fr/r a ratio for arious win sees an irections, an ue to the number of ata oints require to construct confient emirical relationshis; howeer, using such methos on re-existing flux ata from sites with comlicate sensor configurations is erhas still more beneficial than rejecting large erios in the measurement year. 5.3 Future research The ata in the ongoing exeriments were collecte rimarily uring col conitions with few regression oints aboe 15C, an with reominantly low fluxes. This force us to aly the emirical relationshi between irectly measure Ts an Ta eeloe in to the Fc ata from With more ata to be collecte uring the warm season of 26, further tests will inclue comarison of fluxes between the close ath LI-6262 an oen ath LI-75 for a wie range of irectly measure surface temeratures. There are seeral other stes that may hel to further stuy an manage the effects of instrument heating. Haing thermocoules on the surface near the lower an uer winows woul hel to secify the roose correction without the ifficulty of haing to estimate Ts (e.g., Eqs ). Further exeriments on the best sensor orientation might also be one to minimize the effects of instrument surface heating. Particularly, usie-own sensor ositioning might reuce the effect of the heating on the flux in all irections. Further research is also neee to ealuate, test an fine-tune the roose concet for arious enironments, an for ifferent manufacturers of oen-ath CO2 analyzers. 6. Summary an conclusion The influence of surface heating on the CO 2 flux measure with the LI-75 oen ath gas analyzer is being examine in one laboratory stuy an two fiel exeriments. Exeriment in a conectie enironmental chamber showe that the sensor hea electronics substantially contribute to surface heating, an that the surface of the oenath analyzer coul be 5-7 C warmer than the ambient air ue to electrical energy issiation in the instrument hea. The ifference was largest at low ambient temeratures an ecrease linearly with an increase in the ambient temerature. fiel exeriment at the LI-COR facility, which is still in rogress, has shown that large temerature graients (u to 2 C er 1mm) exist just aboe the lower winow of the instrument, imlying strong sensible heat fluxes in the first few mm aboe the instrument surface. The thermal

13 13 bounary layer was thicker uring calm conitions, as execte. During the ay, surface heating occurre ue to solar raiation een when the electronics were not owere. Our measurements also suggest that emission of long-wae raiation cause surface cooling at night to temeratures below ambient. n exeriment is in rogress near Mea, NE also in which we are comaring Fc measure with an oen-ath LI-75 to Fc measure by a close ath LI-6262, which is not affecte by surface heating an is iewe as stanar. To ate, this exeriment has emonstrate, that instrument surface temerature is irectly an strongly relate to air temerature (Eq. 19), but lack of ata uring warm season i not allow for confient extraolation of the relationshi beyon +2 C. n aitional term ue to instrument surface heating was introuce as a art of the Webb- Pearman-Leuning correction. The temerature relationshi eeloe in the Mea exeriment in (Eq. 18) was alie to LI-75 ata from 22-3 in maize an soybean, an was comare against LI-6262 ata from the same erio. Een with artial results from the exeriments aailable, a signifit reuction in aarent CO 2 utake uring off-season erios was obsere as a result of alying the aitional correction, while fluxes uring the growing season hae not been noticeably affecte. The correction also resulte in elimination of most of the wrong signs from off-season oen-ath CO 2 measurements, a signifit reuction in ariability of the ata, elimination of the offset between measurements mae with the LI-6262 an the LI- 75, an in a signifit imroement in offseason integrations of CO 2 exchange. framework was create to eelo a site-secific ractical correction ue to instrument heating. The roose aitional term in the Webb-Pearman-Leuning correction may roe useful for future year-roun CO 2 flux research. The technique to estimate the roose correction for re-existing ata may hel in better interretation of reiously collecte information. cknowlegements ata an facilities for the Mea exeriment, an for numerous aluable hel an aice. We also thank Drs. rian miro, Robert Clement, Dai Cook, Marc-nré Giasson, Ryuichi Hirata, Ray Leuning, ill Massman, Walt Oechel, Howar Skinner, an Jon Welles for their interest in eeloing the toic. We acknowlege hel roie by Jerry Ficke of LI-COR on LI-75 s electronics an other roerties, an for assistance with the enironmental chamber exeriment. We also thank Fluxnet, meriflux, Fluxnet-Canaa, an Carbo- Euroe networks for roiing oortunity for the iscussions an for access to flux ata obtaine with LI-75 an LI-6262 instruments. References: meriflux, 26. meriflux Sites Information." <htt://ublic.ornl.go/ameriflux/site-select.cfm> urba G. G., nerson D. J., Xu L., an D. K. McDermitt. 25a. Soling the off-season utake roblem: correcting fluxes measure with the LI-75 for the effects of instrument surface heating. Progress reort of the ongoing stuy. PRT I: THEORY. Poster resentation. meriflux 25 nnual Meeting, ouler, Colorao. urba G. G., nerson D. J., Xu L., an D. K. McDermitt. 25b. Soling the off-season utake roblem: correcting fluxes measure with the LI-75 for the effects of instrument surface heating. Progress reort of the ongoing stuy. PRT II: RESULTS. Poster Presentation. meriflux 25 nnual Meeting, ouler, Colorao. Cambell, G. S., an J. M. Norman, Enironmental iohysics. Secon Eition. Sringer Science, New York: 286. Gates, D. M iohysical Ecology. Sringer-Verlag, NY. Hirata R., Hirano T., Mogami J., Fujinuma Y., Inukai K., Saigusa N., an S. Yamamoto, 25. CO 2 flux measure by an oen-ath system oer a larch forest uring snowcoere season. Phyton, 45: Monteith, J. L., Gas exchange in lant communities. Enironmental Control of Plant Growth. L. T. Eans, e. caemic Press, New York: Suyker,.E., S.. Verma, G.G. urba, an T.J. rkebauer. 25. Gross rimary rouction an ecosystem resiration of irrigate maize an irrigate soybean uring a growing season. gricultural an Forest Meteorology, 131: Verma, S.., an. E. Suyker. Personal Communications. Setember-December, 25. Webb, E. K., Pearman, G. I., an R. Leuning, 198. Correction of flux measurements for ensity effects ue to heat an water aour transfer. Quart. J. R. Met. Soc., 16: 85-1 We woul like to acknowlege a number of scientists who contribute to this stuy ia aluable iscussions, aice, ata sets, an other ersonal communications. Particularly we want to thank Drs. Shashi Verma an nrew Suyker for roiing

[ K] McAlpine, hmwk #5, ATMS 611. J.D. McAlpine ATMS 611 HMWK #5

[ K] McAlpine, hmwk #5, ATMS 611. J.D. McAlpine ATMS 611 HMWK #5 McAline, hmwk 5, ATMS 6 J.D. McAline ATMS 6 HMWK 5 3.43 Rainro evaorating (is 2 C) : 8 C, calculate mixing ratio of - sat. mixing ratio at 2 C is 8.7 g/kg (assume wet-bulb temerature) - L e = 2.25*0 6

More information

The thermal wind 1. v g

The thermal wind 1. v g The thermal win The thermal win Introuction The geostrohic win is etermine by the graient of the isobars (on a horizontal surface) or isohyses (on a ressure surface). On a ressure surface the graient of

More information

The Effect of a Finite Measurement Volume on Power Spectra from a Burst Type LDA

The Effect of a Finite Measurement Volume on Power Spectra from a Burst Type LDA The Effect of a Finite Measurement Volume on Power Sectra from a Burst Tye LDA Preben Buchhave 1,*, Clara M. Velte, an William K. George 3 1. Intarsia Otics, Birkerø, Denmark. Technical University of Denmark,

More information

Part 3. Atmospheric Thermodynamics The Gas Laws

Part 3. Atmospheric Thermodynamics The Gas Laws Part 3. Atmosheric Thermoynamics The Gas Laws Eq. of state V = mrt mass gas constant for 1 kg of a gas For ry air ρ = m/ = ρrt or α = RT where α = 1/ ρ = ρ R T, where R a R M * = uniersal gas constant

More information

ANALYSIS OF PNEUMATIC FINE PARTICLE PEENING PROCESS BY USING A HIGH-SPEED-CAMERA

ANALYSIS OF PNEUMATIC FINE PARTICLE PEENING PROCESS BY USING A HIGH-SPEED-CAMERA International Journal of Moern Physics B Vol. 24, Nos. 15 & 16 (21) 347 352 Worl Scientific Publishing Comany DOI: 1.1142/S217979216669 ANALYSIS OF PNEUMATIC FINE PARTICLE PEENING PROCESS BY USING A HIGH-SPEED-CAMERA

More information

By completing this chapter, the reader will be able to:

By completing this chapter, the reader will be able to: hater 4. Mechanics of Particles Particle mechanics governs many rinciles of article measurement instrumentation an air cleaning technologies. Therefore, this chater rovies the funamentals of article mechanics.

More information

Skiba without unstable equlibrium in a linear quadratic framework

Skiba without unstable equlibrium in a linear quadratic framework Skiba without unstable equlibrium in a linear quaratic framework Richar F. Hartl Institute of Management, University of Vienna, Brünnerstraße 8, -0 Vienna, ustria Richar.Hartl@univie.ac.at Tel. +43--477-3809

More information

Numerical study of the movement of fine particle in sound wave field

Numerical study of the movement of fine particle in sound wave field Available online at www.scienceirect.com ScienceDirect Energy Proceia 75 (2015 ) 2415 2420 The 7 th International Conference on Alie Energy ICAE2015 Numerical stuy of the movement of fine article in soun

More information

5. THERMAL CONVERSION OF SOLAR RADIATION. Content

5. THERMAL CONVERSION OF SOLAR RADIATION. Content 5. Introuction 5. THEMAL CONVESION OF SOLA ADIATION Content 5. Introuction 5. Collectors without concentration 5.. Otical efficiency of the flat collector 5.. Thermal efficiency of the flat collector 5..3

More information

Submitted to the Journal of Hydraulic Engineering, ASCE, January, 2006 NOTE ON THE ANALYSIS OF PLUNGING OF DENSITY FLOWS

Submitted to the Journal of Hydraulic Engineering, ASCE, January, 2006 NOTE ON THE ANALYSIS OF PLUNGING OF DENSITY FLOWS Submitte to the Journal of Hyraulic Engineering, ASCE, January, 006 NOTE ON THE ANALYSIS OF PLUNGING OF DENSITY FLOWS Gary Parker 1, Member, ASCE an Horacio Toniolo ABSTRACT This note is evote to the correction

More information

Discharge initiation and plasma column formation in aspect ratio A=2 tokamak.

Discharge initiation and plasma column formation in aspect ratio A=2 tokamak. Discharge initiation an lasma column formation in asect ratio A toama... Khayrutinov E.A. Azizov, A.D. Baralov, G.G.Glaush, I.L.Taibaeva, Ph.W.West 3 Troits, Moscow eg., ussia NNC, K 3 General Atomics,

More information

Lecture 6 : Dimensionality Reduction

Lecture 6 : Dimensionality Reduction CPS290: Algorithmic Founations of Data Science February 3, 207 Lecture 6 : Dimensionality Reuction Lecturer: Kamesh Munagala Scribe: Kamesh Munagala In this lecture, we will consier the roblem of maing

More information

Free Vibration Analysis of a Model Structure with New Tuned Cradle Mass Damper

Free Vibration Analysis of a Model Structure with New Tuned Cradle Mass Damper Proc. Schl. Eng. Tokai Tokai Univ., Univ., Ser. ESer. E 37 (0) 37(0)3-8 3-8 Free Vibration Analysis of a Moel Structure with New Tune Crale Mass Damer by Jitjinakun AMONPHAN * an Yoji SHIMAZAKI * (Receive

More information

6B.7 EXPERIMENTAL VALIDATION OF THE WEBB CORRECTION FOR CO2 FLUX WITH AN OPEN-PATH GAS ANALYZER

6B.7 EXPERIMENTAL VALIDATION OF THE WEBB CORRECTION FOR CO2 FLUX WITH AN OPEN-PATH GAS ANALYZER 6B.7 EXPERIMENAL VALIDAION OF HE WEBB CORRECION FOR CO2 FLUX WIH AN OPEN-PAH GAS ANALYZER Osamu sukamoto* an Fumiyoshi Kono Okayama University, Okayama, Japan 1. INRODUCION urbulent flux by the ey covariance

More information

1. Radiative transfer with emission. Atmosphere and surfaces emit infrared and microwave radiation. Recall the Beer-Bouguer-Lambert law for emission

1. Radiative transfer with emission. Atmosphere and surfaces emit infrared and microwave radiation. Recall the Beer-Bouguer-Lambert law for emission Lecture 8. Princiles of assive remote sensing using emission an alications: Remote sensing of atmosheric ath-integrate quantities clou liqui water content an reciitable water vaor.. Raiative transfer with

More information

Combining Time Series and Cross-sectional Data for Current Employment Statistics Estimates 1

Combining Time Series and Cross-sectional Data for Current Employment Statistics Estimates 1 JSM015 - Surey Research Methos Section Combining Time Series an Cross-sectional Data for Current Employment Statistics Estimates 1 Julie Gershunskaya U.S. Bureau of Labor Statistics, Massachusetts Ae NE,

More information

Fundamental Laws of Motion for Particles, Material Volumes, and Control Volumes

Fundamental Laws of Motion for Particles, Material Volumes, and Control Volumes Funamental Laws of Motion for Particles, Material Volumes, an Control Volumes Ain A. Sonin Department of Mechanical Engineering Massachusetts Institute of Technology Cambrige, MA 02139, USA August 2001

More information

On Using FASTEM2 for the Special Sensor Microwave Imager (SSM/I) March 15, Godelieve Deblonde Meteorological Service of Canada

On Using FASTEM2 for the Special Sensor Microwave Imager (SSM/I) March 15, Godelieve Deblonde Meteorological Service of Canada On Using FASTEM2 for the Secial Sensor Microwave Imager (SSM/I) March 15, 2001 Godelieve Deblonde Meteorological Service of Canada 1 1. Introduction Fastem2 is a fast model (multile-linear regression model)

More information

A MECHANISM FOR THE TEMPERATURE DEPENDENCE OF THE SURFACE RELAXATION RATE IN CARBONATES

A MECHANISM FOR THE TEMPERATURE DEPENDENCE OF THE SURFACE RELAXATION RATE IN CARBONATES SCA2002-27 1/12 A MECHANISM FOR THE TEMPERATURE DEPENDENCE OF THE SURFACE RELAXATION RATE IN CARBONATES Christian Straley, Schlumberger-Doll Research, Rigefiel, Connecticut ABSTRACT Unexlaine temerature

More information

Robust Control of Robot Manipulators Using Difference Equations as Universal Approximator

Robust Control of Robot Manipulators Using Difference Equations as Universal Approximator Proceeings of the 5 th International Conference of Control, Dynamic Systems, an Robotics (CDSR'18) Niagara Falls, Canaa June 7 9, 218 Paer No. 139 DOI: 1.11159/csr18.139 Robust Control of Robot Maniulators

More information

Course Outline. Section Instructor Time Atmospheric Variables Cotton 1 class. vd Heever ~4 weeks

Course Outline. Section Instructor Time Atmospheric Variables Cotton 1 class. vd Heever ~4 weeks Course Outline 1 2 3 4 5 6 7 8 9 Section Instructor ime Atmosheric Variables Cotton 1 class Communication with the atmos Jirak 1 class Atmosheric hermoynamics Atmosheric Dynamics Heeer ~4 weeks Heeer ~4

More information

Soil Plasticity and the Structured Cam Clay Model

Soil Plasticity and the Structured Cam Clay Model Soil Plasticity an the Structure Cam Clay Moel Dr. Martin D. Liu School of Ciil Engineering Uniersity of New South Wales Syney, NSW 2052, Australia Tel.: 2-9385-5474 Fax: 2-9385-639 e-mail: martin.liu@unsw.eu.au

More information

P1D.6 IMPACTS OF THE OCEAN SURFACE VELOCITY ON WIND STRESS COEFFICIENT AND WIND STRESS OVER GLOBAL OCEAN DURING

P1D.6 IMPACTS OF THE OCEAN SURFACE VELOCITY ON WIND STRESS COEFFICIENT AND WIND STRESS OVER GLOBAL OCEAN DURING P1D.6 IMPATS OF THE OEAN SURFAE VELOITY ON WIND STRESS OEFFIIENT AND WIND STRESS OVER GLOBAL OEAN DURING 1958-001 Zengan Deng 1* Lian Xie Ting Yu 1 an Kejian Wu 1 1. Physical Oceanography Laboratory, Ocean

More information

6. Friction and viscosity in gasses

6. Friction and viscosity in gasses IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner

More information

Effect of the variable porosity on the heat transfer process in solar air receiver

Effect of the variable porosity on the heat transfer process in solar air receiver Engineering Conerences International ECI Digital Archives Sixth International Conerence on Porous Meia an Its Alications in Science, Engineering an Inustry Proceeings 7-5-2016 Eect o the variable orosity

More information

kg m kg kg m =1 slope

kg m kg kg m =1 slope (5) Win loa Wen structure blocks te flow of win, te win's kinetic energy is converte into otential energy of ressure, wic causes a win loaing. ensity an velocity of air te angle of incience of te win 3

More information

A Critical State Sand Model with Elastic-Plastic Coupling

A Critical State Sand Model with Elastic-Plastic Coupling A Critical State Sand Model with Elastic-Plastic Couling Ali Lashkari * and Ali Golchin Deartment of Ciil & Enironmental Engineering Shiraz Uniersity of Technology, Shiraz, Iran lashkari@sutech.ac.ir,

More information

An Improved Calibration Method for a Chopped Pyrgeometer

An Improved Calibration Method for a Chopped Pyrgeometer 96 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 17 An Imroved Calibration Method for a Choed Pyrgeometer FRIEDRICH FERGG OtoLab, Ingenieurbüro, Munich, Germany PETER WENDLING Deutsches Forschungszentrum

More information

Fundamental Laws of Motion for Particles, Material Volumes, and Control Volumes

Fundamental Laws of Motion for Particles, Material Volumes, and Control Volumes 1 Funamental Laws of Motion for Particles, Material Volumes, an Control Volumes Ain A. Sonin Department of Mechanical Engineering Massachusetts Institute of Technology Cambrige, MA 02139, USA March 2003

More information

Multidisciplinary System Design Optimization (MSDO)

Multidisciplinary System Design Optimization (MSDO) Multiiscilinary System Design Otimization (MSDO) Graient Calculation an Sensitivity Analysis Lecture 9 Olivier e Weck Karen Willco Massachusetts Institute of Technology - Prof. e Weck an Prof. Willco Toay

More information

COMPUTER SIMULATION OF A LABORATORY HYDRAULIC SYSTEM WITH MATLAB-SIMULINK

COMPUTER SIMULATION OF A LABORATORY HYDRAULIC SYSTEM WITH MATLAB-SIMULINK DVNCED ENGINEERING 4(20101, ISSN 1846-5900 COMPUTER SIMULTION OF LORTORY HYDRULIC SYSTEM WITH MTL-SIMULINK Grego, G. & Siminiati, D. bstract: The article resents some selected roblems related to modeling

More information

Application of Measurement System R&R Analysis in Ultrasonic Testing

Application of Measurement System R&R Analysis in Ultrasonic Testing 17th Worl Conference on Nonestructive Testing, 5-8 Oct 8, Shanghai, China Alication of Measurement System & Analysis in Ultrasonic Testing iao-hai ZHANG, Bing-ya CHEN, Yi ZHU Deartment of Testing an Control

More information

Water vapour transport through perforated foils

Water vapour transport through perforated foils Water apour transport through perforate foils Wim an er Spoel TU Delft, Faculty of Ciil Engineering an Geosciences, Section Builing Engineering, Steinweg, 68 CN Delft, The Netherlans E-mail: w.h.anerspoel@citg.tuelft.nl

More information

He s Homotopy Perturbation Method for solving Linear and Non-Linear Fredholm Integro-Differential Equations

He s Homotopy Perturbation Method for solving Linear and Non-Linear Fredholm Integro-Differential Equations nternational Journal of Theoretical an Alie Mathematics 2017; 3(6): 174-181 htt://www.scienceublishinggrou.com/j/ijtam oi: 10.11648/j.ijtam.20170306.11 SSN: 2575-5072 (Print); SSN: 2575-5080 (Online) He

More information

JJMIE Jordan Journal of Mechanical and Industrial Engineering

JJMIE Jordan Journal of Mechanical and Industrial Engineering JJMIE Jordan Journal of Mechanical and Industrial Engineering Volume, Number, Jun. 8 ISSN 995-6665 Pages 7-75 Efficiency of Atkinson Engine at Maximum Power Density using emerature Deendent Secific Heats

More information

Wolfgang POESSNECKER and Ulrich GROSS*

Wolfgang POESSNECKER and Ulrich GROSS* Proceedings of the Asian Thermohysical Proerties onference -4 August, 007, Fukuoka, Jaan Paer No. 0 A QUASI-STEADY YLINDER METHOD FOR THE SIMULTANEOUS DETERMINATION OF HEAT APAITY, THERMAL ONDUTIVITY AND

More information

Improving the performance of geothermal power plant by substituting water steam power plant for the binary power plant

Improving the performance of geothermal power plant by substituting water steam power plant for the binary power plant oak, orsukieicz-ozur an Stachel Imroving the erformance of thermal oer lant by substituting ater steam oer lant for the binary oer lant Wlaysla oak, Aleksanra orsukieicz-ozur an Aleksaner A. Stachel Technical

More information

Consistency and asymptotic normality

Consistency and asymptotic normality Consistency an asymtotic normality Class notes for Econ 842 Robert e Jong March 2006 1 Stochastic convergence The asymtotic theory of minimization estimators relies on various theorems from mathematical

More information

University of North Carolina-Charlotte Department of Electrical and Computer Engineering ECGR 4143/5195 Electrical Machinery Fall 2009

University of North Carolina-Charlotte Department of Electrical and Computer Engineering ECGR 4143/5195 Electrical Machinery Fall 2009 University of North Carolina-Charlotte Deartment of Electrical and Comuter Engineering ECG 4143/5195 Electrical Machinery Fall 9 Problem Set 5 Part Due: Friday October 3 Problem 3: Modeling the exerimental

More information

Mod p 3 analogues of theorems of Gauss and Jacobi on binomial coefficients

Mod p 3 analogues of theorems of Gauss and Jacobi on binomial coefficients ACTA ARITHMETICA 2.2 (200 Mo 3 analogues of theorems of Gauss an Jacobi on binomial coefficients by John B. Cosgrave (Dublin an Karl Dilcher (Halifax. Introuction. One of the most remarkable congruences

More information

Bivariate distributions characterized by one family of conditionals and conditional percentile or mode functions

Bivariate distributions characterized by one family of conditionals and conditional percentile or mode functions Journal of Multivariate Analysis 99 2008) 1383 1392 www.elsevier.com/locate/jmva Bivariate istributions characterize by one family of conitionals an conitional ercentile or moe functions Barry C. Arnol

More information

Finite-Sample Bias Propagation in the Yule-Walker Method of Autoregressive Estimation

Finite-Sample Bias Propagation in the Yule-Walker Method of Autoregressive Estimation Proceedings of the 7th World Congress The International Federation of Automatic Control Seoul, Korea, July 6-, 008 Finite-Samle Bias Proagation in the Yule-Walker Method of Autoregressie Estimation Piet

More information

A secure approach for embedding message text on an elliptic curve defined over prime fields, and building 'EC-RSA-ELGamal' Cryptographic System

A secure approach for embedding message text on an elliptic curve defined over prime fields, and building 'EC-RSA-ELGamal' Cryptographic System International Journal of Comuter Science an Information Security (IJCSIS), Vol. 5, No. 6, June 7 A secure aroach for embeing message tet on an ellitic curve efine over rime fiels, an builing 'EC-RSA-ELGamal'

More information

Colin Cameron: Asymptotic Theory for OLS

Colin Cameron: Asymptotic Theory for OLS Colin Cameron: Asymtotic Theory for OLS. OLS Estimator Proerties an Samling Schemes.. A Roama Consier the OLS moel with just one regressor y i = βx i + u i. The OLS estimator b β = ³ P P i= x iy i canbewrittenas

More information

Integrating Lidar and Atmospheric Boundary Layer Measurements to Determine Fluxes and Dynamics of Particulate Emissions from an Agriculture Facility

Integrating Lidar and Atmospheric Boundary Layer Measurements to Determine Fluxes and Dynamics of Particulate Emissions from an Agriculture Facility Integrating Lidar and Atmosheric Boundary Layer Measurements to Determine Fluxes and Dynamics of Particulate Emissions from an Agriculture Facility Lawrence His 1, John H. Prueger 2, Jerry Hatfield 2,

More information

Homogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers

Homogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers Excert from the roceedings of the COMSOL Conference 1 aris Homogeneous and Inhomogeneous Model for Flow and Heat ransfer in orous Materials as High emerature Solar Air Receivers Olena Smirnova 1 *, homas

More information

ELECTRON DIFFRACTION

ELECTRON DIFFRACTION ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the

More information

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity AP Physics Multiple Choice Practice Electrostatics 1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity. A soli conucting sphere is given a positive charge Q.

More information

Lenny Jones Department of Mathematics, Shippensburg University, Shippensburg, Pennsylvania Daniel White

Lenny Jones Department of Mathematics, Shippensburg University, Shippensburg, Pennsylvania Daniel White #A10 INTEGERS 1A (01): John Selfrige Memorial Issue SIERPIŃSKI NUMBERS IN IMAGINARY QUADRATIC FIELDS Lenny Jones Deartment of Mathematics, Shiensburg University, Shiensburg, Pennsylvania lkjone@shi.eu

More information

7. Introduction to Large Sample Theory

7. Introduction to Large Sample Theory 7. Introuction to Large Samle Theory Hayashi. 88-97/109-133 Avance Econometrics I, Autumn 2010, Large-Samle Theory 1 Introuction We looke at finite-samle roerties of the OLS estimator an its associate

More information

for which the garbage will fall within the can. If the range is large, the shot is easy; if the range is small, the target is hard to hit.

for which the garbage will fall within the can. If the range is large, the shot is easy; if the range is small, the target is hard to hit. EN40: Dynamics an Vibrations Homework : Soling equations of motion for particles Due Friay Feb 18 011 School of Engineering Brown Uniersity 1. Suppose that you wish to throw a piece of garbage into a garbage

More information

Colin Cameron: Brief Asymptotic Theory for 240A

Colin Cameron: Brief Asymptotic Theory for 240A Colin Cameron: Brief Asymtotic Theory for 240A For 240A we o not go in to great etail. Key OLS results are in Section an 4. The theorems cite in sections 2 an 3 are those from Aenix A of Cameron an Trivei

More information

Hotelling s Two- Sample T 2

Hotelling s Two- Sample T 2 Chater 600 Hotelling s Two- Samle T Introduction This module calculates ower for the Hotelling s two-grou, T-squared (T) test statistic. Hotelling s T is an extension of the univariate two-samle t-test

More information

A Simple Exchange Economy with Complex Dynamics

A Simple Exchange Economy with Complex Dynamics FH-Kiel Universitf Alie Sciences Prof. Dr. Anreas Thiemer, 00 e-mail: anreas.thiemer@fh-kiel.e A Simle Exchange Economy with Comlex Dynamics (Revision: Aril 00) Summary: Mukherji (999) shows that a stanar

More information

To understand how scrubbers work, we must first define some terms.

To understand how scrubbers work, we must first define some terms. SRUBBERS FOR PARTIE OETION Backgroun To unerstan how scrubbers work, we must first efine some terms. Single roplet efficiency, η, is similar to single fiber efficiency. It is the fraction of particles

More information

Node-voltage method using virtual current sources technique for special cases

Node-voltage method using virtual current sources technique for special cases Node-oltage method using irtual current sources technique for secial cases George E. Chatzarakis and Marina D. Tortoreli Electrical and Electronics Engineering Deartments, School of Pedagogical and Technological

More information

CONTROL AND PERFORMANCE OF A NINE PHASE SYNCHRONOUS RELUCTANCE DRIVE

CONTROL AND PERFORMANCE OF A NINE PHASE SYNCHRONOUS RELUCTANCE DRIVE CONTRO AN PERFORMANCE OF A NINE PHASE SYNCHRONOUS REUCTANCE RIVE Abstract C.E. Coates*,. Platt** an V.J. Gosbell ** * epartment of Electrical an Computer Engineering Uniersity of Newcastle ** epartment

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

Keywords: pile, liquefaction, lateral spreading, analysis ABSTRACT

Keywords: pile, liquefaction, lateral spreading, analysis ABSTRACT Key arameters in seudo-static analysis of iles in liquefying sand Misko Cubrinovski Deartment of Civil Engineering, University of Canterbury, Christchurch 814, New Zealand Keywords: ile, liquefaction,

More information

3. High Temperature Gases FPK1 2009/MZ 1/23

3. High Temperature Gases FPK1 2009/MZ 1/23 3. High Temerature Gases FP 9/MZ /3 Terms and Concets Dissociation Diatomic element gases Bond energy, dissociation energy (enthali) Flood s dissociation diagram Vaorization, eaoration, boiling, sublimation

More information

Convergence Analysis of Terminal ILC in the z Domain

Convergence Analysis of Terminal ILC in the z Domain 25 American Control Conference June 8-, 25 Portlan, OR, USA WeA63 Convergence Analysis of erminal LC in the Domain Guy Gauthier, an Benoit Boulet, Member, EEE Abstract his aer shows how we can aly -transform

More information

9.3 Laminar Flat-Plate Boundary Layer: Exact Solution w-19

9.3 Laminar Flat-Plate Boundary Layer: Exact Solution w-19 9.3 Laminar Flat-Plate Boundary Layer: Exact Solution w-19 Laminar Flat-Plate Boundary Layer: Exact Solution The solution for the laminar boundary layer on a horizontal flat late was obtained by Prtl s

More information

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments

TMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary

More information

dn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential

dn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential Chem 467 Sulement to Lectures 33 Phase Equilibrium Chemical Potential Revisited We introduced the chemical otential as the conjugate variable to amount. Briefly reviewing, the total Gibbs energy of a system

More information

Analysis of Pressure Transient Response for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia

Analysis of Pressure Transient Response for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia roceedings World Geothermal Congress 00 Bali, Indonesia, 5-9 Aril 00 Analysis of ressure Transient Resonse for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia Jorge A.

More information

MATH , 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208

MATH , 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208 MATH 321-03, 06 Differential Equations Section 03: MWF 1:00pm-1:50pm McLaury 306 Section 06: MWF 3:00pm-3:50pm EEP 208 Instructor: Brent Deschamp Email: brent.eschamp@ssmt.eu Office: McLaury 316B Phone:

More information

Novel Algorithm for Sparse Solutions to Linear Inverse. Problems with Multiple Measurements

Novel Algorithm for Sparse Solutions to Linear Inverse. Problems with Multiple Measurements Novel Algorithm for Sarse Solutions to Linear Inverse Problems with Multile Measurements Lianlin Li, Fang Li Institute of Electronics, Chinese Acaemy of Sciences, Beijing, China Lianlinli1980@gmail.com

More information

State Estimation with ARMarkov Models

State Estimation with ARMarkov Models Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

The Nottingham eprints service makes this work by researchers of the University of Nottingham available open access under the following conditions.

The Nottingham eprints service makes this work by researchers of the University of Nottingham available open access under the following conditions. Li, Xia (2016) Internal structure quantification for granular constitutie modeling. Journal of Engineering Mechanics. C4016001/1-C4016001/16. ISSN 1943-7889 Access from the Uniersity of Nottingham reository:

More information

Setting up the Mathematical Model Review of Heat & Material Balances

Setting up the Mathematical Model Review of Heat & Material Balances Setting u the Mathematical Model Review of Heat & Material Balances Toic Summary... Introduction... Conservation Equations... 3 Use of Intrinsic Variables... 4 Well-Mixed Systems... 4 Conservation of Total

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow outflow Part I. Regular tasks for MAE598/494 Task 1 MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the

More information

arxiv: v1 [physics.data-an] 26 Oct 2012

arxiv: v1 [physics.data-an] 26 Oct 2012 Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch

More information

Probabilistic Learning

Probabilistic Learning Statistical Machine Learning Notes 11 Instructor: Justin Domke Probabilistic Learning Contents 1 Introuction 2 2 Maximum Likelihoo 2 3 Examles of Maximum Likelihoo 3 3.1 Binomial......................................

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

More information

PHOTOINDUCED SUPERCONDUCTIVITY IN CUPRATES

PHOTOINDUCED SUPERCONDUCTIVITY IN CUPRATES Int. J. hem. Sci.: 8(3,, 473-5 PHOOIDUED SUPERODUIIY I UPRAES AUJ UWAL, S. L. KAKAI *a an M. L. KALRA Deartment of Physics, University ollege of Science, M. L. Sukhaia University, UDAIPUR 33 (RAJ. IDIA

More information

Why are we interested in aerosols?

Why are we interested in aerosols? Chemistry an Aerosols Bernhar Vogel Chemistry an Aerosol bernhar.vogel@kit.eu Why are we intereste in aerosols? Aerosols have an imact on human health, BUT they have also an imact on climate an weather.

More information

An effective permeability model to predict field-dependent modulus of magnetorheological elastomers

An effective permeability model to predict field-dependent modulus of magnetorheological elastomers Available online at www.scienceirect.com Communications in Nonlinear Science an Numerical Simulation 13 (2008) 1910 1916 www.elsevier.com/locate/cnsns An effective ermeability moel to reict fiel-eenent

More information

The effect of nonvertical shear on turbulence in a stably stratified medium

The effect of nonvertical shear on turbulence in a stably stratified medium The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:

More information

INVESTIGATION OF THE DEISGN AND PERFORMANCE OF REPEATING SPACE TRACK CONSTELLATIONS

INVESTIGATION OF THE DEISGN AND PERFORMANCE OF REPEATING SPACE TRACK CONSTELLATIONS INVESTIGATION OF THE DEISGN AND PEFOMANCE OF EPEATING SPACE TACK CONSTELLATIONS Michelle K. Perez Avisor: Dr. Christoher Hall Virginia Polytechnic Institute an State University, Blacsburg, VA, 24060 Abstract

More information

J/psi and Upsilon production cross-sections at the ATLAS experiment

J/psi and Upsilon production cross-sections at the ATLAS experiment J/si an Usilon rouction crosssections at the ALAS exeriment Lancaster University, UK Email: susan.cheatham@cern.ch New ALAS measurements of romt an nonromt J/ψ an of Usilon rouction at a centreof mass

More information

3.4 Design Methods for Fractional Delay Allpass Filters

3.4 Design Methods for Fractional Delay Allpass Filters Chater 3. Fractional Delay Filters 15 3.4 Design Methods for Fractional Delay Allass Filters Above we have studied the design of FIR filters for fractional delay aroximation. ow we show how recursive or

More information

High speed wind tunnels 2.0 Definition of high speed. 2.1 Types of high speed wind tunnels

High speed wind tunnels 2.0 Definition of high speed. 2.1 Types of high speed wind tunnels Module Lectures 6 to 1 High Seed Wind Tunnels Keywords: Blow down wind tunnels, Indraft wind tunnels, suersonic wind tunnels, c-d nozzles, second throat diffuser, shocks, condensation in wind tunnels,

More information

The extreme case of the anisothermal calorimeter when there is no heat exchange is the adiabatic calorimeter.

The extreme case of the anisothermal calorimeter when there is no heat exchange is the adiabatic calorimeter. .4. Determination of the enthaly of solution of anhydrous and hydrous sodium acetate by anisothermal calorimeter, and the enthaly of melting of ice by isothermal heat flow calorimeter Theoretical background

More information

Model for Longitudinal Uncertainty during Controlled Time of Arrival Operations

Model for Longitudinal Uncertainty during Controlled Time of Arrival Operations Eleenth USA/Europe Air Traffic Management Research an Deelopment Sear (ATM015) Moel for Longituinal Uncertainty uring Controlle Time of Arrial Operations Dai De Smet Naigation an CNS Research Unit EUROCONTROL

More information

CHS GUSSET PLATE CONNECTIONS ANALYSES Theoretical and Experimental Approaches

CHS GUSSET PLATE CONNECTIONS ANALYSES Theoretical and Experimental Approaches EUROSTEEL 8, 3-5 Setember 8, Graz, Austria 561 CHS GUSSET PLATE CONNECTIONS ANALYSES Theoretical an Exerimental Aroaches Arlene M. S. Freitas a, Daniela G. V. Minchillo b, João A. V. Requena c, Afonso

More information

Reconstructions for some coupled-physics inverse problems

Reconstructions for some coupled-physics inverse problems Reconstructions for some couple-physics inerse problems Guillaume al Gunther Uhlmann March 9, 01 Abstract This letter announces an summarizes results obtaine in [8] an consiers seeral natural extensions.

More information

Q ABS (x,t s ) = S o /4 S(x)a p (x,t s ).

Q ABS (x,t s ) = S o /4 S(x)a p (x,t s ). Lecture 14 ICE The feedback of exanding and contracting ice sheets has often been offered as a lausible exlanation for how the contrasting climates of the glacialinterglacial times can be both (relatively)

More information

Competition in the health care market: a two-sided. approach

Competition in the health care market: a two-sided. approach Cometition in the health care market: a two-sie aroach Mario Pezzino Giacomo Pignataro Abstract Two ientical hositals comete for atients an octors choosing locations (e.g. secializations) on a Hotelling

More information

Digitally delicate primes

Digitally delicate primes Digitally elicate rimes Jackson Hoer Paul Pollack Deartment of Mathematics University of Georgia Athens, Georgia 30602 Tao has shown that in any fixe base, a ositive roortion of rime numbers cannot have

More information

MATH 2710: NOTES FOR ANALYSIS

MATH 2710: NOTES FOR ANALYSIS MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite

More information

Chapter 6. Thermodynamics and the Equations of Motion

Chapter 6. Thermodynamics and the Equations of Motion Chater 6 hermodynamics and the Equations of Motion 6.1 he first law of thermodynamics for a fluid and the equation of state. We noted in chater 4 that the full formulation of the equations of motion required

More information

INVESTIGATION OF AERATED SIPHON

INVESTIGATION OF AERATED SIPHON INVETIGATION OF AERATED IPHON Detlef Aigner, Hans-B. Horlacher Institute for Hydraulic Engineering and Alied Hydromechanics, Dresden Uniersity of Technology, 006 Dresden, Germany, hone: +49 35 46334397;

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

DEVELOPMENT OF A BRANCH AND PRICE APPROACH INVOLVING VERTEX CLONING TO SOLVE THE MAXIMUM WEIGHTED INDEPENDENT SET PROBLEM

DEVELOPMENT OF A BRANCH AND PRICE APPROACH INVOLVING VERTEX CLONING TO SOLVE THE MAXIMUM WEIGHTED INDEPENDENT SET PROBLEM DEVELOPMENT OF A BRANCH AND PRICE APPROACH INVOLVING VERTEX CLONING TO SOLVE THE MAXIMUM WEIGHTED INDEPENDENT SET PROBLEM A Thesis by SANDEEP SACHDEVA Submitted to the Office of Graduate Studies of Texas

More information

CREATION AND DEVELOPMENT OF THE NATIONAL HUMIDITY STANDARD OF GASES

CREATION AND DEVELOPMENT OF THE NATIONAL HUMIDITY STANDARD OF GASES CREATION AND DEVELOPMENT OF THE NATIONAL HUMIDITY STANDARD OF GASES Efimia LUCHIAN, Constantin BORDIANU, Grigore BUZUC, Victoria CEBAN National Institute of Metrology, Chisinau, Republic of Molova tel.:

More information

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2 Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z

More information

Sources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences

Sources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences Sources an Sinks of Available Potential Energy in a Moist Atmosphere Olivier Pauluis 1 Courant Institute of Mathematical Sciences New York University Submitte to the Journal of the Atmospheric Sciences

More information