The dataset supports the findings of the manuscript entitled ‘Exploring the Temporal Dynamics of Methane Ebullition in a Subtropical Freshwater Reservoir’. The results of the manuscript are based on continuous in-situ measurements conducted at Passaúna Reservoir, located in the southern part of Brazil (South America). The monitoring was carried out from 2017 to 2020, during which a comprehensive set of environmental variables was obtained from various studies. The primary objective of the manuscript was to comprehend the temporal dynamics of ebullition flux in the reservoir. Therefore, time series data of ebullition and relevant environmental variables were analyzed across different time scales, ranging from minutes to daily resolutions. Statistical and data driven models were tested to predict ebullition at varying time scales, considering the influence of environmental variables. The findings are discussed in the manuscript. Here, xlsx files and Matlab scripts are provided. The xlsx files contain time series data for environmental variables ('Environmental_Variables_TimeSeries') and ebullition flux ('Ebullition_Flux_TimeSeries'). Separate sheets were utilized for different time intervals, namely 5 minutes (dt=5min), 10 minutes (dt=10min), 1 hour (dt=1hr), and 1 day (dt=1d). Explanations and units are provided within the column labels, while 'NaN' denotes missing data. Matlab scripts for all the empirical models that were tested, as outlined in the manuscript's supporting information (Tables S1 and S2), have been made available (detailed in the following table). These scripts were developed using MatLab R2023a. File (.m) Description Main scripts TableS1_Empirical_Models_Literature Empirical models from the literature tested. (Table S1 in the manuscript) TableS2_FirstPart_Reffited_Models Empirical models from the literature refitted. (Table S2 in the manuscript) TableS2_SecondPart_New_Models New empirical models implemented. (Table S2 in the manuscript) TableS2_SecondPart_New_Models_GAM_Timescales Generalized additive models (GAM) applied to predict ebullition at different timescales. (Table S2 in the manuscript). PredictedR2_GAM_models Calculate the predicted R-squared for GAM. Functions (need for the main scripts) load_TimeSeries Import time series of environmental variables from excel sheets into MatLab. load_EbullitionTS Import time series of ebullition from excel sheets into MatLab. units_description Description containing units of the variables. Bin_xdata Creates data binning of X and Y based on X data. RelativeError Calculate the relative error of accumulated flux between measured and simulated. ANN5Neurons_Retrained_Deshmukh2014 Trained Artificial Neural Network based on the ANN proposed by Deshmukh et al. (2014). ANN_Passauna_20Neurons Trained Artificial Neural Network with the addition of more environmental variables. Financial support The field measurements were financed by the German Federal Ministry of Education and Research (BMBF, Grant 02WGR1431A), in the framework of the research project MuDak-WRM (https://www.mudak-wrm.kit.edu). This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brasil (CAPES) – Finance Code 001. Tobias Bleninger received productivity stipend from the National Council for Scientific and Technological Development (CNPq, grant no. 312211/2020-1, call no. 09/2020). Michael Männich received productivity stipend from the National Council for Scientific and Technological Development (CNPq, grant no. 308744/2021-7, call no. 04/2021). Andreas Lorke received financial support from the German Research Foundation (DFG, grant number LO1150/16-1).