The biosphere plays a significant role in climate mitigation. In the course of photosynthesis, plants sequester approximately a third of the carbon dioxide (CO2) that humans release into the atmosphere each year. The protection of terrestrial ecosystems is therefore essential for the global carbon, water, and energy cycles, as well as for offsetting atmospheric carbon dioxide emissions. In order to better preserve and restore global ecosystems for endangered biodiversity, species, and human well-being, it is imperative to accurately measure the ecosystem-atmosphere mass (e.g., carbon) and energy exchanges - i.e., fluxes. The limitations of current observational and/or estimation techniques, however, make it difficult to quantify fluxes for terrestrial ecosystems, particularly those under management. Accordingly, it remains unclear what the exact values of global ecosystem carbon fluxes are and how they vary over time. Eddy covariance (EC) is a micrometeorological method that provides continuous and direct flux measurements. However, these measurements are only representative of the tower footprint, which ranges in radius from hundreds to several kilometres. It is common to upscale or extrapolate eddy covariance fluxes by integrating satellite remote sensing data in order to determine fluxes over a larger area. It should be noted, however, that many current flux upscaling studies focus on gross primary productivity rather than the net exchange, which necessitates deducting ecosystem respiration from the estimates. In addition, present flux upscaling estimations have a relatively coarse spatial resolution for the purpose of monitoring field-scale fluxes and assessing management activities. The goal of this thesis is to fill the knowledge gaps and develop continuous field-scale flux estimates for ecosystems under management practices such as livestock grazing. In order to accomplish this goal, the thesis will answer three major questions: 1) What factors determine the performance of flux upscaling; 2) How can flux gaps of different lengths be filled reliably; 3) Can field-scale estimates accurately reflect livestock grazing's effects on ecosystem fluxes? Chapter 2 examines the feasibility of consistently estimating global carbon, water, and energy fluxes with the flux upscaling technique. This 'consistent upscaling' means using the same machine learning algorithm as well as satellite and meteorological data over the course of the flux upscaling. According to the results, the consistent flux upscaling explains 70 % of the daily variability of carbon, water, and energy fluxes. Furthermore, this chapter analyses the effects of algorithms, predictor sets, and eddy covariance itself on flux upscaling performance, concluding that the spatiotemporal sampling density of eddy covariance towers is critical to the flux upscaling accuracy. The eXtreme Gradient Boosting (Xgboost) outperforms other machine learning algorithms with at least 14 % higher accuracy and an average of 90 % shorter run time. Using different predictors results in a difference in flux upscaling accuracy of less than 4 %. The spatial sampling density of eddy covariance towers, however, can affect the flux upscaling accuracy to a degree larger than 50 %. This chapter lays the technical foundation for the entire thesis. Chapter 3 improves the method for filling in eddy covariance data gaps since incomplete data may compromise the accuracy of both eddy covariance and flux upscaling. In this chapter, a study at more than a hundred eddy covariance towers around the world demonstrates that the random forest algorithm can fill gaps even longer than a month. Depending on how many environmental drivers are used, random forest has a 15 % to 30 % improvement in gap-filling accuracy over the research standard method. Random forest exhibits greater stability in filling long gaps, as the gap length increases from one day to a month, its accuracy drop is 81 % less than the standard research method. In Chapter 4, the gap-filling study from the previous chapter is extended to methane fluxes and ecosystems under various management practices (e.g., livestock grazing). The random forest method is 29 % to 54 % more accurate than the standard method when it comes to gap-filling methane fluxes. When using modelled environmental drivers, random forest maintains 90 % of gap-filling accuracy compared to using measured drivers. By doing so, it is possible to promote the use of eddy covariance applications when measured environmental drivers are not available. When gaps exceed three months, the R2 of random forest gap-filling remains 0.7. In addition, the random forest is found to preserve environment-flux responses well in the filled gaps. Chapters 3 and 4 contribute to the aim of this thesis - i.e., reliable field-scale flux estimation - by providing complete eddy covariance time series. Chapter 5 examines how satellite spatial resolution affects the extrapolation of eddy covariance fluxes from tower footprints to the entire field. It also evaluates the proposed field-scale flux estimates in light of the performance of reproducing flux spatial variability and recognising grazing periods. According to the results, satellite spatial resolution has little impact on the estimated flux sum of the whole field (R2 difference ≪ 0.1), but it has a great impact on the estimated flux spatial variability. It is observed that the flux estimates are highly consistent with measurements (R2 = 0.7 and annual bias < 2 Mg ha-1 yr-1 for carbon fluxes), and they clearly distinguish between fluxes occurring during grazing periods and those occurring during non-grazing periods. Since Sentinel-2's imagery textiles are of greater clarity, it is superior to other satellite platforms in extrapolating tower-based fluxes (in southwest England). In conclusion, this thesis accomplishes its goal of filling the knowledge gap in reliably estimating ecosystem fluxes with multi-source observation and machine learning. The results confirm that eddy-covariance carbon, water, and energy fluxes can be scaled up in a consistent manner. Specifically, this research emphasises the importance of spatially well-distributed towers for upscaling eddy covariance fluxes, particularly in tropical regions and/or ecosystem types such as evergreen broadleaf forests. This research also identifies the validity of the random forest algorithm in gap-filling - random forest is effective in filling varying-length gaps while maintaining relationships between fluxes and environmental variables. Furthermore, this research highlights the importance of satellite platform selection in quantifying the spatial variability of field-scale flux estimates. As a final note, it validates machine learning algorithms for obtaining fluxes that are resistant to disturbances caused by management activities such as livestock grazing. There are important implications of this research for climate mitigation efforts that are supported by the management and protection of terrestrial ecosystems. Findings from this research are beneficial for pinpointing carbon uptake capacity of global land ecosystems, for determining where to install future eddy covariance towers, and for evaluating the effects of other management activities on ecosystem fluxes.