The present study proposes a multi-objective framework for structure selection of nonlinear systems which are represented by polynomial NARX models. This framework integrates the key components of Multi-Criteria Decision Making (MCDM) which include preference handling, Multi-Objective Evolutionary Algorithms (MOEAs) and a posteriori selection. To this end, three well-known MOEAs such as NSGA-II, SPEA-II and MOEA/D are thoroughly investigated to determine if there exists any significant difference in their search performance. The sensitivity of all these MOEAs to various qualitative and quantitative parameters, such as the choice of recombination mechanism, crossover and mutation probabilities, is also studied. These issues are critically analyzed considering seven discretetime and a continuous-time benchmark nonlinear system as well as a practical case study of non-linear wave-force modeling. The results of this investigation demonstrate that MOEAs can be tailored to determine the correct structure of nonlinear systems. Further, it has been established through frequency domain analysis that it is possible to identify multiple valid discrete-time models for continuous-time systems. A rigorous statistical analysis of MOEAs via performance sweet spots in the parameter space convincingly demonstrates that these algorithms are robust over a wide range of control parameters.of correct terms which captures the system dynamics is one of the major challenges of the system identification [6,7], and, therefore, is the main focus of this investigation.The exhaustive search to solve the structure selection problem is intractable even for moderate number of terms 'n', as this would require an evaluation of 2 n term subsets/structures. The structure selection problem is often 'NP-Hard' [13] and therefore requires an efficient search strategy. Since the search for system structure is primarily dependent on the limited number of observed data, it involves bias-variance trade-off [2, 3, 7], i.e., too sparse a structure may yield a high bias error (under-fitting) whereas a very complex structure may yield a high variance error (over-fitting). Further, earlier investigations [14,15] show that the over-fitted structures tend to induce 'ghost' dynamics which are not present in the actual system. Pursuing these arguments, it is clear that determination of the size of the identified structure (i.e., number of terms or 'cardinality') is a critical aspect of the structure selection process. However, this issue has received relatively less attention in the existing investigations.Among the existing structure selection approaches, the Orthogonal Forward Regression (OFR) proposed by Billings and Korenberg [16] has extensively been studied. It is a sequential model building approach in which the most significant term is added in each step, which is determined on the basis of a statistical criterion. Over the years, this notion of the original OFR approach has been refined to further improve the search performance. These include OFR approaches...